diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 6a2335f790d64449bccf9a71a1f01fc09317ffc6..1e5d2d20abbdb16a9a4acae04f6e685c50e376ab 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -125,7 +125,7 @@ int main (int argc, char *argv[]) try
     //  Construct the scalar function space basis corresponding to the GFE space
     //////////////////////////////////////////////////////////////////////////////////
 
-    typedef P2NodalBasis<typename GridType::LeafGridView,double> FEBasis;
+    typedef P1NodalBasis<typename GridType::LeafGridView,double> FEBasis;
     FEBasis feBasis(grid->leafGridView());
 
     SolutionType x(feBasis.size());