Meshing Columnar Bug using -clratio

[RichardPA]

Hi Romain,

I was testing out the -clratio option to build a denser mesh on the side of a columnar model and ran into the bug below. I think this may be due to using -clratio and other options because I haven't run into this bug before. I have no questions, I just wanted to report this. I'll attach the columnar images in case they're of use.

Cheers,

Richard

[rquey]

Thanks @RichardPA for the bug report. Could you copy/paste your -T -and -M commands? Thanks.

[RichardPA]

Yes, I will paste them below.

$ neper -T -n 50 -id 1 -domain "cube(5,5,5)" -morpho "diameq:lognormal(.5,0.005),1-sphericity:lognormal(.13,0.04),columnar(z)" -oricrysym hexagonal -reg 1 -mloop 8 -o simulation
$ neper -M simulation.tess -order 2 -part 1 -meshqualmin 0.9 -rcl 1.5 -rcl 1 -clratio 1:1:2 -mesh2dmaxtime 1000000000 -mesh2diter 2147483647 -mesh3dmaxtime 1000000000 -mesh3diter 2147483647 -format msh,ori,bcs,vtk -o simulation

Cheers,

Richard

[rquey]

I cannot reproduce this bug:

$ neper -T -n 50 -id 1 -domain "cube(5,5,5)" -morpho "diameq:lognormal(.5,0.005),1-sphericity:lognormal(.13,0.04),columnar(z)" -oricrysym hexagonal -reg 1 -mloop 8 -o simulation

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.1.2                              
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.           
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2020, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n 50 -id 1 -domain cube(5,5,5) -morpho          
         diameq:lognormal(.5,0.005),1-sphericity:lognormal(.13,0.04),columnar(z)
         -oricrysym hexagonal -reg 1 -mloop 8 -o simulation          
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   =97.757957106
Info   :     > Iteration  97833: fmin=56.062852164 f=56.062852165                                        
Info   :     > Reached `eps' criterion.      
Info   : Regularizing tessellation...            
Info   :   - loop 8/8: 100% del=10           
Info   : Writing results...                                                                              
Info   :     [o] Writing file `simulation.tess'...
Info   :     [o] Wrote file `simulation.tess'.                                                           
Info   : Elapsed time: 85.460 secs.          
========================================================================     

$ neper -M simulation.tess -order 2 -part 1 -meshqualmin 0.9 -rcl 1.5 -rcl 1 -clratio 1:1:2 -mesh2dmaxtime 1000000000 -mesh2diter 2147483647 -mesh3dmaxtime 1000000000 -mesh3diter 2147483647 -format msh,ori,bcs,vtk -o simulation                                                       

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.1.2
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2020, and GNU GPL'd, by Romain Quey.
Info   : Ignoring initialization file.
Info   : ---------------------------------------------------------------
Info   : MODULE  -M loaded with arguments:
Info   : [ini file] (none) 
Info   : [com line] simulation.tess -order 2 -part 1 -meshqualmin 0.9
         -rcl 1.5 -rcl 1 -clratio 1:1:2 -mesh2dmaxtime 1000000000
         -mesh2diter 2147483647 -mesh3dmaxtime 1000000000 -mesh3diter
         2147483647 -format msh,ori,bcs,vtk -o simulation
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   : Loading input data...
Info   :   - Loading tessellation...
Info   :     [i] Parsing file `simulation.tess'...
Info   :     [i] Parsed file `simulation.tess'.
Info   : Meshing...
Info   :   - Preparing... (cl = 0.6786) 100%
Info   :   - 0D meshing... 100%
Info   :   - 1D meshing... 100%
Info   :   - 2D meshing... 100% (0.61|0.9/86%|12%| 2%)  
Info   :   - Fixing 2D-mesh pinches...
Info   :   - 3D meshing... 100% (0.78|0.91/94%| 6%| 0%) 
Info   :   - Switching mesh to order 2...
Info   : Partitioning mesh...
Info   :   - Initializing...
Info   :   - Preparing data...
Info   :   - Partitioning nodes...
Info   :   - Partitioning elements...
Info   :   - Partitioning properties:
Info   :     > Nb of nodes per partition: min = 8240, max = 8240
Info   :     > Nb of elts  per partition: min = 5276, max = 5276
Info   :     > Nb of communications: 0
Info   : Searching nsets and fasets...
Info   : Writing mesh results...
Info   :   - Preparing mesh...
Info   :   - Mesh properties:
Info   :     > Node number:     8240
Info   :     > Elt  number:     5276
Info   :     > Mesh volume:  125.000
Info   :   - Writing mesh...
Info   :     [o] Writing file `simulation.msh'...
Info   :     [o] Wrote file `simulation.msh'.
Info   :     [o] Writing file `simulation.ori'...
Info   :     [o] Wrote file `simulation.ori'.
Info   :     [o] Writing file `simulation.bcs'...
Info   :     [o] Wrote file `simulation.bcs'.
Info   :     [o] Writing file `simulation.vtk'...
Info   :     [o] Wrote file `simulation.vtk'.
Info   : Elapsed time: 5.738 secs.
========================================================================