| 328 |
|
|
| 329 |
|
void SlabStatistics::processFrame(int istep) { |
| 330 |
|
RealType z; |
| 331 |
+ |
|
| 332 |
|
hmat_ = currentSnapshot_->getHmat(); |
| 333 |
|
for (int i = 0; i < nBins_; i++) { |
| 334 |
|
z = (((RealType)i + 0.5) / (RealType)nBins_) * hmat_(2,2); |
| 354 |
|
SpatialStatistics(info, filename, sele, nbins), coordinateOrigin_(V3Zero) { |
| 355 |
|
|
| 356 |
|
binWidth_ = 1.0; |
| 357 |
+ |
|
| 358 |
+ |
Globals* simParams = info->getSimParams(); |
| 359 |
+ |
RNEMDParameters* rnemdParams = simParams->getRNEMDParameters(); |
| 360 |
+ |
bool hasCoordinateOrigin = rnemdParams->haveCoordinateOrigin(); |
| 361 |
|
|
| 362 |
+ |
if (hasCoordinateOrigin) { |
| 363 |
+ |
coordinateOrigin_ = rnemdParams->getCoordinateOrigin(); |
| 364 |
+ |
} else { |
| 365 |
+ |
coordinateOrigin_ = V3Zero; |
| 366 |
+ |
} |
| 367 |
+ |
|
| 368 |
|
r_ = new OutputData; |
| 369 |
|
r_->units = "Angstroms"; |
| 370 |
|
r_->title = "R"; |
| 384 |
|
ShellStatistics::~ShellStatistics() { |
| 385 |
|
} |
| 386 |
|
|
| 387 |
< |
int ShellStatistics::getBin(Vector3d pos) { |
| 387 |
> |
int ShellStatistics::getBin(Vector3d pos) { |
| 388 |
|
Vector3d rPos = pos - coordinateOrigin_; |
| 389 |
|
return int(rPos.length() / binWidth_); |
| 390 |
|
} |