diff --git a/include/Eisvogel/Kernels.hh b/include/Eisvogel/Kernels.hh
index 2aecf0a62..7df8b2b1a 100644
--- a/include/Eisvogel/Kernels.hh
+++ b/include/Eisvogel/Kernels.hh
@@ -5,13 +5,6 @@
 #include <cmath>
 #include "Common.hh"
 
-struct Kernel {
-  
-  virtual std::size_t Support() const = 0;
-  virtual scalar_t operator()(scalar_t arg) const = 0;
-  virtual scalar_t CDF(int arg) const = 0;
-};
-
 struct KeysCubicInterpolationKernelNew {
 
   const static int Support = 2;
@@ -52,11 +45,4 @@ struct KeysCubicInterpolationKernelNew {
   }
 };
 
-struct KeysCubicInterpolationKernel : public Kernel {
-
-  std::size_t Support() const;
-  scalar_t operator()(scalar_t arg) const;
-  scalar_t CDF(int arg) const;
-};
-
 #endif
diff --git a/src/Kernels.cxx b/src/Kernels.cxx
deleted file mode 100644
index 4c3770620..000000000
--- a/src/Kernels.cxx
+++ /dev/null
@@ -1,41 +0,0 @@
-#include "Eisvogel/Kernels.hh"
-#include <cmath>
-
-std::size_t KeysCubicInterpolationKernel::Support() const {
-  return 2;
-}
-
-scalar_t KeysCubicInterpolationKernel::operator()(scalar_t arg) const {
-  scalar_t abs_arg = std::fabs(arg);
-
-  if(abs_arg < 1.0) {
-    return 1.0 + abs_arg * abs_arg * (-2.5 + 1.5 * abs_arg);
-  }
-  else if(abs_arg < 2.0) {
-    return 2.0 + abs_arg * (-4.0 + (2.5 - 0.5 * abs_arg) * abs_arg);
-  }
-  else {
-    return 0.0;
-  }
-}
-
-scalar_t KeysCubicInterpolationKernel::CDF(int arg) const {
-
-  if(arg <= -2) {
-    return 0.0;
-  }
-  else if(arg == -1) {
-    return -0.0416667;
-  }
-  else if(arg == 0) {
-    return 0.5;
-  }
-  else if(arg == 1) {
-    return 1.04167;
-  }
-  else if(arg >= 2) {
-    return 1.0;
-  }
-
-  return 1.0;
-}