Fix comprod! to use mul_prod instead of add_sum (#26905)

base/accumulate.jl

@@ -125,7 +125,7 @@ Cumulative product of `A` along the dimension `dims`, storing the result in `B`.
 See also [`cumprod`](@ref).
 """
 cumprod!(B::AbstractArray{T}, A; dims::Integer) where {T} =
- accumulate!(add_sum, B, A, dims=dims)
+ accumulate!(mul_prod, B, A, dims=dims)
 
 """
  cumprod!(y::AbstractVector, x::AbstractVector)

test/arrayops.jl

@@ -2071,6 +2071,16 @@ end
 @test cumsum([1 2; 3 4], dims=2) == [1 3; 3 7]
 @test cumsum([1 2; 3 4], dims=3) == [1 2; 3 4]
 
+@test cumprod!(Vector{Int}(undef, 1), [5], dims=2) == [5]
+@test cumprod!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=3) == [1 2; 3 4]
+@test cumprod!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=1) == [1 2; 3 8]
+@test cumprod!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=2) == [1 2; 3 12]
+
+@test cumsum!(Vector{Int}(undef, 1), [5], dims=2) == [5]
+@test cumsum!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=1) == [1 2; 4 6]
+@test cumsum!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=2) == [1 3; 3 7]
+@test cumsum!(Matrix{Int}(undef, 2, 2), [1 2; 3 4], dims=3) == [1 2; 3 4]
+
 # issue #18363
 @test_throws DimensionMismatch cumsum!([0,0], 1:4)
 @test cumsum(Any[])::Vector{Any} == Any[]