neperfepx / neper

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

-transform rmcell is buggy #569

Open rquey opened 2 years ago

rquey commented 2 years ago

-transform rmcell very often generates errors. For instance:

neper -T -n 2 -transform "rmcell(id>1)"

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.5.0
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] -tesrformat ascii
Info   : [com line] -n 2 -transform rmcell(id>1)
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :   - Setting seeds... 
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   : Transforming tessellation...
Info   :   - Removing cells...
Info   :     > Removed 1 cells...
Info   : Writing results...
Info   :     [o] Writing file `n2-id1.tess'...
Segmentation fault (core dumped)

Of note, a possible workaround is to use -meshpoly in -M (providing the complementary expression as argument).