From f5b2d5436a16166410d9e58c68b9bf5d6f8db389 Mon Sep 17 00:00:00 2001 From: Francois Fleuret Date: Fri, 10 Jul 2009 16:25:03 +0200 Subject: [PATCH] *** empty log message *** --- mash.cc | 82 ++++++++++++++++++++++++++++++++++++++++++++++----------- 1 file changed, 67 insertions(+), 15 deletions(-) diff --git a/mash.cc b/mash.cc index 6eafd6c..ed85d2c 100644 --- a/mash.cc +++ b/mash.cc @@ -15,7 +15,7 @@ using namespace std; #include "universe.h" int main(int argc, char **argv) { - int nb_heur = 100; + int nb_heur = 75; int w = 1000, h = 1000; Universe universe(nb_heur + 1, w, h); @@ -42,31 +42,83 @@ int main(int argc, char **argv) { universe.add_polygon(p); } - scalar_t xs[] = { -30, 30, 30, -30 }; - scalar_t ys[] = { -10, -10, 10, 10 }; - Polygon *blocks[nb_heur]; - for(int h = 0; h < nb_heur; h++) { - blocks[h] = new Polygon(1.0, 1.0, 1.0, 0.0, xs, ys, 4); + for(int i = 0; i < nb_heur; i++) { + + int nb_max_edges = 10; + int xt[nb_max_edges * nb_max_edges], yt[nb_max_edges * nb_max_edges]; + int used[nb_max_edges * nb_max_edges]; + int fail, x, y, n; do { - blocks[h]->set_position(drand48() * w, drand48() * h/10, M_PI * drand48()); - blocks[h]->set_speed(0, 0, 0); - } while(universe.collide(blocks[h])); - universe.initialize(blocks[h]); - universe.add_polygon(blocks[h]); + fail = 0; + for(int k = 0; k < nb_max_edges * nb_max_edges; k++) { + used[k] = 0; + } + x = 0; y = 0; n = 0; + while(!used[x + nb_max_edges * y]) { + used[x + nb_max_edges * y] = 1; + xt[n] = x; yt[n] = y; + n++; + switch(int(drand48() * 4)) { + case 0: x++; break; + case 1: x--; break; + case 2: y++; break; + case 3: y--; break; + default: abort(); + } + if(x < 0) x = 0; + else if(x >= nb_max_edges) x = nb_max_edges-1; + if(y < 0) y = 0; + else if(y >= nb_max_edges) y = nb_max_edges-1; + } + fail = (n < 8) || (x != 0) || (y != 0); + } while(fail); + + scalar_t xs[nb_max_edges], ys[nb_max_edges]; + int nb_edges = 0; + for(int l = 0; l < n; l++) { + int pl = (l + n - 1)%n, nl = (l + 1)%n; + if(xt[nl] - xt[l] != xt[l] - xt[pl] || yt[nl] - yt[l] != yt[l] - yt[pl]) { + xs[nb_edges] = scalar_t(xt[l]) * 15; + ys[nb_edges] = scalar_t(yt[l]) * 15; + nb_edges++; + } + } + + cout << "n = " << n << " nb_edges = " << nb_edges << endl; + + blocks[i] = new Polygon(1.0, 1.0, 1.0, 0.0, xs, ys, nb_edges); + do { + blocks[i]->set_position(drand48() * w/2 + w/4, drand48() * h/2 + h/4, M_PI * drand48()); + blocks[i]->set_speed(0, 0, 0); + // universe.initialize(blocks[i]); + for(int j = 0; j < i + 1; j++) { + universe.initialize(blocks[j]); + } + } while(universe.collide(blocks[i])); + universe.add_polygon(blocks[i]); + cout << "added " << i+1 << " blocks." << endl; } - scalar_t dt = 0.01; + scalar_t dt = 0.1; - for(int n = 0; n < 1000; n++) { - for(int h = 0; h < nb_heur; h++) { - blocks[h]->apply_force(dt, blocks[h]->_center_x, blocks[h]->_center_y, 0, 10); + for(int n = 0; n < 100000; n++) { + for(int i = 0; i < nb_heur; i++) { + blocks[i]->apply_force(dt, blocks[i]->_center_x, blocks[i]->_center_y, 0, 0.25); } universe.update(dt); + if(n%1000 == 0) { + char buffer[1024]; + sprintf(buffer, "/tmp/mash_%06d.fig", n); + ofstream os(buffer); + universe.print_fig(os); + cout << "Wrote " << buffer << endl; + } } ofstream os("/tmp/mash.fig"); universe.print_fig(os); + exit(0); } -- 2.39.5