diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index c2105b09c90018f5fbd7dab45bf18ed3c7d1a331..3fb31da5cab1f4d50b8d939de1d8970b34b842ff 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -1,8 +1,6 @@
 #ifndef HARMONIC_ENERGY_LOCAL_STIFFNESS_HH
 #define HARMONIC_ENERGY_LOCAL_STIFFNESS_HH
 
-//#define HARMONIC_ENERGY_FD_GRADIENT
-
 #include <dune/common/fmatrix.hh>
 #include <dune/grid/common/quadraturerules.hh>
 
diff --git a/harmonicmaps-eoc.cc b/harmonicmaps-eoc.cc
index 5953cae044b61a0069d065008b127881b97ca1de..8755788a4b92659b0fbb3e561b1ec2ec4c63a858 100644
--- a/harmonicmaps-eoc.cc
+++ b/harmonicmaps-eoc.cc
@@ -1,5 +1,7 @@
 #include <config.h>
 
+#define HARMONIC_ENERGY_FD_GRADIENT
+
 #include <dune/common/bitsetvector.hh>
 #include <dune/common/configparser.hh>
 
diff --git a/harmonicmaps.cc b/harmonicmaps.cc
index 97baa03c842c012e9037fdc24ee2a0befcf4477b..80627f7c0287fd2b64201e7b40c10c61fb3e1625 100644
--- a/harmonicmaps.cc
+++ b/harmonicmaps.cc
@@ -2,6 +2,8 @@
 
 #include <fenv.h>
 
+//#define HARMONIC_ENERGY_FD_GRADIENT
+
 #define UNITVECTOR2
 //#define UNITVECTOR3
 //#define ROTATION2