neperfepx / neper

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

2D cube cut #877

Open ruibmpinto opened 4 months ago

ruibmpinto commented 4 months ago

Maybe I was overly optimistic, but I was trying to generate a slightly complex 2D tesselation, something like:

example_2d_geom

I tried to follow roughly what was done here: https://github.com/neperfepx/neper/discussions/826

The command I ran and its output can be found below (I have also tried running with a cube with a finite corner radius and I get the same output). Is there a way to create such a geometry?

Thank you

$ neper -T -n 1000 -dim 2 -domain "square(0.1,0.1125)" -transform "cut(cylinder(0.028,0.0095,0,0,0,1,0.003),cylinder(0.028,0.0845,0,0,0,1,0.003),cylinder(0.028,0.0405,0,0,0,1,0.003),cube(0.025,0.1,0,0.0095,0,0.1,0.00))"

======================== N e p e r ======================= Info : A software package for polycrystal generation and meshing. Info : Version 4.8.2-9 Info : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full) Info : Running on 20 threads. Info : https://neper.info Info : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey. Info : No initialization file found (`/home/ruibmpinto/.neperrc'). Info : --------------------------------------------------------------- Info : MODULE -T loaded with arguments: Info : [ini file] (none) Info : [com line] -n 1000 -dim 2 -domain square(0.1,0.1125) -transform cut(cylinder(0.028,0.0095,0,0,0,1,0.003),cylinder(0.028,0.0845,0,0,0,1,0.003),cylinder(0.028,0.0405,0,0,0,1,0.003),cube(0.025,0.1,0,0.0095,0,0.1,0.00)) Info : --------------------------------------------------------------- Info : Reading input data... Info : Creating domain... Info : Creating tessellation... Info : - Setting seeds... Info : - Running tessellation... Info : Generating crystal orientations... Info : Transforming tessellation... Info : - Cutting (experimental)... Info : > cylinder(0.028,0.0095,0,0,0,1,0.003)... Info : > cylinder(0.028,0.0845,0,0,0,1,0.003)... Info : > cylinder(0.028,0.0405,0,0,0,1,0.003)... Info : > cube(0.025,0.1,0,0.0095,0,0.1,0.00)... Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you.