In one of the PEPSKit testruns, we are running into an issue where the within the rrule for svdsolve
, the eigensolver is failing to converge all requested values. In this case, the output vals
and vecs
are not of the expected length, leading to an out of bounds error.
The stacktrace is found here
The offending line is this line .
There is already a warning being thrown, but I guess we have to limit n
to the size of the computed vals and vecs.