From f0b22bf5db6fca3c715e44632a0a04d697e3e1fb Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 25 Mar 2015 18:56:30 +0100
Subject: [PATCH] [Problem] Make L2 error relative to surface area
---
src/sand-wedge.cc | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/sand-wedge.cc b/src/sand-wedge.cc
index c8795737..6f7dff09 100644
--- a/src/sand-wedge.cc
+++ b/src/sand-wedge.cc
@@ -190,11 +190,12 @@ int main(int argc, char *argv[]) {
matrices.damping);
myAssembler.assembleMass(body.getDensityField(), matrices.mass);
- ScalarMatrix frictionalBoundaryMass;
+ ScalarMatrix relativeFrictionalBoundaryMass;
myAssembler.assembleFrictionalBoundaryMass(frictionalBoundary,
- frictionalBoundaryMass);
+ relativeFrictionalBoundaryMass);
+ relativeFrictionalBoundaryMass /= frictionalBoundary.area();
EnergyNorm<ScalarMatrix, ScalarVector> const stateEnergyNorm(
- frictionalBoundaryMass);
+ relativeFrictionalBoundaryMass);
// Assemble forces
Vector gravityFunctional;
--
GitLab