19 #define A(m,n) BLKADDR(A, double, (m), (n))
20 #define T(m,n) BLKADDR(T, double, (m), (n))
21 #define T2(m,n) BLKADDR(T, double, (m), (n)+A.nt)
33 int ldaM, ldam, ldaMRD;
34 int tempkmin, tempkn, tempMm, tempnn, tempmm, tempMRDm;
44 for (k = 0; k <
min(A.
mt, A.
nt); k++) {
45 tempkn = k == A.
nt-1 ? A.
n-k*A.
nb : A.
nb;
46 for (M = k; M < A.
mt; M += BS) {
47 tempMm = M == A.
mt-1 ? A.
m-M*A.
mb : A.
mb;
48 tempkmin =
min(tempMm, tempkn);
51 plasma->
quark, &task_flags,
52 tempMm, tempkn, ib, T.
nb,
55 for (n = k+1; n < A.
nt; n++) {
56 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
58 plasma->
quark, &task_flags,
60 tempMm, tempnn, tempkmin, ib, T.
nb,
65 for (m = M+1; m <
min(M+BS, A.
mt); m++) {
66 tempmm = m == A.
mt-1 ? A.
m-m*A.
mb : A.
mb;
69 plasma->
quark, &task_flags,
70 tempmm, tempkn, ib, T.
nb,
75 for (n = k+1; n < A.
nt; n++) {
76 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
78 plasma->
quark, &task_flags,
80 A.
nb, tempnn, tempmm, tempnn, A.
nb, ib, T.
nb,
88 for (RD = BS; RD < A.
mt-k; RD *= 2) {
89 for (M = k; M+RD < A.
mt; M += 2*RD) {
90 tempMRDm = M+RD == A.
mt-1 ? A.
m-(M+RD)*A.
mb : A.
mb;
94 plasma->
quark, &task_flags,
95 tempMRDm, tempkn, ib, T.
nb,
100 for (n = k+1; n < A.
nt; n++) {
101 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
103 plasma->
quark, &task_flags,
105 A.
nb, tempnn, tempMRDm, tempnn, A.
nb, ib, T.
nb,