36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
|
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
39 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
39 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
> |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
42 |
|
|
43 |
|
/** |
61 |
|
#include "io/ForceFieldOptions.hpp" |
62 |
|
#include "UseTheForce/ForceField.hpp" |
63 |
|
#include "nonbonded/SwitchingFunction.hpp" |
64 |
+ |
#ifdef IS_MPI |
65 |
+ |
#include <mpi.h> |
66 |
+ |
#endif |
67 |
|
|
68 |
|
using namespace std; |
69 |
|
namespace OpenMD { |
1194 |
|
|
1195 |
|
det = intTensor.determinant(); |
1196 |
|
sysconstants = geomCnst/(RealType)nGlobalIntegrableObjects_; |
1197 |
< |
volume = 4.0/3.0*NumericConstant::PI*pow(sysconstants,3.0/2.0)*sqrt(det); |
1197 |
> |
volume = 4.0/3.0*NumericConstant::PI*pow(sysconstants,geomCnst)*sqrt(det); |
1198 |
|
return; |
1199 |
|
} |
1200 |
|
|
1210 |
|
|
1211 |
|
detI = intTensor.determinant(); |
1212 |
|
sysconstants = geomCnst/(RealType)nGlobalIntegrableObjects_; |
1213 |
< |
volume = 4.0/3.0*NumericConstant::PI*pow(sysconstants,3.0/2.0)*sqrt(detI); |
1213 |
> |
volume = 4.0/3.0*NumericConstant::PI*pow(sysconstants,geomCnst)*sqrt(detI); |
1214 |
|
return; |
1215 |
|
} |
1216 |
|
/* |