Fix Matrix(TransposeMap{CustomMap}) (#159)

Author: dkarrasch

Project.toml

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

src/conversion.jl

@@ -5,7 +5,8 @@ function Base.Matrix{T}(A::LinearMap) where {T}
     v = fill(zero(T), N)
     @inbounds for i in 1:N
         v[i] = one(T)
-        _unsafe_mul!(view(mat, :, i), A, v)
+        # need mul!, e.g., for TransposeMap{<:CustomMap}
+        mul!(view(mat, :, i), A, v)
         v[i] = zero(T)
     end
     return mat

test/linearmaps.jl

@@ -117,4 +117,5 @@ end
         @test mul!(ones(3), A', x, α, β) == fill(β, 3) + fill(15α, 3)
         @test mul!(ones(3, 2), A', [x x], α, β) == fill(β, 3, 2) + fill(15α, 3, 2)
     end
+    @test Matrix(A') == Matrix(A)'
 end