ExpmV.jl
ExpmV.jl copied to clipboard
marcusps
Julia package to compute the result of expm(t*A)*v when A is a sparse matrix, without computing expm(t*A).
Any chance of an update to be compatible with Julia 1.x? I'd be willing to help with this; #10 already updates to Julia 0.6 if that can be merged, and...
Julia 0.6 warns that `ceil()` and `abs()` without the broadcasting dot are deprecated.
[LinearMaps.jl](https://github.com/Jutho/LinearMaps.jl) can map a matrix vector multiplication to a linear operator. Is it possible to support `expmv` fallback for a general linear operator? I need it very much in my...
It looks like may be using `norm(matrix)`. In Julia 0.7, this will compute the Frobenius norm (`vecnorm` in Julia 0.6), due to JuliaLang/julia#27401. If you want the induced/operator norm as...
The code currently computes the 1 norm using `norm(A^m,1)` instead of estimating more efficiently. This ends up being the bottleneck in the calculation of not so sparse matrices. The MATLAB...
I was recently using ExmpV.jl for solving a Lindblad master equation for a time-independent system of coupled harmonic oscillators, i.e. a very sparse `A` that is constant over many `exmpv`...
Register
Hey, would you consider finishing up the last pieces and registering? What exactly needs to be done for that? Where would you need help? I am considering using this as...
My problem of interest is imaginary time evolution of a Schrodinger-like equation, so I need to discretize time and re-normalize the wavefunction after each time step in order to avoid...
