int main(int argc, char **argv)
{
#ifdef USE_WGS84
plan_tests(9 + 36);
#else
plan_tests(9 + 36 + 18);
#endif
const GeoPoint a(Angle::Degrees(7.7061111111111114),
Angle::Degrees(51.051944444444445));
const GeoPoint b(Angle::Degrees(7.599444444444444),
Angle::Degrees(51.099444444444444));
const GeoPoint c(Angle::Degrees(4.599444444444444),
Angle::Degrees(47.099444444444444));
fixed distance = Distance(a, b);
#ifdef USE_WGS84
ok1(distance > fixed(9150) && distance < fixed(9160));
#else
ok1(distance > fixed(9130) && distance < fixed(9140));
#endif
Angle bearing = Bearing(a, b);
ok1(bearing.Degrees() > fixed(304));
ok1(bearing.Degrees() < fixed(306));
bearing = Bearing(b, a);
ok1(bearing.Degrees() > fixed(124));
ok1(bearing.Degrees() < fixed(126));
distance = ProjectedDistance(a, b, a);
ok1(is_zero(distance));
distance = ProjectedDistance(a, b, b);
#ifdef USE_WGS84
ok1(distance > fixed(9150) && distance < fixed(9180));
#else
ok1(distance > fixed(9120) && distance < fixed(9140));
#endif
const GeoPoint middle(a.longitude.Fraction(b.longitude, fixed(0.5)),
a.latitude.Fraction(b.latitude, fixed(0.5)));
distance = ProjectedDistance(a, b, middle);
#ifdef USE_WGS84
ok1(distance > fixed(9150/2) && distance < fixed(9180/2));
#else
ok1(distance > fixed(9100/2) && distance < fixed(9140/2));
#endif
fixed big_distance = Distance(a, c);
ok1(big_distance > fixed(494000) && big_distance < fixed(495000));
TestLinearDistance();
return exit_status();
}
int main(int argc, char **argv)
{
plan_tests(9 + 36 + 18);
const GeoPoint a(Angle::degrees(fixed(7.7061111111111114)),
Angle::degrees(fixed(51.051944444444445)));
const GeoPoint b(Angle::degrees(fixed(7.599444444444444)),
Angle::degrees(fixed(51.099444444444444)));
const GeoPoint c(Angle::degrees(fixed(4.599444444444444)),
Angle::degrees(fixed(47.099444444444444)));
fixed distance = Distance(a, b);
ok1(distance > fixed(9130) && distance < fixed(9140));
Angle bearing = Bearing(a, b);
ok1(bearing.value_degrees() > fixed(304));
ok1(bearing.value_degrees() < fixed(306));
bearing = Bearing(b, a);
ok1(bearing.value_degrees() > fixed(124));
ok1(bearing.value_degrees() < fixed(126));
distance = ProjectedDistance(a, b, a);
ok1(is_zero(distance));
distance = ProjectedDistance(a, b, b);
ok1(distance > fixed(9120) && distance < fixed(9140));
const GeoPoint middle(a.Longitude.Fraction(b.Longitude, fixed_half),
a.Latitude.Fraction(b.Latitude, fixed_half));
distance = ProjectedDistance(a, b, middle);
ok1(distance > fixed(9100/2) && distance < fixed(9140/2));
fixed big_distance = Distance(a, c);
ok1(big_distance > fixed(494000) && big_distance < fixed(495000));
TestLinearDistance();
return exit_status();
}
