From c8a0c2c8153ec4f0ca00f07b4e41066129e875ad Mon Sep 17 00:00:00 2001
From: Bernd Hahnebach <bernd@bimstatik.org>
Date: Thu, 19 Sep 2019 09:06:42 +0200
Subject: [PATCH] FEM: gmsh tools, make it easier to use an own working dir

---
 src/Mod/Fem/femmesh/gmshtools.py | 51 ++++++++++++++++++++++++++++----
 1 file changed, 45 insertions(+), 6 deletions(-)

diff --git a/src/Mod/Fem/femmesh/gmshtools.py b/src/Mod/Fem/femmesh/gmshtools.py
index 9a5f08083d..70e039e765 100644
--- a/src/Mod/Fem/femmesh/gmshtools.py
+++ b/src/Mod/Fem/femmesh/gmshtools.py
@@ -139,16 +139,22 @@ class GmshTools():
         self.gmsh_bin = ""
         self.error = False
 
-    def create_mesh(self):
+    def update_mesh_data(self):
         self.start_logs()
         self.get_dimension()
-        self.get_tmp_file_paths()
-        self.get_gmsh_command()
         self.get_group_data()
         self.get_region_data()
         self.get_boundary_layer_data()
+
+    def write_gmsh_input_files(self):
         self.write_part_file()
         self.write_geo()
+
+    def create_mesh(self):
+        self.update_mesh_data()
+        self.get_tmp_file_paths()
+        self.get_gmsh_command()
+        self.write_gmsh_input_files()
         error = self.run_gmsh_with_geo()
         self.read_and_set_new_mesh()
         return error
@@ -822,17 +828,16 @@ class GmshTools():
 
 """
 # simple example how to use the class GmshTools
-
 import Part, ObjectsFem
-
 doc = App.ActiveDocument
+
 box_obj = doc.addObject("Part::Box", "Box")
 doc.recompute()
+box_obj.ViewObject.Visibility = False
 
 femmesh_obj = ObjectsFem.makeMeshGmsh(doc, box_obj.Name + "_Mesh")
 femmesh_obj.Part = box_obj
 doc.recompute()
-box_obj.ViewObject.Visibility = False
 
 from femmesh.gmshtools import GmshTools as gt
 gmsh_mesh = gt(femmesh_obj)
@@ -842,6 +847,40 @@ doc.recompute()
 
 """
 
+
+"""
+# more sophisticated example which changes the mesh size
+import Part, ObjectsFem
+doc = App.ActiveDocument
+
+box_obj = doc.addObject("Part::Box", "Box")
+doc.recompute()
+box_obj.ViewObject.Visibility = False
+
+from femmesh.gmshtools import GmshTools
+max_mesh_sizes = [0.5, 1, 2, 3, 5, 10]
+for len in max_mesh_sizes:
+    quantity_len = "{}".format(len)
+    print("\n\n Start length = {}".format(quantity_len))
+    femmesh_obj = ObjectsFem.makeMeshGmsh(doc, box_obj.Name + "_Mesh")
+    femmesh_obj.Part = box_obj
+    femmesh_obj.CharacteristicLengthMax = "{}".format(quantity_len)
+    femmesh_obj.CharacteristicLengthMin = "{}".format(quantity_len)
+    doc.recompute()
+    gm = GmshTools(femmesh_obj)
+    gm.update_mesh_data()
+    # set the tmp file path to some user path including the length
+    gm.get_tmp_file_paths("/tmp/fcgm_" + str(len), True)
+    gm.get_gmsh_command()
+    gm.write_gmsh_input_files()
+    error = gm.run_gmsh_with_geo()
+    print(error)
+    gm.read_and_set_new_mesh()
+    doc.recompute()
+    print("Done length = {}".format(quantity_len))
+
+"""
+
 """
 TODO
 class GmshTools should be splittet in two classes