35 |
|
* |
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). |
38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
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 |
+ |
#ifdef IS_MPI |
44 |
+ |
#include <mpi.h> |
45 |
+ |
#endif |
46 |
+ |
|
47 |
|
#include "config.h" |
48 |
|
#include <cmath> |
49 |
|
|
56 |
|
#include "utils/StringUtils.hpp" |
57 |
|
#include "selection/SelectionEvaluator.hpp" |
58 |
|
#include "selection/SelectionManager.hpp" |
55 |
– |
#ifdef IS_MPI |
56 |
– |
#include <mpi.h> |
57 |
– |
#endif |
59 |
|
|
59 |
– |
|
60 |
|
namespace OpenMD { |
61 |
|
|
62 |
|
RestraintForceManager::RestraintForceManager(SimInfo* info): ForceManager(info) { |
101 |
|
if (myType.compare("MOLECULAR")==0){ |
102 |
|
|
103 |
|
int molIndex; |
104 |
– |
std::vector<Vector3d> ref; |
104 |
|
Vector3d refCom; |
105 |
|
|
106 |
|
if (!stamp[i]->haveMolIndex()) { |