-currently: working to address the instability of the Voronoi code relating to the radius of input sphere data (i.e., 2.0 works, 1.9999 doesn't; 87.0 doesn't)
-the orientation of the Delaunay triangle facet normals should be checked more robustly now-- I normalized Delaunay triangle centroid coordinates and this seems to help
-minor fix to the way that the sphere centroid array is assigned
-currently: trying to replace some of the arc cosine distance calculations with more robust haversine calculations to avoid math domain errors (but this is causing accuracy to drop for some surface area calculations)
-wrote a simple unit test for the new haversine code (gets the correct distance to several decimal places)
-wrote a haversine unit test for a sphere of larger radius