Skip to content

num::kernel::matvec_real_complex

Defined in header <kernel/complex.hpp>

template <std::floating_point T>
void matvec_real_complex(std::complex<T> *x, const T *Q, const std::complex<T> *y,
                         idx m, idx n) noexcept;

Mixed real-matrix, complex-vector product \(x = Q y\).

Arises when a real orthogonal basis is applied to a complex Krylov coordinate vector, as in projecting a resolvent solution back from the Krylov subspace. Kept separate from matvec because the scalar types differ on the two sides.