37 #ifndef TEMPLATE_BLAS_SYR2_HEADER 38 #define TEMPLATE_BLAS_SYR2_HEADER 43 const Treal *x,
const integer *incx,
const Treal *y,
const integer *incy,
47 integer a_dim1, a_offset, i__1, i__2;
53 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 125 a_offset = 1 + a_dim1 * 1;
128 jx = jy = kx = ky = 0;
135 }
else if (*incx == 0) {
137 }
else if (*incy == 0) {
147 if (*n == 0 || *alpha == 0.) {
152 if (*incx != 1 || *incy != 1) {
156 kx = 1 - (*n - 1) * *incx;
161 ky = 1 - (*n - 1) * *incy;
171 if (*incx == 1 && *incy == 1) {
173 for (j = 1; j <= i__1; ++j) {
174 if (x[j] != 0. || y[j] != 0.) {
175 temp1 = *alpha * y[j];
176 temp2 = *alpha * x[j];
178 for (i__ = 1; i__ <= i__2; ++i__) {
179 a_ref(i__, j) =
a_ref(i__, j) + x[i__] * temp1 + y[
188 for (j = 1; j <= i__1; ++j) {
189 if (x[jx] != 0. || y[jy] != 0.) {
190 temp1 = *alpha * y[jy];
191 temp2 = *alpha * x[jx];
195 for (i__ = 1; i__ <= i__2; ++i__) {
196 a_ref(i__, j) =
a_ref(i__, j) + x[ix] * temp1 + y[iy]
210 if (*incx == 1 && *incy == 1) {
212 for (j = 1; j <= i__1; ++j) {
213 if (x[j] != 0. || y[j] != 0.) {
214 temp1 = *alpha * y[j];
215 temp2 = *alpha * x[j];
217 for (i__ = j; i__ <= i__2; ++i__) {
218 a_ref(i__, j) =
a_ref(i__, j) + x[i__] * temp1 + y[
227 for (j = 1; j <= i__1; ++j) {
228 if (x[jx] != 0. || y[jy] != 0.) {
229 temp1 = *alpha * y[jy];
230 temp2 = *alpha * x[jx];
234 for (i__ = j; i__ <= i__2; ++i__) {
235 a_ref(i__, j) =
a_ref(i__, j) + x[ix] * temp1 + y[iy]
int integer
Definition: template_blas_common.h:40
int template_blas_syr2(const char *uplo, 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_syr2.h:42
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46