In contrast to what is written in the documentation, propagating backwards in imaginary time is not unstable; it will generally converge to the largest eigenstate with a nrom growing to infinity, but that is a feature, not a bug.
I am not quite sure what propagation in negative imaginary time is good for (increasing temperature? Spin states?), but it seems wrong to restrict the propagators that way for no good reason. While not really urgent, the cost is rather low.
As an additional piece of work, removed OdeRelax. There is little point in this class, because relaxation with Chebychev polynomials should always be faster.
Diff:
Diff: