From e5bd1a8c789cc0b87d72d6d91dce2352902c243a Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sat, 25 Feb 2012 19:33:29 +0100
Subject: [PATCH] Print evolution of state
---
src/one-body-sample.cc | 12 ++++++++++++
src/one-body-sample.parset | 1 +
2 files changed, 13 insertions(+)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index dd0b31d8..b753f0f9 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -452,6 +452,18 @@ int main(int argc, char *argv[]) {
u4 += u4_diff;
s4_old = *s4_new;
+ if (parset.get<bool>("printEvolution")) {
+ // Print displacement on frictional boundary
+ boost::format const formatter("s[%02d] = %+3e, "
+ "%|40t|u[%02d] = %+3e");
+ for (size_t i = 0; i < frictionalNodes.size(); ++i)
+ if (frictionalNodes[i][0]) {
+ std::cout << boost::format(formatter) % i % (*s4_new)[i] % i % u4[i]
+ << std::endl;
+ break; // Only print the first element
+ }
+ }
+
// Compute von Mises stress and write everything to a file
if (parset.get<bool>("writeVTK")) {
auto const displacement =
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index 9629b741..b91d6eda 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -2,6 +2,7 @@
timesteps = 30
verbose = false
+printEvolution = true
printFrictionalBoundary = false
printProgress = false
printDifference = false
--
GitLab