1 files changed, 43 insertions, 0 deletions
diff --git a/src/Math/Plane.c b/src/Math/Plane.c
new file mode 100644
index 0000000..23711c5
--- /dev/null
+++ b/src/Math/Plane.c
@@ -0,0 +1,43 @@
+#include <BH/Math.h>
+#include <string.h>
+
+
+#define EPSILON 0.00001f
+#define PI      3.14159265358979323846f
+
+
+int BH_PlaneFromPoints(const float *a,
+                       const float *b,
+                       const float *c,
+                       float *out)
+{
+    float tmp1[3], tmp2[3];
+
+    BH_Vec3fSub(b, a, tmp1);
+    BH_Vec3fSub(c, a, tmp2);
+    BH_Vec3fCross(tmp2, tmp1, tmp1);
+    if (BH_Vec3fNormalEx(tmp1, tmp1) < EPSILON)
+        return BH_ERROR;
+
+    out[3] = BH_Vec3fDot(a, tmp1);
+    memcpy(out, tmp1, sizeof(tmp1));
+    return BH_OK;
+}
+
+
+float BH_PlaneDistance(const float *plane,
+                       const float *point)
+{
+    return BH_Vec3fDot(plane, point) - plane[3];
+}
+
+
+void BH_PlaneClosestPoint(const float *plane,
+                          const float *point,
+                          float *out)
+{
+    float tmp[3];
+
+    BH_Vec3fScale(plane, BH_PlaneDistance(plane, point), tmp);
+    BH_Vec3fSub(point, tmp, out);
+}