| 164 |
|
|
| 165 |
|
if (gayBerneData != NULL) { |
| 166 |
|
GayBerneParam gayBerneParam = gayBerneData->getData(); |
| 167 |
< |
currShape = new Ellipsoid(V3Zero, |
| 167 |
> |
currShape = new Ellipsoid(V3Zero, |
| 168 |
|
gayBerneParam.GB_l / 2.0, |
| 169 |
|
gayBerneParam.GB_d / 2.0, |
| 170 |
|
Mat3x3d::identity()); |
| 226 |
|
simError(); |
| 227 |
|
} |
| 228 |
|
|
| 229 |
– |
|
| 229 |
|
HydroProp* currHydroProp = currShape->getHydroProp(simParams->getViscosity(),simParams->getTargetTemp()); |
| 230 |
+ |
|
| 231 |
|
std::map<std::string, HydroProp*>::iterator iter = hydroPropMap.find(integrableObject->getType()); |
| 232 |
|
if (iter != hydroPropMap.end()) |
| 233 |
|
hydroProps_.push_back(iter->second); |
| 259 |
|
return props; |
| 260 |
|
} |
| 261 |
|
|
| 262 |
< |
void LDForceManager::postCalculation(bool needStress){ |
| 262 |
> |
void LDForceManager::postCalculation(){ |
| 263 |
|
SimInfo::MoleculeIterator i; |
| 264 |
|
Molecule::IntegrableObjectIterator j; |
| 265 |
|
Molecule* mol; |
| 323 |
|
genRandomForceAndTorque(randomForceBody, randomTorqueBody, index, variance_); |
| 324 |
|
Vector3d randomForceLab = Atrans * randomForceBody; |
| 325 |
|
Vector3d randomTorqueLab = Atrans * randomTorqueBody; |
| 326 |
+ |
|
| 327 |
|
integrableObject->addFrc(randomForceLab); |
| 328 |
|
integrableObject->addTrq(randomTorqueLab + cross(rcrLab, randomForceLab )); |
| 329 |
|
|
| 469 |
|
if(!simParams->getUsePeriodicBoundaryConditions()) |
| 470 |
|
veloMunge->removeAngularDrift(); |
| 471 |
|
|
| 472 |
< |
ForceManager::postCalculation(needStress); |
| 472 |
> |
ForceManager::postCalculation(); |
| 473 |
|
} |
| 474 |
|
|
| 475 |
|
void LDForceManager::genRandomForceAndTorque(Vector3d& force, Vector3d& torque, unsigned int index, RealType variance) { |