diff --git a/dune/gfe/localgfetestfunctionbasis.hh b/dune/gfe/localgfetestfunctionbasis.hh
index 6e786c8c36044811385cdb969a919db5beef1f7f..6cd71292ea4b043e6c2d6767a6043779f35c6749 100644
--- a/dune/gfe/localgfetestfunctionbasis.hh
+++ b/dune/gfe/localgfetestfunctionbasis.hh
@@ -160,7 +160,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateFunction
         /** \todo This call internally keeps computing the value of the gfe function at 'local'.
          * This is expensive.  Eventually we should precompute it once and reused the result. */
         localGFEFunction_.evaluateDerivativeOfValueWRTCoefficient (local, i, derivative);
-        Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficients_[i].orthonormalFrame();
+        Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficient(i).orthonormalFrame();
         
         for (int j=0; j<spaceDim; j++)
             derivative.mv(basisVectors[j], out[i][j]);
@@ -182,7 +182,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateJacobian
         Tensor3< double, embeddedDim, embeddedDim, dim > derivative;
         localGFEFunction_.evaluateDerivativeOfGradientWRTCoefficient (local, i, derivative);
 
-        Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficients_[i].orthonormalFrame();
+        Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficient(i).orthonormalFrame();
         
         for (int j=0; j<spaceDim; j++) {