37 #ifndef TEMPLATE_BLAS_GER_HEADER 38 #define TEMPLATE_BLAS_GER_HEADER 43 const Treal *x,
const integer *incx,
const Treal *y,
const integer *incy,
47 integer a_dim1, a_offset, i__1, i__2;
52 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 110 a_offset = 1 + a_dim1 * 1;
118 }
else if (*incx == 0) {
120 }
else if (*incy == 0) {
130 if (*m == 0 || *n == 0 || *alpha == 0.) {
138 jy = 1 - (*n - 1) * *incy;
142 for (j = 1; j <= i__1; ++j) {
144 temp = *alpha * y[jy];
146 for (i__ = 1; i__ <= i__2; ++i__) {
147 a_ref(i__, j) =
a_ref(i__, j) + x[i__] * temp;
158 kx = 1 - (*m - 1) * *incx;
161 for (j = 1; j <= i__1; ++j) {
163 temp = *alpha * y[jy];
166 for (i__ = 1; i__ <= i__2; ++i__) {
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_ger(const integer *m, const integer *n, const Treal *alpha, const Treal *x, const integer *incx, const Treal *y, const integer *incy, Treal *a, const integer *lda)
Definition: template_blas_ger.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146