[Gmsh] Questions gmsh + bug report

Christophe Geuzaine geuzaine at acm.caltech.edu
Wed Mar 3 23:20:17 CET 2004


Xavier Dardenne wrote:

> Bonjour,
> 
> J'utilise gmsh 1.51 sous RedHat 9.0 installé à partir du package RPM.
> 
> 1)
> 
> Je voudrais tout d'abord signaler un bug : 
> Systématiquement, si la dernière ligne du fichier .geo est commentée
> (//), gmsh plante à l'ouverture de ce fichier.

Good catch! I fixed it.

> 
> 
> 2)
> 
> Depuis que j'ai mis à jour ma version 1.44 avec la 1.50 puis la 1.51,
> j'ai constaté ceci :
> Lorsqu'une surface est recombinée en quadrilatères (Recombine Surface),
> et qu'ensuite, cette surface est copiée (ex:
> Translate{.,.,.}{Duplicata{Surface{.};}}), dans la version 1.44, la
> copie est aussi recombinée et pour les versions plus récentes, je suis
> obligé de recombiner chaque copie. Ce qui pour moi est assez gênant vu
> que je ne connaît pas à priori les numéros de surface des copies.
> Existe-t-il un moyen de remédier à ceci ?

We actually decided to stop copying the meshing constraints for
"duplicated" surfaces because it led to too many problems. For example
it made it impossible to create certain extruded geometries, as there
was no way to remove the mesh constraints a posteriori.

The solution to your problem is probably to use the information returned
by the transformation commands, e.g.

s[] = Translate {1,0,0} { Duplicata{ Surface{1}; } } ;
Recombine Surface {1,s[0]};


> 
> 3)
> 
> Lors du post-processing, lorsqu l'on utilise des Scalar Triangles, il
> suffit de donner une valeur a chaque sommet et gmsh fait une
> interpolation linéaire à l'intérieur du triangle. Y aurait-il un moyen
> rapide de changer cette interpolation linéaire en quadratique ? Par
> exemple en changeant légèrement le code source (appel de fct. openGL)?
> 

Yes, it would be possible. But since you also need more data to
interpolate to higher order, it means that we need to introduce new
(or, at least modify the existing) primitives, change the file format, etc.

The best way to deal with this at the moment is to subdivide your
primitives *before* writing the post-processing files. (And that's
actually what we would need to do in Gmsh too, if we wanted to treat
higher order primitives: OpenGL does not know anything about
non-straight lines or non-flat surfaces...).

Thanks for the bug report!

Take care,

Christophe

-- 
Christophe A. Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org