void compute_vertex_areas(const Vertices vertices, int section);
void compute_edge_lengths(const Vertices vertices, int section);
+void compute_rim_twist_angles(const Vertices vertices, int section);
-extern double vertex_areas[N], vertex_mean_edge_lengths[N], edge_lengths[N][V6];
+//extern double vertex_areas[N], edge_lengths[N][V6];
+extern double vertex_mean_edge_lengths[N];
extern const double edge_angle_cost_circcircrat;
double edge_length_variation_cost(const Vertices vertices, int section);
double prop_edge_length_variation_cost(const Vertices vertices, int section);
double rim_proximity_cost(const Vertices vertices, int section);
+double rim_twist_cost(const Vertices vertices, int section);
double edge_angle_cost(const Vertices vertices, int section);
double small_triangles_cost(const Vertices vertices, int section);
double nonequilateral_triangles_cost(const Vertices vertices, int section);