3 #include "../common/monolish_common.hpp" 
 1262 [[nodiscard]] 
double 
 1986 [[nodiscard]] 
double 
 2031 [[nodiscard]] 
double 
void asum(const vector< double > &x, double &ans)
vector absolute sum
 
void axpy(const double alpha, const vector< double > &x, vector< double > &y)
axpy: y = ax + y
 
void axpyz(const double alpha, const vector< double > &x, const vector< double > &y, vector< double > &z)
axpyz: z = ax + y
 
void dot(const vector< double > &x, const vector< double > &y, double &ans)
inner product (dot)
 
void copy(const matrix::Dense< double > &A, matrix::Dense< double > &C)
Dense matrix copy (C=A)
 
void nrm1(const vector< double > &x, double &ans)
nrm1: sum(abs(x[0:N]))
 
void nrm2(const vector< double > &x, double &ans)
nrm2: ||x||_2
 
void scal(const double alpha, vector< double > &x)
scal: x = alpha * x
 
void sum(const vector< double > &x, double &ans)
vector sum
 
void times(const double alpha, const matrix::Dense< double > &A, matrix::Dense< double > &C)
Dense matrix times: C = alpha * A.
 
void vecadd(const vector< double > &a, const vector< double > &b, vector< double > &y)
element by element addition of vector a and vector b.
 
void vecsub(const vector< double > &a, const vector< double > &b, vector< double > &y)
element by element subtract of vector a and vector b.
 
void xpay(const double alpha, const vector< double > &x, vector< double > &y)
xpay: y = x + ay