X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?a=blobdiff_plain;f=energy.c;h=aa2dd47fa4364eed4fc5c616ca2336ba9573eee9;hb=d045b4e5e53daf45808a65f8811a1c6fbf7ecc98;hp=2f8a854cde076656b02c3c0799cc2fbdd5634cf1;hpb=c343db120c3480b7b907cfc22cdd59edd3f01388;p=moebius2.git diff --git a/energy.c b/energy.c index 2f8a854..aa2dd47 100644 --- a/energy.c +++ b/energy.c @@ -516,7 +516,7 @@ void compute_rim_twist_angles(const Vertices vertices, int section) { /* we are interested in the angle subtended at the rim, from the * rim's point of view. */ K distance[k]= vertices[v][k] - oncircle[k]; - double distance_positive_z= distance[3]; + double distance_positive_z= distance[2]; double distance_radial_outwards= dotprod(distance, oncircle); rim_vertex_angles[v]= atan2(distance_positive_z, distance_radial_outwards); }