Reached unknown criterion

[nippyawn]

Hello,

I want to improve the accuracy of tessellations.
I have imposed strict epsilon (eps) conditions.
However, when I reach an unknown criterion, the generation stops with the default eps conditions.
On the other hand, when I test an easy model, it reaches the eps criterion and the generation is completed successfully.

Please advise me on why the "Reached unknown criterion" occurs.

neper -T -n from_morpho -domain "cube(269,127,371)" -morpho "tesr:file(test_0527-c.tesr)" -transform "grow,autocrop,resetorigin,renumber" -morphooptialgo praxis -morphooptistop "eps<1e-10" -morphooptialgomaxiter "max(varnb,1000000000)" -morphooptiobj "tesr:pts(region=all,res=10)+val(bounddist)" -o test_0527_model2

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.8.2
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 16 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : No initialization file found (`/home/ishikawa/.neperrc').
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n from_morpho -domain cube(269,127,371) -morpho
         tesr:file(test_0527-c.tesr) -transform
         grow,autocrop,resetorigin,renumber -morphooptialgo praxis
         -morphooptistop eps<1e-10 -morphooptialgomaxiter
         max(varnb,1000000000) -morphooptiobj
         tesr:pts(region=all,res=10)+val(bounddist) -o test_0527_model2
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :     [i] Parsing file `test_0527-c.tesr'...
Info   :     [i] Parsed file `test_0527-c.tesr'.
Info   :       . Number of cells: 4804
Info   :   - Setting seeds... 100%
100%
Info   :       . Listing cell voxels... 100.0%
Info   :       . Sampling factor along x = 1.714389.
Info   :       . Sampling factor along y = 1.714389.
Info   :       . Sampling factor along z = 1.714389.
Info   :       . Number of voxels per cell: min = 1, avg = 551, max = 7045.
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =15.544134879
Info   :     > Iter     1241345: fmin==15.393094072 f==15.401476854 ((25 min) 26 void cells)
Info   :     > Final solution  : f   =15.393094072 (1241345 iterations)
Info   :     > Reached unknown criterion.
Info   : Generating crystal orientations...
Info   : Transforming tessellation...
Warning:     > Skipping `grow'...
Warning:     > Skipping `autocrop'...
Warning:     > Skipping `resetorigin'...
Warning:     > Skipping `renumber'...
Info   : Writing results...
Info   :     [o] Writing file `test_0527_model2.tess'...
Info   :     [o] Wrote file `test_0527_model2.tess'.
Info   : Elapsed time: 13201.712 secs.
========================================================================

neper -T -n from_morpho -morpho "tesr:file(n50-id1.tesr)" -transform "grow,autocrop,resetorigin,renumber" -morphooptialgo praxis -morphooptistop "eps<1e-3" -morphooptialgomaxiter "max(varnb,100000)" -morphooptiobj "tesr:pts(region=all,res=10)+val(bounddist)" -o test

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.8.2
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 16 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : No initialization file found (`/home/ishikawa/.neperrc').
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n from_morpho -morpho tesr:file(n50-id1.tesr)
         -transform grow,autocrop,resetorigin,renumber -morphooptialgo
         praxis -morphooptistop eps<1e-3 -morphooptialgomaxiter
         max(varnb,100000) -morphooptiobj
         tesr:pts(region=all,res=10)+val(bounddist) -o test
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :     [i] Parsing file `n50-id1.tesr'...
Info   :     [i] Parsed file `n50-id1.tesr'.
Info   :       . Number of cells: 50
Info   :   - Setting seeds... 100%
100%
Info   :       . Listing cell voxels... 100.0%
Warning:       . Sampling factor along x = 0.673556.  Setting to 1...
Warning:       . Sampling factor along y = 0.673556.  Setting to 1...
Warning:       . Sampling factor along z = 0.673556.  Setting to 1...
Info   :       . Number of voxels per cell: min = 70, avg = 160, max = 314.
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =0.004279264
Info   :     > Iteration    349: fmin=0.002549979 f=0.002549979
Info   :     > Final solution  : f   =0.002549979 (349 iterations)
Info   :     > Reached `eps' criterion.
Info   : Generating crystal orientations...
Info   : Transforming tessellation...
Warning:     > Skipping `grow'...
Warning:     > Skipping `autocrop'...
Warning:     > Skipping `resetorigin'...
Warning:     > Skipping `renumber'...
Info   : Writing results...
Info   :     [o] Writing file `test.tess'...
Info   :     [o] Wrote file `test.tess'.
Info   : Elapsed time: 0.556 secs.
========================================================================

Thank you

[rquey]

The unknown criterion generally appears when the problem is not properly formulated and optimization fails to converge. This is confirmed by the large f and fmin values in your first case.