diff --git a/src/timestepping.cc b/src/timestepping.cc
index 0003fd227d924231da6a8ce7ce8e347519ce8e07..af9edb2a07ea6e2b0df86024c6da1ab563f5c253 100644
--- a/src/timestepping.cc
+++ b/src/timestepping.cc
@@ -3,7 +3,6 @@
 #endif
 
 #include <dune/istl/matrixindexset.hh>
-#include <dune/fufem/arithmetic.hh>
 
 #include "timestepping.hh"
 
diff --git a/src/timestepping/backward_euler.cc b/src/timestepping/backward_euler.cc
index a39f0e5415a0ec0faf028badcd489a1444e51477..35cabff66d7ebe59b6296d6d2eceaf9626546c11 100644
--- a/src/timestepping/backward_euler.cc
+++ b/src/timestepping/backward_euler.cc
@@ -1,3 +1,5 @@
+#include <dune/solvers/common/arithmetic.hh>
+
 template <class Vector, class Matrix, class Function, size_t dim>
 BackwardEuler<Vector, Matrix, Function, dim>::BackwardEuler(
     Matrix const &_A, Matrix const &_M, Matrix const &_C,
diff --git a/src/timestepping/newmark.cc b/src/timestepping/newmark.cc
index 51eeb7a36f79283e4fd9ed928203b976499c7bf4..12e39abcac299a3ec2f6455b6da211b23a0b6495 100644
--- a/src/timestepping/newmark.cc
+++ b/src/timestepping/newmark.cc
@@ -1,3 +1,5 @@
+#include <dune/solvers/common/arithmetic.hh>
+
 template <class Vector, class Matrix, class Function, size_t dim>
 Newmark<Vector, Matrix, Function, dim>::Newmark(
     Matrix const &_A, Matrix const &_M, Matrix const &_C,