diff --git a/ql/experimental/math/laplaceinterpolation.cpp b/ql/experimental/math/laplaceinterpolation.cpp
index 278ee534757..93748bd9a65 100644
--- a/ql/experimental/math/laplaceinterpolation.cpp
+++ b/ql/experimental/math/laplaceinterpolation.cpp
@@ -209,8 +209,7 @@ namespace QuantLib {
         }
     }
 
-    template <class M>
-    void laplaceInterpolation(M& A,
+    void laplaceInterpolation(Matrix& A,
                               const std::vector<Real>& x,
                               const std::vector<Real>& y,
                               Real relTol) {
@@ -243,11 +242,4 @@ namespace QuantLib {
         }
     }
 
-    // template instantiations for matrix classes we want to support
-
-    template void laplaceInterpolation(Matrix& A,
-                                       const std::vector<Real>& x,
-                                       const std::vector<Real>& y,
-                                       Real relTol);
-
 } // namespace QuantLib
diff --git a/ql/experimental/math/laplaceinterpolation.hpp b/ql/experimental/math/laplaceinterpolation.hpp
index 6ce1d8cc8d3..ee0745802eb 100644
--- a/ql/experimental/math/laplaceinterpolation.hpp
+++ b/ql/experimental/math/laplaceinterpolation.hpp
@@ -25,6 +25,7 @@
 #define quantlib_laplace_interpolation
 
 #include <ql/math/array.hpp>
+#include <ql/math/matrix.hpp>
 #include <ql/shared_ptr.hpp>
 #include <ql/types.hpp>
 
@@ -64,8 +65,7 @@ namespace QuantLib {
     /*! Convenience function that Laplace-interpolates null values in a given matrix.
         If the x or y grid or both are not given, an equidistant grid is assumed. */
 
-    template <class M>
-    void laplaceInterpolation(M& A,
+    void laplaceInterpolation(Matrix& A,
                               const std::vector<Real>& x = {},
                               const std::vector<Real>& y = {},
                               Real relTol = 1E-6);