3-scale tessellation w/ lamellar not working?

[xshang93]

I use the following code to perform a test 3-scale tessellation, with the third scale being lamellar. However, The code always just stuck at the beginning of the first cell in the third scale, while consuming 99% of my CPU. When I use the same configuration to do 2-scale with lamellar, or 3-scale without lamellar, it works fine.

Command I used:

neper -T -n 3::2::from_morpho -id 1::1 -morpho "gg::gg::lamellar(w=file(lamwidth),v=(1,1,0))" -ori 'random' -format tess -o Polycrystal_laminated_3scale

lamwidth file:

1 	 10.0
2 	 0.4:0.1
3 	 10.0
4 	 0.4:0.1
5 	 10.0
6 	 10.0

What I get is:

(base) Xiaos-Mac:Geometry_test Xiao$ ./*NRE.sh

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.2.1-7
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   : No initialization file found (/Users/Xiao/.neperrc').
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n 3::2::from_morpho -id 1::1 -morpho
         gg::gg::lamellar(w=file(lamwidth),v=(1,1,0)) -ori random -format
         tess -o Polycrystal_laminated_3scale
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellations...
Info   :   [level] Running in multiscale mode (3 scales)...
Info   :   [    1]  1/1 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =10.231945625
Info   :     > Iteration    729: fmin=1.182046613 f=1.182046995  
Info   :     > Final solution  : f   =1.182046613 (729 iterations)
Info   :     > Reached eps' criterion.
Info   :   [    2]  1/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =51.000716122
Info   :     > Iteration    164: fmin=0.339765872 f=0.339766865  
Info   :     > Final solution  : f   =0.339765872 (164 iterations)
Info   :     > Reached `eps' criterion.
Info   :   [    2]  2/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =1.538774059
Info   :     > Iteration    341: fmin=1.286550132 f=1.286550134
Info   :     > Final solution  : f   =1.286550132 (341 iterations)
Info   :     > Reached eps' criterion.
Info   :   [    2]  3/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =53.849127976
Info   :     > Iteration    268: fmin=0.339364417 f=0.339364417  
Info   :     > Final solution  : f   =0.339364417 (268 iterations)
Info   :     > Reached eps' criterion.
Info   :   [    3]  1/6 Info   :   - Running tessellation...            

Any thoughts?

[rquey]

When you are at scale 3, the file refers to cells at scale 2 and so the multiscale cell ids (that come first on the line) should be of the form scale1_id::scale2_id:

lamwidth:

1::1 	 10.0
1::2 	 0.4:0.1
2::1 	 10.0
2::2 	 0.4:0.1
3::1 	 10.0
3::2 	 10.0

neper -T -n 3::2::from_morpho -id 1::1 -morpho "gg::gg::lamellar(w=file(lamwidt
h),v=(1,1,0))" -ori 'random' -format tess -o Polycrystal_laminated_3scale
                                                                                                         
========================    N   e   p   e   r    ======================= 
Info   : A software package for polycrystal generation and meshing.      
Info   : Version 4.2.1-15                           
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 `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] -tesrformat ascii
Info   : [com line] -n 3::2::from_morpho -id 1::1 -morpho
         gg::gg::lamellar(w=file(lamwidth),v=(1,1,0)) -ori random -format
         tess -o Polycrystal_laminated_3scale
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellations...
Info   :   [level] Running in multiscale mode (3 scales)...
Info   :   [    1]  1/1 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =10.231945625
Info   :     > Iteration    522: fmin=1.182061699 f=1.182061930  
Info   :     > Final solution  : f   =1.182061699 (522 iterations)
Info   :     > Reached `eps' criterion.
Info   :   [    2]  1/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =53.257306190
Info   :     > Iteration    246: fmin=0.338960255 f=0.338960263  
Info   :     > Final solution  : f   =0.338960255 (246 iterations)
Info   :     > Reached `eps' criterion.
Info   :   [    2]  2/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =1.465868429
Info   :     > Iteration    292: fmin=1.152645669 f=1.152645761
Info   :     > Final solution  : f   =1.152645669 (292 iterations)
Info   :     > Reached `eps' criterion.
Info   :   [    2]  3/3 Info   :   - Setting seeds... 100%              
Info   :   - Generating crystal orientations...
Info   :   - Running tessellation...
Info   :     > Initial solution: f   =53.540422811
Info   :     > Iteration    252: fmin=0.339120432 f=0.339120454  
Info   :     > Final solution  : f   =0.339120432 (252 iterations)
Info   :     > Reached `eps' criterion.
Info   :   [    3]  1/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   [    3]  2/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   [    3]  3/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   [    3]  4/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   [    3]  5/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   [    3]  6/6 Info   :   - Running tessellation...            
Info   :   - Generating crystal orientations...
Info   :   - Generating crystal orientations...
Info   : Flattening 3-scale tessellation...
Info   :   [level] Flattening 3-scale tessellation...
Info   :   [    2]  3/3: id=4                                            
Info   :   [    1]  1/1: id=1                                            
Info   : Writing results...
Info   :     [o] Writing file `Polycrystal_laminated_3scale.tess'...
Info   :     [o] Wrote file `Polycrystal_laminated_3scale.tess'.
Info   : Elapsed time: 1.434 secs.
========================================================================

The infinite loop issue was fixed in version 4.2.1-15.

[xshang93]

This is amazing. Thank you for the reply. So essentially this is similar to what has been answered in #325 , where the -ori (or -morpho in this case) files do need to be (or at least indicate that it is) a multiscale file.

Thank you very much for your help.
XS

Looking forward to the new neper version btw!

[rquey]

Versions x.y.z-patch_number are the continuous development versions available from the github repo:
https://github.com/rquey/neper
https://github.com/rquey/neper/commits/main