-nogui doesn't save tetgen generated mesh

Clearly defined bug reports and their fixes

-nogui doesn't save tetgen generated mesh

Postby mmcker » 21 Sep 2017, 00:12

Hi,

I would like to be able to get an elmer mesh using tetgen from commandline.

ElmerGUI -i cube.poly -o elmer -nogui -e

In GUI mode, it works fine - I can mesh easily.

However, below is a transcript from the command line attempt.

Note the line "Unable to open elmer/mesh.header"

(the subdirectory elmer has been created).

I haven't fully gone through code, but I suspect that the mesh file handle hasn't been opened in this mode before the attempt is made to write to it.

Regards,
Mark



Load /usr/local/bin/../share/ElmerGUI/edf/egini.xml... done
Load tetgen plugin... done
Constructing ElmergridAPI... done
libGL error: No matching fbConfigs or visuals found
libGL error: failed to load driver: swrast
Load /usr/local/bin/../share/ElmerGUI/edf/edf.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/heatequation.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/resultoutput.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/electrostatics.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/helmholtz.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/navier-stokes.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/meshdeform.xml... done
Load /usr/local/bin/../share/ElmerGUI/edf/linearelasticity.xml... done
Reading input file cube.poly
Selected tetlib for smesh/poly-format
Opening /home/elm/tg_tst/cube.poly.
Sending start request to mesh generator...
Starting meshing thread with low priority
tetlib: control string: nnJApq1.414V
tetlib: input points: 16
Constructing Delaunay tetrahedralization.
Creating initial tetrahedralization.
Incrementally inserting points.
Mesh generator started
41 Flips (T23 14, T32 4, T22 5, T44 18)
Delaunay seconds: 0.004014
Creating surface mesh.
Constructing mapping from indices to points.
Constructing mapping from points to tetrahedra.
Unifying segments.
Constructing mapping from points to subfaces.
Merging coplanar facets.
Marking acute vertices.
Constructing mapping from points to segments.
Perturbing vertices.
Removing degenerate subfaces.
16 break points.
Delaunizing segments.
Constructing mapping from points to tetrahedra.
Queuing missing segments.
0 protect points.
R1: 0, R2: 0, R3: 0.
Constraining facets.
Constructing mapping from points to tetrahedra.
The biggest cavity: 0 faces, 0 vertices
Enlarged 0 times
Segment and facet seconds: 0.006873
Removing unwanted tetrahedra.
Marking holes for elimination.
Marking concavities for elimination.
Marking neighbors of marked tetrahedra.
Deleting marked tetrahedra.
Hole seconds: 2.9e-05
Repairing mesh.
Repair seconds: 3.2e-05
Adding Steiner points to enforce quality.
Marking sharp segments.
Deciding feature-point sizes.
Constructing mapping from points to segments.
0 feature points.
Splitting encroached subsegments.
1 split points.
Splitting encroached subfaces.
72 split points.
Splitting bad tetrahedra.
379 refinement points.
Totally added 452 points.
Quality seconds: 0.019868
Optimizing mesh.
12 edges are flipped.
Optimize seconds: 0.0013

Writing nodes.
Writing elements.
Writing faces.
Writing neighbors.

Output seconds: 0.000121
Total running seconds: 0.032299

Statistics:

Input points: 16
Input facets: 12
Input segments: 24
Input holes: 1
Input regions: 0

Mesh points: 484
Mesh tetrahedra: 2278
Mesh triangles: 4780
Mesh subfaces: 448
Mesh subsegments: 70

Mesh quality statistics:

Smallest volume: 7.0458e-05 | Largest volume: 25850
Shortest edge: 0.099926 | Longest edge: 78.991
Smallest facangle: 20.914 | Largest facangle: 135.9782
Smallest dihedral: 7.4639 | Largest dihedral: 164.9027

Aspect ratio histogram:
< 1.5 : 43 | 6 - 10 : 42
1.5 - 2 : 455 | 10 - 15 : 1
2 - 2.5 : 734 | 15 - 25 : 0
2.5 - 3 : 460 | 25 - 50 : 0
3 - 4 : 419 | 50 - 100 : 0
4 - 6 : 124 | 100 - : 0
(A tetrahedron's aspect ratio is its longest edge length divided by its
smallest side height)

Face angle histogram:
0 - 10 degrees: 0 | 90 - 100 degrees: 682
10 - 20 degrees: 0 | 100 - 110 degrees: 336
20 - 30 degrees: 523 | 110 - 120 degrees: 141
30 - 40 degrees: 1618 | 120 - 130 degrees: 31
40 - 50 degrees: 1746 | 130 - 140 degrees: 7
50 - 60 degrees: 893 | 140 - 150 degrees: 0
60 - 70 degrees: 876 | 150 - 160 degrees: 0
70 - 80 degrees: 1557 | 160 - 170 degrees: 0
80 - 90 degrees: 1150 | 170 - 180 degrees: 0
Minimum input face angle is 90 (degree).

Dihedral angle histogram:
0 - 5 degrees: 0 | 80 - 110 degrees: 1399
5 - 10 degrees: 3 | 110 - 120 degrees: 345
10 - 20 degrees: 90 | 120 - 130 degrees: 233
20 - 30 degrees: 291 | 130 - 140 degrees: 132
30 - 40 degrees: 660 | 140 - 150 degrees: 79
40 - 50 degrees: 741 | 150 - 160 degrees: 47
50 - 60 degrees: 425 | 160 - 170 degrees: 10
60 - 70 degrees: 68 | 170 - 175 degrees: 0
70 - 80 degrees: 33 | 175 - 180 degrees: 0
Minimum input facet dihedral angle is 90 (degree).

Memory allocation statistics:

Maximum number of vertices: 484
Maximum number of tetrahedra: 2305
Maximum number of subfaces: 452
Maximum number of segments: 70
Approximate heap memory used by the mesh (K bytes): 371.391

tetlib: nodes: 484
tetlib: elements: 2278
tetlib: boundary elements: 448
Mesh generation ready
Old mesh cleared
Found 672 edges on boundary
Bcs/materials on surface elements: 3
Bcs/materials on edge elements: 0
Bcs/materials on point elements: 0
Generating lists to display
Initialize GL
Vendor: Mesa project: www.mesa3d.org
Renderer: Mesa GLX Indirect
Version: 1.2 (1.5 Mesa 6.4.1)
Generated 8 lists
Input file processed
Apply 0 operations
Saving mesh files
Saving elmer mesh files
Unable to open elmer/mesh.header
Ready
Exiting
mmcker
 
Posts: 7
Joined: 21 Sep 2017, 00:00

Re: -nogui doesn't save tetgen generated mesh

Postby mmcker » 21 Sep 2017, 09:55

Please ignore.

Putting full path in -o filename fixes the problem.

-o /home/elm/tg_tst/elmer

I'm a little surprised. Not even ./elmer worked
mmcker
 
Posts: 7
Joined: 21 Sep 2017, 00:00


Return to Bug reports

Who is online

Users browsing this forum: No registered users and 1 guest