37 #ifndef TEMPLATE_BLAS_SYR2K_HEADER 38 #define TEMPLATE_BLAS_SYR2K_HEADER 43 const integer *k,
const Treal *alpha,
const Treal *a,
45 const Treal *beta, Treal *c__,
const integer *ldc)
48 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
56 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 57 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1] 58 #define c___ref(a_1,a_2) c__[(a_2)*c_dim1 + a_1] 158 a_offset = 1 + a_dim1 * 1;
161 b_offset = 1 + b_dim1 * 1;
164 c_offset = 1 + c_dim1 * 1;
183 }
else if (*lda <
maxMACRO(1,nrowa)) {
185 }
else if (*ldb <
maxMACRO(1,nrowa)) {
195 if (*n == 0 || ( (*alpha == 0. || *k == 0) && *beta == 1. ) ) {
203 for (j = 1; j <= i__1; ++j) {
205 for (i__ = 1; i__ <= i__2; ++i__) {
213 for (j = 1; j <= i__1; ++j) {
215 for (i__ = 1; i__ <= i__2; ++i__) {
225 for (j = 1; j <= i__1; ++j) {
227 for (i__ = j; i__ <= i__2; ++i__) {
235 for (j = 1; j <= i__1; ++j) {
237 for (i__ = j; i__ <= i__2; ++i__) {
252 for (j = 1; j <= i__1; ++j) {
255 for (i__ = 1; i__ <= i__2; ++i__) {
259 }
else if (*beta != 1.) {
261 for (i__ = 1; i__ <= i__2; ++i__) {
267 for (l = 1; l <= i__2; ++l) {
269 temp1 = *alpha *
b_ref(j, l);
270 temp2 = *alpha *
a_ref(j, l);
272 for (i__ = 1; i__ <= i__3; ++i__) {
274 * temp1 +
b_ref(i__, l) * temp2;
284 for (j = 1; j <= i__1; ++j) {
287 for (i__ = j; i__ <= i__2; ++i__) {
291 }
else if (*beta != 1.) {
293 for (i__ = j; i__ <= i__2; ++i__) {
299 for (l = 1; l <= i__2; ++l) {
301 temp1 = *alpha *
b_ref(j, l);
302 temp2 = *alpha *
a_ref(j, l);
304 for (i__ = j; i__ <= i__3; ++i__) {
306 * temp1 +
b_ref(i__, l) * temp2;
319 for (j = 1; j <= i__1; ++j) {
321 for (i__ = 1; i__ <= i__2; ++i__) {
325 for (l = 1; l <= i__3; ++l) {
331 c___ref(i__, j) = *alpha * temp1 + *alpha * temp2;
334 temp1 + *alpha * temp2;
342 for (j = 1; j <= i__1; ++j) {
344 for (i__ = j; i__ <= i__2; ++i__) {
348 for (l = 1; l <= i__3; ++l) {
354 c___ref(i__, j) = *alpha * temp1 + *alpha * temp2;
357 temp1 + *alpha * temp2;
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
#define c___ref(a_1, a_2)
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
bool logical
Definition: template_blas_common.h:41
int template_blas_syr2k(const char *uplo, const char *trans, const integer *n, const integer *k, const Treal *alpha, const Treal *a, const integer *lda, const Treal *b, const integer *ldb, const Treal *beta, Treal *c__, const integer *ldc)
Definition: template_blas_syr2k.h:42
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46