From 0aca93090a203817b44fefb3a8cc34df93ed2ef8 Mon Sep 17 00:00:00 2001
From: Thomas Witkowski <thomas.witkowski@gmx.de>
Date: Fri, 29 Aug 2008 08:24:31 +0000
Subject: [PATCH] * Nothing important
---
AMDiS/src/Estimator.cc | 1 -
AMDiS/src/ProblemVec.cc | 3 +++
2 files changed, 3 insertions(+), 1 deletion(-)
diff --git a/AMDiS/src/Estimator.cc b/AMDiS/src/Estimator.cc
index e50da03d..5550158d 100755
--- a/AMDiS/src/Estimator.cc
+++ b/AMDiS/src/Estimator.cc
@@ -21,7 +21,6 @@ namespace AMDiS {
TraverseStack stack;
ElInfo *elInfo = stack.traverseFirst(mesh, -1, traverseFlag);
while (elInfo) {
- std::cout << elInfo->getElement()->getIndex() << std::endl;
estimateElement(elInfo);
elInfo = stack.traverseNext(elInfo);
}
diff --git a/AMDiS/src/ProblemVec.cc b/AMDiS/src/ProblemVec.cc
index 398f78bf..f0be8cca 100644
--- a/AMDiS/src/ProblemVec.cc
+++ b/AMDiS/src/ProblemVec.cc
@@ -552,9 +552,12 @@ namespace AMDiS {
for (int i = 0; i < nComponents; i++) {
DOFVector<double> *tmp = NEW DOFVector<double>(componentSpaces[i], "tmp");
tmp->interpol(exactSolutionFcts[i]);
+ double t = tmp->max();
*tmp -= *(solution_->getDOFVector(i));
double l2Error = tmp->L2Norm();
+ double maxError = tmp->max() / t;
MSG("L2 error = %.8e\n", l2Error);
+ MSG("Max error = %.8e\n", maxError);
DELETE tmp;
}
} else {
--
GitLab