diff --git a/AMDiS/src/Estimator.cc b/AMDiS/src/Estimator.cc index e50da03d12791363677f8f2c2c0ba0c685da9078..5550158d397422adc7e1f4798f66c0093ef6ad76 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 398f78bfcccf5ac027c646dbe73ecdc0f51b64a5..f0be8cca313e6b6c24c0aa6776300df959857074 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 {