Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.



pde toolbox
Posted:
Mar 14, 2013 2:41 PM


Hi,
I create a geometry using the pde toolbox GUI. I then export the decomposed geometry matrix g into the command space. Now I try to use the [p,e,t] = initmesh(g) and it works fine.
However, if I copy the same g matrix from the command space to an .m file and try to use initmesh I receive the following error:  ??? Error using ==> pdevoron at 32 Geometry error.
Error in ==> initmesh at 208 [p,t,c,h]=pdevoron(p,t,[],h,x,y,tol,Hmax,Hgrad);
Error in ==> mhp at 63 [p,e,t] = initmesh(g);  I examined my geometry and found that there was a gap. Theoretically the two points must overlap since they are the same. Is this a visualization problem or is there really a gap in my boundary? How can I close this gap? Why does it work in the command space but not when I use the matrix in a .m file. Thanks. Below is my geometry matrix:
g =[
2 2 2 2 2 2 2 4 1 0 0.2242 0.7154 0 0.0950 0.4698 1 1 0.2242 0.7154 1 0.0950 0.4698 0.8447 0.8447 0 0 0 0 0 0.0950 0.4698 0.4698 0.4698 0 0 0 0.0950 0.4698 0.8447 0.8447 1 1 1 1 0 0 0 1 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0.4698 0 0 0 0 0 0 0 0.4698 0 0 0 0 0 0 0 0.5302 0 0 0 0 0 0 0 0.5302 0 0 0 0 0 0 0 0 ];
Thanks,  Sai



