| 301 |
|
int nSurfaceSDs = surfaceMesh_->getNs(); |
| 302 |
|
|
| 303 |
|
|
| 304 |
< |
std::vector<Triangle*> sMesh = surfaceMesh_->getMesh(); |
| 304 |
> |
std::vector<Triangle> sMesh = surfaceMesh_->getMesh(); |
| 305 |
|
int nTriangles = sMesh.size(); |
| 306 |
|
|
| 307 |
|
|
| 319 |
|
/* Loop over the mesh faces and apply random force to each of the faces*/ |
| 320 |
|
|
| 321 |
|
|
| 322 |
< |
std::vector<Triangle*>::iterator face; |
| 322 |
> |
std::vector<Triangle>::iterator face; |
| 323 |
|
std::vector<StuntDouble*>::iterator vertex; |
| 324 |
|
int thisNumber = 0; |
| 325 |
|
for (face = sMesh.begin(); face != sMesh.end(); ++face){ |
| 326 |
|
|
| 327 |
< |
Triangle* thisTriangle = *face; |
| 328 |
< |
std::vector<StuntDouble*> vertexSDs = thisTriangle->getVertices(); |
| 327 |
> |
Triangle thisTriangle = *face; |
| 328 |
> |
std::vector<StuntDouble*> vertexSDs = thisTriangle.getVertices(); |
| 329 |
|
|
| 330 |
|
/* Get Random Force */ |
| 331 |
< |
Vector3d unitNormal = thisTriangle->getNormal(); |
| 331 |
> |
Vector3d unitNormal = thisTriangle.getNormal(); |
| 332 |
|
unitNormal.normalize(); |
| 333 |
|
Vector3d randomForce = -randNums[thisNumber] * unitNormal; |
| 334 |
< |
Vector3d centroid = thisTriangle->getCentroid(); |
| 334 |
> |
Vector3d centroid = thisTriangle.getCentroid(); |
| 335 |
|
|
| 336 |
< |
Vector3d langevinForce = randomForce - gamma_t_*thisTriangle->getFacetVelocity(); |
| 336 |
> |
Vector3d langevinForce = randomForce - gamma_t_*thisTriangle.getFacetVelocity(); |
| 337 |
|
|
| 338 |
|
for (vertex = vertexSDs.begin(); vertex != vertexSDs.end(); ++vertex){ |
| 339 |
|
if ((*vertex) != NULL){ |
