muladd(matrix, matrix, scalar) gives unexpected result #1162
Open
Description
This seems wrong:
julia> muladd(zeros(2,2), zeros(2,2), 3)
2×2 Matrix{Float64}:
3.0 3.0
3.0 3.0It's not clear that muladd(X, Y, z) should even be defined when z is a scalar, since it is defined as X * Y + z and we stopped doing implicit vector + scalar broadcasting a long time ago.
This was added in JuliaLang/julia#37065 by @mcabbot, but I don't see much discussion of the scalar case in that PR?
Activity