Skip to content
Snippets Groups Projects
Commit 6a1750c0 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Cleanup

parent 579b048f
No related branches found
No related tags found
No related merge requests found
...@@ -238,7 +238,6 @@ int main(int argc, char *argv[]) { ...@@ -238,7 +238,6 @@ int main(int argc, char *argv[]) {
{ {
double const gravity = 9.81; double const gravity = 9.81;
double const density = parset.get<double>("body.density"); double const density = parset.get<double>("body.density");
{ {
MassAssembler<GridType, P1Basis::LocalFiniteElement, MassAssembler<GridType, P1Basis::LocalFiniteElement,
P1Basis::LocalFiniteElement, P1Basis::LocalFiniteElement,
...@@ -246,7 +245,7 @@ int main(int argc, char *argv[]) { ...@@ -246,7 +245,7 @@ int main(int argc, char *argv[]) {
OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis) OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
.assemble(localMass, massMatrix); .assemble(localMass, massMatrix);
massMatrix *= density; massMatrix *= density;
}; }
{ {
double volume = 1.0; double volume = 1.0;
for (int i = 0; i < dims; ++i) for (int i = 0; i < dims; ++i)
...@@ -261,7 +260,7 @@ int main(int argc, char *argv[]) { ...@@ -261,7 +260,7 @@ int main(int argc, char *argv[]) {
// V * g * rho / A = sigma_n // V * g * rho / A = sigma_n
// m^d * N/kg * kg/m^d / m^(d-1) = N/m^(d-1) // m^d * N/kg * kg/m^d / m^(d-1) = N/m^(d-1)
normalStress = volume * gravity * density / area; normalStress = volume * gravity * density / area;
}; }
{ {
SmallVector weightedGravitationalDirection(0); SmallVector weightedGravitationalDirection(0);
weightedGravitationalDirection[1] = -density * gravity; weightedGravitationalDirection[1] = -density * gravity;
...@@ -271,7 +270,7 @@ int main(int argc, char *argv[]) { ...@@ -271,7 +270,7 @@ int main(int argc, char *argv[]) {
gravityFunction); gravityFunction);
FunctionalAssembler<P1Basis>(p1Basis) FunctionalAssembler<P1Basis>(p1Basis)
.assemble(gravityFunctionalAssembler, gravityFunctional, true); .assemble(gravityFunctionalAssembler, gravityFunctional, true);
}; }
} }
SingletonVectorType surfaceNormalStress(finestSize); SingletonVectorType surfaceNormalStress(finestSize);
surfaceNormalStress = normalStress; surfaceNormalStress = normalStress;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment