Changeset 6442 for branches/numpy/anuga/coordinate_transforms/geo_reference.py

Timestamp:

Mar 4, 2009, 8:37:22 AM (14 years ago)

Author:

rwilson

Message:

Added tests for get_relative().

File:

• branches/numpy/anuga/coordinate_transforms/geo_reference.py (modified) (3 diffs)

branches/numpy/anuga/coordinate_transforms/geo_reference.py

@@ -318 +318 @@
                 msg = 'Single point must have two elements'
                 raise ShapeError, msg
-                #points = reshape(points, (1,2))
 
         if not points.shape[1] == 2:
@@ -348 +347 @@
         """
 
-        is_list = False
-        if type(points) == types.ListType:
-            is_list = True
-
-        points = ensure_numeric(points, num.float)
+        # remember if we got a list
+        is_list = isinstance(points, list)
+
+        # convert to a numeric array, force a copy
+        points = ensure_numeric(copy.copy(points), num.float)
+
+        # sanity checks
         if len(points.shape) == 1:
             #One point has been passed
@@ -358 +359 @@
                 msg = 'Single point must have two elements'
                 raise ShapeError, msg
-                #points = reshape(points, (1,2))
 
         if not points.shape[1] == 2: