X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=potentials%2Fharmonic_oscillator.c;h=8addd3a980a9ebb43c79c4a2ed6c994ecd11f2fa;hb=1965279af348fbb9b69459d01516bbcd52b6f240;hp=b203453598ae34a149a85dc1b79827d4f8402b01;hpb=20409ee0c545235c9246edde2d3cda5de0ddabde;p=physik%2Fposic.git

diff --git a/potentials/harmonic_oscillator.c b/potentials/harmonic_oscillator.c
index b203453..8addd3a 100644
--- a/potentials/harmonic_oscillator.c
+++ b/potentials/harmonic_oscillator.c
@@ -19,6 +19,45 @@
 #include "../math/math.h"
 #include "harmonic_oscillator.h"
 
+int harmonic_oscillator_set_params(t_moldyn *moldyn,int element) {
+
+	t_ho_params *p;
+
+	// set cutoff before parameters (actually only necessary for some pots)
+	if(moldyn->cutoff==0.0) {
+		printf("[harmonic oscillator] WARNING: no cutoff!\n");
+		return -1;
+	}
+
+	/* alloc mem for potential parameters */
+	moldyn->pot_params=malloc(sizeof(t_ho_params));
+	if(moldyn->pot_params==NULL) {
+		perror("[harmonic oscillator] pot params alloc");
+		return -1;
+	}
+
+	/* these are now ho parameters */
+	p=moldyn->pot_params;
+
+	switch(element) {
+		case SI:
+			/* type: silicon */
+			p->spring_constant=HO_SC_SI;
+			p->equilibrium_distance=HO_ED_SI;
+			break;
+		case C:
+			/* type: carbon */
+			p->spring_constant=HO_SC_C;
+			p->equilibrium_distance=HO_ED_C;
+			break;
+		default:
+			printf("[harmonic oscillator] WARNING: element\n");
+			return -1;
+	}
+
+	return 0;
+}
+
 int harmonic_oscillator(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
 
 	t_ho_params *params;
@@ -55,3 +94,19 @@ int harmonic_oscillator(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
 	return 0;
 }
 
+int harmonic_oscillator_check_2b_bond(t_moldyn *moldyn,
+                                      t_atom *ai,t_atom *aj,u8 bc) {
+
+	t_3dvec distance;
+	double d;
+
+	v3_sub(&distance,&(aj->r),&(ai->r));
+	if(bc) check_per_bound(moldyn,&distance);
+	d=v3_norm(&distance);
+
+	if(d>moldyn->cutoff)
+		return FALSE;
+
+	return TRUE;
+}
+