neperfepx / neper

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

Bug when using different morphology and a transform sphere cut #584

Open shahmeerbaweja92 opened 2 years ago

shahmeerbaweja92 commented 2 years ago

Describe the bug

I am trying to reproduce an RVE with a centered, small spherical void inside it using cubic morphology instead of using 'voronoi', and then later mesh with c3d8 (hex) elements.

$ neper -T -n from_morpho -transform "cut(sphere(0.5,0.5,0.5,0.106))" -morpho "cube(7)" -format tess -o mgvoidsmall_cube

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.4.2-33
Info   : Built with: gsl|muparser|opengjk|nlopt|libscotch
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : No initialization file found (`/home/sbaweja/.neperrc').
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n from_morpho -transform
         cut(sphere(0.5,0.5,0.5,0.106)) -morpho cube(7) -format tess -o
         mgvoidsmall_cube
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :   - Running tessellation...
Info   :   - Generating crystal orientations...
Info   : Transforming tessellation...
Info   :   - Cutting (experimental)...
Info   :     > sphere(0.5,0.5,0.5,0.106)...
Error  : You have discovered a bug in Neper!  Please file an issue at
         https://github.com/rquey/neper/issues.  Thank you.

Aborted (core dumped)