geoffroyleconte

In-place solve for inverse operators

4

I would be nice to add in-place solves to LinearOperators in the linalg.jl file.

`*` between operators and Numbers is assumed to be a commutative operation, but this could not be the expected behaviour for some specific Numbers (quaternions? ).

convert function for solvers

2

It could be practical to have something like: ```julia convert(::Type{KrylovSolver{T, FC}}, solver::MinresSolver) where {T, FC} = MinresSolver(...) ``` to convert a `solver` from `Float32` to `Float64` for example. However I'm...

It would be interesting to use vectorization macros such as `@simd` or `@avx` to increase the speed of this package and if possible to add GPU support. However this would...

possible method error upper triangle

1

``` n, m = size(A) A = sprand(Float64, 10, 5, 0.25) λ = 1.0e-3 K1 = triu([sparse(1.0I,n,n) A; A' -λ^2 * I]) S1 = ldl(Symmetric(K1, :U)) K2 = [sparse(1.0I,n,n) A;...

The only solver useful for most users is iterations/solvers/augmented/K2LDL.jl. Not precompiling all the other solvers files and preconditioner files might improve a lot the precompilation time of the package.