From f0ad74cb0cc4f85bffc0f986a85588ce5f21f8d8 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 22 Mar 2012 14:33:05 +0100
Subject: [PATCH] Write coefficient of friction to file
---
src/one-body-sample.cc | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index d51ffa4b..fbe1792b 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -208,6 +208,7 @@ int main(int argc, char *argv[]) {
ignoreNodes);
std::fstream octave_writer("data", std::fstream::out);
+ std::fstream coefficient_writer("coefficient", std::fstream::out);
timer.reset();
auto const timesteps = parset.get<size_t>("timesteps");
octave_writer << "# name: A" << std::endl << "# type: matrix" << std::endl
@@ -298,8 +299,9 @@ int main(int argc, char *argv[]) {
auto const eta = parset.get<double>("boundary.friction.eta");
auto const mu = parset.get<double>("boundary.friction.mu");
- std::cout << mu + a *std::log(V * eta) +
- b * (state - std::log(eta * L)) << std::endl;
+ coefficient_writer << mu + a *std::log(V * eta) +
+ b * (state - std::log(eta * L))
+ << std::endl;
break;
}
@@ -328,6 +330,7 @@ int main(int argc, char *argv[]) {
<< timer.elapsed() << "s" << std::endl;
octave_writer.close();
+ coefficient_writer.close();
if (parset.get<bool>("printFrictionalBoundary")) {
// Print displacement on frictional boundary
--
GitLab