Found one bug in the code when calculating mprod. Fixed.

Samo Penic

2021-07-27 3197854e89eb9ffe5ff3137e274d16806c2028a5

 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++){