//inline int virial_calc(t_atom *a,t_3dvec *f,t_3dvec *d)
// __attribute__((always_inline));
int check_per_bound(t_moldyn *moldyn,t_3dvec *a);
//inline int virial_calc(t_atom *a,t_3dvec *f,t_3dvec *d)
// __attribute__((always_inline));
int check_per_bound(t_moldyn *moldyn,t_3dvec *a);