diff --git a/src/operators/common_defaults.jl b/src/operators/common_defaults.jl index 260cee71c..a5f4c9130 100644 --- a/src/operators/common_defaults.jl +++ b/src/operators/common_defaults.jl @@ -16,6 +16,8 @@ for op in (:*, :/, :\), T in (:AbstractArray, :Number) end LinearAlgebra.mul!(Y::AbstractArray, L::AbstractDiffEqLinearOperator, B::AbstractArray) = mul!(Y, convert(AbstractMatrix,L), B) +LinearAlgebra.mul!(Y::AbstractArray, L::AbstractDiffEqLinearOperator, B::AbstractArray, α::Number, β::Number) = + mul!(Y, convert(AbstractMatrix,L), B, α, β) LinearAlgebra.ldiv!(Y::AbstractVecOrMat, L::AbstractDiffEqLinearOperator, B::AbstractVecOrMat) = ldiv!(Y, convert(AbstractMatrix,L), B) for pred in (:isreal, :issymmetric, :ishermitian, :isposdef)