neperfepx / neper

Polycrystal generation and meshing
http://neper.info
GNU General Public License v3.0
214 stars 53 forks source link

Find a bug while doing 2d meshing #937

Open limorit opened 2 weeks ago

limorit commented 2 weeks ago

Describe the bug

I'm meshing a .tesr file whose data come from a EBSD test. There is no bug when I do -elttype quad, but when I change it to -elttype tri, there comes an Error.

Reproduce

Here is an example reaults of -elttype quad :

$ neper -M T4_ori_neper.tesr -elttype quad -rcl 0.07 -dim 2 -order 1 -format msh

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.10.2-1
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 28 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2024, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/xxx/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -M loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] T4_ori_neper.tesr -elttype quad -rcl 0.07 -dim 2
         -order 1 -format msh
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   : Loading input data...
Info   :   - Loading raster tessellation...
Info   :     [i] Parsing file `T4_ori_neper.tesr'...
Info   :     [i] Parsed file `T4_ori_neper.tesr'.
Info   : Meshing...
Info   :   - Preparing...

Info   :   - 2D meshing...
Info   : Searching nsets and fasets...
Info   : Writing mesh results...
Info   :   - Preparing mesh...
Info   :   - Mesh properties:
Info   :     > Node number:    74418
Info   :     > Elt  number:    73868
Info   :     > Mesh area:    0.026
Info   :   - Writing mesh...
Info   :     [o] Writing file `T4_ori_neper.msh'...
Info   :     [o] Wrote file `T4_ori_neper.msh'.
Info   : Elapsed time: 0.145 secs.
========================================================================

And here is the an example reaults of -elttype tri :

$ neper -M T4_ori_neper.tesr -elttype tri -rcl 0.07 -dim 2 -order 1 -format msh

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.10.2-1
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 28 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2024, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/xxx/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -M loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] T4_ori_neper.tesr -elttype tri -rcl 0.07 -dim 2
         -order 2 -format msh
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   : Loading input data...
Info   :   - Loading raster tessellation...
Info   :     [i] Parsing file `T4_ori_neper.tesr'...
Info   :     [i] Parsed file `T4_ori_neper.tesr'.
Info   :   - Reconstructing topology...
Info   :     > Reconstructing 1D mesh... 100%
Info   :     > Reconstructing 0D mesh... 100%
Error  : You have discovered a bug in Neper!  Please file an issue at
         https://github.com/rquey/neper/issues.  Thank you.

Aborted (core dumped)