From 1d93573e28e7698bad81bd77ab10d85a63ba9dfc Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 9 Feb 2010 16:21:56 +0000
Subject: [PATCH] remove code for interpolation via recursive geodesics. It is
invariant under permutation of the vertices
[[Imported from SVN: r5511]]
---
src/localgeodesicfefunction.hh | 21 ---------------------
1 file changed, 21 deletions(-)
diff --git a/src/localgeodesicfefunction.hh b/src/localgeodesicfefunction.hh
index 744d7c12..439bc063 100644
--- a/src/localgeodesicfefunction.hh
+++ b/src/localgeodesicfefunction.hh
@@ -48,27 +48,6 @@ template <int dim, class ctype, class TargetSpace>
TargetSpace LocalGeodesicFEFunction<dim,ctype,TargetSpace>::
evaluate(const Dune::FieldVector<ctype, dim>& local)
{
-#if 0 // Interpolation using recursive geodesic cones, doesn't work :-(
- ctype extraCoord = 1;
- for (int i=0; i<dim; i++)
- extraCoord -= local[i];
-
- ctype normalizingFactor = extraCoord;
-
- TargetSpace result = coefficients_[0];
-
- for (int i=0; i<dim; i++) {
- assert(local[i]>=0);
- normalizingFactor += local[i];
- result = TargetSpace::interpolate(result, coefficients_[i+1],
- (std::abs(normalizingFactor) > 1e-10)
- ? local[i] / normalizingFactor
- : 0);
- }
-
- return result;
-#endif
-
#warning Some code out-commented
#if 0
Dune::FieldVector<ctype, dim+1> barycentricCoordinates;
--
GitLab