Merge pull request #84 from mcrrobinson/master

Reusing 3x3 array

Author: leif81

opendis/RangeCoordinates.py

@@ -112,6 +112,8 @@ class GPS:
 
     wgs84 = WGS84()
 
+    _IDENTITY_3X3 = identity(3)
+
     def ecef2lla(self, ecef, tolerance=1e-9):
         """Convert Earth-centered, Earth-fixed coordinates to lat, lon, alt.
         Input: ecef - (x, y, z) in (m, m, m)
@@ -516,20 +518,16 @@ def rotate_3x3(self, theta, normal_vec):
               the rotation.
         """
 
-        n_v_t = normal_vec.transpose()
-
         n_x = array([[      0  , -normal_vec[2], normal_vec[1] ],
                      [normal_vec[2] ,       0  , -normal_vec[0]],
                      [-normal_vec[1], normal_vec[0] ,   0      ]
                     ]
                    )
 
 
-        I = identity(3)
-
         nnt = normal_vec * self.transpose(normal_vec)
 
-        return (1 - cos(theta)) * nnt + cos(theta) * I + sin(theta) * n_x
+        return (1 - cos(theta)) * nnt + cos(theta) * self._IDENTITY_3X3 + sin(theta) * n_x
 
 
     def transpose(self, arr):