Grain Orientation output

[LiZa111111111]

Hello,
I have implemented the EBSD pixel orientation data in Neper and would need to find out the mean grain orientation from this. Is this possible with Neper?

[rquey]

Yes, it is, thanks to the new -S module.

If you have a .tesr, first export it as a simulation directory (.sim), which is how general post-processing now works:

$ neper -T -loadtesr foo.tesr -for sim

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.6.1-3
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, 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] -loadtesr foo.tesr -for sim
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   : Importing raster tessellation...
Info   :     [i] Parsing file `foo.tesr'...
Info   :     [i] Parsed file `foo.tesr'.
Info   :   - Volume has 20x20x20 (8000) points and 2 cells.
Info   :     > The cells are numbered contiguously from 1.
Info   : Writing results...
Info   :     [o] Writing directory `foo.sim'...
Info   :     [o] Wrote directory `foo.sim'.
Info   : Elapsed time: 0.037 secs.
========================================================================

Then, you can ask for the cell orientations:

$ neper -S foo.sim -rescell ori

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.6.1-3
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -S loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] foo.sim -rescell ori
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   :     > Input files: tesr
Info   :     > Cell number      : 2
Info   :     > Step number      : 0
Info   : Running post-processing...
Info   :   - Loading inputs...
Info   :     > simulation.tesr...
Info   :   - Writing results...
Info   :     [o] Writing directory `foo.sim/results/cells'...
Info   :       . ori (cubic crystal symmetry) ..................... 100%
Info   :     [o] Wrote directory `foo.sim/results/cells'.
Info   : Elapsed time: 0.018 secs.
========================================================================

This populates the foo.sim directory with the new results, as you can see as follows:

$ neper -S foo.sim

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.6.1-3
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -S loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] foo.sim/
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   :     > Input files: tesr
Info   :     > Cell number      : 2
Info   :     > Step number      : 0
Info   :     > Cell results     : ori
Info   : Elapsed time: 0.002 secs.
========================================================================

The structure of the simulation directory is as follows:

$ tree foo.sim
foo.sim
├── inputs
│   └── simulation.tesr
└── results
    └── cells
        └── ori
            └── ori.step0

Now, it you wanted to export, for instance, the cell sizes, you would run

$ neper -S foo.sim -rescell size

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.6.1-3
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -S loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] foo.sim -rescell size
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   :     > Input files: tesr
Info   :     > Cell number      : 2
Info   :     > Step number      : 0
Info   :     > Cell results     : ori
Info   : Running post-processing...
Info   :   - Loading inputs...
Info   :     > simulation.tesr...
Info   :   - Writing results...
Info   :     [o] Writing directory `foo.sim/results/cells'...
Info   :       . size ............................................. 100%
Info   :     [o] Wrote directory `foo.sim/results/cells'.
Info   : Elapsed time: 0.006 secs.
========================================================================

which, again, populates the .sim directory:

$ neper -S foo.sim

========================    N   e   p   e   r    =======================
Info   : A software package for polycrystal generation and meshing.
Info   : Version 4.6.1-3
Info   : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info   : Running on 8 threads.
Info   : <https://neper.info>
Info   : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info   : Loading initialization file `/home/rquey/.neperrc'...
Info   : ---------------------------------------------------------------
Info   : MODULE  -S loaded with arguments:
Info   : [ini file] (none)
Info   : [com line] foo.sim
Info   : ---------------------------------------------------------------
Info   : Reading input data...
Info   :   - Reading arguments...
Info   :     > Input files: tesr
Info   :     > Cell number      : 2
Info   :     > Step number      : 0
Info   :     > Cell results     : ori size
Info   : Elapsed time: 0.001 secs.
========================================================================

$ tree foo.sim
foo.sim
├── inputs
│   └── simulation.tesr
└── results
    └── cells
        ├── ori
        │   └── ori.step0
        └── size
            └── size.step0

[LiZa111111111]

Hello,
Thank you.
I have now tested this.
It works for the generation of cell sizes but not for the orientations. I got the message "Variable 'ori' is not scalar (0 value)".
In the tesr-file, the pixel orientations are defined in Euler-bunge degrees. Do you think the issue is within the definition of angles?

[LiZa111111111]

I managed to solve this issue by adding the *oridef part in my tesr-file.
Now it works fine! Thanks again.