From de94c6e4b1613c628f78863150f2a494d749cba4 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 29 Nov 2011 11:10:53 +0100
Subject: [PATCH] Use a time-dependent neumann condition
---
src/one-body-sample.cc | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 94530eb8..9611eb14 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -101,8 +101,9 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis,
int run) { // constant sample function on neumann boundary
BoundaryPatch<GridView> neumannBoundary(gridView, neumannNodes);
LocalVectorType SampleVector(0);
- // FIXME: random values (time-independent)
- SampleVector[0] = 1;
+ // FIXME: random values (time-dependent)
+ SampleVector[0] =
+ run / 10.0; // If we leave this fixed, we get into trouble for some reason
SampleVector[1] = 0;
ConstantFunction<LocalVectorType, LocalVectorType> fNeumann(SampleVector);
NeumannBoundaryAssembler<GridType, LocalVectorType> neumannBoundaryAssembler(
--
GitLab