| 233 |
|
|
| 234 |
|
kinetic *= 0.5; |
| 235 |
|
eTemp = (2.0 * kinetic) / |
| 236 |
< |
(info_->getNFluctuatingCharges() * PhysicalConstants::kb ); |
| 236 |
> |
(info_->getNFluctuatingCharges() * PhysicalConstants::kb ); |
| 237 |
|
|
| 238 |
|
snap->setElectronicTemperature(eTemp); |
| 239 |
|
} |
| 920 |
|
} |
| 921 |
|
|
| 922 |
|
RealType Thermo::getHullVolume(){ |
| 923 |
+ |
#ifdef HAVE_QHULL |
| 924 |
|
Snapshot* snap = info_->getSnapshotManager()->getCurrentSnapshot(); |
| 924 |
– |
|
| 925 |
– |
#ifdef HAVE_QHULL |
| 925 |
|
if (!snap->hasHullVolume) { |
| 926 |
|
Hull* surfaceMesh_; |
| 927 |
< |
|
| 927 |
> |
|
| 928 |
|
Globals* simParams = info_->getSimParams(); |
| 929 |
|
const std::string ht = simParams->getHULL_Method(); |
| 930 |
|
|
| 956 |
|
// Compute surface Mesh |
| 957 |
|
surfaceMesh_->computeHull(localSites_); |
| 958 |
|
snap->setHullVolume(surfaceMesh_->getVolume()); |
| 959 |
< |
|
| 959 |
> |
|
| 960 |
|
delete surfaceMesh_; |
| 961 |
|
} |
| 962 |
< |
|
| 962 |
> |
|
| 963 |
|
return snap->getHullVolume(); |
| 964 |
|
#else |
| 965 |
|
return 0.0; |