From 66c1a8d020ff3a4fe4c8028159a51e81f9581c93 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 28 Aug 2007 12:43:47 +0000
Subject: [PATCH] write total convergence rate to file. Read explicitly given
parameter file, if requested
[[Imported from SVN: r1616]]
---
dirneucoupling.cc | 12 ++++++++----
1 file changed, 8 insertions(+), 4 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 6513e8b6..b72c48a3 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -60,7 +60,10 @@ int main (int argc, char *argv[]) try
// parse data file
ConfigParser parameterSet;
- parameterSet.parseFile("dirneucoupling.parset");
+ if (argc==2)
+ parameterSet.parseFile(argv[1]);
+ else
+ parameterSet.parseFile("dirneucoupling.parset");
// read solver settings
const int numLevels = parameterSet.get<int>("numLevels");
@@ -546,6 +549,10 @@ int main (int argc, char *argv[]) try
<< " convRate: " << std::pow(totalConvRate[i+1-backTrace], 1/((double)i+1-backTrace))
<< std::endl;
+ std::ofstream convRateFile("convrate");
+ convRateFile << damping << " " << std::pow(totalConvRate[i+1-backTrace], 1/((double)i+1-backTrace))
+ << std::endl;
+
// //////////////////////////////
// Output result
@@ -556,9 +563,6 @@ int main (int argc, char *argv[]) try
writeRod(rodX, "rod3d.result");
-// for (int i=0; i<rodX.size(); i++)
-// std::cout << rodX[i] << std::endl;
-
} catch (Exception e) {
std::cout << e << std::endl;
--
GitLab