Custom Diameq Distribution Function

[whishei]

Hi,

I am having an issue trying to generate a large grain tessellation when I give an input file for the diameq distribution. The tessellation works for 100 and 125 grain, but fails for anything 150 grains and up. Do you have a sample input file format for the distribution function? Here is an image of the probability density function and I have attached the Prob.txt file for the distribution function.

$ neper -T -n 150 -morpho "diameq:file(Prob.txt)" -statcell diameq -o Test150

========================    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
Info   : Running on 6 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2020, and GNU GPL'd, by Romain Quey.
Info   : No initialization file found (`/home/swhitman/.neperrc').
Info   : ---------------------------------------------------------------
Info   : MODULE  -T loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] -n 150 -morpho diameq:file(Prob.txt) -statcell diameq
         -o Test150
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Creating domain...
Info   : Creating tessellation...
Info   :     [i] Parsing file `Prob.txt'...
Info   :     [i] Parsed file `Prob.txt'.
Aborted (core dumped)

Prob.txt

[rquey]

The current command is

$ neper -T -n 150 -morpho "diameq:custom(Prob.txt)" -statcell diameq -o Test150

With your command, Neper was trying to read individual diameq values in Prob.txt, which is why it failed for -n larger than the number of entries in the file.

2 advices:

• specify a sphericity distribution (e.g. the default) to avoid getting weird grain shapes.
• smooth your distribution (although Neper takes it just fine).

[whishei]

Thank you!