Commit 5ee2abef authored by Mathieu Faverge's avatar Mathieu Faverge

RIP morse

parent 165fe037
......@@ -57,9 +57,9 @@ set(CHAMELEON_CONTROL
../control/descriptor.c
../control/workspace.c
../control/tile.c
../control/morse_f77.c
../control/morse_mf77.c
# ../control/morsewinthread.c
../control/chameleon_f77.c
../control/chameleon_mf77.c
# ../control/chameleonwinthread.c
)
set(flags_to_add "")
......@@ -68,7 +68,7 @@ foreach(_prec ${CHAMELEON_PRECISION})
endforeach()
set_source_files_properties(../control/tile.c PROPERTIES COMPILE_FLAGS "${flags_to_add}")
# Generate the morse sources for all possible precisions
# Generate the chameleon sources for all possible precisions
# ------------------------------------------------------
set(CHAMELEON_SRCS_GENERATED "")
set(ZSRC
......@@ -251,7 +251,7 @@ precisions_rules_py(CHAMELEON_SRCS_GENERATED "${ZSRC}"
set(CONTROL_SRCS_GENERATED "")
set(ZSRC
../control/workspace_z.c
../control/morse_zf77.c
../control/chameleon_zf77.c
)
precisions_rules_py(CONTROL_SRCS_GENERATED "${ZSRC}"
......@@ -264,21 +264,21 @@ set(CHAMELEON_SRCS
${CONTROL_SRCS_GENERATED}
)
# Generate the morse fortran sources for all possible precisions
# Generate the chameleon fortran sources for all possible precisions
# --------------------------------------------------------------
if(HAVE_ISO_C_BINDING)
set(CHAMELEON_SRCS_F_GENERATED "")
set(ZSRCF
../control/morse_zcf90.F90
../control/morse_zf90.F90
../control/morse_zf90_wrappers.F90
../control/chameleon_zcf90.F90
../control/chameleon_zf90.F90
../control/chameleon_zf90_wrappers.F90
)
precisions_rules_py(CHAMELEON_SRCS_F_GENERATED "${ZSRCF}"
PRECISIONS "${CHAMELEON_PRECISION}"
TARGETDIR "control" )
set(CHAMELEON_SRCSF
../control/morse_f90.f90
../control/chameleon_f90.f90
${CHAMELEON_SRCS_F_GENERATED}
)
endif(HAVE_ISO_C_BINDING)
......
......@@ -53,19 +53,19 @@
* Identifies this function call (for exception handling purposes).
*
*/
void morse_pzbuild( cham_uplo_t uplo, CHAM_desc_t *A, void *user_data, void* user_build_callback,
void chameleon_pzbuild( cham_uplo_t uplo, CHAM_desc_t *A, void *user_data, void* user_build_callback,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request )
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
int m, n;
int ldam;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
for (m = 0; m < A->mt; m++) {
ldam = BLKLDD(A, m);
......@@ -81,5 +81,5 @@ void morse_pzbuild( cham_uplo_t uplo, CHAM_desc_t *A, void *user_data, void* use
}
}
RUNTIME_options_finalize( &options, morse);
RUNTIME_options_finalize( &options, chamctxt);
}
......@@ -20,7 +20,7 @@
*/
#include "control/common.h"
void morse_pzgebrd_ge2gb(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
void chameleon_pzgebrd_ge2gb(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
int k;
......@@ -31,34 +31,34 @@ void morse_pzgebrd_ge2gb(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
for (k = 0; k < A->nt; k++) {
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
A1 = morse_desc_submatrix(A, k*A->mb, k*A->nb, A->m-k*A->mb, tempkn);
A2 = morse_desc_submatrix(A, k*A->mb, (k+1)*A->nb, A->m-k*A->mb, A->n-(k+1)*A->nb);
T1 = morse_desc_submatrix(T, k*T->mb, k*T->nb, T->m-k*T->mb, T->nb );
A1 = chameleon_desc_submatrix(A, k*A->mb, k*A->nb, A->m-k*A->mb, tempkn);
A2 = chameleon_desc_submatrix(A, k*A->mb, (k+1)*A->nb, A->m-k*A->mb, A->n-(k+1)*A->nb);
T1 = chameleon_desc_submatrix(T, k*T->mb, k*T->nb, T->m-k*T->mb, T->nb );
if ( D != NULL ) {
D1 = morse_desc_submatrix(D, k*D->mb, k*D->nb, D->m-k*D->mb, tempkn);
D1 = chameleon_desc_submatrix(D, k*D->mb, k*D->nb, D->m-k*D->mb, tempkn);
}
morse_pzgeqrf( A1, T1, D1,
chameleon_pzgeqrf( A1, T1, D1,
sequence, request);
morse_pzunmqr( ChamLeft, ChamConjTrans,
chameleon_pzunmqr( ChamLeft, ChamConjTrans,
A1, A2, T1, D1,
sequence, request);
if (k+1 < A->nt){
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
A1 = morse_desc_submatrix(A, k*A->mb, (k+1)*A->nb, tempkm, A->n-(k+1)*A->nb);
A2 = morse_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
T1 = morse_desc_submatrix(T, k*T->mb, (k+1)*T->nb, T->mb, T->n-(k+1)*T->nb);
A1 = chameleon_desc_submatrix(A, k*A->mb, (k+1)*A->nb, tempkm, A->n-(k+1)*A->nb);
A2 = chameleon_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
T1 = chameleon_desc_submatrix(T, k*T->mb, (k+1)*T->nb, T->mb, T->n-(k+1)*T->nb);
if ( D != NULL ) {
D1 = morse_desc_submatrix(D, k*D->mb, (k+1)*D->nb, tempkm, D->n-(k+1)*D->nb);
D1 = chameleon_desc_submatrix(D, k*D->mb, (k+1)*D->nb, tempkm, D->n-(k+1)*D->nb);
}
morse_pzgelqf( A1, T1, D1,
chameleon_pzgelqf( A1, T1, D1,
sequence, request);
morse_pzunmlq( ChamRight, ChamConjTrans,
chameleon_pzunmlq( ChamRight, ChamConjTrans,
A1, A2, T1, D1,
sequence, request);
}
......@@ -68,33 +68,33 @@ void morse_pzgebrd_ge2gb(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
for (k = 0; k < A->mt; k++) {
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
A1 = morse_desc_submatrix(A, k*A->mb, k*A->nb, tempkm, A->n-k*A->nb);
A2 = morse_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, A->n-k*A->nb);
T1 = morse_desc_submatrix(T, k*T->mb, k*T->nb, T->mb, T->n-k*T->nb);
A1 = chameleon_desc_submatrix(A, k*A->mb, k*A->nb, tempkm, A->n-k*A->nb);
A2 = chameleon_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, A->n-k*A->nb);
T1 = chameleon_desc_submatrix(T, k*T->mb, k*T->nb, T->mb, T->n-k*T->nb);
if ( D != NULL ) {
D1 = morse_desc_submatrix(D, k*D->mb, k*D->nb, tempkm, D->n-k*D->nb);
D1 = chameleon_desc_submatrix(D, k*D->mb, k*D->nb, tempkm, D->n-k*D->nb);
}
morse_pzgelqf( A1, T1, D1,
chameleon_pzgelqf( A1, T1, D1,
sequence, request);
morse_pzunmlq( ChamRight, ChamConjTrans,
chameleon_pzunmlq( ChamRight, ChamConjTrans,
A1, A2, T1, D1,
sequence, request);
if (k+1 < A->mt){
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
A1 = morse_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, tempkn);
A2 = morse_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
T1 = morse_desc_submatrix(T, (k+1)*T->mb, k*T->nb, T->m-(k+1)*T->mb, T->nb );
A1 = chameleon_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, tempkn);
A2 = chameleon_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
T1 = chameleon_desc_submatrix(T, (k+1)*T->mb, k*T->nb, T->m-(k+1)*T->mb, T->nb );
if ( D != NULL ) {
D1 = morse_desc_submatrix(D, (k+1)*D->mb, k*D->nb, D->m-(k+1)*D->mb, tempkn);
D1 = chameleon_desc_submatrix(D, (k+1)*D->mb, k*D->nb, D->m-(k+1)*D->mb, tempkn);
}
morse_pzgeqrf( A1, T1, D1,
chameleon_pzgeqrf( A1, T1, D1,
sequence, request);
morse_pzunmqr( ChamLeft, ChamConjTrans,
chameleon_pzunmqr( ChamLeft, ChamConjTrans,
A1, A2, T1, D1,
sequence, request);
}
......
......@@ -36,10 +36,10 @@
/**
* Parallel tile LQ factorization - dynamic scheduling
*/
void morse_pzgelqf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
void chameleon_pzgelqf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
size_t ws_worker = 0;
size_t ws_host = 0;
......@@ -49,10 +49,10 @@ void morse_pzgelqf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
int tempkm, tempkn, tempmm, tempnn;
int ib, minMNT;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -90,7 +90,7 @@ void morse_pzgelqf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
RUNTIME_options_ws_alloc( &options, ws_worker, ws_host );
for (k = 0; k < minMNT; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
......@@ -169,10 +169,10 @@ void morse_pzgelqf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
A->get_rankof( A, m, k ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
(void)D;
}
......@@ -30,11 +30,11 @@
/*
* Parallel tile LQ factorization (reduction Householder) - dynamic scheduling
*/
void morse_pzgelqf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
void chameleon_pzgelqf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
CHAM_desc_t *TS, CHAM_desc_t *TT, CHAM_desc_t *D,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
CHAM_desc_t *T;
size_t ws_worker = 0;
......@@ -47,10 +47,10 @@ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
int ib;
int *tiles;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -91,7 +91,7 @@ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
/* The number of the factorization */
for (k = 0; k < K; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
ldak = BLKLDD(A, k);
......@@ -199,10 +199,10 @@ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
A->get_rankof( A, m, k ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
free(tiles);
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
}
......@@ -38,10 +38,10 @@
/*
* Parallel tile LQ factorization (reduction Householder) - dynamic scheduling
*/
void morse_pzgelqfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
void chameleon_pzgelqfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
size_t ws_worker = 0;
size_t ws_host = 0;
......@@ -52,10 +52,10 @@ void morse_pzgelqfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
int tempkmin, tempkm, tempNn, tempnn, tempmm, tempNRDn;
int ib;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -88,7 +88,7 @@ void morse_pzgelqfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
/* The number of the factorization */
for (k = 0; k < K; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
ldak = BLKLDD(A, k);
......@@ -210,10 +210,10 @@ void morse_pzgelqfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
A->get_rankof( A, m, k ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
(void)D;
}
......@@ -29,12 +29,12 @@
/**
* Parallel tile matrix-matrix multiplication - dynamic scheduling
*/
void morse_pzgemm(cham_trans_t transA, cham_trans_t transB,
void chameleon_pzgemm(cham_trans_t transA, cham_trans_t transB,
CHAMELEON_Complex64_t alpha, CHAM_desc_t *A, CHAM_desc_t *B,
CHAMELEON_Complex64_t beta, CHAM_desc_t *C,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
int m, n, k;
......@@ -44,10 +44,10 @@ void morse_pzgemm(cham_trans_t transA, cham_trans_t transB,
CHAMELEON_Complex64_t zbeta;
CHAMELEON_Complex64_t zone = (CHAMELEON_Complex64_t)1.0;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
for (m = 0; m < C->mt; m++) {
tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb;
......@@ -141,5 +141,5 @@ void morse_pzgemm(cham_trans_t transA, cham_trans_t transB,
}
}
}
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
}
......@@ -36,10 +36,10 @@
/**
* Parallel tile QR factorization - dynamic scheduling
*/
void morse_pzgeqrf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
void chameleon_pzgeqrf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
size_t ws_worker = 0;
size_t ws_host = 0;
......@@ -50,10 +50,10 @@ void morse_pzgeqrf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
int ib;
int minMNT = chameleon_min(A->mt, A->nt);
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -85,7 +85,7 @@ void morse_pzgeqrf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
RUNTIME_options_ws_alloc( &options, ws_worker, ws_host );
for (k = 0; k < minMNT; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
......@@ -165,10 +165,10 @@ void morse_pzgeqrf(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D,
A->get_rankof( A, k, n ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
(void)D;
}
......@@ -29,11 +29,11 @@
/**
* Parallel tile QR factorization (reduction Householder) - dynamic scheduling
*/
void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
void chameleon_pzgeqrf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
CHAM_desc_t *TS, CHAM_desc_t *TT, CHAM_desc_t *D,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
CHAM_desc_t *T;
size_t ws_worker = 0;
......@@ -46,10 +46,10 @@ void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
int ib;
int *tiles;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -89,7 +89,7 @@ void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
/* The number of the factorization */
for (k = 0; k < K; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
/* The number of geqrt to apply */
......@@ -197,10 +197,10 @@ void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, CHAM_desc_t *A,
A->get_rankof( A, k, n ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
free(tiles);
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
}
......@@ -38,10 +38,10 @@
/**
* Parallel tile QR factorization (reduction Householder) - dynamic scheduling
*/
void morse_pzgeqrfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
void chameleon_pzgeqrfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request)
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
size_t ws_worker = 0;
size_t ws_host = 0;
......@@ -52,10 +52,10 @@ void morse_pzgeqrfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
int tempkmin, tempkn, tempMm, tempnn, tempmm, tempMRDm;
int ib;
morse = morse_context_self();
chamctxt = chameleon_context_self();
if (sequence->status != CHAMELEON_SUCCESS)
return;
RUNTIME_options_init(&options, morse, sequence, request);
RUNTIME_options_init(&options, chamctxt, sequence, request);
ib = CHAMELEON_IB;
......@@ -86,7 +86,7 @@ void morse_pzgeqrfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
K = chameleon_min(A->mt, A->nt);
for (k = 0; k < K; k++) {
RUNTIME_iteration_push(morse, k);
RUNTIME_iteration_push(chamctxt, k);
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
for (M = k; M < A->mt; M += BS) {
......@@ -209,10 +209,10 @@ void morse_pzgeqrfrh(CHAM_desc_t *A, CHAM_desc_t *T, CHAM_desc_t *D, int BS,
A->get_rankof( A, k, n ) );
}
RUNTIME_iteration_pop(morse);
RUNTIME_iteration_pop(chamctxt);
}
RUNTIME_options_ws_free(&options);
RUNTIME_options_finalize(&options, morse);
RUNTIME_options_finalize(&options, chamctxt);
(void)D;
}
......@@ -40,10 +40,10 @@
/**
* Parallel tile LU factorization - dynamic scheduling
*/
void morse_pzgetrf_incpiv( CHAM_desc_t *A, CHAM_desc_t *L, CHAM_desc_t *D, int *IPIV,
void chameleon_pzgetrf_incpiv( CHAM_desc_t *A, CHAM_desc_t *L, CHAM_desc_t *D, int *IPIV,
RUNTIME_sequence_t *sequence, RUNTIME_request_t *request )
{
CHAM_context_t *morse;
CHAM_context_t *chamctxt;
RUNTIME_option_t options;
size_t ws_worker = 0;
size_t ws_host = 0;
......@@ -54,10 +54,10 @@ void morse_pzgetrf_incpiv( CHAM_desc_t *A, CHAM_desc_t *L, CHAM_desc_t *D, int *
int ib;
int minMNT = chameleon_min(A->mt, A->nt);
morse = morse_context_self();