[Gmsh] gmsh bug
Alessandro Mori
mori at ing.unifi.it
Tue May 15 11:17:14 CEST 2001
Dear Dr. Geuzaine,
I report you a bug found in gmsh. It is important to note that including
directly the function
in test1.geo there is no error.
I have downloaded your program a couple of weeks ago, but I have already
found your program
very useful, and very impressive.
Is there any documentation or list of reserved words apart the tutorial?
Gmsh 1.17
-------------------------------------------------------
Operating system : Linux 2.2.16-3
Build date : Tue Apr 17 13:34:12 CEST 2001
Build host : elap21.montefiore.ulg.ac.be
Packager : geuzaine
-------------------------------------------------------
Geometry
Loading 'test1.geo'
Info : Including 'test2.geo'
Fatal Error : Segmentation violation (invalid memory reference)
: ------------------------------------------------------
: You have discovered a bug in Gmsh! You may report it
: by e-mail (together with any helpful data permitting to
: reproduce it) to <gmsh at geuz.org>
--
--------------------------------------------------------------------------
Alessandro Mori Laboratorio Elettromagnetismo e Microonde
Dr.(ph.d.) Ing. Dipartimento Elettronica e Telecomunicazioni
Universita' degli Studi di Firenze
Fax : +39.055.4796441 Address:
Phone : +39.055.4796381 v. Santa Marta, 3
e-mail: mori at ing.unifi.it 50139 Firenze (Italy)
--------------------------------------------------------------------------
-------------- next part --------------
Include "test2.geo";
lc = 0.005;
hp = 0.0315;
wp0 = 0.05;
lp0 = 0.05;
ncr = 2;
wpc[] = {0.05,0.05};
lpc[] = {0.05,0.05};
rpc[] = {0.375,0.65};
npc[] = {8,12};
Call CircPatch;
-------------- next part --------------
Function CircPatch
nptc = newp;
nrgc = newreg;
Point(nptc) = { lp0/2, wp0/2, hp, lc};
Point(nptc+1) = { lp0/2, -wp0/2, hp, lc};
Point(nptc+2) = {-lp0/2, -wp0/2, hp, lc};
Point(nptc+3) = {-lp0/2, wp0/2, hp, lc};
Line(nrgc) = {nptc,nptc+1};
Line(nrgc+1) = {nptc+1,nptc+2};
Line(nrgc+2) = {nptc+2,nptc+3};
Line(nrgc+3) = {nptc+3,nptc};
Line Loop(nrgc) = {nrgc,nrgc+1,nrgc+2,nrgc+3};
Plane Surface(nrgc) = {nrgc};
Return