BLAS Level 1 Extension: raxpy

raxpy (defined in namespace flens::blas) computes \(y \leftarrow y +\frac{1}{\alpha} x\) where \(x, y\) are vectors and \(\alpha\) is a constant scalar. raxpy is short for reciprocal of alpha times x plus y.

We also provide a variant that does the computation for matrices instead of vectors.

Vector Variant

template <typename ALPHA, typename VX, typename VY>
    typename RestrictTo<IsDenseVector<VX>::value
                     && IsDenseVector<VY>::value,
             void>::Type
    raxpy(const ALPHA &alpha, const VX &x, VY &&y);

alpha	`(input) real scalar` Reciprocal scaling factor for the values in \(x\).
x	`(input) real or complex valued DenseVector` Vector \(x\).
y	`(input/output) real or complex valued DenseVector` Vector \(y\).

Matrix Variant (Level 1 Extension)

template <typename ALPHA, typename MA, typename MB>
    typename RestrictTo<IsGeMatrix<MA>::value
                     && IsGeMatrix<MB>::value,
             void>::Type
    raxpy(Transpose trans, const ALPHA &alpha, const MA &A, MB &&B);

trans

(input)
Specifiy the operation, i.e.

NoTrans	\(B \leftarrow B + \frac{1}{\alpha}\,A\)
Trans	\(B \leftarrow B + \frac{1}{\alpha}\,A^T\)
ConjTrans	\(B \leftarrow B + \frac{1}{\alpha}\,A^H\)

alpha

(input) real scalar
Reciprocal scaling factor for the values in \(x\).

(input) real or complex valued GeMatrix
Matrix \(A\).

(input/output) real or complex valued GeMatrix
Matrix \(B\).