diff --git a/dune/gfe/geodesicfeassembler.hh b/dune/gfe/geodesicfeassembler.hh
index 327176c5762917c1ac13f86a5de7a5d2c11c3c60..4b3e76f86a13c577a7985c3ca4f39af7fc51be71 100644
--- a/dune/gfe/geodesicfeassembler.hh
+++ b/dune/gfe/geodesicfeassembler.hh
@@ -78,9 +78,9 @@ getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
 
         const typename Basis::LocalFiniteElement& lfe = basis_.getLocalFiniteElement(*it);
         
-        for (int i=0; i<lfe.localBasis().size(); i++) {
+        for (size_t i=0; i<lfe.localBasis().size(); i++) {
             
-            for (int j=0; j<lfe.localBasis().size(); j++) {
+            for (size_t j=0; j<lfe.localBasis().size(); j++) {
                 
                 int iIdx = basis_.index(*it,i);
                 int jIdx = basis_.index(*it,j);
@@ -204,7 +204,7 @@ computeEnergy(const std::vector<TargetSpace>& sol) const
 
         std::vector<TargetSpace> localSolution(nDofs);
 
-        for (int i=0; i<nDofs; i++)
+        for (size_t i=0; i<nDofs; i++)
             localSolution[i] = sol[basis_.index(*it,i)];
 
         energy += localStiffness_->energy(*it, basis_.getLocalFiniteElement(*it), localSolution);