From cd7bafd51395b699e16d2ff4c141d02c40a47170 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 23 Jan 2011 13:58:06 +0000
Subject: [PATCH] hand the continuum Dirichlet values to the complex, not the
Steklov-Poincare step
[[Imported from SVN: r6836]]
---
dirneucoupling.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 357c2457..9f0d6f49 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -196,7 +196,10 @@ int main (int argc, char *argv[]) try
for (int i=0; i<dirichletNodes.size(); i++)
dirichletNodes[i] = dirichletBoundary.containsVertex(i);
- sampleOnBitField(*complex.continua_["continuum"].grid_, dirichletValues[0], dirichletValues.back(), dirichletNodes);
+ sampleOnBitField(*complex.continua_["continuum"].grid_,
+ dirichletValues[0],
+ complex.continua_["continuum"].dirichletValues_,
+ dirichletNodes);
/////////////////////////////////////////////////////////////////////
// Create the two interface boundary patches
@@ -459,7 +462,6 @@ int main (int argc, char *argv[]) try
&rodLocalStiffness,
&rodSolver,
&stiffnessMatrix3d,
- &dirichletValues.back(),
solver,
&localAssembler);
--
GitLab