1 |
|
#include "GridBuilder.hpp" |
2 |
– |
#include "MatVec3.h" |
2 |
|
#define PI 3.14159265359 |
3 |
|
|
4 |
|
|
24 |
|
double sigTemp, sTemp, epsTemp, sigProbe; |
25 |
|
double minDist = 10.0; //minimum start distance |
26 |
|
|
28 |
– |
sList = sGrid; |
27 |
|
sigList = sigmaGrid; |
28 |
+ |
sList = sGrid; |
29 |
|
epsList = epsGrid; |
30 |
|
forcefield = forceField; |
31 |
|
|
67 |
|
|
68 |
|
for (k =0; k < gridwidth; k++) { |
69 |
|
thetaVal = thetaMin + k*thetaStep; |
71 |
– |
printf("Theta step %i\n", k); |
70 |
|
for (j=0; j < gridwidth; j++) { |
71 |
< |
phiVal = j*phiStep; |
71 |
> |
//s2kit10 is actually taking a grid from -pi/2 to 3pi/3, not 0 to 2pi... |
72 |
> |
phiVal = j*phiStep - 0.5*PI; |
73 |
> |
if (phiVal<0.0) |
74 |
> |
phiVal += 2*PI; |
75 |
|
|
76 |
|
rbMol->setEuler(0.0, thetaVal, phiVal); |
77 |
|
|
224 |
|
epsOut << epsList[k] << "\n0\n"; |
225 |
|
} |
226 |
|
} |
227 |
+ |
|