35 #ifndef TEMPLATE_BLAS_SYRK_HEADER
36 #define TEMPLATE_BLAS_SYRK_HEADER
41 const Treal *alpha,
const Treal *a,
const integer *lda,
const Treal *beta,
45 integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3;
52 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
53 #define c___ref(a_1,a_2) c__[(a_2)*c_dim1 + a_1]
137 a_offset = 1 + a_dim1 * 1;
140 c_offset = 1 + c_dim1 * 1;
159 }
else if (*lda <
maxMACRO(1,nrowa)) {
169 if (*n == 0 || ( (*alpha == 0. || *k == 0) && *beta == 1. ) ) {
177 for (j = 1; j <= i__1; ++j) {
179 for (i__ = 1; i__ <= i__2; ++i__) {
187 for (j = 1; j <= i__1; ++j) {
189 for (i__ = 1; i__ <= i__2; ++i__) {
199 for (j = 1; j <= i__1; ++j) {
201 for (i__ = j; i__ <= i__2; ++i__) {
209 for (j = 1; j <= i__1; ++j) {
211 for (i__ = j; i__ <= i__2; ++i__) {
226 for (j = 1; j <= i__1; ++j) {
229 for (i__ = 1; i__ <= i__2; ++i__) {
233 }
else if (*beta != 1.) {
235 for (i__ = 1; i__ <= i__2; ++i__) {
241 for (l = 1; l <= i__2; ++l) {
242 if (
a_ref(j, l) != 0.) {
243 temp = *alpha *
a_ref(j, l);
245 for (i__ = 1; i__ <= i__3; ++i__) {
257 for (j = 1; j <= i__1; ++j) {
260 for (i__ = j; i__ <= i__2; ++i__) {
264 }
else if (*beta != 1.) {
266 for (i__ = j; i__ <= i__2; ++i__) {
272 for (l = 1; l <= i__2; ++l) {
273 if (
a_ref(j, l) != 0.) {
274 temp = *alpha *
a_ref(j, l);
276 for (i__ = j; i__ <= i__3; ++i__) {
291 for (j = 1; j <= i__1; ++j) {
293 for (i__ = 1; i__ <= i__2; ++i__) {
296 for (l = 1; l <= i__3; ++l) {
301 c___ref(i__, j) = *alpha * temp;
312 for (j = 1; j <= i__1; ++j) {
314 for (i__ = j; i__ <= i__2; ++i__) {
317 for (l = 1; l <= i__3; ++l) {
322 c___ref(i__, j) = *alpha * temp;