aboutsummaryrefslogtreecommitdiff
path: root/src/geod_interface.c
AgeCommit message (Collapse)Author
2015-10-01Take 2 of the fix to dmstor.Charles Karney
I'm making this pull request soon after the release of 4.9.2. It will cause the results that people get out of proj.4 to change very slightly. If there are problems, we'll get a chance to iron them out well before the next release. The important change is to use DEG_TO_RAD for degree to radian conversions in dmstor.c instead of the slightly inaccurate number used earlier. This necessitates a change to geod_interface.c (which previously had to work aroung the previous bad behavior). PJ_aeqd.c now does conversions in a compatible manner. In src/proj_api.h, DEG_TO_RAD and RAD_TO_DEG are both given to 17 significant figures (this is just a cosmetic change). I've "fixed" the testvarious tests so that they still pass (on my system at least). Everyone should be suitably skeptical of my fixes. (1) Minor changes to "Test bug 244" and only ask for nanometer (instead of picometer) accuracy on positions. (2) 2 lon_wrap tests now return 0dE instead of 360dE (now it's tight?) (3) The results for the forward healpix projection of (-180, +90) and (-180, -90) are now different. I have put the new values into tv_out.dist. I'm fairly confident that the new values are OK, since this projection has various cuts which meet at the poles. It would be good if someone who knows about this projection can verify this.
2015-08-16Drop in the latest geodesic library from GeographicLib (version 1.44).Charles Karney
http://geographiclib.sourceforge.net/1.44/C/index.html The changes are: - Improve accuracy of calculations by evaluating trigonometric functions more carefully and replacing the series for the reduced length with one with a smaller truncation error. - The allowed ranges for longitudes and azimuths is now unlimited; it used to be [-540d, 540d). - Enforce the restriction of latitude to [-90d, 90d] by returning NaNs if the latitude is outside this range. - The inverse calculation sets s12 to zero for coincident points at pole (instead of returning a tiny quantity). This commit also includes a work-around for an inaccurate value for pi/180 in dmstor.c (see the definitions of DEG_IN and DEG_OUT in geod_interface.c).
2013-05-10Major upgrade to geodesic support from Charles Karney (#197). Frank Warmerdam
Syncs geodesic routines with GeographicLib. Adds geodesic.3 man page. geod_* api exposed publically. geodesic.h is installed. git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@2333 4e78687f-474d-0410-85f9-8d5e500ac6b2
2012-12-07replacement of geodesic engine with one from Charles Karney (#197)Frank Warmerdam
git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@2297 4e78687f-474d-0410-85f9-8d5e500ac6b2