| 184 |
|
void LangevinHullForceManager::postCalculation(){ |
| 185 |
|
|
| 186 |
|
int nTriangles, thisFacet; |
| 187 |
< |
RealType area, thisArea, thisMass; |
| 187 |
> |
RealType area, thisArea; |
| 188 |
|
vector<Triangle> sMesh; |
| 189 |
|
Triangle thisTriangle; |
| 190 |
|
vector<Triangle>::iterator face; |
| 200 |
|
|
| 201 |
|
// Compute surface Mesh |
| 202 |
|
surfaceMesh_->computeHull(localSites_); |
| 203 |
< |
// Get total area and number of surface stunt doubles |
| 204 |
< |
area = surfaceMesh_->getArea(); |
| 203 |
> |
// Get number of surface stunt doubles |
| 204 |
|
sMesh = surfaceMesh_->getMesh(); |
| 205 |
|
nTriangles = sMesh.size(); |
| 206 |
|
|
| 218 |
|
unitNormal = thisTriangle.getUnitNormal(); |
| 219 |
|
centroid = thisTriangle.getCentroid(); |
| 220 |
|
facetVel = thisTriangle.getFacetVelocity(); |
| 222 |
– |
thisMass = thisTriangle.getFacetMass(); |
| 221 |
|
|
| 222 |
|
langevinForce = V3Zero; |
| 223 |
|
|