From 7cef2174a35fe1fd2965d4f7b285231e0808b42e Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Mon, 11 May 2020 15:22:22 +0200
Subject: [PATCH] Only print Neumann-Values once, print out number of surface
shell vertices
---
src/film-on-substrate.cc | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/film-on-substrate.cc b/src/film-on-substrate.cc
index 8158d7af..812ced29 100644
--- a/src/film-on-substrate.cc
+++ b/src/film-on-substrate.cc
@@ -256,8 +256,10 @@ int main (int argc, char *argv[]) try
auto neumannBoundary = std::make_shared<BoundaryPatch<GridView>>(gridView, neumannVertices);
BoundaryPatch<GridView> surfaceShellBoundary(gridView, surfaceShellVertices);
- if (mpiHelper.rank()==0)
+ if (mpiHelper.rank()==0) {
std::cout << "Neumann boundary has " << neumannBoundary->numFaces() << " faces\n";
+ std::cout << "Shell boundary has " << surfaceShellBoundary.numFaces() << " faces\n";
+ }
BitSetVector<1> dirichletNodes(feBasis.size(), false);
@@ -360,8 +362,7 @@ int main (int argc, char *argv[]) try
#if DUNE_VERSION_LT(DUNE_ELASTICITY, 2, 8)
std::shared_ptr<NeumannFunction> neumannFunction;
- neumannFunction = std::make_shared<NeumannFunction>(parameterSet.get<FieldVector<double,dim> >("neumannValues"),
- homotopyParameter);
+ neumannFunction = std::make_shared<NeumannFunction>(neumannValues, homotopyParameter);
#else
// A constant vector-valued function, for simple Neumann boundary values
std::shared_ptr<std::function<Dune::FieldVector<double,dim>(Dune::FieldVector<double,dim>)>> neumannFunctionPtr;
--
GitLab