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