neperfepx / neper

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

Periodic Tessellation error #889

Open rsooraj opened 4 months ago

rsooraj commented 4 months ago

Describe the bug

Trying to produce periodic tessellation, with periodicity along X-direction, Neper shows the following error after running for some number of iterations.

The commands used and the error displayed are as follows.

$ neper -T -n from_morpho -oricrysym cubic -periodicity x -domain "cube(2,5,5)" -morpho voronoi,"diameq:lognormal(mean=0.11,sigma=0.042),1-sphericity:lognormal(0.145,0.03)" -morphooptistop itermax=1000000 -statcell diameq,body,npolynb,sphericity,convexity,nseednb -o c_2x5x5_per_lognorm_mu_11e-2_sig_42e-3_10L_iters_VOR -oridescriptor euler-bunge -format ply,tess

======================== N e p e r ======================= Info : A software package for polycrystal generation and meshing. Info : Version 4.5.1-4 Info : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full) Info : Running on 64 threads. Info : https://neper.info Info : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey. Info : No initialization file found (`/home/sooraj/.neperrc'). Info : --------------------------------------------------------------- Info : MODULE -T loaded with arguments: Info : [ini file] (none) Info : [com line] -n from_morpho -oricrysym cubic -periodicity x -domain cube(2,5,5) -morpho voronoi,diameq:lognormal(mean=0.11,sigma=0.042),1-sphericity:lognormal(0.145,0.03) -morphooptistop itermax=1000000 -statcell diameq,body,npolynb,sphericity,convexity,nseednb -o c_2x5x5_per_lognorm_mu_11e-2_sig_42e-3_10L_iters_VOR -oridescriptor euler-bunge -format ply,tess Info : --------------------------------------------------------------- Info : Reading input data... Info : Creating domain... Info : Creating tessellation... Info : . Number of cells: 47696 Info : - Setting seeds... 100% Info : - Generating crystal orientations... Info : - Running tessellation... Info : > Initial solution: f =0.358254279 Info : > Iteration 9030: fmin=0.343747920 f=0.343747920Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you.

Aborted (core dumped)