Closed DanielVandH closed 1 month ago
Lowering used to be implemented before we moved the package to decomposition based connection coefficients. You will see a TODO comment in semiclassical_jacobimatrix
about this as well. It's in principle not an issue to implement this using QL or reverse Cholesky, would be nice to have.
Missed that TODO comment. I think for now I'm going to just look at implementing the inv
approach into ldiv
rather than working with QL / reverse Cholesky itself, mimicing what's done in Jacobi
This works now on master
The issue is with
semijacobi_ldiv
:Since
Δb < 0
, none of the conditions in the first block ofif
statements are satisfied, but since it ends in anelseif
it just returnsNothing
. I'll have to try and make it handle negatives one at a time and then inverse at the end, just like is done with e.g.Jacobi
.