public class Main extends MeshViewer {
- public static int verts = 0;
+ public static int verts = 1;
public static final Random random = new Random();
goal.ignorecollision = true;
// rotate to align major axis -- this probably needs to be done by a human.
- goal.transform(new Matrix(new Vec(0, 0, 1), (float)(Math.PI/2)));
+ goal.transform(Matrix.rotate(new Vec(0, 0, 1), (float)(Math.PI/2)));
float goal_width = goal.diagonal().dot(new Vec(1, 0, 0));
float goal_height = goal.diagonal().dot(new Vec(0, 1, 0));
float halfup = 0;
translations = new Matrix[] {
- /*
- new Matrix(new Vec(lshift, depth, halfup)),
- new Matrix(new Vec(rshift, depth, halfup)),
- new Matrix(new Vec(lshift, -depth, halfup)),
- new Matrix(new Vec(rshift, -depth, halfup)),
- */
- new Matrix(new Vec(0, depth, halfup)),
- new Matrix(new Vec(0, -depth, halfup)),
+ Matrix.translate(new Vec(lshift, depth, halfup)),
+ Matrix.translate(new Vec(rshift, depth, halfup)),
+ Matrix.translate(new Vec(lshift, -depth, halfup)),
+ Matrix.translate(new Vec(rshift, -depth, halfup)),
+ /*
+ Matrix.translate(new Vec(0, depth, halfup)),
+ Matrix.translate(new Vec(0, -depth, halfup)),
+ */
- new Matrix(new Vec(lshift, 0, height)),
- new Matrix(new Vec(rshift, 0, height)),
- new Matrix(new Vec(lshift, 0, -height)),
- new Matrix(new Vec(rshift, 0, -height)),
+ Matrix.translate(new Vec(lshift, 0, height)),
+ Matrix.translate(new Vec(rshift, 0, height)),
+ Matrix.translate(new Vec(lshift, 0, -height)),
+ Matrix.translate(new Vec(rshift, 0, -height)),
- new Matrix(new Vec( width, 0, 0)),
- new Matrix(new Vec(-width, 0, 0)),
+ Matrix.translate(new Vec( width, 0, 0)),
+ Matrix.translate(new Vec(-width, 0, 0)),
};
}
}
- //xMesh.Vert mid = lbf.getE(mbn).shatter();
+ //xMesh.Vertex mid = lbf.getE(mbn).shatter();
// rescale to match volume
float factor = (float)Math.pow(tile.volume() / goal.volume(), 1.0/3.0);
- goal.transform(new Matrix(factor));
+ goal.transform(Matrix.scale(factor));
// translate to match centroid
- goal.transform(new Matrix(tile.centroid().minus(goal.centroid())));
+ goal.transform(Matrix.translate(tile.centroid().minus(goal.centroid())));
+ goal.makeVerticesImmutable();
//tx.e2.shatter();
//tx.e3.shatter();
tile.score_against = goal;
goal.score_against = tile;
- tile.tilemesh = true;
}
public synchronized void breakit() {
- if (verts > 800) return;
- //while(verts < 800) {
+ int oldverts = verts;
+ System.out.println("doubling vertices.");
PriorityQueue<Mesh.E> es = new PriorityQueue<Mesh.E>();
- for(Mesh.E e : tile.edges()) es.add(e);
- for(int i=0; i<10; i++) {
+ for(Mesh.T t : tile) {
+ es.add(t.e1());
+ es.add(t.e2());
+ es.add(t.e3());
+ }
+ for(int i=0; i<oldverts; i++) {
Mesh.E e = es.poll();
verts++;
- System.out.println("shatter " + e);
+ //System.out.println("shatter " + e);
e.shatter();
}
tile.rebindPoints();
- //}
}
- public synchronized void rand(float temp, Mesh.Vert p) {
+ public synchronized void rand(float temp, Mesh.Vertex p) {
//p.reComputeError();
p.reComputeErrorAround();
float hits = 0;
float misses = 0;
public void anneal() throws Exception {
- int verts = 0;
float hightemp = 10;
float temp = hightemp;
float last = 10;
while(true) {
synchronized(this) {
double ratio = (hits+misses==0) ? 1 : (hits / (hits+misses));
- System.out.println("temp="+temp + " ratio="+(Math.ceil(ratio*100)));
hits = 0;
misses = 0;
float gamma = 0;
accepts = (int)(Math.ceil(ratio*100));
temps = (int)(Math.ceil(temp*1000));
vertss = tile.size();
- if (breaks) {
- breaks = false;
+ if (breaks > 0) { while (breaks>0) {
+ breaks--;
breakit();
//gamma = 1;
gamma = 1;
//temp = last * 0.8f;
//last = temp;
//temp = hightemp;
- } else
+ } } else
if (acceptance > 0.96) gamma = 0.4f;
else if (acceptance > 0.9) gamma = 0.5f;
else if (acceptance > 0.8) gamma = 0.65f;
gamma = 0.9f;
} else if (acceptance > 0.15) {
gamma = 0.95f;
+ } else if (acceptance > 0.10) {
+ gamma = 0.98f;
} else {
breakit();
//gamma = 1;
temp = temp * gamma;
- HashSet<Mesh.Vert> hs = new HashSet<Mesh.Vert>();
- for(Mesh.Vert p : tile.vertices()) hs.add(p);
- for(int i=0; i<10; i++) {
- repaint();
- for(Mesh.Vert v : hs) {
- if (anneal) rand(temp,v);
- Thread.yield();
- repaint();
+ HashSet<Mesh.Vertex> hs = new HashSet<Mesh.Vertex>();
+ for(Mesh.Vertex p : tile.vertices()) hs.add(p);
+ Mesh.Vertex[] pts = (Mesh.Vertex[])hs.toArray(new Mesh.Vertex[0]);
+
+ int count = 0;
+ long then = System.currentTimeMillis();
+ for(int i=0; i<40; i++) {
+ if (anneal) {
+ count++;
+ Mesh.Vertex v = pts[Math.abs(random.nextInt()) % pts.length];
+ rand(temp,v);
}
+ Thread.yield();
+ repaint();
}
- tile.rebuildPointSet();
- repaint();
- //breakit();
- repaint();
- goal.unApplyQuadricToNeighborAll();
+ System.out.println("temp="+temp + " ratio="+(Math.ceil(ratio*100)) + " " +
+ "points_per_second=" +
+ (count*1000)/((double)(System.currentTimeMillis()-then)));
+
+ for(Mesh.Vertex p : goal.vertices()) p.recomputeFundamentalQuadricIfNeighborChanged();
repaint();
- tile.recomputeAllFundamentalQuadrics();
+ for(Mesh.Vertex p : tile.vertices()) p.recomputeFundamentalQuadric();
repaint();
- goal.applyQuadricToNeighborAll();
+
+ synchronized(safeTriangles) {
+ safeTriangles.clear();
+ for(Mesh.T t : tile) if (t.shouldBeDrawn()) safeTriangles.add(t);
+ }
}
}
}
+
public static void main(String[] s) throws Exception {
StlFile stlf = new StlFile();
stlf.load("fish.stl");
projects / anneal.git / blobdiff