diff --git a/src/cosserat-continuum.cc b/src/cosserat-continuum.cc
index 94075e484c1e8bd02f36ca76cc2eafbf9c83b9fc..e6a29062f94fe283cacba18dbaf04c5c11373253 100644
--- a/src/cosserat-continuum.cc
+++ b/src/cosserat-continuum.cc
@@ -8,17 +8,10 @@
 #include <adolc/drivers/drivers.h>    // use of "Easy to Use" drivers
 #include <adolc/taping.h>
 
+#include <dune/fufem/adolc.hh>
 #include <dune/fufem/utilities/adolcnamespaceinjections.hh>
 
 #include <dune/common/typetraits.hh>
-namespace Dune {
-  template <>
-  struct IsNumber<adouble>
-  {
-    constexpr static bool value = true;
-  };
-}
-
 #include <dune/common/bitsetvector.hh>
 #include <dune/common/parametertree.hh>
 #include <dune/common/parametertreeparser.hh>
diff --git a/src/gradient-flow.cc b/src/gradient-flow.cc
index 875b800cab2d1a6e78608b79116467ac5a8739e1..01ae6d7b14b3270b57cd8975b2dd87e3f615297c 100644
--- a/src/gradient-flow.cc
+++ b/src/gradient-flow.cc
@@ -2,22 +2,15 @@
 // vi: set et ts=4 sw=2 sts=2:
 #include <config.h>
 
+#include <array>
+
 // Includes for the ADOL-C automatic differentiation library
 // Need to come before (almost) all others.
 #include <adolc/adouble.h>
+#include <dune/fufem/adolc.hh>
 #include <dune/fufem/utilities/adolcnamespaceinjections.hh>
 
 #include <dune/common/typetraits.hh>
-namespace Dune {
-  template <>
-  struct IsNumber<adouble>
-  {
-    constexpr static bool value = true;
-  };
-}
-
-#include <array>
-
 #include <dune/common/bitsetvector.hh>
 #include <dune/common/parametertree.hh>
 #include <dune/common/parametertreeparser.hh>
diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 6cbf4956b3ae5b9956d4e0ccaaf7c287ca74bcd6..0792f4ff7a16c02ecb6d6715dd2766a9e7e34519 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -1,22 +1,15 @@
 #include <config.h>
 
 #include <fenv.h>
+#include <array>
 
 // Includes for the ADOL-C automatic differentiation library
 // Need to come before (almost) all others.
 #include <adolc/adouble.h>
+#include <dune/fufem/adolc.hh>
 #include <dune/fufem/utilities/adolcnamespaceinjections.hh>
 
 #include <dune/common/typetraits.hh>
-namespace Dune {
-  template <>
-  struct IsNumber<adouble>
-  {
-    constexpr static bool value = true;
-  };
-}
-
-#include <array>
 
 #include <dune/common/bitsetvector.hh>
 #include <dune/common/parametertree.hh>
diff --git a/src/mixed-cosserat-continuum.cc b/src/mixed-cosserat-continuum.cc
index 9d6ef7178f77f210bc67421a45cba5f95d4bde6b..1a7e2e80b9f57abe84cc72e5ed24ebce8b564bc3 100644
--- a/src/mixed-cosserat-continuum.cc
+++ b/src/mixed-cosserat-continuum.cc
@@ -1,6 +1,7 @@
 #include <config.h>
 
 #include <fenv.h>
+#include <array>
 
 // Includes for the ADOL-C automatic differentiation library
 // Need to come before (almost) all others.
@@ -8,19 +9,11 @@
 #include <adolc/drivers/drivers.h>    // use of "Easy to Use" drivers
 #include <adolc/taping.h>
 
+#include <dune/fufem/adolc.hh>
 #include <dune/fufem/utilities/adolcnamespaceinjections.hh>
 
 
 #include <dune/common/typetraits.hh>
-namespace Dune {
-  template <>
-  struct IsNumber<adouble>
-  {
-    constexpr static bool value = true;
-  };
-}
-#include <array>
-
 #include <dune/common/bitsetvector.hh>
 #include <dune/common/parametertree.hh>
 #include <dune/common/parametertreeparser.hh>