Changeset 6463 for branches/numpy/anuga/abstract_2d_finite_volumes/general_mesh.py

Timestamp:

Mar 5, 2009, 10:57:08 AM (15 years ago)

Author:

rwilson

Message:

Modified General_mesh.get_node(absolute=True) to return a copy if modifying the object node data.

File:

• branches/numpy/anuga/abstract_2d_finite_volumes/general_mesh.py (modified) (3 diffs)

branches/numpy/anuga/abstract_2d_finite_volumes/general_mesh.py

@@ +1 @@
+import copy
 import numpy as num
 
@@ -252 +253 @@
         Default is False as many parts of ANUGA expects relative coordinates.
         (To see which, switch to default absolute=True and run tests).
+
+        Note: This method returns a modified _copy_ of the nodes slice if 
+              absolute is True.  If absolute is False, just return the slice.
+              This is related to the ensure_numeric() returning a copy problem.
         """
 
@@ -257 +262 @@
         if absolute is True:
             if not self.geo_reference.is_absolute():
+                # get a copy so as not to modify the internal self.nodes array
+                V = copy.copy(V)
                 V += num.array([self.geo_reference.get_xllcorner(),
                                 self.geo_reference.get_yllcorner()], num.float)