diff --git a/compute/pzgelqf.c b/compute/pzgelqf.c index 3cec88b5714ad8d1c0ea2e84084be1062348d7db..90e8351dd22a09d4e95b19c5810c939c7ca6899f 100644 --- a/compute/pzgelqf.c +++ b/compute/pzgelqf.c @@ -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, diff --git a/compute/pzgelqf_param.c b/compute/pzgelqf_param.c index 0977e086927cada79e5ee8c62302507d9c126f69..1c93e77899fb4020caeb68223c61a67ee2d836b7 100644 --- a/compute/pzgelqf_param.c +++ b/compute/pzgelqf_param.c @@ -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; } diff --git a/compute/pzgelqfrh.c b/compute/pzgelqfrh.c index 0ff1d84ddb0c2ff61c7d4a4e1d63f319471114a8..22a89797ae94bb876d876f2633514d021aaa0309 100644 --- a/compute/pzgelqfrh.c +++ b/compute/pzgelqfrh.c @@ -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, diff --git a/compute/pzgeqrf.c b/compute/pzgeqrf.c index f5b3b744c50d7b8b64412b0d0b7842557ab1a3bf..7146dbc8375719e5ecb623e1e0cb2f8ad21d7a11 100644 --- a/compute/pzgeqrf.c +++ b/compute/pzgeqrf.c @@ -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, diff --git a/compute/pzgeqrf_param.c b/compute/pzgeqrf_param.c index a0d76c887c159c5fce569aff010338fc64dd9b67..03e90369a5a17dbb46bd1bdf4136c40953720313 100644 --- a/compute/pzgeqrf_param.c +++ b/compute/pzgeqrf_param.c @@ -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; } diff --git a/compute/pzgeqrfrh.c b/compute/pzgeqrfrh.c index 357ca25fc69cd606027726ff01116149718687d0..a82d78315ca28125bd34ffa42a5289c2cce66f94 100644 --- a/compute/pzgeqrfrh.c +++ b/compute/pzgeqrfrh.c @@ -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, diff --git a/compute/pztpgqrt.c b/compute/pztpgqrt.c index a070dc85157fe4e1c0623c823617e4c826fe9ce3..b659d45390eca482ef869a698d37bb6e1a7e1614 100644 --- a/compute/pztpgqrt.c +++ b/compute/pztpgqrt.c @@ -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, diff --git a/compute/pztpqrt.c b/compute/pztpqrt.c index 47415ff055e846d4e690709174ec926e5b347028..6e956cff65f831876d52453ffbe4220616b39b53 100644 --- a/compute/pztpqrt.c +++ b/compute/pztpqrt.c @@ -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, diff --git a/compute/pzunglq.c b/compute/pzunglq.c index de9c1409fe92f64fddc26dc0de39554ac0cc8ee9..dcf53f600f420950660a570ba8c061424bacc770 100644 --- a/compute/pzunglq.c +++ b/compute/pzunglq.c @@ -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, diff --git a/compute/pzunglq_param.c b/compute/pzunglq_param.c index 8346cbd39d92b3917276a6043db912eefeb785de..16af09a3a2a35521f1248c671f743fb77f068d53 100644 --- a/compute/pzunglq_param.c +++ b/compute/pzunglq_param.c @@ -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; } diff --git a/compute/pzunglqrh.c b/compute/pzunglqrh.c index 2ae69ce0328649c2c0386eddd646f95c9976477a..591de3742c7733e98d9c1284638e3f99b0cc6af0 100644 --- a/compute/pzunglqrh.c +++ b/compute/pzunglqrh.c @@ -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, diff --git a/compute/pzungqr.c b/compute/pzungqr.c index e305f329f32e5a67e24fe693314a13b482ea5377..140c0e6b273c92019cd0feac31d62a200c6320e2 100644 --- a/compute/pzungqr.c +++ b/compute/pzungqr.c @@ -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, diff --git a/compute/pzungqr_param.c b/compute/pzungqr_param.c index d78459ae9e963ce8ec4d7c1e8d34b65c8834598f..849909b3f7df4d3f2902d951a8cfe6f907c0bc00 100644 --- a/compute/pzungqr_param.c +++ b/compute/pzungqr_param.c @@ -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; } diff --git a/compute/pzungqrrh.c b/compute/pzungqrrh.c index 89853c6cf6870a815c60a198443816d5a22a7f02..621d75e259ca23c0ea9525e938255207a3cfddc3 100644 --- a/compute/pzungqrrh.c +++ b/compute/pzungqrrh.c @@ -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, diff --git a/compute/pzunmlq.c b/compute/pzunmlq.c index 28def3d1140780a8052aaccb6b7c27cded251117..f2d07a65d472791a7af9ad739ae1376f1f9a1293 100644 --- a/compute/pzunmlq.c +++ b/compute/pzunmlq.c @@ -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, diff --git a/compute/pzunmlq_param.c b/compute/pzunmlq_param.c index 8fce435225b335e7b119261955f57d705497196a..d8ca7f2d77447eb6d53a7df6b920a1d70848e1cc 100644 --- a/compute/pzunmlq_param.c +++ b/compute/pzunmlq_param.c @@ -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; } diff --git a/compute/pzunmlqrh.c b/compute/pzunmlqrh.c index 6dfffc232f00eb15fde9455834b5c23a1636c4eb..f6b38ab833e945da1ffbe59b15be86e84a93d318 100644 --- a/compute/pzunmlqrh.c +++ b/compute/pzunmlqrh.c @@ -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, diff --git a/compute/pzunmqr.c b/compute/pzunmqr.c index 02202f41e0e59ba224ecc50500e6a6bd31b3d50b..140a65fcdbc91510cf16b79bd338aec12622c0e7 100644 --- a/compute/pzunmqr.c +++ b/compute/pzunmqr.c @@ -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, diff --git a/compute/pzunmqr_param.c b/compute/pzunmqr_param.c index d09db8c5ec0531e7e1defcbdb1b724ea5531267f..c9c33d6bb606f1961ce501ea656b1ba8f43b4faf 100644 --- a/compute/pzunmqr_param.c +++ b/compute/pzunmqr_param.c @@ -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; } diff --git a/compute/pzunmqrrh.c b/compute/pzunmqrrh.c index 882d587ac9471b3bb620c9fec4d706f4c253f008..08e7992d28f6514c5cea891695d5e50cf735bb06 100644 --- a/compute/pzunmqrrh.c +++ b/compute/pzunmqrrh.c @@ -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,