PostGIS  3.0.6dev-r@@SVN_REVISION@@

◆ lw_dist3d_tri_tri()

int lw_dist3d_tri_tri ( LWTRIANGLE tri1,
LWTRIANGLE tri2,
DISTPTS3D dl 
)

triangle to triangle calculation

Definition at line 918 of file measures3d.c.

919 {
920  PLANE3D plane1, plane2;
921  int planedef1, planedef2;
922 
923  if (dl->mode == DIST_MAX)
924  return lw_dist3d_ptarray_ptarray(tri1->points, tri2->points, dl);
925 
926  planedef1 = define_plane(tri1->points, &plane1);
927  planedef2 = define_plane(tri2->points, &plane2);
928 
929  if (!planedef1 || !planedef2)
930  {
931  /* Neither define a plane: Return distance line to line */
932  if (!planedef1 && !planedef2)
933  return lw_dist3d_ptarray_ptarray(tri1->points, tri2->points, dl);
934 
935  /* Only tri defines a plane: Return distance from line (tri1) to tri2 */
936  else if (!planedef1)
937  return lw_dist3d_ptarray_tri(tri1->points, tri2, &plane2, dl);
938 
939  /* Only poly defines a plane: Return distance from line (tri2) to tri1 */
940  else
941  return lw_dist3d_ptarray_tri(tri2->points, tri1, &plane1, dl);
942  }
943 
944  /* What we do here is to compare the boundary of one polygon with the other polygon
945  and then take the second boundary comparing with the first polygon */
946  dl->twisted = 1;
947  if (!lw_dist3d_ptarray_tri(tri1->points, tri2, &plane2, dl))
948  return LW_FALSE;
949  if (dl->distance < dl->tolerance) /* Just check if the answer already is given*/
950  return LW_TRUE;
951 
952  dl->twisted = -1; /* because we switch the order of geometries we switch "twisted" to -1 which will give the
953  right order of points in shortest line. */
954  return lw_dist3d_ptarray_tri(tri2->points, tri1, &plane1, dl);
955 }
#define LW_FALSE
Definition: liblwgeom.h:108
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:107
int define_plane(POINTARRAY *pa, PLANE3D *pl)
Definition: measures3d.c:1451
int lw_dist3d_ptarray_tri(POINTARRAY *pa, LWTRIANGLE *tri, PLANE3D *plane, DISTPTS3D *dl)
Computes pointarray to triangle distance.
Definition: measures3d.c:1373
int lw_dist3d_ptarray_ptarray(POINTARRAY *l1, POINTARRAY *l2, DISTPTS3D *dl)
Finds all combinations of segments between two pointarrays.
Definition: measures3d.c:1088
#define DIST_MAX
Definition: measures.h:43
int twisted
Definition: measures3d.h:45
double distance
Definition: measures3d.h:40
int mode
Definition: measures3d.h:43
double tolerance
Definition: measures3d.h:47
POINTARRAY * points
Definition: liblwgeom.h:481

References define_plane(), DIST_MAX, DISTPTS3D::distance, lw_dist3d_ptarray_ptarray(), lw_dist3d_ptarray_tri(), LW_FALSE, LW_TRUE, DISTPTS3D::mode, LWTRIANGLE::points, DISTPTS3D::tolerance, and DISTPTS3D::twisted.

Referenced by lw_dist3d_distribute_bruteforce().

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