diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index cb9734e97751680de83ffb6d00dba537bae8c23f..203595dcc4e42ba61db9097e45306f798229f57a 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -212,6 +212,7 @@ int main (int argc, char *argv[]) try
     typedef TargetSpace::rebind<adouble>::other ATargetSpace;
     typedef LocalGeodesicFEFunction<dim, double, FEBasis::LocalView::Tree::FiniteElement, ATargetSpace> LocalInterpolationRule;
     //typedef GFE::LocalProjectedFEFunction<dim, double, FEBasis::LocalView::Tree::FiniteElement, ATargetSpace> LocalInterpolationRule;
+    std::cout << "Using local interpolation: " << className<LocalInterpolationRule>() << std::endl;
 
     // Assembler using ADOL-C
     std::shared_ptr<LocalGeodesicFEStiffness<FEBasis,ATargetSpace> > localEnergy;