| 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; |
| 135 |
|
// velocity half step |
| 136 |
|
vel += (dt2 /mass * OOPSEConstant::energyConvert) * frc; |
| 137 |
|
|
| 138 |
< |
integrableObject->setVel(vel); |
| 138 |
> |
integrableObject->setVel(vel); |
| 139 |
|
|
| 140 |
|
if (integrableObject->isDirectional()){ |
| 141 |
|
|
| 162 |
|
} |
| 163 |
|
|
| 164 |
|
|
| 165 |
< |
double NVE::calcConservedQuantity() { |
| 165 |
> |
RealType NVE::calcConservedQuantity() { |
| 166 |
|
return thermo.getTotalE(); |
| 167 |
|
} |
| 168 |
|
|