Skip to content
Snippets Groups Projects
Commit 2f6d030c authored by Elias Pipping's avatar Elias Pipping
Browse files

[Cleanup] Rename: createRHS -> computeExternalForces

parent 16da6254
No related branches found
No related tags found
No related merge requests found
...@@ -230,13 +230,13 @@ int main(int argc, char *argv[]) { ...@@ -230,13 +230,13 @@ int main(int argc, char *argv[]) {
myAssembler.assembleBodyForce(gravity, density, zenith, gravityFunctional); myAssembler.assembleBodyForce(gravity, density, zenith, gravityFunctional);
// Problem formulation: right-hand side // Problem formulation: right-hand side
auto const createRHS = [&](double _relativeTime, Vector &_ell) { auto const computeExternalForces = [&](double _relativeTime, Vector &_ell) {
myAssembler.assembleNeumann(neumannBoundary, _ell, neumannFunction, myAssembler.assembleNeumann(neumannBoundary, _ell, neumannFunction,
_relativeTime); _relativeTime);
_ell += gravityFunctional; _ell += gravityFunctional;
}; };
Vector ell(fineVertexCount); Vector ell(fineVertexCount);
createRHS(0.0, ell); computeExternalForces(0.0, ell);
// {{{ Initial conditions // {{{ Initial conditions
ScalarVector alpha_initial(fineVertexCount); ScalarVector alpha_initial(fineVertexCount);
...@@ -379,7 +379,7 @@ int main(int argc, char *argv[]) { ...@@ -379,7 +379,7 @@ int main(int argc, char *argv[]) {
timeSteppingScheme->nextTimeStep(); timeSteppingScheme->nextTimeStep();
auto const relativeTime = double(run) / double(timesteps); auto const relativeTime = double(run) / double(timesteps);
createRHS(relativeTime, ell); computeExternalForces(relativeTime, ell);
Matrix velocityMatrix; Matrix velocityMatrix;
Vector velocityRHS(fineVertexCount); Vector velocityRHS(fineVertexCount);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment