I am trying to evaluate the grid function at points with some highly refined mesh. However, I find that when the mesh is fine enough, the program will crash after evaluating the grid function at several hundred points. You can find example codes showing the problem. Here I have a square with the top boundary replaced by 3 steps of Koch snowflake. And I want to evaluate the grid function on this pre-fractal structure. But the program just crashes when I try to run it.
However, It still works if I replace the mesh with a coarser one.
I wonder if it is that I am not coding it in the right way so it does not work.
Thanks in advance for any help!