Jeff Fessler
Jeff Fessler
I've got it 90% done. Famous saying about the last 10%... Will send something soon.
Related to the original question, v0.6 of https://github.com/JeffFessler/LinearMapsAA.jl now provides a type (`LinearMapAO`, where O is for Operator) that supports linear mappings between `AbstractArrays` of arbitrary dimensions, quite similar to...
A complication is that `false*I` is also of the same type but is mathematically equivalent to `0I`. So a slightly more complicated function would be needed to check the value...
I was curious what other objects do with `* I` so I ran a simple test: `Diagonal(1:5) * I` returns type `Diagonal{Int64,StepRange{Int64,Int64}}` whereas `Diagonal(1:5)` has type `Diagonal{Int64,UnitRange{Int64}}` which is almost...
what about the more general case of, say, `1.0I` ``` function Base.:(*)(A1::LinearMap, A2::UniformScaling) # needs refinement to exclude Bool? if A2.λ == 1 # or true ? return A1 else...
@platawiec fyi I have added the "simply return A" feature to https://github.com/JeffFessler/LinearMapsAA.jl (master branch only, not yet tagged)
This was discussed in #38. Nice to see another vote in favor of `AbstractMatrix` though. See https://github.com/JeffFessler/LinearMapsAA.jl
When I first brought up #38, I also felt pretty strongly that a `LinearMap` should be a subtype of `AbstractMatrix`, which led to `LinearMapsAA`. Over time, however, I found that...
Ok, I did a test of `LinearMapsAA.jl` and it passed, but with only one deprecation warning so maybe that package does not actually use so many `FunctionMap`s. So then I...