diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index b959bfa3eb627956318722c285a85cf9a93c1873..adbfeaeb1b6b2593a518db5d4b73ce423a8d5755 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -238,7 +238,6 @@ int main(int argc, char *argv[]) { { double const gravity = 9.81; double const density = parset.get<double>("body.density"); - { MassAssembler<GridType, P1Basis::LocalFiniteElement, P1Basis::LocalFiniteElement, @@ -246,7 +245,7 @@ int main(int argc, char *argv[]) { OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis) .assemble(localMass, massMatrix); massMatrix *= density; - }; + } { double volume = 1.0; for (int i = 0; i < dims; ++i) @@ -261,7 +260,7 @@ int main(int argc, char *argv[]) { // V * g * rho / A = sigma_n // m^d * N/kg * kg/m^d / m^(d-1) = N/m^(d-1) normalStress = volume * gravity * density / area; - }; + } { SmallVector weightedGravitationalDirection(0); weightedGravitationalDirection[1] = -density * gravity; @@ -271,7 +270,7 @@ int main(int argc, char *argv[]) { gravityFunction); FunctionalAssembler<P1Basis>(p1Basis) .assemble(gravityFunctionalAssembler, gravityFunctional, true); - }; + } } SingletonVectorType surfaceNormalStress(finestSize); surfaceNormalStress = normalStress;