diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index fa4b1f8a80da9b8962df0cdd5192f18c101f62d5..94984e3e44eb6782d2d0ad983eb7fa6be26b5d1a 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -189,7 +189,7 @@ int main (int argc, char *argv[]) try
                                   FEBasis::LocalFiniteElement,
                                   TargetSpace> localGFEADOLCStiffness(localEnergy.get());
 
-    GeodesicFEAssembler<FEBasis,TargetSpace> assembler(grid->leafGridView(), &localGFEADOLCStiffness);
+    GeodesicFEAssembler<FEBasis,TargetSpace> assembler(feBasis, &localGFEADOLCStiffness);
 
     // /////////////////////////////////////////////////
     //   Create a Riemannian trust-region solver