#include "../moldyn.h"
#include "../math/math.h"
-//#include "harmonic_oscillator.h"
+#include "harmonic_oscillator.h"
int harmonic_oscillator(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
double d,f;
double sc,equi_dist;
- params=moldyn->pot2b_params;
+ params=moldyn->pot_params;
sc=params->spring_constant;
equi_dist=params->equilibrium_distance;