PostGIS  2.5.0dev-r@@SVN_REVISION@@

◆ lw_dist3d_pt_pt()

int lw_dist3d_pt_pt ( POINT3DZ p1,
POINT3DZ p2,
DISTPTS3D dl 
)

Compares incoming points and stores the points closest to each other or most far away from each other depending on dl->mode (max or min)

Definition at line 809 of file measures3d.c.

References DISTPTS3D::distance, LW_TRUE, LWDEBUGF, DISTPTS3D::mode, DISTPTS3D::p1, DISTPTS3D::p2, DISTPTS3D::twisted, POINT3DZ::x, POINT3DZ::y, and POINT3DZ::z.

Referenced by lw_dist3d_point_point(), lw_dist3d_pt_poly(), lw_dist3d_pt_seg(), lw_dist3d_ptarray_ptarray(), and lw_dist3d_seg_seg().

810 {
811  double dx = thep2->x - thep1->x;
812  double dy = thep2->y - thep1->y;
813  double dz = thep2->z - thep1->z;
814  double dist = sqrt ( dx*dx + dy*dy + dz*dz);
815  LWDEBUGF(2, "lw_dist3d_pt_pt called (with points: p1.x=%f, p1.y=%f,p1.z=%f,p2.x=%f, p2.y=%f,p2.z=%f)",thep1->x,thep1->y,thep1->z,thep2->x,thep2->y,thep2->z );
816 
817  if (((dl->distance - dist)*(dl->mode))>0) /*multiplication with mode to handle mindistance (mode=1) and maxdistance (mode = (-1)*/
818  {
819  dl->distance = dist;
820 
821  if (dl->twisted>0) /*To get the points in right order. twisted is updated between 1 and (-1) every time the order is changed earlier in the chain*/
822  {
823  dl->p1 = *thep1;
824  dl->p2 = *thep2;
825  }
826  else
827  {
828  dl->p1 = *thep2;
829  dl->p2 = *thep1;
830  }
831  }
832  return LW_TRUE;
833 }
double distance
Definition: measures3d.h:41
POINT3DZ p2
Definition: measures3d.h:43
POINT3DZ p1
Definition: measures3d.h:42
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:75
int twisted
Definition: measures3d.h:45
int mode
Definition: measures3d.h:44
#define LWDEBUGF(level, msg,...)
Definition: lwgeom_log.h:88
Here is the caller graph for this function: