num::kernel::spmv_axpy

Defined in header <kernel/sparse.hpp>

template <std::floating_point T, std::integral Index>
void spmv_axpy(T *y, T alpha, const T *val, const Index *row_ptr, const Index *col_idx,
               const T *x, T beta, std::type_identity_t<Index> m) noexcept;

Fused CSR SpMV and vector accumulation \(\mathbf{y} \leftarrow \alpha A \mathbf{x} + \beta \mathbf{y}\).