Debug Issue #8

#8

pykonal/__version__.py

@@ -1,7 +1,7 @@
 __major_version__  = 0
 __minor_version__  = 3
 __patch__          = 2
-__release__        = "b1"
+__release__        = "b2"
 __version_tuple__ = (
     __major_version__,
     __minor_version__,

pykonal/fields.pyx

@@ -296,7 +296,14 @@ cdef class Field3D(object):
         :rtype: numpy.ndarray(shape=(N0,N1,N2,3), dtype=numpy.float)
         """
         if self.coord_sys == "spherical" and coord_sys.lower() == "spherical":
-            return (transformations.sph2sph(self.nodes, origin))
+            force_phi_positive = self.min_coords[2] >= 0
+            return (
+                transformations.sph2sph(
+                    self.nodes,
+                    origin,
+                    force_phi_positive=force_phi_positive
+                )
+            )
         elif self.coord_sys == "cartesian" and coord_sys.lower() == "spherical":
             return (transformations.xyz2sph(self.nodes, origin))
         elif self.coord_sys == "spherical" and coord_sys.lower() == "cartesian":

pykonal/solver.pyx

@@ -634,7 +634,10 @@ class PointSourceSolver(EikonalSolver):
         origin = (r0, t0, p0)
         # Transform the coordinates of the near-field grid nodes to the
         # far-field coordinate system.
-        nodes = self.near_field.vv.transform_coordinates(self.coord_sys, origin)
+        nodes = self.near_field.vv.transform_coordinates(
+            self.coord_sys,
+            origin
+        )
         # Find the indices of the near-field grid nodes that are inside
         # the far-field grid.
         bool_idx = True

pykonal/transformations.py

@@ -26,7 +26,7 @@ def geo2sph(nodes):
     return (sph)
 
 
-def sph2sph(nodes, origin=(0,0,0)):
+def sph2sph(nodes, origin=(0,0,0), force_phi_positive=False):
     """
     Transform spherical coordinates to new spherical coordinate system.
 
@@ -55,11 +55,13 @@ def sph2sph(nodes, origin=(0,0,0)):
     tt  = np.arccos(xyz[...,2] / rr)
     np.seterr(**old)
     pp  = np.arctan2(xyz[...,1], xyz[...,0])
+    if force_phi_positive:
+        pp = np.mod(pp, 2*np.pi)
     rtp = np.moveaxis(np.stack([rr, tt, pp]), 0, -1)
     return (rtp)
 
 
-def xyz2sph(nodes, origin=(0,0,0)):
+def xyz2sph(nodes, origin=(0,0,0), force_phi_positive=False):
     """
     Transform Cartesian coordinates to new spherical coordinate system.
 
@@ -78,6 +80,8 @@ def xyz2sph(nodes, origin=(0,0,0)):
     tt  = np.arccos(xyz[...,2] / rr)
     np.seterr(**old)
     pp  = np.arctan2(xyz[...,1], xyz[...,0])
+    if force_phi_positive:
+        pp = np.mod(pp, 2*np.pi)
     rtp = np.moveaxis(np.stack([rr, tt, pp]), 0, -1)
     return (rtp)