| 105 |
|
|
| 106 |
|
ji += (dt2 * PhysicalConstants::energyConvert) * Tb; |
| 107 |
|
|
| 108 |
< |
rotAlgo->rotate(integrableObject, ji, dt); |
| 108 |
> |
rotAlgo_->rotate(integrableObject, ji, dt); |
| 109 |
|
|
| 110 |
|
integrableObject->setJ(ji); |
| 111 |
|
} |
| 114 |
|
} |
| 115 |
|
} //end for(mol = info_->beginMolecule(i)) |
| 116 |
|
|
| 117 |
< |
rattle->constraintA(); |
| 118 |
< |
|
| 117 |
> |
flucQ_->moveA(); |
| 118 |
> |
rattle_->constraintA(); |
| 119 |
|
} |
| 120 |
|
|
| 121 |
|
void LangevinDynamics::moveB(){ |
| 161 |
|
} |
| 162 |
|
} //end for(mol = info_->beginMolecule(i)) |
| 163 |
|
|
| 164 |
< |
|
| 165 |
< |
rattle->constraintB(); |
| 166 |
< |
|
| 164 |
> |
flucQ_->moveB(); |
| 165 |
> |
rattle_->constraintB(); |
| 166 |
|
} |
| 167 |
|
|
| 168 |
|
|