37 #ifndef TEMPLATE_BLAS_SYMM_HEADER 38 #define TEMPLATE_BLAS_SYMM_HEADER 43 const Treal *alpha,
const Treal *a,
const integer *lda,
const Treal *b,
44 const integer *ldb,
const Treal *beta, Treal *c__,
const integer *ldc)
47 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
55 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 56 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1] 57 #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)) {
193 if (*m == 0 || *n == 0 || ( *alpha == 0. && *beta == 1. ) ) {
200 for (j = 1; j <= i__1; ++j) {
202 for (i__ = 1; i__ <= i__2; ++i__) {
210 for (j = 1; j <= i__1; ++j) {
212 for (i__ = 1; i__ <= i__2; ++i__) {
226 for (j = 1; j <= i__1; ++j) {
228 for (i__ = 1; i__ <= i__2; ++i__) {
229 temp1 = *alpha *
b_ref(i__, j);
232 for (k = 1; k <= i__3; ++k) {
242 a_ref(i__, i__) + *alpha * temp2;
250 for (j = 1; j <= i__1; ++j) {
251 for (i__ = *m; i__ >= 1; --i__) {
252 temp1 = *alpha *
b_ref(i__, j);
255 for (k = i__ + 1; k <= i__2; ++k) {
265 a_ref(i__, i__) + *alpha * temp2;
275 for (j = 1; j <= i__1; ++j) {
276 temp1 = *alpha *
a_ref(j, j);
279 for (i__ = 1; i__ <= i__2; ++i__) {
285 for (i__ = 1; i__ <= i__2; ++i__) {
292 for (k = 1; k <= i__2; ++k) {
294 temp1 = *alpha *
a_ref(k, j);
296 temp1 = *alpha *
a_ref(j, k);
299 for (i__ = 1; i__ <= i__3; ++i__) {
306 for (k = j + 1; k <= i__2; ++k) {
308 temp1 = *alpha *
a_ref(j, k);
310 temp1 = *alpha *
a_ref(k, j);
313 for (i__ = 1; i__ <= i__3; ++i__) {
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_symm(const char *side, const char *uplo, const integer *m, const integer *n, 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_symm.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
#define c___ref(a_1, a_2)
bool logical
Definition: template_blas_common.h:41
side
Definition: Matrix.h:75
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46