The error is this. If I create a 3-D mesh from scratch, with some uniform refinements, it has a certain number of faces, facets, and edges. I then save the mesh to a .vol file. When I reload the mesh, the number of faces, facets, and edges is *different*. Why?
I have created a sample code below.
import numpy as np
import netgen.gui
from ngsolve import *
from netgen.csg import unit_cube
from netgen.csg import CSGeometry
from netgen.csg import Pnt
from netgen.csg import OrthoBrick
#load_mesh = True
load_mesh = False
mesh_FN = "test_mesh.vol"
if load_mesh:
print("Load mesh from file.")
mesh = Mesh(mesh_FN)
else:
print("Create mesh from scratch.")
# make the mesh
cube = OrthoBrick( Pnt(0,0,0), Pnt(1,1,1) ).bc("Gamma")
geo = CSGeometry()
geo.Add(cube)
mesh = Mesh(geo.GenerateMesh(maxh=1.0))
for r in range(3):
mesh.Refine()
# save the mesh
mesh.ngmesh.Save(mesh_FN)
print("Some info about the mesh:")
print("Number of Elements: " + str(mesh.ne))
print("Number of Faces: " + str(mesh.nface))
print("Number of Facets: " + str(mesh.nfacet))
print("Number of Edges: " + str(mesh.nedge))
print("Number of Vertices: " + str(mesh.nv))
When you run it the first time, you should get:
>From scratch:
>Number of Elements: 3072
>Number of Faces: 7530
>Number of Facets: 7530
>Number of Edges: 4905
>Number of Vertices: 729
Next, change "load_mesh" to True and rerun. The output will be:
>When loading from file:
>Number of Elements: 3072
>Number of Faces: 6528
>Number of Facets: 6528
>Number of Edges: 4184
>Number of Vertices: 729
Why is it different?
