array: add .fill()

src/array.hpp

@@ -419,6 +419,15 @@ public:
             static_assert(Shape::size == detail::ShapeType<Dims...>::size, "incompatible shapes");
             return {data(), Base::size(), shape};
         }
+
+    /** Fill with value */
+    void fill(const Scalar& value)
+        {
+            const size_t n = Base::size();
+
+            for (size_t i = 0; i < n; ++i)
+                data_[i] = value;
+        }
 };
 
 /**