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 |
> |
phiVal = j*phiStep + 0.5*PI; |
72 |
|
|
73 |
|
rbMol->setEuler(0.0, thetaVal, phiVal); |
74 |
|
|
221 |
|
epsOut << epsList[k] << "\n0\n"; |
222 |
|
} |
223 |
|
} |
224 |
+ |
|