diff options
Diffstat (limited to 'src/PJ_molodensky.c')
| -rw-r--r-- | src/PJ_molodensky.c | 80 |
1 files changed, 0 insertions, 80 deletions
diff --git a/src/PJ_molodensky.c b/src/PJ_molodensky.c index 73d0e5c2..1b0eb3a1 100644 --- a/src/PJ_molodensky.c +++ b/src/PJ_molodensky.c @@ -312,83 +312,3 @@ PJ *TRANSFORMATION(molodensky,1) { return P; } - -#ifndef PJ_SELFTEST -int pj_molodensky_selftest (void) {return 0;} -#else -int pj_molodensky_selftest (void) { - - PJ_COORD in, ni, res, exp; - PJ *P; - - /* Test the abridged Molodensky first. Example from appendix 3 of Deakin (2004). */ - P = proj_create(PJ_DEFAULT_CTX, - "+proj=molodensky +a=6378160 +rf=298.25 " - "+da=-23 +df=-8.120449e-8 +dx=-134 +dy=-48 +dz=149 " - "+abridged " - ); - if (0==P) - return 10; - - in.lpz.lam = PJ_TORAD(144.9667); - in.lpz.phi = PJ_TORAD(-37.8); - in.lpz.z = 50.0; - - exp.lpz.lam = PJ_TORAD(144.968); - exp.lpz.phi = PJ_TORAD(-37.79848); - exp.lpz.z = 46.378; - - res = proj_trans(P, PJ_FWD, in); - - if (proj_lp_dist(P, res.lp, exp.lp) > 2 ) { /* we don't expect much accurecy here... */ - proj_destroy(P); - return 11; - } - - /* let's try a roundtrip */ - ni = in; - if (proj_roundtrip(P, PJ_FWD, 100, &ni) > 1) { - proj_destroy(P); - return 12; - } - - if (res.lpz.z - exp.lpz.z > 1e-3) { - proj_destroy(P); - return 13; - } - - proj_destroy(P); - - /* Test the abridged Molodensky first. Example from appendix 3 of Deaking (2004). */ - - P = proj_create(PJ_DEFAULT_CTX, - "+proj=molodensky +a=6378160 +rf=298.25 " - "+da=-23 +df=-8.120449e-8 +dx=-134 +dy=-48 +dz=149 " - ); - if (0==P) - return 20; - - res = proj_trans(P, PJ_FWD, in); - - if (proj_lp_dist(P, res.lp, exp.lp) > 2 ) { /* we don't expect much accurecy here... */ - proj_destroy(P); - return 21; - } - - /* let's try a roundtrip */ - ni = in; - if (proj_roundtrip(P, PJ_FWD, 100, &ni) > 1) { - proj_destroy(P); - return 22; - } - - if (res.lpz.z - exp.lpz.z > 1e-3) { - proj_destroy(P); - return 23; - } - - proj_destroy(P); - return 0; -} - -#endif |
