diff --git a/src/numerics2.hpp b/src/numerics2.hpp
index ed8db7e0e..07d611091 100644
--- a/src/numerics2.hpp
+++ b/src/numerics2.hpp
@@ -19,6 +19,7 @@
 #ifndef NUMERICS_HPP
 #define NUMERICS_HPP
 
+#include <array>
 #include <cmath>
 #include <limits>
 #include <cstddef>
@@ -64,6 +65,12 @@ bool is_finite(const double v[N])
    return is_finite;
 }
 
+template <typename T, std::size_t N>
+bool is_finite(const std::array<T, N>& v)
+{
+   return is_finite<N>(&v[0]);
+}
+
 } // namespace flexiblesusy
 
 #endif