19 #define A(m,n) BLKADDR(A, double, m, n)
20 #define L(m,n) BLKADDR(L, double, m, n)
21 #define IPIV(m,n) &(IPIV[(int64_t)A.mb*((int64_t)(m)+(int64_t)A.mt*(int64_t)(n))])
39 int tempkn, tempkm, tempmm, tempnn;
65 while (next_n >= A.
nt && next_k <
min(A.
mt, A.
nt)) {
67 next_n = next_n-A.
nt+next_k;
72 tempmm = m == A.
mt-1 ? A.
m-m*A.
mb : A.
mb;
73 tempkm = k == A.
mt-1 ? A.
m-k*A.
mb : A.
mb;
74 tempkn = k == A.
nt-1 ? A.
n-k*A.
nb : A.
nb;
75 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
87 if (info != 0 && m == A.
mt-1) {
96 tempmm, tempkn, ib, A.
nb,
102 if (info != 0 && m == A.
mt-1) {
114 tempkm, tempnn, tempkm, ib,
123 A.
nb, tempnn, tempmm, tempnn, A.
nb, ib,
151 int tempkm, tempkn, tempmm, tempnn;
160 for (k = 0; k <
min(A.
mt, A.
nt); k++) {
161 tempkm = k == A.
mt-1 ? A.
m-k*A.
mb : A.
mb;
162 tempkn = k == A.
nt-1 ? A.
n-k*A.
nb : A.
nb;
165 plasma->
quark, &task_flags,
166 tempkm, tempkn, ib, L.
nb,
167 A(k, k), ldak,
IPIV(k, k),
169 k == A.
mt-1, A.
nb*k);
171 for (n = k+1; n < A.
nt; n++) {
172 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
174 plasma->
quark, &task_flags,
175 tempkm, tempnn, tempkm, ib, L.
nb,
180 for (m = k+1; m < A.
mt; m++) {
181 tempmm = m == A.
mt-1 ? A.
m-m*A.
mb : A.
mb;
184 plasma->
quark, &task_flags,
185 tempmm, tempkn, ib, L.
nb,
191 m == A.
mt-1, A.
nb*k);
193 for (n = k+1; n < A.
nt; n++) {
194 tempnn = n == A.
nt-1 ? A.
n-n*A.
nb : A.
nb;
196 plasma->
quark, &task_flags,
197 A.
nb, tempnn, tempmm, tempnn, A.
nb, ib, L.
nb,