Changeset 838


Ignore:
Timestamp:
Feb 7, 2005, 3:12:40 PM (20 years ago)
Author:
duncan
Message:

pmesh/mesh - checking in speed-up changes I was half way thru

Location:
inundation/ga/storm_surge/pmesh
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • inundation/ga/storm_surge/pmesh/mesh.py

    r829 r838  
    13061306        #print [xmin, ymin, xmax, ymax]
    13071307     
     1308    def boxsizeVerts(self):
     1309        """
     1310        Returns a list of verts denoting a box or triangle that contains verts on the xmin, ymin, xmax and ymax axis.
     1311        Structure: list of verts
     1312        """
     1313        # FIXME dsg!!! large is a hack
     1314        #You want the kinds package, part of Numeric:
     1315        #In [2]: import kinds
     1316       
     1317        #In [3]: kinds.default_float_kind.M
     1318        #kinds.default_float_kind.MAX         kinds.default_float_kind.MIN
     1319    #kinds.default_float_kind.MAX_10_EXP  kinds.default_float_kind.MIN_10_EXP
     1320        #kinds.default_float_kind.MAX_EXP     kinds.default_float_kind.MIN_EXP
     1321
     1322        #In [3]: kinds.default_float_kind.MIN
     1323        #Out[3]: 2.2250738585072014e-308
     1324
     1325        large = 1e100
     1326        xmin= large
     1327        xmax=-large
     1328        ymin= large
     1329        ymax=-large
     1330        for vertex in self.userVertices:
     1331            if vertex.x < xmin:
     1332                xmin = vertex.x
     1333                xminVert = vertex
     1334            if vertex.x > xmax:
     1335                xmax = vertex.x
     1336                xmaxVert = vertex
     1337               
     1338            if vertex.y < ymin:
     1339                ymin = vertex.y
     1340                yminVert = vertex
     1341            if vertex.y > ymax:
     1342                ymax = vertex.y
     1343                ymaxVert = vertex
     1344        verts, count = self.removeDuplicatedVertices([xminVert,xmaxVert,yminVert,ymaxVert])
     1345         
     1346        return verts
     1347   
    13081348    def boxsize(self):
    13091349        """
  • inundation/ga/storm_surge/pmesh/pmesh.py

    r830 r838  
    11181118            self.clearMesh()
    11191119            self.mesh = newmesh
     1120
     1121            #FIXME - to speed things up, don't visualise the mesh
     1122            # use ResizeToFitWrapper
    11201123            self.visualiseMesh(self.mesh)
    11211124            self.ResizeToFit()
     
    11301133            showerror('File error',
    11311134                      'file ' + ofile + ' could not be found.')
     1135
    11321136        except RuntimeError:
    11331137            showerror('File error',
Note: See TracChangeset for help on using the changeset viewer.