--- trunk/src/integrators/LangevinDynamics.cpp 2006/04/25 02:09:01 945 +++ trunk/src/integrators/LangevinDynamics.cpp 2008/04/18 16:55:15 1237 @@ -56,6 +56,10 @@ namespace oopse { LangevinDynamics::LangevinDynamics(SimInfo* info) : VelocityVerletIntegrator(info){ setForceManager(new LDForceManager(info)); + + // Langevin Dynamics Force Manager needs to know about the half-time step + // size to get convergence on the friction forces: + dynamic_cast(forceMan_)->setDt2(dt2); } void LangevinDynamics::moveA(){ @@ -68,7 +72,7 @@ namespace oopse { Vector3d frc; Vector3d Tb; Vector3d ji; - double mass; + RealType mass; for (mol = info_->beginMolecule(i); mol != NULL; mol = info_->nextMolecule(i)) { for (integrableObject = mol->beginIntegrableObject(j); integrableObject != NULL; @@ -122,7 +126,7 @@ namespace oopse { Vector3d frc; Vector3d Tb; Vector3d ji; - double mass; + RealType mass; for (mol = info_->beginMolecule(i); mol != NULL; mol = info_->nextMolecule(i)) { for (integrableObject = mol->beginIntegrableObject(j); integrableObject != NULL; @@ -162,7 +166,7 @@ namespace oopse { } - double LangevinDynamics::calcConservedQuantity() { + RealType LangevinDynamics::calcConservedQuantity() { return 0.0; }