COST(1000.0, edgewise_vertex_displacement_cost(vertices));
COST(1.0, graph_layout_cost(vertices,vertex_areas));
- COST(1e3, noncircular_rim_cost(vertices));
+ COST(1e-30, noncircular_rim_cost(vertices));
printf("| total %# e |", energy);
if (energy < best_energy) {
printf(" BEST");
best_f= fopen(BEST_F ".new","wb"); if (!best_f) diee("fopen new best");
- r= fwrite(vertices,sizeof(vertices),1,best_f); if (r!=1) diee("fwrite");
+ r= fwrite(vertices,sizeof(Vertices),1,best_f); if (r!=1) diee("fwrite");
if (fclose(best_f)) diee("fclose new best");
if (rename(BEST_F ".new", BEST_F)) diee("rename install new best");
+
+ best_energy= energy;
}
putchar('\n');
flushoutput();
initial_gsl.owner= 0;
step_size_gsl= initial_gsl;
- initial_gsl.data= (double*)initial;
- step_size_gsl.data= (double*)step_size;
+ initial_gsl.data= &initial[0][0];
+ step_size_gsl.data= &step_size[0][0];
FOR_VERTEX(v)
- K step_size[v][k]= 1e-3;
+ K step_size[v][k]= 1e-4;
FOR_RIM_VERTEX(vx,vy,v)
step_size[v][3] *= 0.1;