| 37 |
|
for (i=0; i<bandwidth; i++){ |
| 38 |
|
for (j=0; j<bandwidth; j++){ |
| 39 |
|
releaseProbe(startDist); |
| 40 |
+ |
|
| 41 |
+ |
sigmaGrid.push_back(sigDist); |
| 42 |
+ |
sGrid.push_back(sDist); |
| 43 |
+ |
epsGrid.push_back(epsVal); |
| 44 |
+ |
|
| 45 |
|
stepPhi(phiStep); |
| 46 |
|
} |
| 47 |
|
stepTheta(thetaStep); |
| 61 |
|
double interpFrac; |
| 62 |
|
|
| 63 |
|
probeCoor = farPos; |
| 64 |
+ |
potProgress.clear(); |
| 65 |
+ |
distProgress.clear(); |
| 66 |
|
tooClose = 0; |
| 67 |
|
epsVal = 0; |
| 68 |
|
rhoStep = 0.1; //the distance the probe atom moves between steps |
| 69 |
|
|
| 70 |
+ |
|
| 71 |
|
while (!tooClose){ |
| 72 |
|
calcEnergy(); |
| 73 |
|
potProgress.push_back(potEnergy); |