| 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 | 
> | 
      if (phiVal>=2*PI) | 
| 73 | 
> | 
        phiVal -= 2*PI; | 
| 74 | 
  | 
 | 
| 75 | 
  | 
      rbMol->setEuler(0.0, thetaVal, phiVal); | 
| 76 | 
  | 
 | 
| 223 | 
  | 
    epsOut << epsList[k] << "\n0\n"; | 
| 224 | 
  | 
  } | 
| 225 | 
  | 
} | 
| 226 | 
+ | 
 |