I am building a microstructure where some of the grains have thin lamella of a different phase. To do so, I create a basic polycrystal without the lamella to determine with a volumetric calculation which grains will have the lamella. Then, with this information, I create two files with the lamella width relationship and the normal vectors. Finally, I run the tessellation and the mesh commands. When I was running the tessellation, the following message appeared: "Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you".
neper -T -n $number_grains::from_morpho -domain "cube(0.15,0.15,0.15)" -id 1::1 -morpho "diameq:lognormal(1,0.35),1-sphericity:lognormal(0.145,0.03)::lamellar(w=file(lamwidth.txt),v=file(lamnormal.txt),pos=optimal)" -reg 1 -sel 0.08 -oridescriptor e -group "w<0.006?1:2" -statcell "scaleid(1)" -format tess,ori -o Polycrystal_laminated
======================== N e p e r =======================
Info : A software package for polycrystal generation and meshing.
Info : Version 4.4.2-16
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 : Ignoring initialization file.
Info : ---------------------------------------------------------------
Info : MODULE -T loaded with arguments:
Info : [ini file] (none)
Info : [com line] -n 300::from_morpho -domain cube(0.15,0.15,0.15) -id
1::1 -morpho
diameq:lognormal(1,0.35),1-sphericity:lognormal(0.145,0.03)::lamellar(w=file(lamwidth.txt),v=file(lamnormal.txt),pos=optimal)
-reg 1 -sel 0.08 -oridescriptor e -group w<0.006?1:2 -statcell
scaleid(1) -format tess,ori -o Polycrystal_laminated
Info : ---------------------------------------------------------------
Info : Reading input data...
Info : Creating domain...
Info : Creating tessellations...
Info : [level] Running in multiscale mode (2 scales)...
Info : [ 1] 1/1 Info : - Setting seeds... 100%
Info : - Generating crystal orientations...
Info : - Running tessellation...
Info : > Initial solution: f =2.087236984
Info : > Iteration 94026: fmin=0.000196552 f=0.000196552
Info : > Final solution : f =0.000196552 (94026 iterations)
Info : > Reached `eps' criterion.
Info : [ 2] 1/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 2/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 3/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 4/300 Info : - Running tessellation...
...
Info : - Generating crystal orientations...
Info : [ 2] 296/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 297/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 298/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 299/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : [ 2] 300/300 Info : - Running tessellation...
Info : - Generating crystal orientations...
Info : Flattening 2-scale tessellation...
Info : [level] Flattening 2-scale tessellation...
Info : [ 1] 1/1: id=1
Error : You have discovered a bug in Neper! Please file an issue at
https://github.com/rquey/neper/issues. Thank you.
Aborted
The first 10 lines of the lamwidth.txt and lamnormal.txt files are the following (I also attached the files):
Hi,
I am building a microstructure where some of the grains have thin lamella of a different phase. To do so, I create a basic polycrystal without the lamella to determine with a volumetric calculation which grains will have the lamella. Then, with this information, I create two files with the lamella width relationship and the normal vectors. Finally, I run the tessellation and the mesh commands. When I was running the tessellation, the following message appeared: "Error : You have discovered a bug in Neper! Please file an issue at https://github.com/rquey/neper/issues. Thank you".
The first 10 lines of the lamwidth.txt and lamnormal.txt files are the following (I also attached the files):
lamwidth.txt 1 10.0:9.0 2 10.0:9.0 3 0.0225:0.002625 4 0.0225:0.002625 5 0.0225:0.002625 6 0.0225:0.002625 7 10.0:9.0 8 10.0:9.0 9 10.0:9.0 10 0.0225:0.002625
lamnormal.txt 1 1 0 0 2 1 0 0 3 0.98 -0.15 -0.17 4 -0.83 0.07 0.59 5 -0.94 0.28 -0.24 6 -0.13 0.35 0.24 7 1 0 0 8 1 0 0 9 1 0 0 10 -0.76 -0.94 -0.38
lamnormal.txt lamwidth.txt
Thank you for your continuous support.