diff --git a/src/language/modules/SchemeModule.cpp b/src/language/modules/SchemeModule.cpp
index f2570a3c30b911990e9d8d2c12f7bc2de673f4a6..b9b789cfa8e3d7bc5328062c6085217dbe142120 100644
--- a/src/language/modules/SchemeModule.cpp
+++ b/src/language/modules/SchemeModule.cpp
@@ -65,9 +65,6 @@
 #include <utils/checkpointing/WriteIQuadratureDescriptor.hpp>
 #include <utils/checkpointing/WriteVariableBCDescriptor.hpp>
 
-#warning REMOVE
-#include <scheme/test_reconstruction.hpp>
-
 #include <memory>
 
 SchemeModule::SchemeModule()
@@ -759,12 +756,6 @@ SchemeModule::SchemeModule()
 
                                               ));
 
-  this->_addBuiltinFunction("test_reconstruction", std::function{[](std::shared_ptr<const DiscreteFunctionVariant> df,
-                                                                    size_t degree, size_t number) -> void {
-                              //
-                              test_reconstruction(df, degree, number);
-                            }});
-
   MathFunctionRegisterForVh{*this};
 }
 
diff --git a/src/scheme/CMakeLists.txt b/src/scheme/CMakeLists.txt
index 5a555c426c7f1bda6d1859bd15d49a410a7dabf2..852a6cd751d092a7ab1c59437c1e549ae2f8094f 100644
--- a/src/scheme/CMakeLists.txt
+++ b/src/scheme/CMakeLists.txt
@@ -14,8 +14,6 @@ add_library(
   HyperelasticSolver.cpp
   LoadBalancer.cpp
   PolynomialReconstruction.cpp
-
-  test_reconstruction.cpp
 )
 
 target_link_libraries(
diff --git a/src/scheme/test_reconstruction.cpp b/src/scheme/test_reconstruction.cpp
deleted file mode 100644
index 1353ab46149f333ca4ba2cdf8e72f2011932919d..0000000000000000000000000000000000000000
--- a/src/scheme/test_reconstruction.cpp
+++ /dev/null
@@ -1,46 +0,0 @@
-#include <scheme/test_reconstruction.hpp>
-
-#include <scheme/PolynomialReconstruction.hpp>
-#include <utils/Timer.hpp>
-
-void
-test_reconstruction(std::shared_ptr<const DiscreteFunctionVariant> df, size_t degree, size_t number)
-{
-  if (degree == 1) {
-    PolynomialReconstructionDescriptor descriptor{IntegrationMethodType::cell_center, degree};
-    PolynomialReconstruction rec_builder{descriptor};
-
-    Timer t;
-
-    for (size_t i = 0; i < number; ++i) {
-      auto rec = rec_builder.build(df);
-    }
-
-    std::cout << "*** Elapsed time for " << number << " cell center reconstructions: " << t.seconds() << "s ***\n";
-  }
-
-  {
-    PolynomialReconstructionDescriptor descriptor{IntegrationMethodType::boundary, degree};
-    PolynomialReconstruction rec_builder{descriptor};
-
-    Timer t;
-
-    for (size_t i = 0; i < number; ++i) {
-      auto rec = rec_builder.build(df);
-    }
-
-    std::cout << "*** Elapsed time for " << number << " boundary reconstructions: " << t.seconds() << "s ***\n";
-  }
-  {
-    PolynomialReconstructionDescriptor descriptor{IntegrationMethodType::element, degree};
-    PolynomialReconstruction rec_builder{descriptor};
-
-    Timer t;
-
-    for (size_t i = 0; i < number; ++i) {
-      auto rec = rec_builder.build(df);
-    }
-
-    std::cout << "*** Elapsed time for " << number << " element reconstructions: " << t.seconds() << "s ***\n";
-  }
-}
diff --git a/src/scheme/test_reconstruction.hpp b/src/scheme/test_reconstruction.hpp
deleted file mode 100644
index 36ed481ec14f5edc34a677e3fc16d838e1bb6cf1..0000000000000000000000000000000000000000
--- a/src/scheme/test_reconstruction.hpp
+++ /dev/null
@@ -1,8 +0,0 @@
-#ifndef TEST_RECONSTRUCTION_HPP
-#define TEST_RECONSTRUCTION_HPP
-
-#include <scheme/DiscreteFunctionVariant.hpp>
-#warning REMOVE
-void test_reconstruction(std::shared_ptr<const DiscreteFunctionVariant> df, size_t degree, size_t number);
-
-#endif   // TEST_RECONSTRUCTION_HPP