Make array isapprox more generic by avoiding zip/explicit iteration (…

…#44893)

* Avoid explicit array iteration in isapprox

This is more generic and works for, e.g., gpu arrays

stdlib/LinearAlgebra/src/generic.jl

@@ -1774,7 +1774,8 @@ function isapprox(x::AbstractArray, y::AbstractArray;
         return d <= max(atol, rtol*max(norm(x), norm(y)))
     else
         # Fall back to a component-wise approximate comparison
-        return all(ab -> isapprox(ab[1], ab[2]; rtol=rtol, atol=atol, nans=nans), zip(x, y))
+        # (mapreduce instead of all for greater generality [#44893])
+        return mapreduce((a, b) -> isapprox(a, b; rtol=rtol, atol=atol, nans=nans), &, x, y)
     end
 end