Fix performance issue when RHS is a Matrix (#137)

Author: rasmushenningsson

Project.toml

@@ -1,6 +1,6 @@
 name = "LinearMaps"
 uuid = "7a12625a-238d-50fd-b39a-03d52299707e"
-version = "3.2.1"
+version = "3.2.2"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/LinearMaps.jl

@@ -213,11 +213,11 @@ end
 # the following is of interest in, e.g., subspace-iteration methods
 function mul!(Y::AbstractMatrix, A::LinearMap, X::AbstractMatrix)
     check_dim_mul(Y, A, X)
-    return _generic_mapmat_mul!(Y, A, X)
+    return _unsafe_mul!(Y, A, X)
 end
 function mul!(Y::AbstractMatrix, A::LinearMap, X::AbstractMatrix, α::Number, β::Number)
     check_dim_mul(Y, A, X)
-    return _generic_mapmat_mul!(Y, A, X, α, β)
+    return _unsafe_mul!(Y, A, X, α, β)
 end
 
 function _generic_mapmat_mul!(Y, A, X, α=true, β=false)
@@ -236,6 +236,9 @@ _unsafe_mul!(y, A::AbstractMatrix, x, α, β) = mul!(y, A, x, α, β)
 function _unsafe_mul!(y::AbstractVecOrMat, A::LinearMap, x::AbstractVector, α, β)
     return _generic_mapvec_mul!(y, A, x, α, β)
 end
+function _unsafe_mul!(y::AbstractMatrix, A::LinearMap, x::AbstractMatrix)
+    return _generic_mapmat_mul!(y, A, x)
+end
 function _unsafe_mul!(y::AbstractMatrix, A::LinearMap, x::AbstractMatrix, α, β)
     return _generic_mapmat_mul!(y, A, x, α, β)
 end