[Gmsh] unstable behaviour with 2nd order elements

Christophe Geuzaine cgeuzaine at ulg.ac.be
Wed Nov 19 18:53:28 CET 2008


Francisco Gilabert Villegas wrote:
> Dear all,
> 
> I have the following problem:
> 
> Using the version 2.2.5, gmsh is abruptly closed
> if the simple code listed below try to be processed
> to create a 2D mesh using 2nd order elements:
> 
> //--------------------------------------------
> mm = 0.001;
> Lx = 1000*mm;
> Ly = 1000*mm;
> Lz = 1000*mm;
> lc = 100*mm;
> p1=newp; Point(p1)={-Lx/2,-Ly/2,-Lz/2,lc};
> p2=newp; Point(p2)={-Lx/2,+Ly/2,-Lz/2,lc};
> p3=newp; Point(p3)={+Lx/2,+Ly/2,-Lz/2,lc};
> p4=newp; Point(p4)={+Lx/2,-Ly/2,-Lz/2,lc};
> l1=newl; Line(l1) = {p1, p2};
> l2=newl; Line(l2) = {p2, p3};
> l3=newl; Line(l3) = {p3, p4};
> l4=newl; Line(l4) = {p4, p1};
> ll1=newll; Line Loop(ll1) = {l1, l2, l3, l4};
> s1=news; Plane Surface(s1) = {ll1};
> Recombine Surface{s1};
> //--------------------------------------------
> 
> The older version 1.65 works well and it successfully
> creates 2D and 3D meshes with 2nd order elements.
> 
> My current OS is OpenSuse 11.0. Just for cheking,
> I executed the same piece of code using the same
> version 2.2.5 under Windows Vista Business. Unfortunately,
> the result was the same: suddenly, gmsh is aborted without
> emitting any message or track to follow.
> 
> It seems to me that some extremely severe error occurs.
> Do you know any reason which is provoking this situation ?
> Am I doing something wrong ?
> 

No, it seems that we shipped 2.2.5 with a severe bug in the high order 
code on Windows and Linux (on Mac everything seems OK).

Could you try the latest nightly snapshot and tell us if you still 
observe the problem?

Thanks,

Christophe







> Thank you in advance
> 
> 
> Francisco
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> Dr. Francisco A. Gilabert Villegas
> Mecánica de Medios Continuos y Modelización
> 
> Instituto de Tecnología Cerámica
> Campus Universitario Riu Sec | Avda. Vicent Sos Baynat, s/n
> 12006 Castellón (Spain)
> Tel. 34 964 34 24 24
> Fax 34 964 34 24 25
> 
> CONFIDENCIALIDAD
> La información contenida en este mensaje y/o archivo(s) adjunto(s) es 
> confidencial/privilegiada y está destinada
> a ser leída sólo por la(s) persona(s) a la(s) que va dirigida. Si usted 
> lee este mensaje y no es el destinatario señalado,
> el empleado o el agente responsable de entregar el mensaje al 
> destinatario, o ha recibido esta comunicación por
> error le informamos que está totalmente prohibida, y puede ser ilegal, 
> cualquier divulgación, distribución o
> reproducción de esta comunicación, y le rogamos que nos lo notifique 
> inmediatamente y nos devuelva el mensaje
> original al e-mail: itc at itc.uji.es. Gracias.
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine