Changeset 4808 for anuga_core/source/anuga/utilities/quad.py

Timestamp:

Nov 12, 2007, 9:53:05 AM (17 years ago)

Author:

duncan

Message:

bug fix

File:

• anuga_core/source/anuga/utilities/quad.py (modified) (3 diffs)

anuga_core/source/anuga/utilities/quad.py

@@ -143 +143 @@
             point_id = int(args[1])
             x, y = self.mesh.get_node(point_id, absolute=True)
-
-            #print point_id, x, y
         elif len(args) == 3:
             x = float(args[1])
@@ -445 +443 @@
     #print mesh.coordinates
 
-    
-    nodes = mesh.get_nodes(absolute=True)
-    xmin = min(nodes[:,0])
-    xmax = max(nodes[:,0])
-    ymin = min(nodes[:,1])
-    ymax = max(nodes[:,1])
-
-    # Don't know why this didn't work
-    #xmin, xmax, ymin, ymax = mesh.get_extent(absolute=True)
+    xmin, xmax, ymin, ymax = mesh.get_extent(absolute=True)
     
     # Ensure boundary points are fully contained in region
@@ -473 +463 @@
     
     #FIXME: Use mesh.filename if it exists
+    # why?
     root = Cell(ymin, ymax, xmin, xmax,mesh,
-                #name = .... 
                 max_points_per_cell = max_points_per_cell)