diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index d870f934b7c6de24b623d97834d271d3764ffdfd..ce4cdb901499c8ccbea46b6b9b549439c04d2b51 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -119,7 +119,7 @@ private:
     Dune::Matrix<Dune::FieldMatrix<ctype,1,1> > computeDFdw(const TargetSpace& q) const
     {
         Dune::Matrix<Dune::FieldMatrix<ctype,1,1> > dFdw(embeddedDim,localFiniteElement_.localBasis().size());
-        for (int i=0; i<localFiniteElement_.localBasis().size(); i++) {
+        for (size_t i=0; i<localFiniteElement_.localBasis().size(); i++) {
             Dune::FieldVector<ctype,embeddedDim> tmp = TargetSpace::derivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], q);
             for (int j=0; j<embeddedDim; j++)
                 dFdw[j][i] = tmp[j];
@@ -131,7 +131,7 @@ private:
     {
         Tensor3<ctype,embeddedDim,embeddedDim,embeddedDim> result;
         result = 0;
-        for (int i=0; i<w.size(); i++)
+        for (size_t i=0; i<w.size(); i++)
             result.axpy(w[i][0], TargetSpace::thirdDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i],q));
         return result;
     }
diff --git a/dune/gfe/localgfetestfunction.hh b/dune/gfe/localgfetestfunction.hh
index 3d7fa65f764a25bee533e1e0f5c02f98faffe809..c4a6a00df7913ca6d65d24a96eaea04e703fe60f 100644
--- a/dune/gfe/localgfetestfunction.hh
+++ b/dune/gfe/localgfetestfunction.hh
@@ -107,8 +107,8 @@ void LocalGFETestFunction<dim,ctype,LocalFiniteElement,TargetSpace>::evaluateJac
         
             // Contract the second index of the derivative with the tangent vector at the i-th Lagrange point.
             // Add that to the result.
-            for (size_t k=0; k<embeddedDim; k++)
-                for (size_t l=0; l<embeddedDim; l++)
+            for (int k=0; k<embeddedDim; k++)
+                for (int l=0; l<embeddedDim; l++)
                     for (size_t m=0; m<dim; m++)
                         out[i*spaceDim+j][k][m] += derivative[k][l][m] * basisVectors[j][l];