From 1029fce483b50650b756ed6ddfe5d4d8d01a5c50 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 27 Apr 2010 14:30:18 +0000
Subject: [PATCH] out-comment the use of method computeGeodesicDifference. It
currently doesn't work for all target spaces
[[Imported from SVN: r5959]]
---
harmonicmaps.cc | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/harmonicmaps.cc b/harmonicmaps.cc
index 2f7111e8..c7915287 100644
--- a/harmonicmaps.cc
+++ b/harmonicmaps.cc
@@ -2,8 +2,8 @@
#include <fenv.h>
-//#define UNITVECTOR2
-#define UNITVECTOR3
+#define UNITVECTOR2
+//#define UNITVECTOR3
//#define ROTATION2
//#define REALTUPLE1
@@ -278,7 +278,8 @@ int main (int argc, char *argv[]) try
// Compute error of the initial iterate
typedef BlockVector<FieldVector<double,blocksize> > DifferenceType;
- DifferenceType geodesicDifference = computeGeodesicDifference(exactSolution, initialIterate);
+#warning computeGeodesicDifference outcommented
+ DifferenceType geodesicDifference = DifferenceType(0);//computeGeodesicDifference(exactSolution, initialIterate);
double oldError = std::sqrt(EnergyNorm<BCRSMatrix<FieldMatrix<double, blocksize, blocksize> >, BlockVector<FieldVector<double,blocksize> > >::normSquared(geodesicDifference, hessian));
int i;
@@ -303,7 +304,8 @@ int main (int argc, char *argv[]) try
// Compute error
// /////////////////////////////////////////////////////
- geodesicDifference = computeGeodesicDifference(exactSolution, intermediateSolution);
+#warning computeGeodesicDifference outcommented
+ geodesicDifference = DifferenceType(0);//computeGeodesicDifference(exactSolution, intermediateSolution);
error = std::sqrt(EnergyNorm<BCRSMatrix<FieldMatrix<double, blocksize, blocksize> >, BlockVector<FieldVector<double,blocksize> > >::normSquared(geodesicDifference, hessian));
--
GitLab