num::blas::add

Defined in header <blas/vector_ops.hpp>

void add(const vec<real> &x, const vec<real> &y, vec<real> &z) noexcept;

z <- x + y. No single BLAS call does an out-of-place add; a plain loop is exactly as fast as dcopy+daxpy and skips the extra pass.