checkpoint i think it works!

[anneal.git] / src / edu / berkeley / qfat / Mesh.java

diff --git a/src/edu/berkeley/qfat/Mesh.java b/src/edu/berkeley/qfat/Mesh.java
index b6ff386..64724a6 100644 (file)
--- a/src/edu/berkeley/qfat/Mesh.java
+++ b/src/edu/berkeley/qfat/Mesh.java
@@ -305,11 +305,6 @@ public class Mesh implements Iterable<Mesh.T> {
                 e.p2.bindTo(Matrix.ONE, ebound.p1);
             }
         }
-        public void shatter(BindingGroup bg1, BindingGroup bg2, boolean triangles) {
-            for(E e : set) {
-                e.shatter(e.midpoint(), bg1, bg2, triangles);
-            }
-        }
     }
 
     /** [UNIQUE] an edge */
@@ -327,16 +322,20 @@ public class Mesh implements Iterable<Mesh.T> {
         public boolean intersects(T t) { return t.intersects(p1.p, p2.p); }
 
         public void bindingGroupChanged(edu.berkeley.qfat.geom.BindingGroup newBindingGroup_) {
+
             edu.berkeley.qfat.geom.BindingGroup<E> newBindingGroup =
                 (edu.berkeley.qfat.geom.BindingGroup<E>)newBindingGroup_;
             if (newBindingGroup==null) return;
             if (this==newBindingGroup.getMaster()) return;
-            /*
             for(E eother : (Iterable<E>)newBindingGroup) {
-                this.next.bindTo(newBindingGroup.getMatrix(eother), eother.next);
-                this.prev.bindTo(newBindingGroup.getMatrix(eother), eother.prev);
+                if (next==null || prev==null) continue;
+                if (eother.next==null || eother.prev==null) continue;
+                if (next.isBoundTo(eother.pair.prev.pair) && !prev.isBoundTo(eother.pair.next.pair))
+                    prev.bindTo(next.getBindingMatrix(eother.pair.prev.pair), eother.pair.next.pair);
+                if (!next.isBoundTo(eother.pair.prev.pair) && prev.isBoundTo(eother.pair.next.pair))
+                    next.bindTo(prev.getBindingMatrix(eother.pair.next.pair), eother.pair.prev.pair);
             }
-            */
+
         }
 
         public float stretchRatio() {
@@ -394,40 +393,46 @@ public class Mesh implements Iterable<Mesh.T> {
             }
         }
 
-        public Point shatter() { return shatter(true); }
-        public Point shatter(boolean triangles) { return shatter(midpoint(), null, null, triangles); }
-        public Point shatter(Point mid, BindingGroup bg1, BindingGroup bg2, boolean triangles) {
-            return shatter(mid, bg1, bg2, triangles, false);
-        }
-        public Point shatter(Point mid, BindingGroup bg1, BindingGroup bg2, boolean triangles, boolean leader) {
-            if (shattered || destroyed) return mid;
+        public Point shatter() {
+            if (shattered || destroyed) return null;
             shattered = true;
-
-            Vertex r = next.p2;
-            E next = this.next;
-            E prev = this.prev;
-
-            int old_colorclass = t==null ? 0 : t.colorclass;
-            if (bg1==null) bg1 = new BindingGroup();
-            if (bg2==null) bg2 = new BindingGroup();
-            BindingGroup old_bind_to = bind_to;
-            bind_peers.shatter(bg1, bg2, triangles);
-            old_bind_to.shatter(bg2.other(), bg1.other(), triangles);
-            if (!triangles) {
-                next.shatter(false);
-                prev.shatter(false);
+            E first = null;
+            E firste = null;
+            E firstx = null;
+            E firstq = null;
+            for(E e : (Iterable<E>)getBoundPeers()) {
+                E enext = e.next;
+                E eprev = e.prev;
+                E pnext = e.pair.next;
+                E pprev = e.pair.prev;
+                Point mid = e.midpoint();
+                Vertex r = e.next.p2;
+                Vertex l = e.pair.next.p2;
+                e.destroy();
+                e.pair.destroy();
+                newT(r.p, e.p1.p, mid,    null, 0);
+                newT(r.p, mid,    e.p2.p, null, 0);
+                newT(l.p, mid,    e.p1.p, null, 0);
+                newT(l.p, e.p2.p, mid,    null, 0);
+                /*
+                */
             }
-            pair.shatter();
-            destroy();
-
-            if (triangles) {
-                newT(r.p, p1.p, mid, null, old_colorclass);
-                newT(r.p, mid, p2.p, null, old_colorclass);
-                bg1.add(p1.getE(mid));
-                bg2.add(p2.getE(mid).pair);
-                if (leader) p1.getE(mid).shatter();
+            for(E e : (Iterable<E>)getBoundPeers()) {
+                Point mid = e.midpoint();
+                if (first==null) {
+                    first = e.p1.getE(mid);
+                    firste = e;
+                    firstx = e.pair;
+                    firstq = e.p2.getE(mid).pair;
+                    continue;
+                }
+                e.p1.getE(mid).          bindTo(e.getBindingMatrix(firste), first);
+                e.p1.getE(mid).pair.     bindTo(e.getBindingMatrix(firste), first.pair);
+                e.p2.getE(mid).pair.     bindTo(e.getBindingMatrix(firste), firstq);
+                e.p2.getE(mid).pair.pair.bindTo(e.getBindingMatrix(firste), firstq.pair);
+                //e.p1.getE(mid).pair.prev.bindTo(e.pair.getBindingMatrix(firstx), first.pair.prev);
             }
-            return mid;
+            return null;
         }
 
         public boolean destroyed = false;
@@ -646,14 +651,9 @@ public class Mesh implements Iterable<Mesh.T> {
         public void reinsert() { triangles.remove(this); triangles.add(this); }
 
         public boolean shouldBeDrawn() {
-            /*
-            if (e1().bind_to==null) return false;
-            if (e2().bind_to==null) return false;
-            if (e3().bind_to==null) return false;
-            if (e1().bind_to.set.size() == 0) return false;
-            if (e2().bind_to.set.size() == 0) return false;
-            if (e3().bind_to.set.size() == 0) return false;
-            */
+            if (e1().bindingGroupSize() <= 1) return false;
+            if (e2().bindingGroupSize() <= 1) return false;
+            if (e3().bindingGroupSize() <= 1) return false;
             return true;
         }