Add output to Array. Now Vector uses Array output to print data.

src/algebra/Vector.hpp

@@ -2,7 +2,6 @@
 #define VECTOR_HPP
 
 #include <utils/Array.hpp>
-#include <utils/NaNHelper.hpp>
 #include <utils/PugsAssert.hpp>
 #include <utils/PugsMacros.hpp>
 #include <utils/PugsUtils.hpp>
@@ -29,13 +28,7 @@ class Vector   // LCOV_EXCL_LINE
   friend std::ostream&
   operator<<(std::ostream& os, const Vector& x)
   {
-    if (x.size() > 0) {
-      os << 0 << ':' << NaNHelper(x[0]);
-    }
-    for (size_t i = 1; i < x.size(); ++i) {
-      os << ' ' << i << ':' << NaNHelper(x[i]);
-    }
-    return os;
+    return os << x.m_values;
   }
 
   friend Vector<std::remove_const_t<DataType>>

src/utils/Array.hpp

@@ -2,6 +2,7 @@
 #define ARRAY_HPP
 
 #include <utils/InvalidData.hpp>
+#include <utils/NaNHelper.hpp>
 #include <utils/PugsAssert.hpp>
 #include <utils/PugsMacros.hpp>
 #include <utils/PugsUtils.hpp>
@@ -23,6 +24,17 @@ class [[nodiscard]] Array
   friend Array<std::add_const_t<DataType>>;
 
  public:
+  friend std::ostream& operator<<(std::ostream& os, const Array& x)
+  {
+    if (x.size() > 0) {
+      os << 0 << ':' << NaNHelper(x[0]);
+    }
+    for (size_t i = 1; i < x.size(); ++i) {
+      os << ' ' << i << ':' << NaNHelper(x[i]);
+    }
+    return os;
+  }
+
   PUGS_INLINE size_t size() const noexcept
   {
     return m_values.extent(0);