PostGIS  2.1.10dev-r@@SVN_REVISION@@
int lw_dist3d_pt_pt ( POINT3DZ p1,
POINT3DZ p2,
DISTPTS3D dl 
)

Compares incomming 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 598 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().

599 {
600  double dx = thep2->x - thep1->x;
601  double dy = thep2->y - thep1->y;
602  double dz = thep2->z - thep1->z;
603  double dist = sqrt ( dx*dx + dy*dy + dz*dz);
604  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 );
605 
606  if (((dl->distance - dist)*(dl->mode))>0) /*multiplication with mode to handle mindistance (mode=1) and maxdistance (mode = (-1)*/
607  {
608  dl->distance = dist;
609 
610  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*/
611  {
612  dl->p1 = *thep1;
613  dl->p2 = *thep2;
614  }
615  else
616  {
617  dl->p1 = *thep2;
618  dl->p2 = *thep1;
619  }
620  }
621  return LW_TRUE;
622 }
double distance
Definition: measures3d.h:26
POINT3DZ p2
Definition: measures3d.h:28
POINT3DZ p1
Definition: measures3d.h:27
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:51
int twisted
Definition: measures3d.h:30
int mode
Definition: measures3d.h:29
#define LWDEBUGF(level, msg,...)
Definition: lwgeom_log.h:55

Here is the caller graph for this function: