Trisurf Monte Carlo simulator
Samo Penic
2021-04-19 efb615c065452b58799e32735647b0559698163f
src/bond.c
@@ -1,3 +1,4 @@
/* vim: set ts=4 sts=4 sw=4 noet : */
#include<stdlib.h>
#include "general.h"
#include "vertex.h"
@@ -38,14 +39,27 @@
ts_bool bond_vector(ts_bond *bond){
   
   bond->x=bond->vtx1->x-bond->vtx2->x;
   bond->y=bond->vtx1->y-bond->vtx2->y;
   bond->z=bond->vtx1->z-bond->vtx2->z;
   bond->x = bond->vtx1->x - bond->vtx2->x;
   bond->y = bond->vtx1->y - bond->vtx2->y;
   bond->z = bond->vtx1->z - bond->vtx2->z;
   return TS_SUCCESS;   
}
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++){