@@ -1239,9 +1239,10 @@ def computeGeometricInfo(self):
12391239
12401240 def buildMatlabMeshDataStructures (self ,meshFileBase = 'meshMatlab' ,writeToFile = True ):
12411241 """
1242- build array data structures for matlab finite element mesh representation
1243- and write to a file to view and play with in matlatb. The current matlab support
1244- is mostly for 2d, but this will return basic arrays for 1d and 3d too
1242+ build array data structures for matlab finite element mesh
1243+ representation and write to a file to view and play with in
1244+ matlatb. The current matlab support is mostly for 2d, but this
1245+ will return basic arrays for 1d and 3d too
12451246
12461247 in matlab can then print mesh with
12471248
@@ -1262,6 +1263,7 @@ def buildMatlabMeshDataStructures(self,meshFileBase='meshMatlab',writeToFile=Tru
12621263 row 1 = x coord,
12631264 row 2 = y coord for nodes in mesh
12641265 row 3 = z coord for nodes in mesh ...
1266+
12651267 edge matrix is [2*nd+3 x num faces]
12661268 format:
12671269 row 1 = start vertex number
@@ -1280,10 +1282,11 @@ def buildMatlabMeshDataStructures(self,meshFileBase='meshMatlab',writeToFile=Tru
12801282 ...
12811283 row nd+1 = vertex 3 global number
12821284 row 4 = triangle subdomain number
1285+
12831286 where 1,2,3 is a local counter clockwise numbering of vertices in
12841287 triangle
12851288
1286- """
1289+ """
12871290 matlabBase = 1
12881291 nd = self .nNodes_element - 1
12891292 p = np .zeros ((nd ,self .nNodes_global ),'d' )
@@ -6085,10 +6088,6 @@ def getMeshIntersections(mesh, toPolyhedron, endpoints):
60856088 intersections .update (((tuple (elementIntersections [0 ]), tuple (elementIntersections [1 ])),),)
60866089 return intersections
60876090
6088-
6089- from proteus import default_n as dn
6090- from proteus import default_p as dp
6091-
60926091class MeshOptions :
60936092 """
60946093 Mesh options for the domain
@@ -6101,17 +6100,17 @@ def __init__(self, domain):
61016100 self .Domain = domain
61026101 self .he = 1.
61036102 self .use_gmsh = False
6104- self .genMesh = dp . genMesh
6103+ self .genMesh = True
61056104 self .outputFiles_name = 'mesh'
61066105 self .outputFiles = {'poly' : True ,
61076106 'ply' : False ,
61086107 'asymptote' : False ,
61096108 'geo' : False }
6110- self .restrictFineSolutionToAllMeshes = dn . restrictFineSolutionToAllMeshes
6111- self .parallelPartitioningType = dn . parallelPartitioningType
6112- self .nLayersOfOverlapForParallel = dn . parallelPartitioningType
6113- self .triangleOptions = dn . triangleOptions # defined when setTriangleOptions called
6114- self .nLevels = dn . nLevels
6109+ self .restrictFineSolutionToAllMeshes = False
6110+ self .parallelPartitioningType = MeshParallelPartitioningTypes . node
6111+ self .nLayersOfOverlapForParallel = 1
6112+ self .triangleOptions = "q30DenA" # defined when setTriangleOptions called
6113+ self .nLevels = 1
61156114 if domain is not None :
61166115 self .nd = domain .nd
61176116 if self .nd == 2 :
0 commit comments