src/bond.c
@@ -47,6 +47,19 @@ } ts_bool bond_get_edge_vector(ts_double *vector, ts_bond *bond, ts_vertex *vtx){ if(vtx==bond->vtx2){ vector[0]=bond->x; vector[1]=bond->y; vector[2]=bond->z; } else { vector[0]=-bond->x; vector[1]=-bond->y; vector[2]=-bond->z; } return TS_SUCCESS; } ts_bool bond_list_free(ts_bond_list *blist){ ts_uint i; for(i=0;i<blist->n;i++){