Implement missing distance3d_pt_pt function

git-svn-id: http://svn.osgeo.org/postgis/trunk@13627 b70326c6-7e19-0410-871a-916f4a2858ee

diff --git a/liblwgeom/measures3d.c b/liblwgeom/measures3d.c
index a14205e4ee975385e5f9da8013c3e47b387b482b..f7b022a73a1f8ac146e5b20488e89adbad392a71 100644 (file)
--- a/liblwgeom/measures3d.c
+++ b/liblwgeom/measures3d.c
@@ -620,6 +620,15 @@ lw_dist3d_pt_seg(POINT3DZ *p, POINT3DZ *A, POINT3DZ *B, DISTPTS3D *dl)
        return lw_dist3d_pt_pt(p,&c,dl);
 }
 
+double
+distance3d_pt_pt(const POINT3D *p1, const POINT3D *p2)
+{
+  double dx = p2->x - p1->x;
+  double dy = p2->y - p1->y;
+  double dz = p2->z - p1->z;
+  return sqrt ( dx*dx + dy*dy + dz*dz);
+}
+
 
 /**