matlab_FEM_dismesh

  DistMesh giving a singular FEM matrix?(Hi Anyone here with some experience of using DistMesh in finite element code? I m solving the scalar Helmoltz equation in an annular region in 2D: rho_1 < rho < rho_2 where rho_1 is the radius of a Perfect Electric Conducting cylinder and rho_2 is where the mesh is truncated using an Absorbing Boundary condition. Now, I ve written code to mesh the region myself (by dividing the region into annular rings and picking a fixed number of points on each ring) and I ve written another program which uses DistMesh to mesh the region. I get a nice solution using my own meshing code but the FEM matrix becomes singular when I use DistMesh. This is the DistMesh Code I used: Circle with hole rABC=1.5 rCyl=0.5 fdstring=sprintf( ddiff(dcircle(p,0,0, f),dcircle(p,0,0, f)) ,rABC,rCyl) fd=inline(fdstring, p ) box=[-2,-2 2,2] [p,N]=distmesh2d(fd,@huniform,0.04,box,[]) I tried changing box and the 0.04 value (initial edge length). But every time, the matrix becomes sin)

boundedges.m,622,2006-02-14
circumcenter.m,469,2006-02-14
COPYRIGHT.TXT,1078,2010-06-05
dblock.m,215,2006-02-14
dcircle.m,153,2006-02-14
ddiff.m,121,2006-02-14
dellipse.cpp,1841,2006-02-14
dellipse.dll,7680,2004-06-30
dellipse.mexglx,7291,2004-07-01
dellipsoid.cpp,2529,2006-02-14
dellipsoid.dll,8704,2004-06-30
dellipsoid.mexglx,7922,2004-07-01
dexpr.m,1046,2006-02-14
dintersect.m,121,2006-02-14
distmesh2d.m,3969,2006-02-14
distmeshnd.m,3349,2006-02-14
dmatrix.m,169,2006-02-14
dmatrix3d.m,184,2006-02-14
dpoly.m,781,2006-02-14
drectangle.m,182,2006-02-14
drectangle0.m,421,2006-02-14
dsegment.cpp,1309,2006-02-14
dsegment.dll,5632,2004-06-30
dsegment.mexglx,6124,2004-07-28
dsphere.m,173,2006-02-14
dunion.m,121,2006-02-14
fixmesh.m,443,2006-02-14
hmatrix.m,172,2006-02-14
hmatrix3d.m,187,2006-02-14
huniform.m,137,2006-02-14
meshdemo2d.m,3774,2006-02-14
meshdemond.m,1038,2006-02-14
mkt2t.m,960,2006-02-14
protate.m,154,2006-02-14
pshift.m,142,2006-02-14
simpplot.m,868,2006-02-14
simpqual.m,2186,2006-02-14
simpvol.m,714,2006-02-14
surftri.m,694,2006-02-14
uniformity.m,213,2006-02-14