Loading...
Searching...
No Matches
test_loncxlatc.c
This contains examples on yac_loncxlatc_vec
// Copyright (c) 2024 The YAC Authors
//
// SPDX-License-Identifier: BSD-3-Clause
#include "test_cxc.h"
static void test_loncxlatc(double lon_a, double lat_a, double lon_b, double lat_b,
double lon_c, double lat_c, double lon_d, double lat_d,
double lon_ref_p, double lat_ref_p,
double lon_ref_q, double lat_ref_q, int ref_ret_val);
unsigned const p_between_ab = 1 << 0;
unsigned const q_between_ab = 1 << 1;
unsigned const p_between_cd = 1 << 2;
// circles that intersect
test_loncxlatc( 0, 10, // point a
0, -10, // point b
10, 0, // point c
-10, 0, // point d
0, 0, // reference point p
180, 0, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect
test_loncxlatc(65, 35, // point a
65, 45, // point b
60, 40, // point c
70, 40, // point d
65, 40, // reference point p
65+180, 40, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect but the edges not
test_loncxlatc(-30, -60, // point a
-30, -70, // point b
-31, -65, // point c
-40, -65, // point d
-30, -65, // reference point p
-30+180, -65, // reference point q
p_between_ab); // reference return value
// circles that intersect and the edges touch
test_loncxlatc(-30, -60, // point a
-30, -70, // point b
-30, -65, // point c
-40, -65, // point d
-30, -65, // reference point p
-30+180, -65, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect but the edges not
test_loncxlatc(-45, 80, // point a
-45, 90, // point b
-40, 70, // point c
-50, 70, // point d
-45, 70, // reference point p
-45+180, 70, // reference point q
p_between_cd); // reference return value
// circles that intersect and the edges touch
test_loncxlatc(-45, 70, // point a
-45, 90, // point b
-40, 70, // point c
-50, 70, // point d
-45, 70, // reference point p
-45+180, 70, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect at the north pole and the lon
// edge crosses the north pole
test_loncxlatc(-45, 80, // point a
-45+180, 80, // point b
20, 90, // point c
30, 90, // point d
0, 90, // reference point p
0, -90, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect at the south pole and the lon
// edge crosses the south pole
test_loncxlatc(-45, -80, // point a
-45+180, -80, // point b
20, -90, // point c
30, -90, // point d
0, -90, // reference point p
0, 90, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect at the north pole and the lon
// edge does not cross the north pole
test_loncxlatc(-45, 80, // point a
-45, 85, // point b
20, 90, // point c
30, 90, // point d
0, 90, // reference point p
0, -90, // reference point q
p_between_cd); // reference return value
// circles that intersect at the south pole and the lon
// edge does not cross the south pole
test_loncxlatc(-45, -80, // point a
-45, -85, // point b
20, -90, // point c
30, -90, // point d
0, -90, // reference point p
0, 90, // reference point q
p_between_cd); // reference return value
// circles that intersect at the south pole and the lon
// edge touches the south pole
test_loncxlatc(-45, -80, // point a
-45, -90, // point b
20, -90, // point c
30, -90, // point d
0, -90, // reference point p
0, 90, // reference point q
p_between_ab + p_between_cd); // reference return value
// circles that intersect close to the south pole and the lon
// edge intersects the lat circle twice
test_loncxlatc( 0, -80, // point a
180, -80, // point b
50, -89, // point c
60, -89, // point d
0, -89, // reference point p
180, -89, // reference point q
p_between_ab + q_between_ab); // reference return value
test_loncxlatc( 0, -80, // point a
180, -80, // point b
-10, -89, // point c
10, -89, // point d
0, -89, // reference point p
180, -89, // reference point q
p_between_ab + p_between_cd + q_between_ab); // reference return value
// longitude circle consists of two identical points and is not on
// the latitude circle
test_loncxlatc( 0, 0, // point a
0, 0, // point b
5, 10, // point c
5, 10, // point d
-1, -1, // reference point p
-1, -1, // reference point q
-1); // reference return value
test_loncxlatc( 0, 0, // point a
0, 0, // point b
5, 10, // point c
10, 10, // point d
-1, -1, // reference point p
-1, -1, // reference point q
-1); // reference return value
// longitude circle consists of two identical points and is on the
// latitude circle
test_loncxlatc( 0, 0, // point a
0, 0, // point b
5, 0, // point c
-5, 0, // point d
0, 0, // reference point p
180, 0, // reference point q
p_between_ab + p_between_cd); // reference return value
// latitude circle consists of two identical points and is not on
// the longitude circle
test_loncxlatc( 0, -5, // point a
0, 5, // point b
5, 0, // point c
5, 0, // point d
-1, -1, // reference point p
-1, -1, // reference point q
-1); // reference return value
test_loncxlatc( 0, -5, // point a
0, 5, // point b
5, 10, // point c
5, 10, // point d
-1, -1, // reference point p
-1, -1, // reference point q
-1); // reference return value
// longitude circle consists of two identical points and is on the
// latitude circle
test_loncxlatc( 0, -5, // point a
0, 5, // point b
0, 0, // point c
0, 0, // point d
0, 0, // reference point p
180, 0, // reference point q
p_between_ab + p_between_cd); // reference return value
test_loncxlatc( 0, 90, // point a
140, 60, // point b
90, -45, // point c
39, -45, // point d
140, -45, // reference point p
-40, -45, // reference point q
0); // reference return value
// both edges have length zero and are on a pole
for (int i = 0; i < 4; ++i) {
double lat_coords[2] = {90, -90};
int ret_val[2] = {-1, p_between_ab + p_between_cd};
test_loncxlatc( 20, lat_coords[i&1], // point a
200, lat_coords[i&1], // point b
10, lat_coords[i>>1], // point c
20, lat_coords[i>>1], // point d
0, lat_coords[i>>1], // reference point p
0, -lat_coords[i>>1], // reference point q
ret_val[(i&1)==(i>>1)]); // reference return value
}
{
double a[3] = {-0.01728256081754178,
-0.0024289055230921597,
0.99984769515639127};
double b[3] = {-0.034559857199638548,
-0.0048570711780326408,
0.99939082701909576};
double c[3] = {-0.017298289573378096,
-0.0023142316840916847,
0.99984769515639127};
double d[3] = { 0.01732231884620395,
0.0021269133133512588,
0.99984769515639127};
double a_lon, a_lat, b_lon, b_lat, c_lon, c_lat, d_lon, d_lat;
XYZtoLL(a, &a_lon, &a_lat);
XYZtoLL(b, &b_lon, &b_lat);
XYZtoLL(c, &c_lon, &c_lat);
XYZtoLL(d, &d_lon, &d_lat);
p_between_ab + p_between_cd); // reference return value
}
return TEST_EXIT_CODE;
}
static void test_loncxlatc(double lon_a, double lat_a, double lon_b, double lat_b,
double lon_c, double lat_c, double lon_d, double lat_d,
double lon_ref_p, double lat_ref_p,
double lon_ref_q, double lat_ref_q, int ref_ret_val) {
test_cxc(lon_a, lat_a, lon_b, lat_b, lon_c, lat_c, lon_d, lat_d,
lon_ref_p, lat_ref_p, lon_ref_q, lat_ref_q, ref_ret_val);
}
Generated on Wed Dec 4 2024 11:31:42 for YetAnotherCoupler by
1.11.0