Changeset 7317 for anuga_core/source/anuga/pmesh
- Timestamp:
- Jul 22, 2009, 9:22:11 AM (16 years ago)
- Location:
- anuga_core/source/anuga/pmesh
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/pmesh/graphical_mesh_generator.py
r5188 r7317 17 17 import load_mesh.loadASCII 18 18 from anuga.alpha_shape.alpha_shape import AlphaError 19 import anuga.utilities.log as log 19 20 20 21 # CONSTANTS … … 413 414 dialog = AddVertexDialog(self.canvas) 414 415 if dialog.xyValuesOk: 415 print dialog.x416 print dialog.y416 log.critical(str(dialog.x)) 417 log.critical(str(dialog.y)) 417 418 self.drawVertex(dialog.x*self.SCALE,dialog.y*self.SCALE,None) 418 419 #Since the new vertex may be off screen 419 420 self.ResizeToFit() 420 421 else: 421 print "bad values"422 log.critical("bad values") 422 423 423 424 def windowDefault (self, parent): … … 510 511 smooth_indents=dialog.smooth_indents.get(), 511 512 expand_pinch=dialog.expand_pinch.get()) 512 #print "newsegs",newsegs513 #print "ObjectsToVisuallyDelete",ObjectsToVisuallyDelete514 513 515 514 for drawOb in ObjectsToVisuallyDelete: … … 563 562 564 563 if dialog.ValuesOk: 565 print dialog.minAngle566 print dialog.maxArea564 log.critical(str(dialog.minAngle)) 565 log.critical(str(dialog.maxArea)) 567 566 568 567 self.clearSelections() … … 586 585 else: 587 586 pass 588 print "userMeshChanged = False"587 log.critical("userMeshChanged = False") 589 588 self.UserMeshChanged = False 590 589 self.visualiseMesh(self.mesh) 591 print "Mesh Generation finished"590 log.critical("Mesh Generation finished") 592 591 593 592 def MeshGenAreaAngle (self, minAngle, maxArea, mesh): … … 621 620 pass 622 621 meshArea = 0 623 #print "tempMesh.getTriangulation()", tempMesh.getTriangulation()624 622 meshArea = tempMesh.tri_mesh.calc_mesh_area() 625 623 maxArea = meshArea/numTriangles … … 639 637 #The screen canvas has y 'flipped'. -1* unflips it 640 638 self.lasty = self.starty = -1*self.canvas.canvasy(event.y) 641 print "----------------------"639 log.critical "----------------------") 642 640 self.mouseDownCurFunc( self.lastx, 643 641 self.lasty,event) #!!! remove the event? … … 651 649 if event.widget.find_withtag(CURRENT): # if there's a widget with a tag 652 650 [tag,string] = self.canvas.gettags(CURRENT) # get a list of them 653 print "tag",tag#tags ('M*1008', 'current')651 log.critical("tag %s" % str(tag)) #tags ('M*1008', 'current') 654 652 if tag[:2] == 'M*': #!!! this can be removed when there are 655 653 # only mesh objects on screen 656 654 #key, value = string.split(tag, '*') 657 655 objectID = tag 658 print "Found!! objectID:", objectID656 log.critical("Found!! objectID: %s" % str(objectID)) 659 657 660 658 meshObjects = self.getAllUserMeshObjects() … … 663 661 selMeshObject = meshObjects.getMeshObject(objectID) 664 662 found = True 665 print "Found! selMeshObject", selMeshObject 663 log.critical("Found! selMeshObject: %s" 664 % str(selMeshObject)) 666 665 #Only select one object at a time 667 666 if self.selMeshObject: … … 728 727 #x_scaled = self.SCALE*x 729 728 #y_scaled = -1*self.SCALE*y 730 #print "x y:", x,y731 729 vert = self.Vertices.draw(x-self.mesh.geo_reference.get_xllcorner, 732 730 y-self.mesh.geo_reference.get_yllcorner, … … 749 747 #x_scaled = self.SCALE*x 750 748 #y_scaled = -1*self.SCALE*y 751 #print "x y:", x,y752 749 vert = self.Vertices.draw(x,y,self.mesh,self.uniqueID,self.SCALE,self.canvas,event) 753 750 self.UserMeshChanged = True … … 784 781 if event.widget.find_withtag(CURRENT): # if there's a widget with a tag 785 782 [tag,string] = self.canvas.gettags(CURRENT) # get a list of them 786 print "tag",tag#tags ('M*1008', 'current')783 log.critical("tag %s" % str(tag)) #tags ('M*1008', 'current') 787 784 objectID = tag 788 #print "Found!! objectID:", objectID789 785 if self.Vertices.hasKey(objectID): #isinstance(self.meshObjects[objectID],mesh.Vertex): 790 786 vertex = self.Vertices.getMeshObject(objectID) 791 787 found = True 792 print "Found! vertex", vertex788 log.critical("Found! vertex: %s" % str(vertex)) 793 789 794 790 if found and self.selVertex == vertex: 795 print "The selected vertex has already been selected"791 log.critical("The selected vertex has already been selected") 796 792 #The selected vertex has already been selected 797 793 # therefore deselect it … … 808 804 self.selectVertex(vertex,objectID) 809 805 else: 810 print "vertex is the 1st vertex"806 log.critical("vertex is the 1st vertex") 811 807 #vertex is the 1st vertex 812 808 self.selectVertex(vertex,objectID) 813 809 else: 814 print " There are no widgets. This happen's too much"810 log.critical(" There are no widgets. This happen's too much") 815 811 816 812 … … 844 840 def printGeoReference(self): 845 841 try: 846 print "geo reference", self.mesh.geo_reference842 log.critical("geo reference %s" % str(self.mesh.geo_reference)) 847 843 except: 848 print "no geo reference"844 log.critical("no geo reference") 849 845 850 846 def visualiseMesh(self,mesh): … … 1007 1003 import mesh data from a variety of formats (currently 2!) 1008 1004 """ 1009 print "self.currentPath",self.currentPath1005 log.critical("self.currentPath %s" % str(self.currentPath)) 1010 1006 ofile = tkFileDialog.askopenfilename(initialdir=self.currentPath, 1011 1007 filetypes=[ ("text Mesh", … … 1021 1017 newmesh = mesh.importMeshFromFile(ofile) 1022 1018 self.currentPath, dummy = os.path.split(ofile) 1023 #print "be good self.currentPath",self.currentPath1024 1019 self.currentFilePathName = ofile 1025 1020 self.clearMesh() … … 1041 1036 # Could not get the file name to showup in the title 1042 1037 #appname = ofile + " - " + APPLICATION_NAME 1043 #print "appname",appname1044 1038 1045 1039 except load_mesh.loadASCII.TitleAmountError: … … 1086 1080 Save the current drawing 1087 1081 """ 1088 #print "dsg!!! self.currentFilePathName ",self.currentFilePathName1089 1082 if (self.currentFilePathName[-4:] != ".tsh" or 1090 1083 self.currentFilePathName[-4:] != ".msh"): … … 1120 1113 there is no generated mesh do not prompt. 1121 1114 """ 1122 #print "self.UserMeshChanged",self.UserMeshChanged1123 #print "self.mesh.isTriangulation()",self.mesh.isTriangulation()1124 1115 if (self.UserMeshChanged) and self.mesh.isTriangulation(): 1125 1116 … … 1220 1211 1221 1212 fd.close() 1222 print 'returning m'1213 log.critical('returning m') 1223 1214 return oadtestmesh(ofile) 1224 1215 -
anuga_core/source/anuga/pmesh/mesh.py
r7276 r7317 34 34 ensure_geospatial, ensure_absolute, ensure_numeric 35 35 from anuga.mesh_engine.mesh_engine import generate_mesh 36 import anuga.utilities.log as log 37 36 38 37 39 try: … … 145 147 x = scale*(self.x + xoffset) 146 148 y = -1*scale*(self.y + yoffset) # - since for a canvas - is up 147 #print "draw x:", x148 #print "draw y:", y149 149 cornerOffset= self.VERTEXSQUARESIDELENGTH/2 150 150 … … 187 187 x = scale*(self.x + xoffset) 188 188 y = -1*scale*(self.y + yoffset) # - since for a canvas - is up 189 #print "draw x:", x190 #print "draw y:", y191 189 cornerOffset= self.HOLECORNERLENGTH/2 192 190 return canvas.create_oval(x-cornerOffset, … … 612 610 # Depreciated 613 611 def addRegionEN(self, x,y): 614 print "depreciated, use add_region"612 log.critical("deprecated, use add_region") 615 613 return self.add_region(x,y) 616 614 … … 2266 2264 counter = newmesh.removeDuplicatedUserVertices() 2267 2265 if (counter >0): 2268 print "%i duplicate vertices removed from dataset" % (counter)2266 log.critical("%i duplicate vertices removed from dataset" % counter) 2269 2267 elif (ofile[-4:]== ".tsh" or ofile[-4:]== ".msh"): 2270 2268 dict = import_mesh_file(ofile) … … 2340 2338 region_list[i] = (region_list[i][0],region_list[i][1],i) 2341 2339 else: 2342 print "The region list has a bad size"2340 log.critical("The region list has a bad size") 2343 2341 # raise an error .. 2344 2342 raise Error -
anuga_core/source/anuga/pmesh/mesh_interface.py
r7276 r7317 5 5 import numpy as num 6 6 from anuga.utilities.polygon import inside_polygon 7 7 import anuga.utilities.log as log 8 8 9 9 … … 151 151 #raise Exception, msg 152 152 # Fixme: Use proper Python warning 153 if verbose: print 'WARNING: ', msg153 if verbose: log.critical('WARNING: %s' % msg) 154 154 155 155 … … 187 187 else: 188 188 msg += ' I will ignore it.' 189 print msg189 log.critical(msg) 190 190 191 191 else: … … 313 313 return m 314 314 else: 315 if verbose: print 'Generating mesh to file "%s"' %filename315 if verbose: log.critical("Generating mesh to file '%s'" % filename) 316 316 m.generate_mesh(minimum_triangle_angle=minimum_triangle_angle, 317 317 verbose=verbose)
Note: See TracChangeset
for help on using the changeset viewer.