From 84d01c3ac4bae4ddae787c298f49ac96957d372a Mon Sep 17 00:00:00 2001 From: rodri Date: Sat, 19 Oct 2024 11:43:22 +0000 Subject: med: remove ptin(cone|cylinder) (now part of libgeometry.) --- med.c | 48 ------------------------------------------------ 1 file changed, 48 deletions(-) diff --git a/med.c b/med.c index b2ab44e..ef47356 100644 --- a/med.c +++ b/med.c @@ -122,54 +122,6 @@ Ptquat(Quaternion q, double w) return Pt3(q.i, q.j, q.k, w); } -/* - * p is the point to test - * p0 and p1 are the centers of the circles at each end of the cylinder - * r is the radius of these circles - */ -int -ptincylinder(Point3 p, Point3 p0, Point3 p1, double r) -{ - Point3 p01, p0p, p1p; - double h; - - p01 = subpt3(p1, p0); - p0p = subpt3(p, p0); - p1p = subpt3(p, p1); - h = vec3len(p01); - - if(h == 0) - return 0; - - return dotvec3(p0p, p01) >= 0 && - dotvec3(p1p, p01) <= 0 && - vec3len(crossvec3(p0p, p01))/h <= r; -} - -/* - * p is the point to test - * p0 is the apex - * p1 is the center of the base - * br is the radius of the base - */ -int -ptincone(Point3 p, Point3 p0, Point3 p1, double br) -{ - Point3 p01, p0p; - double h, d, r; - - p01 = subpt3(p1, p0); - p0p = subpt3(p, p0); - h = vec3len(p01); - d = dotvec3(p0p, normvec3(p01)); - - if(h == 0 || d < 0 || d > h) - return 0; - - r = d/h * br; - return vec3len(crossvec3(p0p, p01))/h <= r; -} - static Point randptfromrect(Rectangle *r) { -- cgit v1.2.3