diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
index 7acb4d4d1a92613292df0b07cd9ecf218609f411..293f7c59dd4d13c0dab2b343ed47b6d82974240d 100644
--- a/src/CMakeLists.txt
+++ b/src/CMakeLists.txt
@@ -1,5 +1,2 @@
 add_executable("C1Elements" C1Elements.cc)
 target_link_dune_default_libraries("C1Elements")
-
-add_executable("testMIBasis" testMIBasis.cc)
-target_link_dune_default_libraries("testMIBasis")
\ No newline at end of file
diff --git a/src/testMIBasis.cc b/src/testMIBasis.cc
deleted file mode 100644
index ebd87f6f0c4a9d735dc3f42b56c054a6555eab3b..0000000000000000000000000000000000000000
--- a/src/testMIBasis.cc
+++ /dev/null
@@ -1,49 +0,0 @@
-#include <fstream>
-#include <utility>
-
-#include <dune/geometry/type.hh>
-
-#include <dune/localfunctions/utility/polynomialbasis.hh>
-#include <dune/localfunctions/hermite.hh>
-
-#include <dune/localfunctions/argyris.hh>
-
-int main(int argc, char *argv[])
-{
-  constexpr unsigned int dim = 2;
-  typedef Dune::MultiIndex<dim, double> MI;
-  typedef Dune::MonomialBasis<Dune::GeometryTypes::simplex(dim).toId(), MI> MIBasis;
-
-  MIBasis basis(5);
-  Dune::FieldVector<MI, dim> x;
-  /*
-       * Init MultiIndices
-       * x[0]=(1,0,0) x
-       * x[1]=(0,1,0) y
-       * x[2]=(0,0,1) z
-       */
-  for (unsigned int i = 0; i < dim; ++i)
-    x[i].set(i, 1);
-  std::vector<MI> val(basis.size());
-
-  // val now contains all monomials in $n$ dimensions with degree $\leq order+1$
-  basis.evaluate(x, val);
-  std::cout << "Monomial Basis " << std::endl
-            << val << std::endl;
-
-  Dune::Impl::HermiteLocalBasis<MI, MI, 2> hermite;
-  val.resize(hermite.size());
-  hermite.evaluate(x, val);
-  std::cout << "Hermite Basis 2d " << std::endl << val << std::endl;
-
-  Dune::Impl::ArgyrisLocalBasis<MI, MI> argyris;
-  val.resize(argyris.size());
-  argyris.evaluate(x, val);
-  std::cout << "Argyris Basis " << std::endl << val << std::endl;
-
-  std::vector<decltype(argyris)::Traits::JacobianType> jacs;
-  argyris.evaluateJacobian(x, jacs);
-  std::cout << "Argyris Basis Jacobian " << std::endl << jacs << std::endl;
-
-  return 0;
-}