From 97b0f1132298fcc90f8c0d1b26426e482db72619 Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <nicolas.marsic@gmail.com>
Date: Fri, 9 Nov 2012 14:52:44 +0000
Subject: [PATCH] FunctionSpaceEdge -- Maps at the End
---
FunctionSpace/FunctionSpaceEdge.cpp | 10 ++++------
1 file changed, 4 insertions(+), 6 deletions(-)
diff --git a/FunctionSpace/FunctionSpaceEdge.cpp b/FunctionSpace/FunctionSpaceEdge.cpp
index f7dc3b54ac..dd3bd0dda0 100644
--- a/FunctionSpace/FunctionSpaceEdge.cpp
+++ b/FunctionSpace/FunctionSpaceEdge.cpp
@@ -51,15 +51,14 @@ interpolate(const MElement& element,
for(unsigned int i = 0; i < nFun; i++){
fullVector<double> vi =
- Mapper::grad(Polynomial::at(*fun[i], uvw[0], uvw[1], uvw[2]),
- invJac);
+ Polynomial::at(*fun[i], uvw[0], uvw[1], uvw[2]);
vi.scale(coef[i]);
val.axpy(vi, 1);
}
// Return Interpolated Value //
- return val;
+ return Mapper::grad(val, invJac);
}
fullVector<double> FunctionSpaceEdge::
@@ -88,13 +87,12 @@ interpolateInRefSpace(const MElement& element,
for(unsigned int i = 0; i < nFun; i++){
fullVector<double> vi =
- Mapper::grad(Polynomial::at(*fun[i], uvw(0), uvw(1), uvw(2)),
- invJac);
+ Polynomial::at(*fun[i], uvw(0), uvw(1), uvw(2));
vi.scale(coef[i]);
val.axpy(vi, 1);
}
// Return Interpolated Value //
- return val;
+ return Mapper::grad(val, invJac);
}
--
GitLab