diff --git a/FunctionSpace/FunctionSpaceEdge.cpp b/FunctionSpace/FunctionSpaceEdge.cpp
index f7dc3b54ac08488cb48918b8e963870ba15dd669..dd3bd0dda097cff1a77147e81c8bb49ba2f6c044 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);
 }