From d46adcc8087871aa770db408d787c12555a639ed Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 20 Nov 2007 18:23:32 +0000
Subject: [PATCH] bugfix in getResultantStress: make sure the boundary patch is
on the maxlevel
[[Imported from SVN: r1792]]
---
src/rodassembler.cc | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index a409ccc0..d165e58a 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -1023,6 +1023,10 @@ getResultantForce(const BoundaryPatch<GridType>& boundary,
if (sol.size()!=indexSet.size(gridDim))
DUNE_THROW(Exception, "Solution vector doesn't match the grid!");
+ /** \todo Eigentlich sollte das BoundaryPatch ein LeafBoundaryPatch sein */
+ if (boundary.level() != grid_->maxLevel())
+ DUNE_THROW(Exception, "The boundary patch has to refer to the max level!");
+
FieldVector<double,3> canonicalStress(0);
canonicalTorque = 0;
--
GitLab