@@ -24,20 +24,20 @@ void sinNPhiCosNPhi(int n, float* sinNPhi, float* cosNPhi,
24
24
25
25
template <int N>
26
26
struct SinCos {
27
- static void sinNPhiCosNPhi (double & sinNPhi, double & cosNPhi,
28
- double sinPhi, double cosPhi ) {
27
+ static void sinNPhiCosNPhi (float & sinNPhi, float & cosNPhi,
28
+ float sinPhi, float cosPhi ) {
29
29
float sinNm1Phi;
30
30
float cosNm1Phi;
31
- sinNPhiCosNPhi <N-1 >(sinNm1Phi,cosNm1Phi,sinPhi,cosPhi);
31
+ SinCos <N-1 >:: sinNPhiCosNPhi (sinNm1Phi,cosNm1Phi,sinPhi,cosPhi);
32
32
sinNPhi = cosPhi*sinNm1Phi+sinPhi*cosNm1Phi;
33
33
cosNPhi = cosPhi*cosNm1Phi-sinPhi*sinNm1Phi;
34
34
}
35
35
};
36
36
37
37
template <>
38
38
struct SinCos <1 > {
39
- static void sinNPhiCosNPhi (double & sinNPhi, double & cosNPhi,
40
- double sinPhi, double cosPhi ) {
39
+ static void sinNPhiCosNPhi (float & sinNPhi, float & cosNPhi,
40
+ float sinPhi, float cosPhi ) {
41
41
sinNPhi = sinPhi;
42
42
cosNPhi = cosPhi;
43
43
}
@@ -83,7 +83,7 @@ float dihedral_energy(Dihedral* dihedral_begin, Dihedral* dihedral_end, float* p
83
83
float cosPhase;
84
84
float V;
85
85
float DN;
86
- float IN;
86
+ int IN;
87
87
float x1, y1, z1;
88
88
float x2, y2, z2;
89
89
float x3, y3, z3;
@@ -150,7 +150,7 @@ float old_dihedral_energy(Dihedral* dihedral_begin, Dihedral*dihedral_end, float
150
150
float cosPhase;
151
151
float V;
152
152
float DN;
153
- float IN;
153
+ int IN;
154
154
float x1, y1, z1;
155
155
float x2, y2, z2;
156
156
float x3, y3, z3;
0 commit comments