From b9a80fd28968fb47c4b28cf03592a96d1e9cd945 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 15 May 2012 14:20:42 +0200
Subject: [PATCH] Define time as h * run
---
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 55add0c6..759c0a1e 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -235,6 +235,7 @@ int main(int argc, char *argv[]) {
auto const &neumannFunction = functions.get("neumannCondition");
for (size_t run = 1; run <= timesteps; ++run) {
+ double const time = h * run;
if (parset.get<bool>("printProgress")) {
std::cout << '*';
std::cout.flush();
@@ -242,12 +243,12 @@ int main(int argc, char *argv[]) {
{
assemble_neumann<GridType, GridView, SmallVector, P1Basis>(
- leafView, p1Basis, neumannNodes, b4, neumannFunction, h * run);
+ leafView, p1Basis, neumannNodes, b4, neumannFunction, time);
stiffnessMatrix.mmv(u4, b4);
// Apply Dirichlet condition
for (int i = 0; i < finestSize; ++i)
if (ignoreNodes[i].count() == dim) {
- dirichletFunction.evaluate(h * run, u4_diff[i][0]);
+ dirichletFunction.evaluate(time, u4_diff[i][0]);
u4_diff[i][0] /= refinement_factor;
}
--
GitLab