Xingjian Guo

Full list of exponential integrators to implement

8

Note: many integrators for linear systems also use exponential-like matrix functions and Krylov approximations (Magnus expansion, adiabatic integrators, symplectic exponential methods, etc). Since they require special structure for the intended...

Drop `ishermitian` and `opnorm` requirement for custom operator types.

1

Currently to use a matrix-free operator we need to implement the `LinearAlgebra.ishermitian` and `LinearAlgebra.opnorm` (or provide a custom norm function). These slows down rapid prototyping and for some sparse types...

Add benchmarks

2

From slack: > we should post the benchmarks, but essentially if you want to compute the matrix exponential `exp(t*A)` for a large `t`, timestepping it can help improve the accuracy...

Extending expv using Saad to general phiv?

2

@jagot I have finally went through the paper by Saad. While in theory the algorithm can be extended to arbitrary phi orders and non-Hermitian matrices, I still have some doubts...

Prototype for the first target (#33)

17

Prototype code for #33. The focus is on setting up the matrix-free operators for the generic `L*QB` route and come up with a mock user interface to create the operators,...