[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