ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/integrators/Integrator.cpp
(Generate patch)

Comparing branches/development/src/integrators/Integrator.cpp (file contents):
Revision 1758 by gezelter, Wed Jun 20 13:55:52 2012 UTC vs.
Revision 1870 by gezelter, Tue May 7 19:09:54 2013 UTC

# Line 35 | Line 35
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   */
# Line 43 | Line 43
43   #include "brains/Snapshot.hpp"
44   #include "integrators/Integrator.hpp"
45   #include "integrators/DLM.hpp"
46 + #include "flucq/FluctuatingChargeLangevin.hpp"
47   #include "flucq/FluctuatingChargeNVT.hpp"
48   #include "utils/simError.h"
49  
# Line 52 | Line 53 | namespace OpenMD {
53        needReset(false), velocitizer_(NULL), needVelocityScaling(false),
54        rnemd_(NULL), useRNEMD(false), rotAlgo_(NULL), flucQ_(NULL),
55        rattle_(NULL), dumpWriter(NULL), statWriter(NULL), thermo(info),
56 <      currentSnapshot_(info->getSnapshotManager()->getCurrentSnapshot()) {
56 >      snap(info->getSnapshotManager()->getCurrentSnapshot()) {
57      
58      simParams = info->getSimParams();
59      
# Line 94 | Line 95 | namespace OpenMD {
95      }
96      
97      if (!simParams->getUseInitalTime()) {
98 <      currentSnapshot_->setTime(0.0);
98 >      snap->setTime(0.0);
99      }
100      
101      if (simParams->haveResetTime()) {
# Line 142 | Line 143 | namespace OpenMD {
143      
144      rotAlgo_ = new DLM();
145      rattle_ = new Rattle(info);
146 <    //forceMan_->initialize();
147 <    flucQ_ = new FluctuatingChargeNVT(info, forceMan_);
146 >    if (simParams->getFluctuatingChargeParameters()->havePropagator()) {
147 >      std::string prop = toUpperCopy(simParams->getFluctuatingChargeParameters()->getPropagator());
148 >      if (prop.compare("NVT")==0){
149 >         flucQ_ = new FluctuatingChargeNVT(info);
150 >      } else if (prop.compare("LANGEVIN")==0) {
151 >         flucQ_ = new FluctuatingChargeLangevin(info);
152 >      } else {
153 >        sprintf(painCave.errMsg,
154 >                "Integrator Error: Unknown Fluctuating Charge propagator (%s) requested\n",
155 >                simParams->getFluctuatingChargeParameters()->getPropagator().c_str());
156 >        painCave.isFatal = 1;
157 >      }
158 >    }
159 >    flucQ_->setForceManager(forceMan_);
160    }
161    
162    Integrator::~Integrator(){

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines