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

FreeCAD · Dec 17, 2019 · e7be69e · e7be69e
1 parent 57bfc2e
commit e7be69e
Showing 1 changed file with 11 additions and 10 deletions.
diff --git a/src/Mod/Fem/femmesh/gmshtools.py b/src/Mod/Fem/femmesh/gmshtools.py
@@ -723,6 +723,17 @@ def write_geo(self):
         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 @@ def write_geo(self):
             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