aboutsummaryrefslogtreecommitdiff
path: root/src/PJ_molodensky.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/PJ_molodensky.c')
-rw-r--r--src/PJ_molodensky.c80
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