1 # contains function to compute a mesh in parallel
6 notebook = salome_notebook.NoteBook()
13 from salome.geom import geomBuilder
14 from salome.smesh import smeshBuilder
19 geompy = geomBuilder.New()
27 # First creating all the boxes
32 x_orig = i*(boxsize+offset)
33 y_orig = j*(boxsize+offset)
34 z_orig = k*(boxsize+offset)
36 tmp_box = geompy.MakeBoxDXDYDZ(boxsize, boxsize, boxsize)
38 if not i == j == k == 0:
39 box = geompy.MakeTranslation(tmp_box, x_orig,
44 geompy.addToStudy(box, 'box_{}:{}:{}'.format(i, j, k))
48 # Create fuse of all boxes
49 all_boxes = geompy.MakeCompound(boxes)
50 geompy.addToStudy(all_boxes, 'Compound_1')
52 # Removing duplicates faces and edges
53 all_boxes = geompy.MakeGlueFaces(all_boxes, 1e-07)
54 geompy.addToStudy(all_boxes, 'Glued_Faces_1')
56 rubik_cube = geompy.MakeGlueEdges(all_boxes, 1e-07)
57 geompy.addToStudy(all_boxes, 'rubik_cube')
59 smesh = smeshBuilder.New()
60 print("Creating Parallel Mesh")
61 par_mesh = smesh.ParallelMesh(rubik_cube, name="par_mesh", mesher3D="GMSH")
63 print("Creating hypoehtesis for netgen")
64 NETGEN_3D_Parameters_1 = smesh.CreateHypothesisByAverageLength( 'NETGEN_Parameters',
65 'NETGENEngine', 34.641, 0 )
66 print("Adding hypothesis")
67 par_mesh.AddGlobalHypothesis(NETGEN_3D_Parameters_1, mesher="GMSH" )
69 #Set here particular mesh parameters in 3D
70 for algo3d in par_mesh._algo3d:
71 param3d = algo3d.Parameters()
74 param3d.SetSmouthSteps(2)
75 param3d.SetSizeFactor(1.1)
77 print("Setting parallelism method")
78 par_mesh.SetParallelismMethod(smeshBuilder.MULTITHREAD)
80 print("Setting parallelism options")
81 param = par_mesh.GetParallelismSettings()
84 print("Starting parallel compute")
85 is_done = par_mesh.Compute()
87 raise Exception("Error when computing Mesh")
89 print(" Tetrahedron: ", par_mesh.NbTetras())
90 print(" Triangle: ", par_mesh.NbTriangles())
91 print(" edge: ", par_mesh.NbEdges())
93 assert par_mesh.NbTetras() > 0
94 assert par_mesh.NbTriangles() > 0
95 assert par_mesh.NbEdges() > 0