Minor patch #6

• Some refactoring of code
• Use -1 to represent NULL value instead of np.nan which was causing insidious errors
• Account for 2π periodicity when interpolating in spherical coordinates
• Set EikonalSolver.vvp = EikonalSolver.vv when EikonalSolver.vgrid == EikonalSolver.pgrid

Minor patch #5

• Debug EikonalSolver._get_gradient() which was returning None for Cartesian coordinates.
• Migrate demo notebook to current API

Minor patch #4

Improve stability by dropping imaginary part of solution to quadratic equation. This was previously done only when the discriminant was small negative, but now is done even when it is large negative.

Patch to improve stability

• When the determinant of the quadratic equation in the update scheme is negative but close to zero, an approximate solution to the quadratic will be accepted.

Minor patch to improve stability

• LinearInterpolator3D is returning np.inf and should raise an error instead. This patch checks for np.inf values in the client code, but should be fixed on the LinearInterpolator3D side.

Alpha release implementing ray tracing in spherical coordinates

0.1a1

UPGRADE:: Implement ray-tracing in spherical coordinates

Alpha release including periodic spherical coordinates.

• API is not backwards compatible will 0.0 versions.
• This version implements spherical coordinates.
• Periodicity of spherical coordinate systems is automatically accounted for.
• This version sorted out a bug with the heap sort function.

Beta version of 3D code: minor update 6

• Ray tracing is significantly faster and more stable

Beta version of 3D code: minor update 4

• Debug error during edge case of interpolation method.
• Update examples.
• Update method documentation.

Beta version of 3D code: minor update 3

• Optimized ray-tracer
• Cleaned up some code
• Added detail to setup.py