What and Why
For open quantum systems, a solver that is more efficient than the standard ODE solvers would be nice. One such example is Peter's solver based on Newton polynomials.
TODO
- Understand how the time evolution with Newton polynomials works
- Input parameters
- error bounds and scaling
- Implement a solver that uses Newton polynomials
- Add a demo or convert an existing demo that uses Newton polyomials
- Follwup ticket: Use a SIL solver?
- Used a solver based on Faber polynomials; these looked easier to understand and handle.
- also added some helper functions related to guessing the spectrum and truncating operators
- new demo LinearRamp as an example
Diff:
Diff: