loadmesh option is not working for 2D polycrystals #521

shahul31 · 2022-07-20T15:04:36Z

shahul31
Jul 20, 2022

Hello,

Thank you for the update about loadmesh options, which is really helpful. It is working good in 3D, but this option is not working with 2D case! It fails like this below, it is the same for '-statelset' also in 2D.

shahul@shahul-VirtualBox:~/Neper/FEPX_results/8_Test$ neper -M -loadmesh test50_2D.msh -format inp

======================== N e p e r =======================
Info : A software package for polycrystal generation and meshing.
Info : Version 4.4.2-33
Info : Built with: gsl|muparser|opengjk|openmp|nlopt|libscotch (full)
Info : Running on 4 threads.
Info : https://neper.info
Info : Copyright (C) 2003-2022, and GNU GPL'd, by Romain Quey.
Info : ---------------------------------------------------------------
Info : MODULE -M loaded with arguments:
Info : [ini file] (none)
Info : [com line] -loadmesh test50_2D.msh -format inp
Info : ---------------------------------------------------------------
Info : Reading input data...
Info : - Reading arguments...
Info : Loading input data...
Info : Loading result mesh...
Info : [i] Parsing file test50_2D.msh'... Info : [i] Parsed file test50_2D.msh'.
Info : - Reconstructing mesh...
Error : You have discovered a bug in Neper! Please file an issue at
https://github.com/rquey/neper/issues. Thank you.

Aborted (core dumped)

Answered by rquey

Jul 20, 2022

There is indeed a bug, but it always helps to load the tessellation alongside a mesh, so that Neper does not have to reconstruct the topology:

$ neper -T -n 10 -dim 2

========================    N   e   p   e   r    =======================                                 
Info   : A software package for polycrystal generation and meshing.                                      
Info   : Version 4.4.2-33                           
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 …

View full answer

rquey · 2022-07-20T15:32:57Z

rquey
Jul 20, 2022
Maintainer

There is indeed a bug, but it always helps to load the tessellation alongside a mesh, so that Neper does not have to reconstruct the topology:

$ neper -T -n 10 -dim 2

========================    N   e   p   e   r    =======================                                 
Info   : A software package for polycrystal generation and meshing.                                      
Info   : Version 4.4.2-33                           
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] -n 10 -dim 2                    
Info   : ---------------------------------------------------------------                                 
Info   : Reading input data...                      
Info   : Creating domain...                         
Info   : Creating tessellation...                   
Info   :   - Setting seeds...                       
Info   :   - Generating crystal orientations...                                                          
Info   :   - Running tessellation...                
Info   : Writing results...                         
Info   :     [o] Writing file `n10-id1.tess'...                                                          
Info   :     [o] Wrote file `n10-id1.tess'.                                                              
Info   : Elapsed time: 0.006 secs.                  
========================================================================                                 

$ neper -M n10-id1.tess                                                                       

========================    N   e   p   e   r    =======================                                 
Info   : A software package for polycrystal generation and meshing.                                      
Info   : Version 4.4.2-33                           
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  -M loaded with arguments:                                                               
Info   : [ini file] -tmp /home2/tmp                 
Info   : [com line] n10-id1.tess                    
Info   : ---------------------------------------------------------------                                 
Info   : Reading input data...                      
Info   :   - Reading arguments...                   
Info   : Loading input data...                      
Info   :   - Loading tessellation...                
Info   :     [i] Parsing file `n10-id1.tess'...                                                          
Info   :     [i] Parsed file `n10-id1.tess'.                                                             
Info   : Meshing...       
Info   :   - Preparing... (cl = 0.07906)                                                                 
Info   :   - 0d meshing... 100%                     
Info   :   - 1d meshing... 100%                     
Info   :   - 2d meshing... 100% (0.87|0.9/60%|20%|20%)                                                   
Info   : Searching nsets and fasets...              
Info   : Writing mesh results...                    
Info   :   - Preparing mesh...                      
Info   :   - Mesh properties:                       
Info   :     > Node number:      297                
Info   :     > Elt  number:      535                
Info   :     > Mesh area:    1.000                  
Info   :   - Writing mesh...                        
Info   :     [o] Writing file `n10-id1.msh'...                                                           
Info   :     [o] Wrote file `n10-id1.msh'. 
          Info   : Elapsed time: 0.123 secs.                  
========================================================================                                 

$ neper -M n10-id1.tess -loadmesh n10-id1.msh -format inp                                     

========================    N   e   p   e   r    =======================                                 
Info   : A software package for polycrystal generation and meshing.                                      
Info   : Version 4.4.2-33                           
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  -M loaded with arguments:                                                               
Info   : [ini file] -tmp /home2/tmp                 
Info   : [com line] n10-id1.tess -loadmesh n10-id1.msh -format inp                                       
Info   : ---------------------------------------------------------------                                 
Info   : Reading input data...                      
Info   :   - Reading arguments...                   
Info   : Loading input data...                      
Info   :   - Loading tessellation...                
Info   :     [i] Parsing file `n10-id1.tess'...                                                          
Info   :     [i] Parsed file `n10-id1.tess'.                                                             
Info   : Loading result mesh...                     
Info   :     [i] Parsing file `n10-id1.msh'...                                                           
Info   :     [i] Parsed file `n10-id1.msh'.                                                              
Info   :   - Reconstructing mesh...                 
Info   : Searching nsets and fasets...              
Info   : Writing mesh results...                    
Info   :   - Preparing mesh...                      
Info   :   - Mesh properties:                       
Info   :     > Node number:      297                
Info   :     > Elt  number:      535                
Info   :     > Mesh area:    1.000                  
Info   :   - Writing mesh...                        
Info   :     [o] Writing file `n10-id1.inp'...                                                           
Info   :     [o] Wrote file `n10-id1.inp'.                                                               
Info   : Elapsed time: 0.013 secs.                  
========================================================================

1 reply

shahul31 Jul 20, 2022
Author

Great, thank you very much :) :)

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Uh oh!

loadmesh option is not working for 2D polycrystals #521

Uh oh!

{{title}}

Uh oh!

Replies: 1 comment 1 reply

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{editor}}'s edit

{{editor}}'s edit

Uh oh!

Uh oh!

{{title}}

Uh oh!

Select a reply

Uh oh!

loadmesh option is not working for 2D polycrystals #521

Uh oh!

shahul31 Jul 20, 2022

Replies: 1 comment · 1 reply

Uh oh!

Uh oh!

rquey Jul 20, 2022 Maintainer

Uh oh!

shahul31 Jul 20, 2022 Author

shahul31
Jul 20, 2022

Replies: 1 comment 1 reply

rquey
Jul 20, 2022
Maintainer

shahul31 Jul 20, 2022
Author