neperfepx / neper

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

Bug of generating circle #483

Closed 2020-cai closed 2 years ago

2020-cai commented 2 years ago

======================== N e p e r ======================= Info : A software package for polycrystal generation and meshing. Info : Version 4.0.2 Info : Built with: gsl|muparser|opengjk|nlopt Info : http://neper.info Info : Copyright (C) 2003-2020, and GNU GPL'd, by Romain Quey. Info : --------------------------------------------------------------- Info : MODULE -T loaded with arguments: Info : [ini file] (none) Info : [com line] -n 100 -morpho diameq:0.71lognormal(8,0.5)+0.21lognormal(8,0.5)+0.08*lognormal(8,0.5) -regularization 1 -dim 2 -domain circle(50.8) -group mode -o BC 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 =16.695844393 Info : > Iteration 6451: fmin=0.000004945 f=0.000004953 Info : > Reached `eps' criterion. Info : Regularizing tessellation... Info : - loop 1/2: 65% del=51Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you.

Aborted

rquey commented 2 years ago

Thanks @2020-cai for the bug report.

For a circular domain, regularization generally cannot be pushed as far as for a standard, square domain. Use -rsel 0.5 to stop it earlier (it has failed at 65%, so it can theoretically be pushed to close to 0.65).

2020-cai commented 2 years ago

Thank you

---Original--- From: "Romain @.> Date: Wed, Jun 1, 2022 00:03 AM To: @.>; Cc: @.**@.>; Subject: Re: [neperfepx/neper] Bug of generating circle (Issue #483)

Thanks @2020-cai for the bug report.

For a circular domain, regularization generally cannot be pushed as far as for a standard, square domain. Use -rsel 0.5 to stop it earlier (it has failed at 65%, so it can theoretically be pushed to close to 0.65).

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were mentioned.Message ID: @.***>