From 0857c0dc36267158369f53dca164e38029603652 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 25 Nov 2011 11:24:35 +0000
Subject: [PATCH] print estimated edge size instead of number of vertices (was
buggy anyways)
[[Imported from SVN: r8259]]
---
harmonicmaps-eoc.cc | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/harmonicmaps-eoc.cc b/harmonicmaps-eoc.cc
index 3867b8da..ffd45a25 100644
--- a/harmonicmaps-eoc.cc
+++ b/harmonicmaps-eoc.cc
@@ -241,7 +241,7 @@ int main (int argc, char *argv[]) try
// ///////////////////////////////////////////////////////////
std::ofstream logFile("harmonicmaps-eoc.results");
- logFile << "# vertices max-norm, L2-norm, h1-seminorm" << std::endl;
+ logFile << "# mesh size, max-norm, L2-norm, h1-seminorm" << std::endl;
for (int i=1; i<numLevels; i++) {
@@ -299,7 +299,7 @@ int main (int argc, char *argv[]) try
H1SemiNorm< BlockVector<TargetSpace::CoordinateType> > l2Norm(massMatrix);
// Compute max-norm difference
- std::cout << "Vertices: " << xEmbedded.size() << std::endl;
+ std::cout << "h: " << std::pow(0.5, i-1) << std::endl;
std::cout << "Level: " << i-1
<< ", max-norm error: " << difference.infinity_norm()
<< std::endl;
@@ -312,7 +312,7 @@ int main (int argc, char *argv[]) try
<< ", H1 error: " << h1Norm(difference)
<< std::endl;
- logFile << xEmbedded.size() << " " << difference.infinity_norm()
+ logFile << std::pow(0.5, i-1) << " " << difference.infinity_norm()
<< " " << l2Norm(difference)
<< " " << h1Norm(difference)
<< std::endl;
--
GitLab