From a8bde65de2faae58ed2e37c0d6f91b0be60f0fd0 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 31 Jan 2011 14:03:26 +0000
Subject: [PATCH] fail if no Dirichlet boundary is given at all
[[Imported from SVN: r6910]]
---
dirneucoupling.cc | 3 +++
1 file changed, 3 insertions(+)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index de0edb4a..43bffc2d 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -160,6 +160,9 @@ int main (int argc, char *argv[]) try
rodFactory.create(complex.rods_["rod"].dirichletValues_,
RigidBodyMotion<3>(FieldVector<double,3>(0), Rotation<3,double>::identity()));
BitSetVector<1> rodDNodes(complex.rods_["rod"].dirichletValues_.size(), false);
+
+ // we need at least one Dirichlet side
+ assert(parameterSet.hasKey("dirichletValue0") or parameterSet.hasKey("dirichletValue1"));
if (parameterSet.hasKey("dirichletValue0")){
--
GitLab