checkpoint

[anneal.git] / src / edu / berkeley / qfat / geom / Plane.java

diff --git a/src/edu/berkeley/qfat/geom/Plane.java b/src/edu/berkeley/qfat/geom/Plane.java
index 1f79b06..5f5d990 100644 (file)
--- a/src/edu/berkeley/qfat/geom/Plane.java
+++ b/src/edu/berkeley/qfat/geom/Plane.java
@@ -1,7 +1,10 @@
 package edu.berkeley.qfat.geom;
 import javax.media.opengl.*;
 
-public class Plane {
+public class Plane implements AffineConstraint {
+
+    // FIXME: could be given by
+    // ax+by+cz=d
 
     public final Vec norm;
     public final float dvalue;
@@ -11,14 +14,26 @@ public class Plane {
         this.dvalue = p.x*this.norm.x+p.y*this.norm.y+p.z*this.norm.z;
     }
 
+    /** provided at least one of a,b,c is nonzero, return the Plane representing ax+by+cz=d */
+    public Plane(float a, float b, float c, float d) {
+        throw new RuntimeException("not implemented yet");
+    }
+
     public Point intersect(Plane p1, Plane p2) {
         Plane p3 = this;
         float z = p1.norm.dot(p2.norm.cross(p3.norm));
-        if (Math.abs(z) == 0) return null;  // planes do not intersect at a point
+        if (Math.abs(z) < 0.0001) return null;  // planes do not intersect at a point
         Vec v1 = p2.norm.cross(p3.norm).times(-1 * p1.dvalue);
         Vec v2 = p3.norm.cross(p1.norm).times(-1 * p2.dvalue);
         Vec v3 = p1.norm.cross(p2.norm).times(-1 * p3.dvalue);
-        return Point.ORIGIN.plus(v1.plus(v2).plus(v3).times(1/z));
+        return Point.ZERO.plus(v1.plus(v2).plus(v3).times(1/z));
     }
 
+    public Point getProjection(Point p) {
+        throw new RuntimeException("not implemented yet");
+    }
+
+    public AffineConstraint intersect(AffineConstraint c, float epsilon) {
+        throw new RuntimeException("not implemented yet");
+    }
 }