FEM: gmsh, new interpolation parameter, add missing code (thanks to _UR)

This commit is contained in:
Bernd Hahnebach
2019-12-17 18:12:08 +01:00
parent b09233314b
commit c394eb91ff

View File

@@ -723,6 +723,17 @@ class GmshTools():
geo.write("\n")
geo.write("// mesh order\n")
geo.write("Mesh.ElementOrder = " + self.order + ";\n")
if self.order == "2":
if hasattr(self.mesh_obj, "SecondOrderLinear") and self.mesh_obj.SecondOrderLinear is True:
geo.write(
"Mesh.SecondOrderLinear = 1; // Second order nodes are created "
"by linear interpolation instead by curvilinear\n"
)
else:
geo.write(
"Mesh.SecondOrderLinear = 0; // Second order nodes are created "
"by linear interpolation instead by curvilinear\n"
)
geo.write("\n")
geo.write(
@@ -757,16 +768,6 @@ class GmshTools():
geo.write("Coherence Mesh; // Remove duplicate vertices\n")
else:
geo.write("Mesh " + self.dimension + ";\n")
if hasattr(self.mesh_obj, "SecondOrderLinear") and self.mesh_obj.SecondOrderLinear is True:
geo.write(
"Mesh.SecondOrderLinear = 1; // Second order nodes are created "
"by linear interpolation instead by curvilinear\n"
)
else:
geo.write(
"Mesh.SecondOrderLinear = 0; // Second order nodes are created "
"by linear interpolation instead by curvilinear\n"
)
geo.write("\n")
# save mesh