num::hessenberg_shifted_substitute

Defined in header <linear/factorization/hessenberg.hpp>

template <std::floating_point T, class Index>
void hessenberg_shifted_substitute(std::complex<T> *y, const std::complex<T> *work,
                                   const Index *piv, const std::complex<T> *b,
                                   idx n) noexcept;

Substitute a right-hand side through a factored shifted Hessenberg system.

Separate from the factorization so that many right-hand sides share one factorization at the same shift. y and b may alias.