Muladd! (useful?)

jling · October 30, 2022, 7:57am

JuliaLang/julia/blob/v1.8.2/stdlib/LinearAlgebra/src/matmul.jl#L246


      
          Base.muladd(x::TransposeAbsVec, A::AbstractMatrix, z::Union{Number, AbstractVecOrMat}) =
              transpose(muladd(transpose(A), transpose(x), transpose(z)))
          
          
function Base.muladd(A::StridedMaybeAdjOrTransMat{<:Number}, y::AbstractVector{<:Number}, z::Union{Number, AbstractVector})
              T = promote_type(eltype(A), eltype(y), eltype(z))
              C = similar(A, T, axes(A,1))
              C .= z
              mul!(C, A, y, true, true)
          end
          
          
function Base.muladd(A::StridedMaybeAdjOrTransMat{<:Number}, B::StridedMaybeAdjOrTransMat{<:Number}, z::Union{Number, AbstractVecOrMat})
              T = promote_type(eltype(A), eltype(B), eltype(z))
              C = similar(A, T, axes(A,1), axes(B,2))
              C .= z
              mul!(C, A, B, true, true)
          end
          
          
"""
              mul!(Y, A, B) -> Y
          
          
Calculates the matrix-matrix or matrix-vector product ``AB`` and stores the result in `Y`,

Topic		Replies	Views
In-place matrix multiplication compatible with 0.5 and 0.6 General Usage	3	746	January 10, 2017
Inplace multiplication by a square matrix General Usage	11	7552	February 17, 2017
Are there any 5-argument versions of `mul!` or similar in LinearAlgebra now? Performance	1	944	September 5, 2018
Working with MulAddMul Performance	0	394	March 6, 2020
In-place multiplication methods for UniformScaling type Internals & Design linearalgebra	4	561	October 3, 2018

Muladd! (useful?)

Related topics