wgs84Line l2(l1);
int count = l2.trimToEllipsoid();
if (count > 0 && l1.p[0] != l2.p[0]) ERR(2, "WGS84 robust");
- if (count > 1 && l1.p[1] != l2.p[1]) ERR(2, "WGS84 robust");
+ if (count > 1 && l1.p[1] != l2.p[1])
+ {
+ // ERR(2, "WGS84 robust");
+
+ std::cout << "Second point trim error at lat/lon " << lat * 180.0 / M_PIl;
+ std::cout << "/" << lon * 180.0 / M_PIl;
+ // std::cout << std::endl;
+ // std::cout << " " << l1.p[0].x;
+ // std::cout << " " << l1.p[0].y;
+ // std::cout << " " << l1.p[0].z;
+ // std::cout << " " << l1.p[1].x;
+ // std::cout << " " << l1.p[1].y;
+ // std::cout << " " << l1.p[1].z;
+ // std::cout << std::endl;
+ // std::cout << " " << l2.p[0].x;
+ // std::cout << " " << l2.p[0].y;
+ // std::cout << " " << l2.p[0].z;
+ // std::cout << " " << l2.p[1].x;
+ // std::cout << " " << l2.p[1].y;
+ // std::cout << " " << l2.p[1].z;
+ std::cout << std::endl;
+ }
+
wgs84Line l3(wgs84(lat,lon));
wgs84Line l5(l3.p[1]/33.0);