search

neperfepx / neper

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

Core dumped when creating a mesh image with cohesive elements #377

Closed MiguelVCarvalho closed 2 years ago

MiguelVCarvalho commented 2 years ago

Describe the bug

When attempting to create an image from a .msh file that contains cohesive elements an error occurs:

Info   :     [i] Parsing file `file.msh'...
Error  : Reading of mesh file failed! (Element quantity does not match (leftelts = 23.)
./neper_issue.sh: line 7: 30958 Aborted                 (core dumped) neper --rcfile none -V file.msh -dataelsetcol 'ori:internal' -print file_mesh

The meshless image (using the .tess file) is generated without a problem and, if the option --interface cohesive is removed from the mesh command, the mesh image is also created as expected.

To Reproduce

# Tesselation and Meshing

neper --rcfile none -T -n 5 -id 1 -dim 2 -morpho graingrowth -o file -format tess,ori -oridescriptor e
neper --rcfile none -M -tmp /run/user/1000 file.tess -elttype tri -rcl 1.0 -dim 2 -order 1 -format msh -interface cohesive

# Visualization

neper --rcfile none -V file.tess -datacellcol ori -print file_tess
neper --rcfile none -V file.msh  -dataelsetcol 'ori:internal' -print file_mesh

# Output regarding the last `-V` command:

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.2.1-29
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 12 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2021, and GNU GPL'd, by Romain Quey.
Info   : Ignoring initialization file.
Info   : ---------------------------------------------------------------
Info   : MODULE  -V loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] file.msh -dataelsetcol ori:internal -print file_mesh
Info   : ---------------------------------------------------------------
Info   : Loading mesh...
Info   :     [i] Parsing file `file.msh'...
Error  : Reading of mesh file failed! (Element quantity does not match (leftelts = 23.)
./neper_issue.sh: line 7: 30958 Aborted                 (core dumped) neper --rcfile none -V file.msh -dataelsetcol 'ori:internal' -print file_mesh
rquey commented 2 years ago

Thanks for the bug report (apart from the unnecessary options ;-)). Now fixed.