| 53 |
|
#include "integrators/LDForceManager.hpp" |
| 54 |
|
namespace oopse { |
| 55 |
|
|
| 56 |
< |
|
| 56 |
> |
|
| 57 |
|
LangevinDynamics::LangevinDynamics(SimInfo* info) : VelocityVerletIntegrator(info){ |
| 58 |
|
setForceManager(new LDForceManager(info)); |
| 59 |
|
} |
| 60 |
< |
|
| 60 |
> |
|
| 61 |
|
void LangevinDynamics::moveA(){ |
| 62 |
|
SimInfo::MoleculeIterator i; |
| 63 |
|
Molecule::IntegrableObjectIterator j; |
| 68 |
|
Vector3d frc; |
| 69 |
|
Vector3d Tb; |
| 70 |
|
Vector3d ji; |
| 71 |
< |
double mass; |
| 71 |
> |
RealType mass; |
| 72 |
|
|
| 73 |
|
for (mol = info_->beginMolecule(i); mol != NULL; mol = info_->nextMolecule(i)) { |
| 74 |
|
for (integrableObject = mol->beginIntegrableObject(j); integrableObject != NULL; |
| 122 |
|
Vector3d frc; |
| 123 |
|
Vector3d Tb; |
| 124 |
|
Vector3d ji; |
| 125 |
< |
double mass; |
| 125 |
> |
RealType mass; |
| 126 |
|
|
| 127 |
|
for (mol = info_->beginMolecule(i); mol != NULL; mol = info_->nextMolecule(i)) { |
| 128 |
|
for (integrableObject = mol->beginIntegrableObject(j); integrableObject != NULL; |
| 162 |
|
} |
| 163 |
|
|
| 164 |
|
|
| 165 |
< |
double LangevinDynamics::calcConservedQuantity() { |
| 165 |
> |
RealType LangevinDynamics::calcConservedQuantity() { |
| 166 |
|
return 0.0; |
| 167 |
|
} |
| 168 |
|
|