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; |