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 |
|
|