PostGIS  3.1.6dev-r@@SVN_REVISION@@

◆ test_lw_dist2d_arc_arc()

static void test_lw_dist2d_arc_arc ( void  )
static

Definition at line 894 of file cu_measures.c.

895 {
896  /* lw_dist2d_arc_arc(const POINT2D *A1, const POINT2D *A2, const POINT2D *A3,
897  const POINT2D *B1, const POINT2D *B2, const POINT2D *B3,
898  DISTPTS *dl) */
899  DISTPTS dl;
900  POINT2D A1, A2, A3, B1, B2, B3;
901  int rv;
902 
903  /* Ticket #4326 */
904  lw_dist2d_distpts_init(&dl, DIST_MIN);
905  A1.x = -1.0; A1.y = 4.0;
906  A2.x = 0.0; A2.y = 5.0;
907  A3.x = 1.0; A3.y = 4.0;
908  B1.x = 1.0; B1.y = 6.0;
909  B2.x = 6.0; B2.y = 1.0;
910  B3.x = 9.0; B3.y = 7.0;
911  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
912  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
913  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.0475666, 0.000001);
914 
915  /* Unit semicircle at 0,0 */
916  B1.x = -1; B1.y = 0;
917  B2.x = 0 ; B2.y = 1;
918  B3.x = 1 ; B3.y = 0;
919 
920  /* Arc above the unit semicircle */
921  lw_dist2d_distpts_init(&dl, DIST_MIN);
922  A1.x = -1; A1.y = 3;
923  A2.x = 0 ; A2.y = 2;
924  A3.x = 1 ; A3.y = 3;
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, 1, 0.000001);
928 
929  /* Arc grazes the unit semicircle */
930  lw_dist2d_distpts_init(&dl, DIST_MIN);
931  A1.x = -1; A1.y = 2;
932  A2.x = 0 ; A2.y = 1;
933  A3.x = 1 ; A3.y = 2;
934  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
935  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
936  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0, 0.000001);
937 
938  /* Circles intersect, but arcs do not */
939  lw_dist2d_distpts_init(&dl, DIST_MIN);
940  A1.x = -1; A1.y = 1;
941  A2.x = 0; A2.y = 2;
942  A3.x = 1; A3.y = 1;
943  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
944  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
945  CU_ASSERT_DOUBLE_EQUAL(dl.distance, sqrt(2)-1, 0.000001);
946 
947  /* Circles and arcs intersect */
948  lw_dist2d_distpts_init(&dl, DIST_MIN);
949  A1.x = -1; A1.y = 1;
950  A2.x = 0; A2.y = 0;
951  A3.x = 1; A3.y = 1;
952  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
953  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
954  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0, 0.000001);
955 
956  /* Concentric: and fully parallel */
957  lw_dist2d_distpts_init(&dl, DIST_MIN);
958  A1.x = -2.0; A1.y = 0.0;
959  A2.x = 0.0; A2.y = 2.0;
960  A3.x = 2.0; A3.y = 0.0;
961  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
962  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
963  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
964 
965  /* Concentric: with A fully included in B's range */
966  lw_dist2d_distpts_init(&dl, DIST_MIN);
967  A1.x = -0.5 / sqrt(2.0); A1.y = 0.5 / sqrt(2.0);
968  A2.x = 0.0; A2.y = 0.5;
969  A3.x = 0.5 / sqrt(2.0); A3.y = 0.5 / sqrt(2.0);
970  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
971  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
972  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);
973 
974  /* Concentric: with A partially included in B's range */
975  lw_dist2d_distpts_init(&dl, DIST_MIN);
976  A1.x = -0.5 / sqrt(2.0); A1.y = -0.5 / sqrt(2.0);
977  A2.x = -0.5; A2.y = 0.0;
978  A3.x = -0.5 / sqrt(2.0); A3.y = 0.5 / sqrt(2.0);
979  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
980  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
981  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);
982 
983  /* Concentric: with A and B without parallel segments */
984  lw_dist2d_distpts_init(&dl, DIST_MIN);
985  A1.x = -0.5 / sqrt(2.0); A1.y = -0.5 / sqrt(2.0);
986  A2.x = 0.0; A2.y = -0.5;
987  A3.x = 0.5 / sqrt(2.0); A3.y = -0.5 / sqrt(2.0);
988  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
989  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
990  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.736813, 0.000001);
991 
992  /* Concentric: Arcs are the same */
993  lw_dist2d_distpts_init(&dl, DIST_MIN);
994  A1.x = -1.0; A1.y = 0.0;
995  A2.x = 0.0; A2.y = 1.0;
996  A3.x = 1.0; A3.y = 0.0;
997  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
998  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
999  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.0, 0.000001);
1000 
1001  /* Check different orientations */
1002  B1.x = -10.0; B1.y = 0.0;
1003  B2.x = 0.0 ; B2.y = 10.0;
1004  B3.x = 10.0 ; B3.y = 0.0;
1005 
1006  lw_dist2d_distpts_init(&dl, DIST_MIN);
1007  A1.x = -22.0; A1.y = 0.0;
1008  A2.x = -17.0; A2.y = -5.0;
1009  A3.x = -12.0; A3.y = 0.0;
1010  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1011  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1012  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 2.0, 0.000001);
1013 
1014  lw_dist2d_distpts_init(&dl, DIST_MIN);
1015  A1.x = -19.0; A1.y = 0.0;
1016  A2.x = -14.0; A2.y = -5.0;
1017  A3.x = - 9.0; A3.y = 0.0;
1018  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1019  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1020  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1021 
1022  lw_dist2d_distpts_init(&dl, DIST_MIN);
1023  A1.x = -9.0; A1.y = 0.0;
1024  A2.x = -4.0; A2.y = -5.0;
1025  A3.x = 1.0; A3.y = 0.0;
1026  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1027  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1028  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1029 
1030  lw_dist2d_distpts_init(&dl, DIST_MIN);
1031  A1.x = -1.0; A1.y = 0.0;
1032  A2.x = 4.0; A2.y = -5.0;
1033  A3.x = 9.0; A3.y = 0.0;
1034  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1035  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1036  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1037 
1038  lw_dist2d_distpts_init(&dl, DIST_MIN);
1039  A1.x = 1.0; A1.y = 0.0;
1040  A2.x = 6.0; A2.y = -5.0;
1041  A3.x = 11.0; A3.y = 0.0;
1042  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1043  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1044  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1045 
1046  lw_dist2d_distpts_init(&dl, DIST_MIN);
1047  A1.x = 11.0; A1.y = 0.0;
1048  A2.x = 16.0; A2.y = -5.0;
1049  A3.x = 21.0; A3.y = 0.0;
1050  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1051  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1052  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1053 
1054 
1055  lw_dist2d_distpts_init(&dl, DIST_MIN);
1056  A1.x = -15.0; A1.y = -6.0;
1057  A2.x = -10.0; A2.y = -1.0;
1058  A3.x = - 5.0; A3.y = -6.0;
1059  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1060  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1061  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1.0, 0.000001);
1062 
1063  lw_dist2d_distpts_init(&dl, DIST_MIN);
1064  A1.x = -5.0; A1.y = 0.0;
1065  A2.x = 0.0; A2.y = 5.0;
1066  A3.x = 5.0; A3.y = 0.0;
1067  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1068  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1069  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 5.0, 0.000001);
1070 
1071  lw_dist2d_distpts_init(&dl, DIST_MIN);
1072  A1.x = -5.0; A1.y = 0.0;
1073  A2.x = 0.0; A2.y = -5.0;
1074  A3.x = 5.0; A3.y = 0.0;
1075  rv = lw_dist2d_arc_arc(&A1, &A2, &A3, &B1, &B2, &B3, &dl);
1076  CU_ASSERT_EQUAL( rv, LW_SUCCESS );
1077  CU_ASSERT_DOUBLE_EQUAL(dl.distance, 5.0, 0.000001);
1078 
1079 }
LW_SUCCESS
#define LW_SUCCESS
Definition: liblwgeom.h:111
lw_dist2d_distpts_init
void lw_dist2d_distpts_init(DISTPTS *dl, int mode)
Definition: measures.c:64
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:1575
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:404
POINT2D::x
double x
Definition: liblwgeom.h:404

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: