PostGIS  3.6.1dev-r@@SVN_REVISION@@

◆ test_lw_dist2d_arc_arc()

static void test_lw_dist2d_arc_arc ( void  )
static

Definition at line 908 of file cu_measures.c.

909 {
910  /* lw_dist2d_arc_arc(const POINT2D *A1, const POINT2D *A2, const POINT2D *A3,
911  const POINT2D *B1, const POINT2D *B2, const POINT2D *B3,
912  DISTPTS *dl) */
913  DISTPTS dl;
914  POINT2D A1, A2, A3, B1, B2, B3;
915  int rv;
916 
917  /* Ticket #4326 */
918  lw_dist2d_distpts_init(&dl, DIST_MIN);
919  A1.x = -1.0; A1.y = 4.0;
920  A2.x = 0.0; A2.y = 5.0;
921  A3.x = 1.0; A3.y = 4.0;
922  B1.x = 1.0; B1.y = 6.0;
923  B2.x = 6.0; B2.y = 1.0;
924  B3.x = 9.0; B3.y = 7.0;
925  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
926  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
927  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.0475666, 0.000001);
928 
929  /* Unit semicircle at 0,0 */
930  B1.x = -1; B1.y = 0;
931  B2.x = 0 ; B2.y = 1;
932  B3.x = 1 ; B3.y = 0;
933 
934  /* Arc inside the semicircle */
935  lw_dist2d_distpts_init(&dl, DIST_MIN);
936  A1.x = 0; A1.y = .5;
937  A2.x = -0.3; A2.y = .5;
938  A3.x = 0; A3.y = .6;
939  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
940  printf("%g\n", dl.distance);
941  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
942  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.271798, 0.000001);
943 
944  /* Arc inside the semicircle */
945  lw_dist2d_distpts_init(&dl, DIST_MIN);
946  A1.x = -0.5; A1.y = .5;
947  A2.x = -0.4; A2.y = .2;
948  A3.x = 0; A3.y = 0;
949  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
950  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
951  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.292893, 0.000001);
952 
953  /* Arc above the unit semicircle */
954  lw_dist2d_distpts_init(&dl, DIST_MIN);
955  A1.x = -1; A1.y = 3;
956  A2.x = 0 ; A2.y = 2;
957  A3.x = 1 ; A3.y = 3;
958  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
959  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
960  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1, 0.000001);
961 
962  /* Arc grazes the unit semicircle */
963  lw_dist2d_distpts_init(&dl, DIST_MIN);
964  A1.x = -1; A1.y = 2;
965  A2.x = 0 ; A2.y = 1;
966  A3.x = 1 ; A3.y = 2;
967  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
968  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
969  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0, 0.000001);
970 
971  /* Circles intersect, but arcs do not */
972  lw_dist2d_distpts_init(&dl, DIST_MIN);
973  A1.x = -1; A1.y = 1;
974  A2.x = 0; A2.y = 2;
975  A3.x = 1; A3.y = 1;
976  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
977  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
978  CU_ASSERT_DOUBLE_EQUAL(dl.distance, sqrt(2)-1, 0.000001);
979 
980  /* Circles and arcs intersect */
981  lw_dist2d_distpts_init(&dl, DIST_MIN);
982  A1.x = -1; A1.y = 1;
983  A2.x = 0; A2.y = 0;
984  A3.x = 1; A3.y = 1;
985  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
986  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
987  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0, 0.000001);
988 
989  /* Closed circle */
990  lw_dist2d_distpts_init(&dl, DIST_MIN);
991  A1.x = -2.0; A1.y = -0.1;
992  A2.x = 1.5; A2.y = -0.1;
993  A3.x = -2.0; A3.y = -0.1;
994  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
995  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
996  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.480742, 0.0001);
997 
998  /* Concentric: and fully parallel */
999  lw_dist2d_distpts_init(&dl, DIST_MIN);
1000  A1.x = -2.0; A1.y = 0.0;
1001  A2.x = 0.0; A2.y = 2.0;
1002  A3.x = 2.0; A3.y = 0.0;
1003  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1004  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1005  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1006 
1007  /* Concentric: with A fully included in B's range */
1008  lw_dist2d_distpts_init(&dl, DIST_MIN);
1009  A1.x = -0.5 / sqrt(2.0); A1.y = 0.5 / sqrt(2.0);
1010  A2.x = 0.0; A2.y = 0.5;
1011  A3.x = 0.5 / sqrt(2.0); A3.y = 0.5 / sqrt(2.0);
1012  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1013  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1014  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);
1015 
1016  /* Concentric: with A partially included in B's range */
1017  lw_dist2d_distpts_init(&dl, DIST_MIN);
1018  A1.x = -0.5 / sqrt(2.0); A1.y = -0.5 / sqrt(2.0);
1019  A2.x = -0.5; A2.y = 0.0;
1020  A3.x = -0.5 / sqrt(2.0); A3.y = 0.5 / sqrt(2.0);
1021  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1022  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1023  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);
1024 
1025  /* Concentric: with A and B without parallel segments */
1026  lw_dist2d_distpts_init(&dl, DIST_MIN);
1027  A1.x = -0.5 / sqrt(2.0); A1.y = -0.5 / sqrt(2.0);
1028  A2.x = 0.0; A2.y = -0.5;
1029  A3.x = 0.5 / sqrt(2.0); A3.y = -0.5 / sqrt(2.0);
1030  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1031  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1032  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.736813, 0.000001);
1033 
1034  /* Concentric: Arcs are the same */
1035  lw_dist2d_distpts_init(&dl, DIST_MIN);
1036  A1.x = -1.0; A1.y = 0.0;
1037  A2.x = 0.0; A2.y = 1.0;
1038  A3.x = 1.0; A3.y = 0.0;
1039  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1040  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1041  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.0, 0.000001);
1042 
1043  /* Check different orientations */
1044  B1.x = -10.0; B1.y = 0.0;
1045  B2.x = 0.0 ; B2.y = 10.0;
1046  B3.x = 10.0 ; B3.y = 0.0;
1047 
1048  lw_dist2d_distpts_init(&dl, DIST_MIN);
1049  A1.x = -22.0; A1.y = 0.0;
1050  A2.x = -17.0; A2.y = -5.0;
1051  A3.x = -12.0; A3.y = 0.0;
1052  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1053  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1054  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 2.0, 0.000001);
1055 
1056  lw_dist2d_distpts_init(&dl, DIST_MIN);
1057  A1.x = -19.0; A1.y = 0.0;
1058  A2.x = -14.0; A2.y = -5.0;
1059  A3.x = - 9.0; A3.y = 0.0;
1060  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1061  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1062  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1063 
1064  lw_dist2d_distpts_init(&dl, DIST_MIN);
1065  A1.x = -9.0; A1.y = 0.0;
1066  A2.x = -4.0; A2.y = -5.0;
1067  A3.x = 1.0; A3.y = 0.0;
1068  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1069  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1070  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1071 
1072  lw_dist2d_distpts_init(&dl, DIST_MIN);
1073  A1.x = -1.0; A1.y = 0.0;
1074  A2.x = 4.0; A2.y = -5.0;
1075  A3.x = 9.0; A3.y = 0.0;
1076  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1077  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1078  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1079 
1080  lw_dist2d_distpts_init(&dl, DIST_MIN);
1081  A1.x = 1.0; A1.y = 0.0;
1082  A2.x = 6.0; A2.y = -5.0;
1083  A3.x = 11.0; A3.y = 0.0;
1084  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1085  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1086  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1087 
1088  lw_dist2d_distpts_init(&dl, DIST_MIN);
1089  A1.x = 11.0; A1.y = 0.0;
1090  A2.x = 16.0; A2.y = -5.0;
1091  A3.x = 21.0; A3.y = 0.0;
1092  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1093  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1094  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1095 
1096 
1097  lw_dist2d_distpts_init(&dl, DIST_MIN);
1098  A1.x = -15.0; A1.y = -6.0;
1099  A2.x = -10.0; A2.y = -1.0;
1100  A3.x = - 5.0; A3.y = -6.0;
1101  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1102  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1103  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1104 
1105  lw_dist2d_distpts_init(&dl, DIST_MIN);
1106  A1.x = -5.0; A1.y = 0.0;
1107  A2.x = 0.0; A2.y = 5.0;
1108  A3.x = 5.0; A3.y = 0.0;
1109  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1110  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1111  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 5.0, 0.000001);
1112 
1113  lw_dist2d_distpts_init(&dl, DIST_MIN);
1114  A1.x = -5.0; A1.y = 0.0;
1115  A2.x = 0.0; A2.y = -5.0;
1116  A3.x = 5.0; A3.y = 0.0;
1117  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1118  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1119  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 5.0, 0.000001);
1120 
1121 }
LW_SUCCESS
#define LW_SUCCESS
Definition: liblwgeom.h:97
lw_dist2d_distpts_init
void lw_dist2d_distpts_init(DISTPTS *dl, int mode)
Definition: measures.c:67
lw_dist2d_arc_arc
int lw_dist2d_arc_arc(const POINT2D *A1, const POINT2D *A2, const POINT2D *A3, const POINT2D *B1, const POINT2D *B2, const POINT2D *B3, DISTPTS *dl)
Definition: measures.c:1677
DIST_MIN
#define DIST_MIN
Definition: measures.h:44
DISTPTS::distance
double distance
Definition: measures.h:51
DISTPTS
Structure used in distance-calculations.
Definition: measures.h:50
POINT2D::y
double y
Definition: liblwgeom.h:390
POINT2D::x
double x
Definition: liblwgeom.h:390

References DIST_MIN, DISTPTS::distance, lw_dist2d_arc_arc(), lw_dist2d_distpts_init(), LW_SUCCESS, POINT2D::x, and POINT2D::y.

Referenced by measures_suite_setup().

Here is the call graph for this function:
Here is the caller graph for this function: