neperfepx / neper

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

Error in tessellation #376

Closed AdityaVuppala closed 2 years ago

AdityaVuppala commented 2 years ago

I use the following tessellation command and I get an unknown error. Any anyone help me figure out the issue.

$ neper -T -n 'from_morpho' -id 1 -dim 2 -domain 'square(0.0051, 0.0014)' -regularization 1 -morpho "diameq:lognormal(0.000245,0.0000001),sphericity:lognormal(0.145,0.03,1-x)"

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.2.0
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2021, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/NAS3/vuppala/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] -reg 1
Info   : [com line] -n from_morpho -id 1 -dim 2 -domain square(0.00508,
         0.0014) -regularization 1 -morpho
         diameq:lognormal(0.000245,0.0000001),sphericity:lognormal(0.145,0.03,1-x)
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :       . Number of cells: 150
Info   :   - Setting seeds... 100%
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =82.702669983
Info   :     > Iteration  11494: fmin=37.258408864 f=37.306633120
Info   :     > Iteration  22838: fmin=35.705562086 f=35.710851523 (2 subplex)
Info   :     > Iteration  27127: fmin=35.432430843 f=35.437606098 (3 subplex)
Info   :     > Iteration  31451: fmin=35.302779513 f=35.306960803 (4 subplex)
Info   :     > Iteration  32926: fmin=35.295377739 f=35.329510631 (5 subplex)
Info   :     > Iteration  34421: fmin=35.289242222 f=35.363198234 (6 subplex)
Info   :     > Iteration  35423: fmin=35.289242222 f=35.326945072 (7 subplex)
Info   :     > Iteration  38482: fmin=34.959298598 f=34.966646296 (8 praxis)
Info   :     > Iteration  39900: fmin=34.958681103 f=35.025680695 (9 subplex)
Info   :     > Iteration  40902: fmin=34.958681103 f=34.995289284 (10 subplex)
Info   :     > Iteration  42804: fmin=34.896056014 f=34.930809132 (11 praxis)
Info   :     > Iteration  43806: fmin=34.896056014 f=34.944855445 (12 subplex)
Info   :     > Iteration  45623: fmin=34.883181334 f=34.971881020 (13 praxis)
Info   :     > Iteration  46625: fmin=34.883181334 f=34.913132060 (14 subplex)
Info   :     > Iteration  47627: fmin=34.883181334 f=34.885049346 (15 praxis)
Info   :     > Iteration  47628: fmin=34.883181334 f=34.885014016 (15 praxis)
Info   :     > Final solution  : f   =34.883181334
Warning:     > Failed due to unknown error.
Info   : Regularizing tessellation...
Info   :   -
Info   : Writing results...
Info   :     [o] Writing file `nfrom_morpho-id1.tess'...
Info   :     [o] Wrote file `nfrom_morpho-id1.tess'.
Info   : Elapsed time: 40.186 secs.
========================================================================
rquey commented 2 years ago

Change sphericity:lognormal(0.145,0.03,1-x) into 1-sphericity:lognormal(0.145,0.03).