Commit 2685add3 authored by Mathieu Faverge's avatar Mathieu Faverge

Add comment to specify if it is a TS or TT kernel

parent bbde4462
......@@ -138,6 +138,7 @@ void morse_pzgelqf(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D,
RUNTIME_data_migrate( sequence, A(k, k),
A->get_rankof( A, k, n ) );
/* TS kernel */
MORSE_TASK_ztplqt(
&options,
tempkm, tempnn, 0, ib, T->nb,
......
......@@ -151,12 +151,13 @@ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, MORSE_desc_t *A,
tempnn = n == A->nt-1 ? A->n-n*A->nb : A->nb;
temppn = p == A->nt-1 ? A->n-p*A->nb : A->nb;
/* Tiles killed is a TS */
if (qrtree->gettype(qrtree, k, n) == 0) {
/* TS kernel */
T = TS;
L = 0;
}
else {
/* TT kernel */
T = TT;
L = tempnn;
}
......
......@@ -138,6 +138,7 @@ void morse_pzgelqfrh(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D, int BS,
RUNTIME_data_migrate( sequence, A(k, N),
A->get_rankof( A, k, n ) );
/* TS kernel */
MORSE_TASK_ztplqt(
&options,
tempkm, tempnn, 0, ib, T->nb,
......@@ -172,6 +173,7 @@ void morse_pzgelqfrh(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D, int BS,
RUNTIME_data_migrate( sequence, A(k, N+RD),
A->get_rankof( A, k, N+RD ) );
/* TT kernel */
MORSE_TASK_ztplqt(
&options,
tempkm, tempNRDn, chameleon_min(tempNRDn, tempkm), ib, T->nb,
......
......@@ -134,6 +134,7 @@ void morse_pzgeqrf(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D,
RUNTIME_data_migrate( sequence, A(k, k),
A->get_rankof( A, m, k ) );
/* TS kernel */
MORSE_TASK_ztpqrt(
&options,
tempmm, tempkn, 0, ib, T->nb,
......@@ -147,6 +148,7 @@ void morse_pzgeqrf(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D,
RUNTIME_data_migrate( sequence, A(k, n),
A->get_rankof( A, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseConjTrans,
......
......@@ -149,12 +149,13 @@ void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, MORSE_desc_t *A,
ldap = BLKLDD(A, p);
ldam = BLKLDD(A, m);
/* Tiles killed is a TS */
if (qrtree->gettype(qrtree, k, m) == 0) {
/* TS kernel */
T = TS;
L = 0;
}
else {
/* TT kernel */
T = TT;
L = tempmm;
}
......
......@@ -138,6 +138,7 @@ void morse_pzgeqrfrh(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D, int BS,
RUNTIME_data_migrate( sequence, A(M, k),
A->get_rankof( A, m, k ) );
/* TS kernel */
MORSE_TASK_ztpqrt(
&options,
tempmm, tempkn, 0, ib, T->nb,
......@@ -173,6 +174,7 @@ void morse_pzgeqrfrh(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D, int BS,
RUNTIME_data_migrate( sequence, A(M+RD, k),
A->get_rankof( A, M+RD, k ) );
/* TT kernel */
MORSE_TASK_ztpqrt(
&options,
tempMRDm, tempkn, chameleon_min( tempMRDm, tempkn ), ib, T->nb,
......
......@@ -106,6 +106,7 @@ void morse_pztpgqrt( int L,
for (n = k; n < Q2->nt; n++) {
tempnn = n == Q2->nt-1 ? Q2->n-n*Q2->nb : Q2->nb;
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseNoTrans,
......@@ -123,6 +124,7 @@ void morse_pztpgqrt( int L,
ldqm = BLKLDD(Q1, m);
for (n = k; n < Q1->nt; n++) {
tempnn = n == Q1->nt-1 ? Q1->n-n*Q1->nb : Q1->nb;
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseNoTrans,
......
......@@ -86,6 +86,7 @@ void morse_pztpqrt( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T,
tempmm = m == B->mt-1 ? B->m-m*B->mb : B->mb;
templm = m == maxmt-1 ? tempmm : 0;
ldbm = BLKLDD(B, m);
/* TT kernel */
MORSE_TASK_ztpqrt(
&options,
tempmm, tempkn, templm, ib, T->nb,
......
......@@ -111,6 +111,7 @@ void morse_pzunglq(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T, MORSE_desc
RUNTIME_data_migrate( sequence, Q(m, k),
Q->get_rankof( Q, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
MorseRight, MorseNoTrans,
......
......@@ -103,12 +103,13 @@ void morse_pzunglq_param(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_des
tempnn = n == Q->nt-1 ? Q->n-n*Q->nb : Q->nb;
/* TS or TT */
if(qrtree->gettype(qrtree, k, n) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempnn;
T = TT;
}
......
......@@ -108,6 +108,7 @@ void morse_pzunglqrh(MORSE_desc_t *A, MORSE_desc_t *Q,
RUNTIME_data_migrate( sequence, Q(m, N+RD),
Q->get_rankof( Q, m, N+RD ) );
/* TT kernel */
MORSE_TASK_ztpmlqt(
&options,
MorseRight, MorseNoTrans,
......@@ -134,6 +135,7 @@ void morse_pzunglqrh(MORSE_desc_t *A, MORSE_desc_t *Q,
RUNTIME_data_migrate( sequence, Q(m, n),
Q->get_rankof( Q, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
MorseRight, MorseNoTrans,
......
......@@ -113,6 +113,7 @@ void morse_pzungqr(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T, MORSE_desc
RUNTIME_data_migrate( sequence, Q(k, n),
Q->get_rankof( Q, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseNoTrans,
......
......@@ -109,12 +109,13 @@ void morse_pzungqr_param(const libhqr_tree_t *qrtree,
ldqm = BLKLDD(Q, m);
ldqp = BLKLDD(Q, p);
/* TT or TS */
if(qrtree->gettype(qrtree, k , m) == 0) {
/* TS kernel */
T = TS;
L = 0;
}
else {
/* TT kernel */
T = TT;
L = tempmm;
}
......
......@@ -111,6 +111,7 @@ void morse_pzungqrrh(MORSE_desc_t *A, MORSE_desc_t *Q,
RUNTIME_data_migrate( sequence, Q(M+RD, n),
Q->get_rankof( Q, M+RD, n ) );
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseNoTrans,
......@@ -140,6 +141,7 @@ void morse_pzungqrrh(MORSE_desc_t *A, MORSE_desc_t *Q,
RUNTIME_data_migrate( sequence, Q(m, n),
Q->get_rankof( Q, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
MorseLeft, MorseNoTrans,
......
......@@ -140,6 +140,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(k, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -180,6 +181,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(k, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -242,6 +244,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, k),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -329,6 +332,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, k),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......
......@@ -144,12 +144,13 @@ void morse_pzunmlq_param(const libhqr_tree_t *qrtree,
ldbp = BLKLDD(B, p);
ldbm = BLKLDD(B, m);
/* TT or TS */
if(qrtree->gettype(qrtree, k, m) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = A->nb;
T = TT;
}
......@@ -202,12 +203,13 @@ void morse_pzunmlq_param(const libhqr_tree_t *qrtree,
ldbp = BLKLDD(B, p);
ldbm = BLKLDD(B, m);
/* TT or TS */
if(qrtree->gettype(qrtree, k, m) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = A->nb;
T = TT;
}
......@@ -292,12 +294,13 @@ void morse_pzunmlq_param(const libhqr_tree_t *qrtree,
tempnn = n == B->nt-1 ? B->n-n*B->nb : B->nb;
ldbp = BLKLDD(B, p);
/* TS or TT */
if(qrtree->gettype(qrtree, k, n) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempnn;
T = TT;
}
......@@ -416,10 +419,12 @@ void morse_pzunmlq_param(const libhqr_tree_t *qrtree,
ldbp = BLKLDD(B, p);
if(qrtree->gettype(qrtree, k, n) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempnn;
T = TT;
}
......
......@@ -139,6 +139,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options, side, trans,
tempmm, tempnn, tempkm, 0, ib, T->nb,
......@@ -162,6 +163,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(N+RD, n),
B->get_rankof( B, N+RD, n ) );
/* TT kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -207,6 +209,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(N+RD, n),
B->get_rankof( B, N+RD, n ) );
/* TT kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -233,6 +236,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -302,6 +306,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, N+RD),
B->get_rankof( B, m, N+RD ) );
/* TT kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -327,6 +332,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, m),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -420,6 +426,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......@@ -443,6 +450,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, N+RD),
B->get_rankof( B, m, N+RD ) );
/* TT kernel */
MORSE_TASK_ztpmlqt(
&options,
side, trans,
......
......@@ -141,6 +141,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(k, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
side, trans,
......@@ -182,6 +183,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(k, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
side, trans,
......@@ -246,6 +248,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, k),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
side, trans,
......@@ -334,6 +337,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, k),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options,
side, trans,
......
......@@ -145,12 +145,13 @@ void morse_pzunmqr_param(const libhqr_tree_t *qrtree,
ldbm = BLKLDD(B, m);
ldbp = BLKLDD(B, p);
/* TT or TS */
if(qrtree->gettype(qrtree, k, m) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempmm;
T = TT;
}
......@@ -203,12 +204,13 @@ void morse_pzunmqr_param(const libhqr_tree_t *qrtree,
ldbm = BLKLDD(B, m);
ldbp = BLKLDD(B, p);
/* TT or TS */
if(qrtree->gettype(qrtree, k, m) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempmm;
T = TT;
}
......@@ -294,12 +296,13 @@ void morse_pzunmqr_param(const libhqr_tree_t *qrtree,
ldan = BLKLDD(A, n);
ldbp = BLKLDD(B, p);
/* TS or TT */
if(qrtree->gettype(qrtree, k, n) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempmm;
T = TT;
}
......@@ -419,10 +422,12 @@ void morse_pzunmqr_param(const libhqr_tree_t *qrtree,
ldbp = BLKLDD(B, p);
if(qrtree->gettype(qrtree, k, n) == 0){
/* TS kernel */
L = 0;
T = TS;
}
else {
/* TT kernel */
L = tempmm;
T = TT;
}
......
......@@ -140,6 +140,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempnn, tempkn, 0, ib, T->nb,
......@@ -164,6 +165,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(M+RD, n),
B->get_rankof( B, M+RD, n ) );
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempMRDm, tempnn, tempkn, tempMRDm, ib, T->nb,
......@@ -208,6 +210,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(M+RD, n),
B->get_rankof( B, M+RD, n ) );
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempMRDm, tempnn, tempkn, tempMRDm, ib, T->nb,
......@@ -235,6 +238,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempnn, tempkn, 0, ib, T->nb,
......@@ -301,6 +305,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, M+RD),
B->get_rankof( B, m, M+RD ) );
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempMRDm, tempkn, tempmm, ib, T->nb,
......@@ -328,6 +333,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, m),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempnn, tempkn, 0, ib, T->nb,
......@@ -419,6 +425,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, n),
B->get_rankof( B, m, n ) );
/* TS kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempnn, tempkn, 0, ib, T->nb,
......@@ -442,6 +449,7 @@ void morse_pzunmqrrh( MORSE_enum side, MORSE_enum trans,
RUNTIME_data_migrate( sequence, B(m, M+RD),
B->get_rankof( B, m, M+RD ) );
/* TT kernel */
MORSE_TASK_ztpmqrt(
&options, side, trans,
tempmm, tempMRDm, tempkn, tempmm, ib, T->nb,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment