num::kernel::axpy_norm_sq

Defined in header <kernel/vector.hpp>

template <std::floating_point T>
T axpy_norm_sq(T *y, const T *x, T alpha, idx n) noexcept;

y <- y + alpha*x, returning the new \(\|y\|^2\) from the same traversal.