| 101 |
|
|
| 102 |
|
ji += (dt2 * PhysicalConstants::energyConvert) * Tb; |
| 103 |
|
|
| 104 |
< |
rotAlgo->rotate(integrableObject, ji, dt); |
| 104 |
> |
rotAlgo_->rotate(integrableObject, ji, dt); |
| 105 |
|
|
| 106 |
|
integrableObject->setJ(ji); |
| 107 |
|
} |
| 109 |
|
|
| 110 |
|
} |
| 111 |
|
} //end for(mol = info_->beginMolecule(i)) |
| 112 |
< |
|
| 113 |
< |
rattle->constraintA(); |
| 114 |
< |
|
| 112 |
> |
flucQ_->moveA(); |
| 113 |
> |
rattle_->constraintA(); |
| 114 |
|
} |
| 115 |
|
|
| 116 |
|
void NVE::moveB(){ |
| 156 |
|
} |
| 157 |
|
} //end for(mol = info_->beginMolecule(i)) |
| 158 |
|
|
| 159 |
< |
|
| 160 |
< |
rattle->constraintB(); |
| 162 |
< |
|
| 159 |
> |
flucQ_->moveB(); |
| 160 |
> |
rattle_->constraintB(); |
| 161 |
|
} |
| 162 |
|
|
| 163 |
|
|
| 164 |
|
RealType NVE::calcConservedQuantity() { |
| 165 |
< |
return thermo.getTotalE() ; |
| 165 |
> |
return thermo.getTotalE(); |
| 166 |
|
} |
| 167 |
|
|
| 168 |
|
} //end namespace OpenMD |
