Math Forum Discussions

srao Rao

Posts: 28
Registered: 7/16/10

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

Plain Text

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