From 40e0f44de3ecf66f75567b3cba7d86a2ee2c7ddf Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 14 Aug 2007 13:05:29 +0000
Subject: [PATCH] don't print averaged convergence rate, which doesn't make
much sense for quadratic convergence
[[Imported from SVN: r1508]]
---
rod3d.cc | 8 ++------
1 file changed, 2 insertions(+), 6 deletions(-)
diff --git a/rod3d.cc b/rod3d.cc
index 0d907b4f..bf1f9424 100644
--- a/rod3d.cc
+++ b/rod3d.cc
@@ -10,7 +10,6 @@
#include "../solver/iterativesolver.hh"
-#include "../common/geomestimator.hh"
#include "../common/energynorm.hh"
#include "src/configuration.hh"
@@ -182,7 +181,6 @@ int main (int argc, char *argv[]) try
double error = std::numeric_limits<double>::max();
double oldError = 0;
- double totalConvRate = 1;
SolutionType intermediateSolution(x.size());
@@ -192,7 +190,7 @@ int main (int argc, char *argv[]) try
// from zero anyways
//oldError += computeEnergyNormSquared(exactSol3d, *hessian3d);
- /** \todo Rod error still missing */
+#warning Rod error still missing
oldError = std::sqrt(oldError);
@@ -229,12 +227,10 @@ int main (int argc, char *argv[]) try
double convRate = error / oldError;
- totalConvRate *= convRate;
// Output
std::cout << "Trust-region iteration: " << i << " error : " << error << ", "
- << "convrate " << convRate
- << " total conv rate " << std::pow(totalConvRate, 1/((double)i+1)) << std::endl;
+ << "convrate " << convRate << std::endl;
if (error < 1e-12)
break;
--
GitLab