35 #ifndef TEMPLATE_BLAS_SYR2K_HEADER
36 #define TEMPLATE_BLAS_SYR2K_HEADER
41 const integer *k,
const Treal *alpha,
const Treal *a,
43 const Treal *beta, Treal *c__,
const integer *ldc)
46 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
54 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
55 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
56 #define c___ref(a_1,a_2) c__[(a_2)*c_dim1 + a_1]
156 a_offset = 1 + a_dim1 * 1;
159 b_offset = 1 + b_dim1 * 1;
162 c_offset = 1 + c_dim1 * 1;
181 }
else if (*lda <
maxMACRO(1,nrowa)) {
183 }
else if (*ldb <
maxMACRO(1,nrowa)) {
193 if (*n == 0 || ( (*alpha == 0. || *k == 0) && *beta == 1. ) ) {
201 for (j = 1; j <= i__1; ++j) {
203 for (i__ = 1; i__ <= i__2; ++i__) {
211 for (j = 1; j <= i__1; ++j) {
213 for (i__ = 1; i__ <= i__2; ++i__) {
223 for (j = 1; j <= i__1; ++j) {
225 for (i__ = j; i__ <= i__2; ++i__) {
233 for (j = 1; j <= i__1; ++j) {
235 for (i__ = j; i__ <= i__2; ++i__) {
250 for (j = 1; j <= i__1; ++j) {
253 for (i__ = 1; i__ <= i__2; ++i__) {
257 }
else if (*beta != 1.) {
259 for (i__ = 1; i__ <= i__2; ++i__) {
265 for (l = 1; l <= i__2; ++l) {
267 temp1 = *alpha *
b_ref(j, l);
268 temp2 = *alpha *
a_ref(j, l);
270 for (i__ = 1; i__ <= i__3; ++i__) {
272 * temp1 + b_ref(i__, l) * temp2;
282 for (j = 1; j <= i__1; ++j) {
285 for (i__ = j; i__ <= i__2; ++i__) {
289 }
else if (*beta != 1.) {
291 for (i__ = j; i__ <= i__2; ++i__) {
297 for (l = 1; l <= i__2; ++l) {
299 temp1 = *alpha *
b_ref(j, l);
300 temp2 = *alpha *
a_ref(j, l);
302 for (i__ = j; i__ <= i__3; ++i__) {
304 * temp1 + b_ref(i__, l) * temp2;
317 for (j = 1; j <= i__1; ++j) {
319 for (i__ = 1; i__ <= i__2; ++i__) {
323 for (l = 1; l <= i__3; ++l) {
329 c___ref(i__, j) = *alpha * temp1 + *alpha * temp2;
332 temp1 + *alpha * temp2;
340 for (j = 1; j <= i__1; ++j) {
342 for (i__ = j; i__ <= i__2; ++i__) {
346 for (l = 1; l <= i__3; ++l) {
352 c___ref(i__, j) = *alpha * temp1 + *alpha * temp2;
355 temp1 + *alpha * temp2;