num::kernel::axpby

Defined in header <kernel/vector.hpp>

template <std::floating_point T>
void axpby(T *y, const T *x, T a, T b, idx n) noexcept;

Computes fused scaled vector update \(y_i \leftarrow a x_i + b y_i\) in a single memory pass.