neperfepx / neper

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

Small domain bug #482

Open harshsharma21 opened 2 years ago

harshsharma21 commented 2 years ago

neper -T -n 100 -domain "cube(0.00000001,0.00000001,0.00000001)" -morpho "diameq:lognormal(0.5,0.05),1-sphericity:lognormal(0.145,0.03)" -o 100grains -format 'vtk' -statcell vol,diameq,sphericity

======================== N e p e r ======================= Info : A software package for polycrystal generation and meshing. Info : Version 4.4.2-5 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 : No initialization file found (`/home/m200545me/.neperrc'). Info : --------------------------------------------------------------- Info : MODULE -T loaded with arguments: Info : [ini file] (none) Info : [com line] -n 100 -domain cube(0.00000001,0.00000001,0.00000001) -morpho diameq:lognormal(0.5,0.05),1-sphericity:lognormal(0.145,0.03) -o 100grains -format vtk -statcell vol,diameq,sphericity Info : --------------------------------------------------------------- Info : Reading input data... Info : Creating domain... Info : Creating tessellation... Info : - Setting seeds... 100% Info : - Generating crystal orientations... Info : - Running tessellation... Info : > Initial solution: f =2.815470205 Info : > Iteration 6758: fmin=0.001703365 f=0.001703365Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you.

Aborted (core dumped)

rquey commented 2 years ago

This is a bug... but it comes from your use of a very small domain (10^-8). A workaround is to generate the tessellation in the unit domain and then scale it down using -transform "scale(1e-8,1e-8,1e-8)".

Also, -statcell sphericity is not available for tesrs.

harshsharma21 commented 2 years ago

Thank you