Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Small domain bug #482

Open
harshsharma21 opened this issue May 24, 2022 · 2 comments
Open

Small domain bug #482

harshsharma21 opened this issue May 24, 2022 · 2 comments

Comments

@harshsharma21
Copy link

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
Copy link
Member

rquey commented May 24, 2022

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.

@rquey rquey closed this as completed May 24, 2022
@rquey rquey reopened this May 24, 2022
@harshsharma21
Copy link
Author

Thank you

@rquey rquey changed the title You have discovered a bug in Neper Small domain bug May 31, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants