35 #ifndef TEMPLATE_BLAS_SYMM_HEADER
36 #define TEMPLATE_BLAS_SYMM_HEADER
41 const Treal *alpha,
const Treal *a,
const integer *lda,
const Treal *b,
42 const integer *ldb,
const Treal *beta, Treal *c__,
const integer *ldc)
45 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
53 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
54 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
55 #define c___ref(a_1,a_2) c__[(a_2)*c_dim1 + a_1]
154 a_offset = 1 + a_dim1 * 1;
157 b_offset = 1 + b_dim1 * 1;
160 c_offset = 1 + c_dim1 * 1;
179 }
else if (*lda <
maxMACRO(1,nrowa)) {
191 if (*m == 0 || *n == 0 || ( *alpha == 0. && *beta == 1. ) ) {
198 for (j = 1; j <= i__1; ++j) {
200 for (i__ = 1; i__ <= i__2; ++i__) {
208 for (j = 1; j <= i__1; ++j) {
210 for (i__ = 1; i__ <= i__2; ++i__) {
224 for (j = 1; j <= i__1; ++j) {
226 for (i__ = 1; i__ <= i__2; ++i__) {
227 temp1 = *alpha *
b_ref(i__, j);
230 for (k = 1; k <= i__3; ++k) {
232 temp2 += b_ref(k, j) * a_ref(k, i__);
240 a_ref(i__, i__) + *alpha * temp2;
248 for (j = 1; j <= i__1; ++j) {
249 for (i__ = *m; i__ >= 1; --i__) {
250 temp1 = *alpha *
b_ref(i__, j);
253 for (k = i__ + 1; k <= i__2; ++k) {
255 temp2 += b_ref(k, j) * a_ref(k, i__);
263 a_ref(i__, i__) + *alpha * temp2;
273 for (j = 1; j <= i__1; ++j) {
274 temp1 = *alpha *
a_ref(j, j);
277 for (i__ = 1; i__ <= i__2; ++i__) {
283 for (i__ = 1; i__ <= i__2; ++i__) {
290 for (k = 1; k <= i__2; ++k) {
292 temp1 = *alpha * a_ref(k, j);
294 temp1 = *alpha * a_ref(j, k);
297 for (i__ = 1; i__ <= i__3; ++i__) {
304 for (k = j + 1; k <= i__2; ++k) {
306 temp1 = *alpha * a_ref(j, k);
308 temp1 = *alpha * a_ref(k, j);
311 for (i__ = 1; i__ <= i__3; ++i__) {