Skip to content

num::hessenberg_project

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

template <class Rhs>
vec<cplx> hessenberg_project(const mat<real> &Q, const Rhs &b);

Project a right-hand side onto the Hessenberg basis: \(\tilde b = Q^T b\).

Accepts a real or complex right-hand side. The result is complex either way, since the shift generally is.