37 #ifndef TEMPLATE_BLAS_SYRK_HEADER 38 #define TEMPLATE_BLAS_SYRK_HEADER 43 const Treal *alpha,
const Treal *a,
const integer *lda,
const Treal *beta,
47 integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3;
54 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 55 #define c___ref(a_1,a_2) c__[(a_2)*c_dim1 + a_1] 139 a_offset = 1 + a_dim1 * 1;
142 c_offset = 1 + c_dim1 * 1;
161 }
else if (*lda <
maxMACRO(1,nrowa)) {
171 if (*n == 0 || ( (*alpha == 0. || *k == 0) && *beta == 1. ) ) {
179 for (j = 1; j <= i__1; ++j) {
181 for (i__ = 1; i__ <= i__2; ++i__) {
189 for (j = 1; j <= i__1; ++j) {
191 for (i__ = 1; i__ <= i__2; ++i__) {
201 for (j = 1; j <= i__1; ++j) {
203 for (i__ = j; i__ <= i__2; ++i__) {
211 for (j = 1; j <= i__1; ++j) {
213 for (i__ = j; i__ <= i__2; ++i__) {
228 for (j = 1; j <= i__1; ++j) {
231 for (i__ = 1; i__ <= i__2; ++i__) {
235 }
else if (*beta != 1.) {
237 for (i__ = 1; i__ <= i__2; ++i__) {
243 for (l = 1; l <= i__2; ++l) {
244 if (
a_ref(j, l) != 0.) {
245 temp = *alpha *
a_ref(j, l);
247 for (i__ = 1; i__ <= i__3; ++i__) {
259 for (j = 1; j <= i__1; ++j) {
262 for (i__ = j; i__ <= i__2; ++i__) {
266 }
else if (*beta != 1.) {
268 for (i__ = j; i__ <= i__2; ++i__) {
274 for (l = 1; l <= i__2; ++l) {
275 if (
a_ref(j, l) != 0.) {
276 temp = *alpha *
a_ref(j, l);
278 for (i__ = j; i__ <= i__3; ++i__) {
293 for (j = 1; j <= i__1; ++j) {
295 for (i__ = 1; i__ <= i__2; ++i__) {
298 for (l = 1; l <= i__3; ++l) {
303 c___ref(i__, j) = *alpha * temp;
314 for (j = 1; j <= i__1; ++j) {
316 for (i__ = j; i__ <= i__2; ++i__) {
319 for (l = 1; l <= i__3; ++l) {
324 c___ref(i__, j) = *alpha * temp;
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_syrk(const char *uplo, const char *trans, const integer *n, const integer *k, const Treal *alpha, const Treal *a, const integer *lda, const Treal *beta, Treal *c__, const integer *ldc)
Definition: template_blas_syrk.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
bool logical
Definition: template_blas_common.h:41
#define c___ref(a_1, a_2)
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46