diff --git a/compute/zbuild.c b/compute/zbuild.c index e95654787aa7cc4c6b962bd4cc629e34deff8774..7fe753ef6ae56796ab8041ed1fdf7181727c132c 100644 --- a/compute/zbuild.c +++ b/compute/zbuild.c @@ -124,7 +124,7 @@ int MORSE_zbuild( MORSE_enum uplo, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescOutput, uplo, @@ -137,13 +137,13 @@ int MORSE_zbuild( MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescOutput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -203,13 +203,15 @@ int MORSE_zbuild_Tile( MORSE_enum uplo, MORSE_desc_t *A, morse_fatal_error("MORSE_zbuild_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zbuild_Tile_Async( uplo, A, user_data, user_build_callback, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/zgeadd.c b/compute/zgeadd.c index 4bc293e2f87b2d8524c3ef04bb263d751fe53d42..5478a00aab4624a466d1df0536ed61a8fea76604 100644 --- a/compute/zgeadd.c +++ b/compute/zgeadd.c @@ -91,9 +91,9 @@ * @sa MORSE_sgeadd * ******************************************************************************/ -int MORSE_zgeadd(MORSE_enum trans, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t beta, MORSE_Complex64_t *B, int LDB) +int MORSE_zgeadd( MORSE_enum trans, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t beta, MORSE_Complex64_t *B, int LDB ) { int NB; int Am, An; @@ -152,7 +152,7 @@ int MORSE_zgeadd(MORSE_enum trans, int M, int N, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -169,14 +169,14 @@ int MORSE_zgeadd(MORSE_enum trans, int M, int N, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -228,9 +228,9 @@ int MORSE_zgeadd(MORSE_enum trans, int M, int N, * @sa MORSE_sgeadd_Tile * ******************************************************************************/ -int MORSE_zgeadd_Tile(MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *B) +int MORSE_zgeadd_Tile( MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -242,15 +242,16 @@ int MORSE_zgeadd_Tile(MORSE_enum trans, morse_fatal_error("MORSE_zgeadd_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgeadd_Tile_Async( trans, alpha, A, beta, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -283,10 +284,10 @@ int MORSE_zgeadd_Tile(MORSE_enum trans, * @sa MORSE_sgeadd_Tile_Async * ******************************************************************************/ -int MORSE_zgeadd_Tile_Async(MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgeadd_Tile_Async( MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int M, N; @@ -361,9 +362,11 @@ int MORSE_zgeadd_Tile_Async(MORSE_enum trans, N = B->n; /* Quick return */ - if (M == 0 || N == 0 || - ((alpha == (MORSE_Complex64_t)0.0) && beta == (MORSE_Complex64_t)1.0)) + if ( (M == 0) || (N == 0) || + ((alpha == (MORSE_Complex64_t)0.0) && (beta == (MORSE_Complex64_t)1.0)) ) + { return MORSE_SUCCESS; + } morse_pztradd( MorseUpperLower, trans, alpha, A, beta, B, sequence, request ); diff --git a/compute/zgelqf.c b/compute/zgelqf.c index 3b9a447d3ce13dc440e5030e09cb3d58287dc289..f21a09861d367ee3a418efd78ac55952be647aa4 100644 --- a/compute/zgelqf.c +++ b/compute/zgelqf.c @@ -75,9 +75,9 @@ * @sa MORSE_zgelqs * ******************************************************************************/ -int MORSE_zgelqf(int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT) +int MORSE_zgelqf( int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT ) { int NB; int status; @@ -120,7 +120,7 @@ int MORSE_zgelqf(int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -133,13 +133,13 @@ int MORSE_zgelqf(int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -194,14 +194,16 @@ int MORSE_zgelqf_Tile( MORSE_desc_t *A, MORSE_desc_t *T ) morse_fatal_error("MORSE_zgelqf_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgelqf_Tile_Async( A, T, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -234,8 +236,8 @@ int MORSE_zgelqf_Tile( MORSE_desc_t *A, MORSE_desc_t *T ) * @sa MORSE_zgelqs_Tile_Async * ******************************************************************************/ -int MORSE_zgelqf_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgelqf_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zgelqf_param.c b/compute/zgelqf_param.c index d0b8b3d361537adb9ef65193bd1433fd41f9bb6e..65d497f29b2d934f9e2879b3e88ebeb77edc0981 100644 --- a/compute/zgelqf_param.c +++ b/compute/zgelqf_param.c @@ -72,9 +72,9 @@ * @sa MORSE_zgelqs * ******************************************************************************/ -int MORSE_zgelqf_param(const libhqr_tree_t *qrtree, int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT) +int MORSE_zgelqf_param( const libhqr_tree_t *qrtree, int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT ) { int NB; int status; @@ -117,7 +117,7 @@ int MORSE_zgelqf_param(const libhqr_tree_t *qrtree, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -130,13 +130,13 @@ int MORSE_zgelqf_param(const libhqr_tree_t *qrtree, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -191,14 +191,17 @@ int MORSE_zgelqf_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zgelqf_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgelqf_param_Tile_Async( qrtree, A, TS, TT, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -231,8 +234,8 @@ int MORSE_zgelqf_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE * @sa MORSE_zgelqs_Tile_Async * ******************************************************************************/ -int MORSE_zgelqf_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgelqf_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zgelqs.c b/compute/zgelqs.c index 78acc2a2687085b34795795750fd906c3865d552..30d29ccb35d56124a984251a06c69eb2071b8187 100644 --- a/compute/zgelqs.c +++ b/compute/zgelqs.c @@ -81,10 +81,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zgelqs(int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgelqs( int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -136,7 +136,7 @@ int MORSE_zgelqs(int M, int N, int NRHS, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -153,14 +153,14 @@ int MORSE_zgelqs(int M, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -215,15 +215,17 @@ int MORSE_zgelqs_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B ) morse_fatal_error("MORSE_zgelqs_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgelqs_Tile_Async( A, T, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -257,8 +259,8 @@ int MORSE_zgelqs_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B ) * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zgelqs_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgelqs_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subB; MORSE_desc_t *subA; diff --git a/compute/zgelqs_param.c b/compute/zgelqs_param.c index 28dc028e24b0caaf32ffa2e3b926105186c6d5ab..42ee682dd8c557916821fcd001dbaa2bbd83fab3 100644 --- a/compute/zgelqs_param.c +++ b/compute/zgelqs_param.c @@ -83,10 +83,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zgelqs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgelqs_param( const libhqr_tree_t *qrtree, int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -138,7 +138,7 @@ int MORSE_zgelqs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -155,14 +155,14 @@ int MORSE_zgelqs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -208,7 +208,8 @@ int MORSE_zgelqs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, * @sa MORSE_zgelqf_Tile * ******************************************************************************/ -int MORSE_zgelqs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B ) +int MORSE_zgelqs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, + MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -220,15 +221,18 @@ int MORSE_zgelqs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zgelqs_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgelqs_param_Tile_Async( qrtree, A, TS, TT, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -262,8 +266,8 @@ int MORSE_zgelqs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zgelqs_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgelqs_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subB; MORSE_desc_t *subA; diff --git a/compute/zgels.c b/compute/zgels.c index 5c6fa83135c4b3f52d5b5ce34e9de93a69194f29..9882bcdb435c81ee4bd34b5fdb1a17ee9d6f32bb 100644 --- a/compute/zgels.c +++ b/compute/zgels.c @@ -107,10 +107,10 @@ * @sa MORSE_sgels * ******************************************************************************/ -int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgels( MORSE_enum trans, int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *B, int LDB ) { int i, j; int NB; @@ -170,7 +170,7 @@ int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ if ( M >= N ) { @@ -195,14 +195,14 @@ int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -260,8 +260,8 @@ int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS, * @sa MORSE_sgels_Tile * ******************************************************************************/ -int MORSE_zgels_Tile(MORSE_enum trans, MORSE_desc_t *A, - MORSE_desc_t *T, MORSE_desc_t *B) +int MORSE_zgels_Tile( MORSE_enum trans, MORSE_desc_t *A, + MORSE_desc_t *T, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -273,15 +273,17 @@ int MORSE_zgels_Tile(MORSE_enum trans, MORSE_desc_t *A, morse_fatal_error("MORSE_zgels_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgels_Tile_Async( trans, A, T, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -314,9 +316,9 @@ int MORSE_zgels_Tile(MORSE_enum trans, MORSE_desc_t *A, * @sa MORSE_sgels_Tile_Async * ******************************************************************************/ -int MORSE_zgels_Tile_Async(MORSE_enum trans, MORSE_desc_t *A, - MORSE_desc_t *T, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgels_Tile_Async( MORSE_enum trans, MORSE_desc_t *A, + MORSE_desc_t *T, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subA; MORSE_desc_t *subB; diff --git a/compute/zgels_param.c b/compute/zgels_param.c index d8893a197d22a79e5c7ac34018df3ac850381959..dcfcd3812d8f6a6b3d1085ab1bca86d96e57d992 100644 --- a/compute/zgels_param.c +++ b/compute/zgels_param.c @@ -109,10 +109,10 @@ * @sa MORSE_sgels * ******************************************************************************/ -int MORSE_zgels_param(const libhqr_tree_t *qrtree, MORSE_enum trans, int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgels_param( const libhqr_tree_t *qrtree, MORSE_enum trans, int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *B, int LDB ) { int i, j; int NB; @@ -172,7 +172,7 @@ int MORSE_zgels_param(const libhqr_tree_t *qrtree, MORSE_enum trans, int M, int /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ if ( M >= N ) { @@ -196,14 +196,14 @@ int MORSE_zgels_param(const libhqr_tree_t *qrtree, MORSE_enum trans, int M, int morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -264,8 +264,8 @@ int MORSE_zgels_param(const libhqr_tree_t *qrtree, MORSE_enum trans, int M, int * @sa MORSE_sgels_Tile * ******************************************************************************/ -int MORSE_zgels_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_desc_t *A, - MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B) +int MORSE_zgels_param_Tile( const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_desc_t *A, + MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -277,15 +277,18 @@ int MORSE_zgels_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_ morse_fatal_error("MORSE_zgels_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgels_param_Tile_Async( qrtree, trans, A, TS, TT, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -318,9 +321,9 @@ int MORSE_zgels_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_ * @sa MORSE_sgels_Tile_Async * ******************************************************************************/ -int MORSE_zgels_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_desc_t *A, - MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgels_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_enum trans, MORSE_desc_t *A, + MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subA; MORSE_desc_t *subB; diff --git a/compute/zgemm.c b/compute/zgemm.c index ca9fd933631299f7a71e1c5d932656f530e476a8..92101a69bb9e0f7ea48f6101e7332924e4837b1c 100644 --- a/compute/zgemm.c +++ b/compute/zgemm.c @@ -130,10 +130,10 @@ * @sa MORSE_sgemm * ******************************************************************************/ -int MORSE_zgemm(MORSE_enum transA, MORSE_enum transB, int M, int N, int K, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB, - MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC) +int MORSE_zgemm( MORSE_enum transA, MORSE_enum transB, int M, int N, int K, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB, + MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am, An, Bm, Bn; @@ -210,7 +210,7 @@ int MORSE_zgemm(MORSE_enum transA, MORSE_enum transB, int M, int N, int K, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -231,7 +231,7 @@ int MORSE_zgemm(MORSE_enum transA, MORSE_enum transB, int M, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -239,7 +239,7 @@ int MORSE_zgemm(MORSE_enum transA, MORSE_enum transB, int M, int N, int K, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -300,9 +300,9 @@ int MORSE_zgemm(MORSE_enum transA, MORSE_enum transB, int M, int N, int K, * @sa MORSE_sgemm_Tile * ******************************************************************************/ -int MORSE_zgemm_Tile(MORSE_enum transA, MORSE_enum transB, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C) +int MORSE_zgemm_Tile( MORSE_enum transA, MORSE_enum transB, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -314,16 +314,17 @@ int MORSE_zgemm_Tile(MORSE_enum transA, MORSE_enum transB, morse_fatal_error("MORSE_zgemm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgemm_Tile_Async( transA, transB, alpha, A, B, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -355,10 +356,10 @@ int MORSE_zgemm_Tile(MORSE_enum transA, MORSE_enum transB, * @sa MORSE_sgemm_Tile_Async * ******************************************************************************/ -int MORSE_zgemm_Tile_Async(MORSE_enum transA, MORSE_enum transB, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgemm_Tile_Async( MORSE_enum transA, MORSE_enum transB, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int M, N, K; @@ -459,9 +460,11 @@ int MORSE_zgemm_Tile_Async(MORSE_enum transA, MORSE_enum transB, K = An; /* Quick return */ - if (M == 0 || N == 0 || - ((alpha == (MORSE_Complex64_t)0.0 || K == 0) && beta == (MORSE_Complex64_t)1.0)) + if ( (M == 0) || (N == 0) || + (((alpha == (MORSE_Complex64_t)0.0) || (K == 0)) && (beta == (MORSE_Complex64_t)1.0)) ) + { return MORSE_SUCCESS; + } morse_pzgemm( transA, transB, alpha, A, B, beta, C, sequence, request ); diff --git a/compute/zgeqrf.c b/compute/zgeqrf.c index e2baf0a4fc80e07443ed2dc465a9bf2fd30555f2..abca5cd8c6ffaf6426f136d6e34b885e99a66b8d 100644 --- a/compute/zgeqrf.c +++ b/compute/zgeqrf.c @@ -74,9 +74,9 @@ * @sa MORSE_zgeqrs * ******************************************************************************/ -int MORSE_zgeqrf(int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT) +int MORSE_zgeqrf( int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT ) { int NB; int status; @@ -119,7 +119,7 @@ int MORSE_zgeqrf(int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -132,13 +132,13 @@ int MORSE_zgeqrf(int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -193,14 +193,16 @@ int MORSE_zgeqrf_Tile( MORSE_desc_t *A, MORSE_desc_t *T ) morse_fatal_error("MORSE_zgeqrf_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgeqrf_Tile_Async( A, T, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -233,8 +235,8 @@ int MORSE_zgeqrf_Tile( MORSE_desc_t *A, MORSE_desc_t *T ) * @sa MORSE_zgeqrs_Tile_Async * ******************************************************************************/ -int MORSE_zgeqrf_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgeqrf_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zgeqrf_param.c b/compute/zgeqrf_param.c index 17acd0c69469de7ebfd47ddfaff76810cc7e5f0a..421415824e465e690b3f1cdc55956298aef921d1 100644 --- a/compute/zgeqrf_param.c +++ b/compute/zgeqrf_param.c @@ -77,9 +77,9 @@ * @sa MORSE_zgeqrs * ******************************************************************************/ -int MORSE_zgeqrf_param(const libhqr_tree_t *qrtree, int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT) +int MORSE_zgeqrf_param( const libhqr_tree_t *qrtree, int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT ) { int NB; int status; @@ -122,7 +122,7 @@ int MORSE_zgeqrf_param(const libhqr_tree_t *qrtree, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -135,13 +135,13 @@ int MORSE_zgeqrf_param(const libhqr_tree_t *qrtree, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -200,14 +200,17 @@ int MORSE_zgeqrf_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zgeqrf_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgeqrf_param_Tile_Async( qrtree, A, TS, TT, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -240,8 +243,8 @@ int MORSE_zgeqrf_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE * @sa MORSE_zgeqrs_param_Tile_Async * ******************************************************************************/ -int MORSE_zgeqrf_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgeqrf_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zgeqrs.c b/compute/zgeqrs.c index 31820a5045fe7fc3cb316e9d54631724cd65d664..9f71621df5defb5ae8069a4bca4e8a1abd310f95 100644 --- a/compute/zgeqrs.c +++ b/compute/zgeqrs.c @@ -81,10 +81,10 @@ * @sa MORSE_zgeqrf * *******************************************************************************/ -int MORSE_zgeqrs(int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgeqrs( int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -136,7 +136,7 @@ int MORSE_zgeqrs(int M, int N, int NRHS, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -153,14 +153,14 @@ int MORSE_zgeqrs(int M, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -214,16 +214,17 @@ int MORSE_zgeqrs_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B ) morse_fatal_error("MORSE_zgeqrs_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgeqrs_Tile_Async( A, T, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -257,8 +258,8 @@ int MORSE_zgeqrs_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B ) * @sa MORSE_zgeqrf_Tile_Async * ******************************************************************************/ -int MORSE_zgeqrs_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgeqrs_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subA; MORSE_desc_t *subB; diff --git a/compute/zgeqrs_param.c b/compute/zgeqrs_param.c index ecb20508e12726f99d7c8d80ee28c9858e3622d0..e25b2fd26a56e7bbf258b3eb216ee6a503ab060a 100644 --- a/compute/zgeqrs_param.c +++ b/compute/zgeqrs_param.c @@ -76,10 +76,10 @@ * @sa MORSE_zgeqrf * *******************************************************************************/ -int MORSE_zgeqrs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgeqrs_param( const libhqr_tree_t *qrtree, int M, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -131,7 +131,7 @@ int MORSE_zgeqrs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -148,14 +148,14 @@ int MORSE_zgeqrs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -197,7 +197,8 @@ int MORSE_zgeqrs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS, * @sa MORSE_zgeqrf_Tile * ******************************************************************************/ -int MORSE_zgeqrs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B ) +int MORSE_zgeqrs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, + MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -209,17 +210,18 @@ int MORSE_zgeqrs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zgeqrs_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgeqrs_param_Tile_Async( qrtree, A, TS, TT, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( TS, sequence ); MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -253,9 +255,9 @@ int MORSE_zgeqrs_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE * @sa MORSE_zgeqrf_Tile_Async * ******************************************************************************/ -int MORSE_zgeqrs_param_Tile_Async(const libhqr_tree_t *qrtree, - MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgeqrs_param_Tile_Async( const libhqr_tree_t *qrtree, + MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t *subA; MORSE_desc_t *subB; diff --git a/compute/zgesv_incpiv.c b/compute/zgesv_incpiv.c index b6d28248561fedcc4b3c2310474855c4540c5178..700f6b02e960c8723090f567f33107b467f6dddc 100644 --- a/compute/zgesv_incpiv.c +++ b/compute/zgesv_incpiv.c @@ -86,10 +86,10 @@ * @sa MORSE_sgesv_incpiv * ******************************************************************************/ -int MORSE_zgesv_incpiv(int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descL, int *IPIV, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgesv_incpiv( int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descL, int *IPIV, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -135,7 +135,7 @@ int MORSE_zgesv_incpiv(int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -152,14 +152,14 @@ int MORSE_zgesv_incpiv(int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -219,15 +219,17 @@ int MORSE_zgesv_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_ morse_fatal_error("MORSE_zgesv_incpiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgesv_incpiv_Tile_Async( A, L, IPIV, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( L, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -261,8 +263,8 @@ int MORSE_zgesv_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_ * @sa MORSE_zcgesv_Tile_Async * ******************************************************************************/ -int MORSE_zgesv_incpiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgesv_incpiv_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zgesv_nopiv.c b/compute/zgesv_nopiv.c index e1f5d33c4bceeebbfadf36d1d63fd0c062e11f98..0a5753ad3753d7cef81d1f2a6f91be709d16df02 100644 --- a/compute/zgesv_nopiv.c +++ b/compute/zgesv_nopiv.c @@ -86,9 +86,9 @@ * @sa MORSE_sgesv_nopiv * ******************************************************************************/ -int MORSE_zgesv_nopiv(int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgesv_nopiv( int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -134,7 +134,7 @@ int MORSE_zgesv_nopiv(int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -151,14 +151,14 @@ int MORSE_zgesv_nopiv(int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -211,15 +211,16 @@ int MORSE_zgesv_nopiv_Tile( MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zgesv_nopiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgesv_nopiv_Tile_Async( A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -253,8 +254,8 @@ int MORSE_zgesv_nopiv_Tile( MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_zcgesv_Tile_Async * ******************************************************************************/ -int MORSE_zgesv_nopiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgesv_nopiv_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zgesvd.c b/compute/zgesvd.c index 662794846ccec4d9b03333ea79ba2f254d514298..ba35c4da4e0164dbe0bb319ffea38209cbaa41b1 100644 --- a/compute/zgesvd.c +++ b/compute/zgesvd.c @@ -147,13 +147,13 @@ * @sa MORSE_sgesvd * ******************************************************************************/ -int MORSE_zgesvd(MORSE_enum jobu, MORSE_enum jobvt, - int M, int N, - MORSE_Complex64_t *A, int LDA, - double *S, - MORSE_desc_t *descT, - MORSE_Complex64_t *U, int LDU, - MORSE_Complex64_t *VT, int LDVT) +int MORSE_zgesvd( MORSE_enum jobu, MORSE_enum jobvt, + int M, int N, + MORSE_Complex64_t *A, int LDA, + double *S, + MORSE_desc_t *descT, + MORSE_Complex64_t *U, int LDU, + MORSE_Complex64_t *VT, int LDVT ) { int NB; int status; @@ -212,7 +212,7 @@ int MORSE_zgesvd(MORSE_enum jobu, MORSE_enum jobvt, /* Set MT, NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -225,13 +225,13 @@ int MORSE_zgesvd(MORSE_enum jobu, MORSE_enum jobvt, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -335,16 +335,16 @@ int MORSE_zgesvd(MORSE_enum jobu, MORSE_enum jobvt, * @sa MORSE_sgesvd_Tile * ******************************************************************************/ -int MORSE_zgesvd_Tile(MORSE_enum jobu, MORSE_enum jobvt, - MORSE_desc_t *A, - double *S, - MORSE_desc_t *T, - MORSE_Complex64_t *U, int LDU, - MORSE_Complex64_t *VT, int LDVT) +int MORSE_zgesvd_Tile( MORSE_enum jobu, MORSE_enum jobvt, + MORSE_desc_t *A, + double *S, + MORSE_desc_t *T, + MORSE_Complex64_t *U, int LDU, + MORSE_Complex64_t *VT, int LDVT ) { - MORSE_context_t *morse; + MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; - MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + MORSE_request_t request = MORSE_REQUEST_INITIALIZER; int status; morse = morse_context_self(); @@ -352,13 +352,16 @@ int MORSE_zgesvd_Tile(MORSE_enum jobu, MORSE_enum jobvt, morse_fatal_error("MORSE_zgesvd_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgesvd_Tile_Async( jobu, jobvt, A, S, T, U, LDU, VT, LDVT, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -392,13 +395,13 @@ int MORSE_zgesvd_Tile(MORSE_enum jobu, MORSE_enum jobvt, * @sa MORSE_sgesvd_Tile_Async * ******************************************************************************/ -int MORSE_zgesvd_Tile_Async(MORSE_enum jobu, MORSE_enum jobvt, - MORSE_desc_t *A, - double *S, - MORSE_desc_t *T, - MORSE_Complex64_t *U, int LDU, - MORSE_Complex64_t *VT, int LDVT, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgesvd_Tile_Async( MORSE_enum jobu, MORSE_enum jobvt, + MORSE_desc_t *A, + double *S, + MORSE_desc_t *T, + MORSE_Complex64_t *U, int LDU, + MORSE_Complex64_t *VT, int LDVT, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_desc_t descA; MORSE_desc_t descT; @@ -535,7 +538,7 @@ int MORSE_zgesvd_Tile_Async(MORSE_enum jobu, MORSE_enum jobvt, } } - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); #if !defined(CHAMELEON_SIMULATION) info = LAPACKE_zgbbrd( LAPACK_COL_MAJOR, @@ -561,7 +564,7 @@ int MORSE_zgesvd_Tile_Async(MORSE_enum jobu, MORSE_enum jobvt, //morse_zooplap2tile( descVT, VT, NB, NB, LDVT, N, 0, 0, N, N, sequence, request, morse_desc_mat_free(&(descVT)) ); } - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); subA = NULL; subT = NULL; @@ -606,7 +609,7 @@ int MORSE_zgesvd_Tile_Async(MORSE_enum jobu, MORSE_enum jobvt, //morse_zooptile2lap( descVT, VT, NB, NB, LDVT, N, sequence, request ); } - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); if (subA) { free(subA); free(subUVT); free(subT); diff --git a/compute/zgetrf_incpiv.c b/compute/zgetrf_incpiv.c index f8b17559da71b7a9f391415b81c7836c81b3f2fa..a86915d081d5b92b07269f534348dc985c023f7d 100644 --- a/compute/zgetrf_incpiv.c +++ b/compute/zgetrf_incpiv.c @@ -78,9 +78,9 @@ * @sa MORSE_zgetrs_incpiv * ******************************************************************************/ -int MORSE_zgetrf_incpiv(int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descL, int *IPIV) +int MORSE_zgetrf_incpiv( int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descL, int *IPIV ) { int NB; int status; @@ -121,7 +121,7 @@ int MORSE_zgetrf_incpiv(int M, int N, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -134,13 +134,13 @@ int MORSE_zgetrf_incpiv(int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -198,14 +198,16 @@ int MORSE_zgetrf_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV ) morse_fatal_error("MORSE_zgetrf_incpiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgetrf_incpiv_Tile_Async( A, L, IPIV, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( L, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -238,8 +240,8 @@ int MORSE_zgetrf_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV ) * @sa MORSE_zgetrs_incpiv_Tile_Async * ******************************************************************************/ -int MORSE_zgetrf_incpiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgetrf_incpiv_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zgetrf_nopiv.c b/compute/zgetrf_nopiv.c index 689ca3d2e39c1c4db102efef2df90f0ed330acd2..cf44857a8df9ac437e6b9bd4ef6b5c84de571d94 100644 --- a/compute/zgetrf_nopiv.c +++ b/compute/zgetrf_nopiv.c @@ -72,8 +72,8 @@ * @sa MORSE_sgetrf_nopiv * ******************************************************************************/ -int MORSE_zgetrf_nopiv(int M, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_zgetrf_nopiv( int M, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -114,7 +114,7 @@ int MORSE_zgetrf_nopiv(int M, int N, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, @@ -127,13 +127,13 @@ int MORSE_zgetrf_nopiv(int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -184,14 +184,15 @@ int MORSE_zgetrf_nopiv_Tile( MORSE_desc_t *A ) morse_fatal_error("MORSE_zgetrf_nopiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgetrf_nopiv_Tile_Async( A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -228,9 +229,9 @@ int MORSE_zgetrf_nopiv_Tile( MORSE_desc_t *A ) * @sa MORSE_zgetrs_Tile_Async * ******************************************************************************/ -int MORSE_zgetrf_nopiv_Tile_Async(MORSE_desc_t *A, - MORSE_sequence_t *sequence, - MORSE_request_t *request) +int MORSE_zgetrf_nopiv_Tile_Async( MORSE_desc_t *A, + MORSE_sequence_t *sequence, + MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zgetrs_incpiv.c b/compute/zgetrs_incpiv.c index d8a786e068077b903cec99f0d93097dd31b40a2b..f877c8aa38b8d1ab39c86809842803384e53ba85 100644 --- a/compute/zgetrs_incpiv.c +++ b/compute/zgetrs_incpiv.c @@ -88,10 +88,10 @@ * @sa MORSE_zgetrf_incpiv * ******************************************************************************/ -int MORSE_zgetrs_incpiv(MORSE_enum trans, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descL, int *IPIV, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgetrs_incpiv( MORSE_enum trans, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descL, int *IPIV, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -141,7 +141,7 @@ int MORSE_zgetrs_incpiv(MORSE_enum trans, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -158,14 +158,14 @@ int MORSE_zgetrs_incpiv(MORSE_enum trans, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -223,15 +223,17 @@ int MORSE_zgetrs_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE morse_fatal_error("MORSE_zgetrs_incpiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgetrs_incpiv_Tile_Async( A, L, IPIV, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( L, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -265,8 +267,8 @@ int MORSE_zgetrs_incpiv_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE * @sa MORSE_zgetrf_incpiv_Tile_Async * ******************************************************************************/ -int MORSE_zgetrs_incpiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgetrs_incpiv_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zgetrs_nopiv.c b/compute/zgetrs_nopiv.c index e1f5d857b22dbb414c62ccfd2c7a7efb82d670c6..7833e3eafad5da98d755e1af57c94b89dd3c94d0 100644 --- a/compute/zgetrs_nopiv.c +++ b/compute/zgetrs_nopiv.c @@ -83,9 +83,9 @@ * @sa MORSE_zgetrf_nopiv * ******************************************************************************/ -int MORSE_zgetrs_nopiv(MORSE_enum trans, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zgetrs_nopiv( MORSE_enum trans, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -135,7 +135,7 @@ int MORSE_zgetrs_nopiv(MORSE_enum trans, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -152,14 +152,14 @@ int MORSE_zgetrs_nopiv(MORSE_enum trans, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -211,15 +211,16 @@ int MORSE_zgetrs_nopiv_Tile( MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zgetrs_nopiv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zgetrs_nopiv_Tile_Async( A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -253,8 +254,8 @@ int MORSE_zgetrs_nopiv_Tile( MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_zgetrf_nopiv_Tile_Async * ******************************************************************************/ -int MORSE_zgetrs_nopiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zgetrs_nopiv_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zheevd.c b/compute/zheevd.c index 200d35460d231a979644756310d702d6157a8a78..7346739297368b2bdc181beeaf0a99385fe418ee 100644 --- a/compute/zheevd.c +++ b/compute/zheevd.c @@ -100,10 +100,10 @@ * @sa MORSE_ssyev * ******************************************************************************/ -int MORSE_zheevd(MORSE_enum jobz, MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA, - double *W, - MORSE_desc_t *descT) +int MORSE_zheevd( MORSE_enum jobz, MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA, + double *W, + MORSE_desc_t *descT ) { int NB; int status; @@ -150,7 +150,7 @@ int MORSE_zheevd(MORSE_enum jobz, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -163,13 +163,13 @@ int MORSE_zheevd(MORSE_enum jobz, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -240,12 +240,12 @@ int MORSE_zheevd(MORSE_enum jobz, MORSE_enum uplo, int N, * @sa MORSE_ssyev * ******************************************************************************/ -int MORSE_zheevd_Tile(MORSE_enum jobz, MORSE_enum uplo, - MORSE_desc_t *A, double *W, MORSE_desc_t *T) +int MORSE_zheevd_Tile( MORSE_enum jobz, MORSE_enum uplo, + MORSE_desc_t *A, double *W, MORSE_desc_t *T ) { - MORSE_context_t *morse; + MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; - MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + MORSE_request_t request = MORSE_REQUEST_INITIALIZER; int status; morse = morse_context_self(); @@ -253,15 +253,16 @@ int MORSE_zheevd_Tile(MORSE_enum jobz, MORSE_enum uplo, morse_fatal_error("MORSE_zheevd_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zheevd_Tile_Async( jobz, uplo, A, W, T, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( T, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -481,7 +482,7 @@ int MORSE_zheevd_Tile_Async( MORSE_enum jobz, MORSE_enum uplo, } #endif subA = morse_desc_submatrix(&descA, descA.mb, 0, descA.m -descA.mb, descA.n-descA.nb); - subQ = morse_desc_submatrix(&descQ2, descQ2.mb, 0, descQ2.m-descQ2.mb, descQ2.n ); + subQ = morse_desc_submatrix(&descQ2, descQ2.mb, 0, descQ2.m-descQ2.mb, descQ2.n ); subT = morse_desc_submatrix(&descT, descT.mb, 0, descT.m -descT.mb, descT.n-descT.nb); /* Compute Q2 = Q1 * Q2 */ @@ -505,7 +506,7 @@ int MORSE_zheevd_Tile_Async( MORSE_enum jobz, MORSE_enum uplo, } #endif subA = morse_desc_submatrix(&descA, 0, descA.nb, descA.m -descA.mb, descA.n -descA.nb ); - subQ = morse_desc_submatrix(&descQ2, descQ2.mb, 0, descQ2.m-descQ2.mb, descQ2.n ); + subQ = morse_desc_submatrix(&descQ2, descQ2.mb, 0, descQ2.m-descQ2.mb, descQ2.n ); subT = morse_desc_submatrix(&descT, 0, descT.nb, descT.m -descT.mb, descT.n -descT.nb ); /* Compute Q2 = Q1^h * Q2 */ @@ -520,7 +521,7 @@ int MORSE_zheevd_Tile_Async( MORSE_enum jobz, MORSE_enum uplo, sequence, request ); } - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); free(subA); free(subQ); free(subT); morse_desc_mat_free( &descQ2 ); diff --git a/compute/zhemm.c b/compute/zhemm.c index 37ca9b64483058de28a11ac6048f4d3ac6a97c2f..ce06a71dab415bba74aa660e34c21b270d79bfab 100644 --- a/compute/zhemm.c +++ b/compute/zhemm.c @@ -106,10 +106,10 @@ * @sa MORSE_shemm * ******************************************************************************/ -int MORSE_zhemm(MORSE_enum side, MORSE_enum uplo, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB, - MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC) +int MORSE_zhemm( MORSE_enum side, MORSE_enum uplo, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB, + MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am; @@ -173,7 +173,7 @@ int MORSE_zhemm(MORSE_enum side, MORSE_enum uplo, int M, int N, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -194,7 +194,7 @@ int MORSE_zhemm(MORSE_enum side, MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -202,7 +202,7 @@ int MORSE_zhemm(MORSE_enum side, MORSE_enum uplo, int M, int N, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -265,9 +265,9 @@ int MORSE_zhemm(MORSE_enum side, MORSE_enum uplo, int M, int N, * @sa MORSE_shemm_Tile * ******************************************************************************/ -int MORSE_zhemm_Tile(MORSE_enum side, MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C) +int MORSE_zhemm_Tile( MORSE_enum side, MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -279,16 +279,17 @@ int MORSE_zhemm_Tile(MORSE_enum side, MORSE_enum uplo, morse_fatal_error("MORSE_zhemm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zhemm_Tile_Async(side, uplo, alpha, A, B, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -320,10 +321,10 @@ int MORSE_zhemm_Tile(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_shemm_Tile_Async * ******************************************************************************/ -int MORSE_zhemm_Tile_Async(MORSE_enum side, MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zhemm_Tile_Async( MORSE_enum side, MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; @@ -414,9 +415,11 @@ int MORSE_zhemm_Tile_Async(MORSE_enum side, MORSE_enum uplo, } /* Quick return */ - if (C->m == 0 || C->n == 0 || - ( (alpha == (MORSE_Complex64_t)0.0) && (beta == (MORSE_Complex64_t)1.0) )) + if ( (C->m == 0) || (C->n == 0) || + ( (alpha == (MORSE_Complex64_t)0.0) && (beta == (MORSE_Complex64_t)1.0) ) ) + { return MORSE_SUCCESS; + } morse_pzhemm( side, uplo, alpha, A, B, beta, C, sequence, request ); diff --git a/compute/zher2k.c b/compute/zher2k.c index cd021b2fa47571d7649202b63872712bae384bac..71608495f399ad90c52a5b4dd4c00a9366a5f3b2 100644 --- a/compute/zher2k.c +++ b/compute/zher2k.c @@ -108,9 +108,9 @@ * @sa MORSE_sher2k * ******************************************************************************/ -int MORSE_zher2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, +int MORSE_zher2k( MORSE_enum uplo, MORSE_enum trans, int N, int K, MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, MORSE_Complex64_t *B, int LDB, - double beta, MORSE_Complex64_t *C, int LDC) + double beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am, An; @@ -178,7 +178,7 @@ int MORSE_zher2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -199,7 +199,7 @@ int MORSE_zher2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -207,7 +207,7 @@ int MORSE_zher2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -265,9 +265,9 @@ int MORSE_zher2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, * @sa MORSE_sher2k * ******************************************************************************/ -int MORSE_zher2k_Tile(MORSE_enum uplo, MORSE_enum trans, +int MORSE_zher2k_Tile( MORSE_enum uplo, MORSE_enum trans, MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - double beta, MORSE_desc_t *C) + double beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -279,16 +279,17 @@ int MORSE_zher2k_Tile(MORSE_enum uplo, MORSE_enum trans, morse_fatal_error("MORSE_zher2k_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zher2k_Tile_Async( uplo, trans, alpha, A, B, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -320,10 +321,10 @@ int MORSE_zher2k_Tile(MORSE_enum uplo, MORSE_enum trans, * @sa MORSE_sher2k_Tile_Async * ******************************************************************************/ -int MORSE_zher2k_Tile_Async(MORSE_enum uplo, MORSE_enum trans, +int MORSE_zher2k_Tile_Async( MORSE_enum uplo, MORSE_enum trans, MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, double beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int N, K; @@ -404,9 +405,11 @@ int MORSE_zher2k_Tile_Async(MORSE_enum uplo, MORSE_enum trans, K = An; /* Quick return */ - if ( N == 0 || - ((alpha == (MORSE_Complex64_t)0.0 || K == 0) && beta == (double)1.0)) + if ( (N == 0) || + (((alpha == (MORSE_Complex64_t)0.0) || (K == 0)) && (beta == (double)1.0)) ) + { return MORSE_SUCCESS; + } morse_pzher2k( uplo, trans, alpha, A, B, beta, C, sequence, request ); diff --git a/compute/zherk.c b/compute/zherk.c index abb9a359499f06d55e2fe9dc8dcb8e25e70be98a..b99dd19884401cb400d435abc60dea745f9e333f 100644 --- a/compute/zherk.c +++ b/compute/zherk.c @@ -98,9 +98,9 @@ * @sa MORSE_sherk * ******************************************************************************/ -int MORSE_zherk(MORSE_enum uplo, MORSE_enum trans, int N, int K, - double alpha, MORSE_Complex64_t *A, int LDA, - double beta, MORSE_Complex64_t *C, int LDC) +int MORSE_zherk( MORSE_enum uplo, MORSE_enum trans, int N, int K, + double alpha, MORSE_Complex64_t *A, int LDA, + double beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am, An; @@ -163,7 +163,7 @@ int MORSE_zherk(MORSE_enum uplo, MORSE_enum trans, int N, int K, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -180,14 +180,14 @@ int MORSE_zherk(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -241,9 +241,9 @@ int MORSE_zherk(MORSE_enum uplo, MORSE_enum trans, int N, int K, * @sa MORSE_sherk * ******************************************************************************/ -int MORSE_zherk_Tile(MORSE_enum uplo, MORSE_enum trans, - double alpha, MORSE_desc_t *A, - double beta, MORSE_desc_t *C) +int MORSE_zherk_Tile( MORSE_enum uplo, MORSE_enum trans, + double alpha, MORSE_desc_t *A, + double beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -255,15 +255,16 @@ int MORSE_zherk_Tile(MORSE_enum uplo, MORSE_enum trans, morse_fatal_error("MORSE_zherk_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zherk_Tile_Async( uplo, trans, alpha, A, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -295,10 +296,10 @@ int MORSE_zherk_Tile(MORSE_enum uplo, MORSE_enum trans, * @sa MORSE_sherk_Tile_Async * ******************************************************************************/ -int MORSE_zherk_Tile_Async(MORSE_enum uplo, MORSE_enum trans, - double alpha, MORSE_desc_t *A, - double beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zherk_Tile_Async( MORSE_enum uplo, MORSE_enum trans, + double alpha, MORSE_desc_t *A, + double beta, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int N, K; @@ -375,9 +376,11 @@ int MORSE_zherk_Tile_Async(MORSE_enum uplo, MORSE_enum trans, K = An; /* Quick return */ - if ( N == 0 || - ((alpha == (double)0.0 || K == 0) && beta == (double)1.0)) + if ( (N == 0) || + (((alpha == (double)0.0) || (K == 0)) && (beta == (double)1.0)) ) + { return MORSE_SUCCESS; + } morse_pzherk( uplo, trans, alpha, A, beta, C, sequence, request ); diff --git a/compute/zhetrd.c b/compute/zhetrd.c index 8d74ec7735878aacba7bcded227ec03a4c044bcf..5593d45a4b18c68d3403bdf7bee88be86f2449eb 100644 --- a/compute/zhetrd.c +++ b/compute/zhetrd.c @@ -113,12 +113,12 @@ * @sa MORSE_ssytrd * ******************************************************************************/ -int MORSE_zhetrd(MORSE_enum jobz, MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA, - double *D, - double *E, - MORSE_desc_t *descT, - MORSE_Complex64_t *Q, int LDQ) +int MORSE_zhetrd( MORSE_enum jobz, MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA, + double *D, + double *E, + MORSE_desc_t *descT, + MORSE_Complex64_t *Q, int LDQ ) { int NB; int status; @@ -164,7 +164,7 @@ int MORSE_zhetrd(MORSE_enum jobz, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -177,13 +177,13 @@ int MORSE_zhetrd(MORSE_enum jobz, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -267,9 +267,9 @@ int MORSE_zhetrd(MORSE_enum jobz, MORSE_enum uplo, int N, * @sa MORSE_zhetrd_Tile * ******************************************************************************/ -int MORSE_zhetrd_Tile(MORSE_enum jobz, MORSE_enum uplo, - MORSE_desc_t *A, double *D, double *E, - MORSE_desc_t *T, MORSE_Complex64_t *Q, int LDQ) +int MORSE_zhetrd_Tile( MORSE_enum jobz, MORSE_enum uplo, + MORSE_desc_t *A, double *D, double *E, + MORSE_desc_t *T, MORSE_Complex64_t *Q, int LDQ ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -281,13 +281,16 @@ int MORSE_zhetrd_Tile(MORSE_enum jobz, MORSE_enum uplo, morse_fatal_error("MORSE_zhetrd_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zhetrd_Tile_Async( jobz, uplo, A, D, E, T, Q, LDQ, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -323,14 +326,14 @@ int MORSE_zhetrd_Tile(MORSE_enum jobz, MORSE_enum uplo, * @sa MORSE_ssytrd_Tile_Async * ******************************************************************************/ -int MORSE_zhetrd_Tile_Async(MORSE_enum jobz, - MORSE_enum uplo, - MORSE_desc_t *A, - double *W, - double *E, - MORSE_desc_t *T, - MORSE_Complex64_t *Q, int LDQ, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zhetrd_Tile_Async( MORSE_enum jobz, + MORSE_enum uplo, + MORSE_desc_t *A, + double *W, + double *E, + MORSE_desc_t *T, + MORSE_Complex64_t *Q, int LDQ, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t descA; @@ -417,7 +420,7 @@ int MORSE_zhetrd_Tile_Async(MORSE_enum jobz, morse_pztile2band( uplo, A, &descAB, sequence, request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Reduce band matrix to tridiagonal matrix */ #if !defined(CHAMELEON_SIMULATION) diff --git a/compute/zlacpy.c b/compute/zlacpy.c index fbcc00388442ce54a7b1a2281f2b9f399c11a240..d47417ac4d7bdc868db32dae4d94bdac3fad0ae3 100644 --- a/compute/zlacpy.c +++ b/compute/zlacpy.c @@ -74,9 +74,9 @@ * @sa MORSE_slacpy * ******************************************************************************/ -int MORSE_zlacpy(MORSE_enum uplo, int M, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zlacpy( MORSE_enum uplo, int M, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -129,7 +129,7 @@ int MORSE_zlacpy(MORSE_enum uplo, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -146,13 +146,13 @@ int MORSE_zlacpy(MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return MORSE_SUCCESS; } @@ -202,19 +202,24 @@ int MORSE_zlacpy_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; morse = morse_context_self(); if (morse == NULL) { morse_fatal_error("MORSE_zlacpy_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlacpy_Tile_Async( uplo, A, B, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( B, sequence ); - morse_sequence_destroy(morse, sequence); - return MORSE_SUCCESS; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return status; } /** @@ -244,8 +249,8 @@ int MORSE_zlacpy_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_slacpy_Tile_Async * ******************************************************************************/ -int MORSE_zlacpy_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlacpy_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlange.c b/compute/zlange.c index a523fad21dd5b2e9283217435a4c7014b749284d..e351bb5c7672c141f0e357cb53213f06c39fc6dc 100644 --- a/compute/zlange.c +++ b/compute/zlange.c @@ -84,7 +84,7 @@ * ******************************************************************************/ double MORSE_zlange(MORSE_enum norm, int M, int N, - MORSE_Complex64_t *A, int LDA) + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -132,7 +132,7 @@ double MORSE_zlange(MORSE_enum norm, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -145,12 +145,12 @@ double MORSE_zlange(MORSE_enum norm, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return value; } @@ -193,11 +193,12 @@ double MORSE_zlange(MORSE_enum norm, int M, int N, * @sa MORSE_slange_Tile * ******************************************************************************/ -double MORSE_zlange_Tile(MORSE_enum norm, MORSE_desc_t *A) +double MORSE_zlange_Tile(MORSE_enum norm, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; double value; morse = morse_context_self(); @@ -205,13 +206,16 @@ double MORSE_zlange_Tile(MORSE_enum norm, MORSE_desc_t *A) morse_fatal_error("MORSE_zlange_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlange_Tile_Async( norm, A, &value, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); - morse_sequence_destroy(morse, sequence); - return value; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return ( status == MORSE_SUCCESS ) ? value : (double)status; } /** @@ -241,8 +245,8 @@ double MORSE_zlange_Tile(MORSE_enum norm, MORSE_desc_t *A) * @sa MORSE_slange_Tile_Async * ******************************************************************************/ -int MORSE_zlange_Tile_Async(MORSE_enum norm, MORSE_desc_t *A, double *value, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlange_Tile_Async( MORSE_enum norm, MORSE_desc_t *A, double *value, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlanhe.c b/compute/zlanhe.c index 9bfa983315dd0ad8ba391f064aa295cbe7246573..e43a8111d29f62922cd24050badb702d29875b32 100644 --- a/compute/zlanhe.c +++ b/compute/zlanhe.c @@ -84,7 +84,7 @@ * ******************************************************************************/ double MORSE_zlanhe(MORSE_enum norm, MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -132,7 +132,7 @@ double MORSE_zlanhe(MORSE_enum norm, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -145,12 +145,12 @@ double MORSE_zlanhe(MORSE_enum norm, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return value; } @@ -197,11 +197,12 @@ double MORSE_zlanhe(MORSE_enum norm, MORSE_enum uplo, int N, * @sa MORSE_slanhe_Tile * ******************************************************************************/ -double MORSE_zlanhe_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) +double MORSE_zlanhe_Tile( MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; double value; morse = morse_context_self(); @@ -209,13 +210,16 @@ double MORSE_zlanhe_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) morse_fatal_error("MORSE_zlanhe_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlanhe_Tile_Async( norm, uplo, A, &value, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); - morse_sequence_destroy(morse, sequence); - return value; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return ( status == MORSE_SUCCESS ) ? value : (double)status; } /** @@ -245,8 +249,8 @@ double MORSE_zlanhe_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) * @sa MORSE_slanhe_Tile_Async * ******************************************************************************/ -int MORSE_zlanhe_Tile_Async(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A, double *value, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlanhe_Tile_Async( MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A, double *value, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlansy.c b/compute/zlansy.c index 1cff553583b51b71f73ed243894555cdf2f1a6dd..cb97a551881bb1cdb44580e6f3ca7d6631500964 100644 --- a/compute/zlansy.c +++ b/compute/zlansy.c @@ -84,7 +84,7 @@ * ******************************************************************************/ double MORSE_zlansy(MORSE_enum norm, MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -132,7 +132,7 @@ double MORSE_zlansy(MORSE_enum norm, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -145,12 +145,12 @@ double MORSE_zlansy(MORSE_enum norm, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return value; } @@ -197,11 +197,12 @@ double MORSE_zlansy(MORSE_enum norm, MORSE_enum uplo, int N, * @sa MORSE_slansy_Tile * ******************************************************************************/ -double MORSE_zlansy_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) +double MORSE_zlansy_Tile( MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; double value; morse = morse_context_self(); @@ -209,13 +210,16 @@ double MORSE_zlansy_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) morse_fatal_error("MORSE_zlansy_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlansy_Tile_Async( norm, uplo, A, &value, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); - morse_sequence_destroy(morse, sequence); - return value; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return ( status == MORSE_SUCCESS ) ? value : (double)status; } /** @@ -245,8 +249,8 @@ double MORSE_zlansy_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A) * @sa MORSE_slansy_Tile_Async * ******************************************************************************/ -int MORSE_zlansy_Tile_Async(MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A, double *value, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlansy_Tile_Async( MORSE_enum norm, MORSE_enum uplo, MORSE_desc_t *A, double *value, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlantr.c b/compute/zlantr.c index 079e2f92d7de495c43baec2ee0015aa4626c2377..f803abd84474cbe35a720c2bd5e98d20bd727a5d 100644 --- a/compute/zlantr.c +++ b/compute/zlantr.c @@ -96,7 +96,7 @@ * ******************************************************************************/ double MORSE_zlantr(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, - int M, int N, MORSE_Complex64_t *A, int LDA) + int M, int N, MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -152,7 +152,7 @@ double MORSE_zlantr(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -165,12 +165,12 @@ double MORSE_zlantr(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return value; } @@ -219,11 +219,12 @@ double MORSE_zlantr(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, * @sa MORSE_slantr_Tile * ******************************************************************************/ -double MORSE_zlantr_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A) +double MORSE_zlantr_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; double value; morse = morse_context_self(); @@ -231,13 +232,16 @@ double MORSE_zlantr_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, MORS morse_fatal_error("MORSE_zlantr_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlantr_Tile_Async( norm, uplo, diag, A, &value, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); - morse_sequence_destroy(morse, sequence); - return value; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return ( status == MORSE_SUCCESS ) ? value : (double)status; } /** @@ -267,9 +271,9 @@ double MORSE_zlantr_Tile(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, MORS * @sa MORSE_slantr_Tile_Async * ******************************************************************************/ -int MORSE_zlantr_Tile_Async(MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, - MORSE_desc_t *A, double *value, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlantr_Tile_Async( MORSE_enum norm, MORSE_enum uplo, MORSE_enum diag, + MORSE_desc_t *A, double *value, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlascal.c b/compute/zlascal.c index 68b84cdaa4e04a38a29324861977cc85433c0147..b37d8f61881884d7b7dcd943f7c7b5b44c55f0be 100644 --- a/compute/zlascal.c +++ b/compute/zlascal.c @@ -62,8 +62,8 @@ * @sa MORSE_slascal * ******************************************************************************/ -int MORSE_zlascal(MORSE_enum uplo, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA) +int MORSE_zlascal( MORSE_enum uplo, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -111,7 +111,7 @@ int MORSE_zlascal(MORSE_enum uplo, int M, int N, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -124,13 +124,13 @@ int MORSE_zlascal(MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -174,8 +174,8 @@ int MORSE_zlascal(MORSE_enum uplo, int M, int N, * @sa MORSE_slascal_Tile * ******************************************************************************/ -int MORSE_zlascal_Tile(MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A) +int MORSE_zlascal_Tile( MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -187,11 +187,15 @@ int MORSE_zlascal_Tile(MORSE_enum uplo, morse_fatal_error("MORSE_zlascal_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlascal_Tile_Async( uplo, alpha, A, sequence, &request ); - morse_sequence_wait(morse, sequence); + + MORSE_Desc_Flush( A, sequence ); + + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -224,9 +228,9 @@ int MORSE_zlascal_Tile(MORSE_enum uplo, * @sa MORSE_slascal_Tile_Async * ******************************************************************************/ -int MORSE_zlascal_Tile_Async(MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlascal_Tile_Async( MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t descA; @@ -273,7 +277,9 @@ int MORSE_zlascal_Tile_Async(MORSE_enum uplo, /* Quick return */ if ( (descA.m == 0) || (descA.n == 0) || (alpha == (MORSE_Complex64_t)1.0) ) + { return MORSE_SUCCESS; + } morse_pzlascal( uplo, alpha, A, sequence, request ); diff --git a/compute/zlaset.c b/compute/zlaset.c index 9b0d5572071991363916daa62671b47e5e6f0d96..41093bc7d0468aa7cf7e5b0c7df96d7da21b69e8 100644 --- a/compute/zlaset.c +++ b/compute/zlaset.c @@ -75,9 +75,9 @@ * @sa MORSE_slaset * ******************************************************************************/ -int MORSE_zlaset(MORSE_enum uplo, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t beta, - MORSE_Complex64_t *A, int LDA) +int MORSE_zlaset( MORSE_enum uplo, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t beta, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -125,7 +125,7 @@ int MORSE_zlaset(MORSE_enum uplo, int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -138,12 +138,12 @@ int MORSE_zlaset(MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return MORSE_SUCCESS; } @@ -184,26 +184,30 @@ int MORSE_zlaset(MORSE_enum uplo, int M, int N, * @sa MORSE_slaset_Tile * ******************************************************************************/ -int MORSE_zlaset_Tile(MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_Complex64_t beta, - MORSE_desc_t *A) +int MORSE_zlaset_Tile( MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_Complex64_t beta, + MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; MORSE_request_t request = MORSE_REQUEST_INITIALIZER; + int status; morse = morse_context_self(); if (morse == NULL) { morse_fatal_error("MORSE_zlaset_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlaset_Tile_Async( uplo, alpha, beta, A, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); - morse_sequence_destroy(morse, sequence); - return MORSE_SUCCESS; + morse_sequence_wait( morse, sequence ); + status = sequence->status; + morse_sequence_destroy( morse, sequence ); + return status; } /** @@ -233,10 +237,10 @@ int MORSE_zlaset_Tile(MORSE_enum uplo, * @sa MORSE_slaset_Tile_Async * ******************************************************************************/ -int MORSE_zlaset_Tile_Async(MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_Complex64_t beta, - MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlaset_Tile_Async( MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_Complex64_t beta, + MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zlauum.c b/compute/zlauum.c index 963ad6e18089c378be28a4ab56858fd7bb02a19b..7c29cf3908488ddd67fa090718b3daaab92556dd 100644 --- a/compute/zlauum.c +++ b/compute/zlauum.c @@ -77,8 +77,8 @@ * @sa MORSE_zpotri * ******************************************************************************/ -int MORSE_zlauum(MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_zlauum( MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -119,7 +119,7 @@ int MORSE_zlauum(MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -132,13 +132,13 @@ int MORSE_zlauum(MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -195,14 +195,15 @@ int MORSE_zlauum_Tile( MORSE_enum uplo, MORSE_desc_t *A ) morse_fatal_error("MORSE_zlauum_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zlauum_Tile_Async( uplo, A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -237,8 +238,8 @@ int MORSE_zlauum_Tile( MORSE_enum uplo, MORSE_desc_t *A ) * @sa MORSE_zpotri_Tile_Async * ******************************************************************************/ -int MORSE_zlauum_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zlauum_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zplghe.c b/compute/zplghe.c index 5836630efb52ec812e62b8fbd22b726ff1fb6c03..5f8deeee1f8174ba6eb659b97d140243beefb537 100644 --- a/compute/zplghe.c +++ b/compute/zplghe.c @@ -113,7 +113,7 @@ int MORSE_zplghe( double bump, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescOutput, uplo, @@ -126,13 +126,13 @@ int MORSE_zplghe( double bump, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescOutput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -191,13 +191,15 @@ int MORSE_zplghe_Tile( double bump, MORSE_enum uplo, MORSE_desc_t *A, morse_fatal_error("MORSE_zplghe_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zplghe_Tile_Async( bump, uplo, A, seed, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/zplgsy.c b/compute/zplgsy.c index 832a5a785a93f196f8c5c1543ee9140674c68a69..31ec6fdc714d98fd2d8e0ca6a5fd56c2c13be896 100644 --- a/compute/zplgsy.c +++ b/compute/zplgsy.c @@ -113,7 +113,7 @@ int MORSE_zplgsy( MORSE_Complex64_t bump, MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescOutput, uplo, @@ -126,13 +126,13 @@ int MORSE_zplgsy( MORSE_Complex64_t bump, MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescOutput, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -192,13 +192,15 @@ int MORSE_zplgsy_Tile( MORSE_Complex64_t bump, MORSE_enum uplo, morse_fatal_error("MORSE_zplgsy_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zplgsy_Tile_Async( bump, uplo, A, seed, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/zplrnt.c b/compute/zplrnt.c index 9494b1a3c334de9a58f32b6b24dbd7990448ff2d..8035b30ea27c4e6386cff37c1e4ef878e237a221 100644 --- a/compute/zplrnt.c +++ b/compute/zplrnt.c @@ -111,7 +111,7 @@ int MORSE_zplrnt( int M, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescOutput, MorseUpperLower, @@ -124,13 +124,13 @@ int MORSE_zplrnt( int M, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescOutput, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -183,13 +183,15 @@ int MORSE_zplrnt_Tile( MORSE_desc_t *A, morse_fatal_error("MORSE_zplrnt_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zplrnt_Tile_Async( A, seed, sequence, &request ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -226,7 +228,7 @@ int MORSE_zplrnt_Tile( MORSE_desc_t *A, int MORSE_zplrnt_Tile_Async( MORSE_desc_t *A, unsigned long long int seed, MORSE_sequence_t *sequence, - MORSE_request_t *request) + MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zposv.c b/compute/zposv.c index 69797306995f40ed74afbe0e800082c959513ff3..6f4615ce69fbe0b6ae145d5bbe7ef5ba3ec94941 100644 --- a/compute/zposv.c +++ b/compute/zposv.c @@ -95,9 +95,9 @@ * @sa MORSE_sposv * ******************************************************************************/ -int MORSE_zposv(MORSE_enum uplo, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zposv( MORSE_enum uplo, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -148,7 +148,7 @@ int MORSE_zposv(MORSE_enum uplo, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -165,14 +165,14 @@ int MORSE_zposv(MORSE_enum uplo, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -237,15 +237,16 @@ int MORSE_zposv_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zposv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zposv_Tile_Async( uplo, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -278,8 +279,8 @@ int MORSE_zposv_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_sposv_Tile_Async * ******************************************************************************/ -int MORSE_zposv_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zposv_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zpotrf.c b/compute/zpotrf.c index 4a718629ef2e682c8f5ef8d185c262bd751bca04..f37ee389b55fdc2ab7d11c0c02560f4b91b9eda4 100644 --- a/compute/zpotrf.c +++ b/compute/zpotrf.c @@ -83,8 +83,8 @@ * @sa MORSE_zpotrs * ******************************************************************************/ -int MORSE_zpotrf(MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_zpotrf( MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -125,7 +125,7 @@ int MORSE_zpotrf(MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -138,13 +138,13 @@ int MORSE_zpotrf(MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -207,14 +207,15 @@ int MORSE_zpotrf_Tile( MORSE_enum uplo, MORSE_desc_t *A ) morse_fatal_error("MORSE_zpotrf_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zpotrf_Tile_Async( uplo, A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -248,8 +249,8 @@ int MORSE_zpotrf_Tile( MORSE_enum uplo, MORSE_desc_t *A ) * @sa MORSE_zpotrs_Tile_Async * ******************************************************************************/ -int MORSE_zpotrf_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zpotrf_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zpotri.c b/compute/zpotri.c index 7378823a484ad738c9328233726a73fc8d4f3ad4..6662a1c661608b63269d537b8def7ed4c54c4f77 100644 --- a/compute/zpotri.c +++ b/compute/zpotri.c @@ -74,8 +74,8 @@ * @sa MORSE_zpotrf * ******************************************************************************/ -int MORSE_zpotri(MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_zpotri( MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -116,7 +116,7 @@ int MORSE_zpotri(MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -129,13 +129,13 @@ int MORSE_zpotri(MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -195,14 +195,15 @@ int MORSE_zpotri_Tile( MORSE_enum uplo, MORSE_desc_t *A ) morse_fatal_error("MORSE_zpotri_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zpotri_Tile_Async( uplo, A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -237,8 +238,8 @@ int MORSE_zpotri_Tile( MORSE_enum uplo, MORSE_desc_t *A ) * @sa MORSE_zpotrf_Tile_Async * ******************************************************************************/ -int MORSE_zpotri_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zpotri_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zpotrimm.c b/compute/zpotrimm.c index 2a97aa4a3c409448e3ceb6f9b4e521c0e530604a..ec2f8759a6a6dd690545751c165cbfd71044c36b 100644 --- a/compute/zpotrimm.c +++ b/compute/zpotrimm.c @@ -74,10 +74,10 @@ * @sa MORSE_zpotrf * ******************************************************************************/ -int MORSE_zpotrimm(MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB, - MORSE_Complex64_t *C, int LDC) +int MORSE_zpotrimm( MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB, + MORSE_Complex64_t *C, int LDC ) { int NB; int status; @@ -128,7 +128,7 @@ int MORSE_zpotrimm(MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -149,7 +149,7 @@ int MORSE_zpotrimm(MORSE_enum uplo, int N, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -157,7 +157,7 @@ int MORSE_zpotrimm(MORSE_enum uplo, int N, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -217,14 +217,17 @@ int MORSE_zpotrimm_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORS morse_fatal_error("MORSE_zpotrimm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zpotrimm_Tile_Async( uplo, A, B, C, sequence, &request ); - MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( B, sequence ); + MORSE_Desc_Flush( C, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -259,8 +262,8 @@ int MORSE_zpotrimm_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORS * @sa MORSE_zpotrf_Tile_Async * ******************************************************************************/ -int MORSE_zpotrimm_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zpotrimm_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zpotrs.c b/compute/zpotrs.c index a78834af34e9dcffe0bdedb7db3d28ad49a5face..7ff43c5698b9f6eda87fe6342edf3abdca8ed95d 100644 --- a/compute/zpotrs.c +++ b/compute/zpotrs.c @@ -80,9 +80,9 @@ * @sa MORSE_zpotrf * ******************************************************************************/ -int MORSE_zpotrs(MORSE_enum uplo, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zpotrs( MORSE_enum uplo, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -132,7 +132,7 @@ int MORSE_zpotrs(MORSE_enum uplo, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -149,14 +149,14 @@ int MORSE_zpotrs(MORSE_enum uplo, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -213,15 +213,16 @@ int MORSE_zpotrs_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zpotrs_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zpotrs_Tile_Async( uplo, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -255,8 +256,8 @@ int MORSE_zpotrs_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_zpotrf_Tile_Async * ******************************************************************************/ -int MORSE_zpotrs_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zpotrs_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zsymm.c b/compute/zsymm.c index 45603d2b366fcc1179614e1904910f14be93b111..a6662f3e9adeb3b6c2468e4ee859f325c94395af 100644 --- a/compute/zsymm.c +++ b/compute/zsymm.c @@ -106,10 +106,10 @@ * @sa MORSE_ssymm * ******************************************************************************/ -int MORSE_zsymm(MORSE_enum side, MORSE_enum uplo, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB, - MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC) +int MORSE_zsymm( MORSE_enum side, MORSE_enum uplo, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB, + MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am; @@ -173,7 +173,7 @@ int MORSE_zsymm(MORSE_enum side, MORSE_enum uplo, int M, int N, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -194,7 +194,7 @@ int MORSE_zsymm(MORSE_enum side, MORSE_enum uplo, int M, int N, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -202,7 +202,7 @@ int MORSE_zsymm(MORSE_enum side, MORSE_enum uplo, int M, int N, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -265,9 +265,9 @@ int MORSE_zsymm(MORSE_enum side, MORSE_enum uplo, int M, int N, * @sa MORSE_ssymm_Tile * ******************************************************************************/ -int MORSE_zsymm_Tile(MORSE_enum side, MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C) +int MORSE_zsymm_Tile( MORSE_enum side, MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -279,16 +279,17 @@ int MORSE_zsymm_Tile(MORSE_enum side, MORSE_enum uplo, morse_fatal_error("MORSE_zsymm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsymm_Tile_Async(side, uplo, alpha, A, B, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -320,10 +321,10 @@ int MORSE_zsymm_Tile(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_ssymm_Tile_Async * ******************************************************************************/ -int MORSE_zsymm_Tile_Async(MORSE_enum side, MORSE_enum uplo, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zsymm_Tile_Async( MORSE_enum side, MORSE_enum uplo, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_Complex64_t beta, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; @@ -414,9 +415,11 @@ int MORSE_zsymm_Tile_Async(MORSE_enum side, MORSE_enum uplo, } /* Quick return */ - if (C->m == 0 || C->n == 0 || - ( (alpha == (MORSE_Complex64_t)0.0) && (beta == (MORSE_Complex64_t)1.0) )) + if ( (C->m == 0) || (C->n == 0) || + ( (alpha == (MORSE_Complex64_t)0.0) && (beta == (MORSE_Complex64_t)1.0) ) ) + { return MORSE_SUCCESS; + } morse_pzsymm( side, uplo, alpha, A, B, beta, C, sequence, request ); diff --git a/compute/zsyr2k.c b/compute/zsyr2k.c index 2b3ae43b86c16cc700790f6e22e0c58ce66db332..84a501a4808176c987358ebbca2a9d9a0a5a2d27 100644 --- a/compute/zsyr2k.c +++ b/compute/zsyr2k.c @@ -108,9 +108,9 @@ * @sa MORSE_ssyr2k * ******************************************************************************/ -int MORSE_zsyr2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, +int MORSE_zsyr2k( MORSE_enum uplo, MORSE_enum trans, int N, int K, MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, MORSE_Complex64_t *B, int LDB, - MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC) + MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am, An; @@ -178,7 +178,7 @@ int MORSE_zsyr2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -199,7 +199,7 @@ int MORSE_zsyr2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); @@ -207,7 +207,7 @@ int MORSE_zsyr2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -265,9 +265,9 @@ int MORSE_zsyr2k(MORSE_enum uplo, MORSE_enum trans, int N, int K, * @sa MORSE_ssyr2k * ******************************************************************************/ -int MORSE_zsyr2k_Tile(MORSE_enum uplo, MORSE_enum trans, +int MORSE_zsyr2k_Tile( MORSE_enum uplo, MORSE_enum trans, MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_Complex64_t beta, MORSE_desc_t *C) + MORSE_Complex64_t beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -279,16 +279,17 @@ int MORSE_zsyr2k_Tile(MORSE_enum uplo, MORSE_enum trans, morse_fatal_error("MORSE_zsyr2k_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsyr2k_Tile_Async( uplo, trans, alpha, A, B, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -320,10 +321,10 @@ int MORSE_zsyr2k_Tile(MORSE_enum uplo, MORSE_enum trans, * @sa MORSE_ssyr2k_Tile_Async * ******************************************************************************/ -int MORSE_zsyr2k_Tile_Async(MORSE_enum uplo, MORSE_enum trans, +int MORSE_zsyr2k_Tile_Async( MORSE_enum uplo, MORSE_enum trans, MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_Complex64_t beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int N, K; diff --git a/compute/zsyrk.c b/compute/zsyrk.c index f9ed92792fd8ac5ef4400a6da6a386a31a8864d6..2dd3e6ac15d9a225cbc45c5b79fde39396386a61 100644 --- a/compute/zsyrk.c +++ b/compute/zsyrk.c @@ -98,9 +98,9 @@ * @sa MORSE_ssyrk * ******************************************************************************/ -int MORSE_zsyrk(MORSE_enum uplo, MORSE_enum trans, int N, int K, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC) +int MORSE_zsyrk( MORSE_enum uplo, MORSE_enum trans, int N, int K, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t beta, MORSE_Complex64_t *C, int LDC ) { int NB; int Am, An; @@ -163,7 +163,7 @@ int MORSE_zsyrk(MORSE_enum uplo, MORSE_enum trans, int N, int K, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpperLower, @@ -180,14 +180,14 @@ int MORSE_zsyrk(MORSE_enum uplo, MORSE_enum trans, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -241,9 +241,9 @@ int MORSE_zsyrk(MORSE_enum uplo, MORSE_enum trans, int N, int K, * @sa MORSE_ssyrk * ******************************************************************************/ -int MORSE_zsyrk_Tile(MORSE_enum uplo, MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *C) +int MORSE_zsyrk_Tile( MORSE_enum uplo, MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -255,15 +255,16 @@ int MORSE_zsyrk_Tile(MORSE_enum uplo, MORSE_enum trans, morse_fatal_error("MORSE_zsyrk_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsyrk_Tile_Async( uplo, trans, alpha, A, beta, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -295,10 +296,10 @@ int MORSE_zsyrk_Tile(MORSE_enum uplo, MORSE_enum trans, * @sa MORSE_ssyrk_Tile_Async * ******************************************************************************/ -int MORSE_zsyrk_Tile_Async(MORSE_enum uplo, MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zsyrk_Tile_Async( MORSE_enum uplo, MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int N, K; diff --git a/compute/zsysv.c b/compute/zsysv.c index d787cd332741816d24c472056b9195b8a73725db..a53b8cb57eb575da39f9b05387dec10f8483b09f 100644 --- a/compute/zsysv.c +++ b/compute/zsysv.c @@ -92,9 +92,9 @@ * @sa MORSE_csysv * ******************************************************************************/ -int MORSE_zsysv(MORSE_enum uplo, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zsysv( MORSE_enum uplo, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -145,7 +145,7 @@ int MORSE_zsysv(MORSE_enum uplo, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -162,14 +162,14 @@ int MORSE_zsysv(MORSE_enum uplo, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -230,15 +230,16 @@ int MORSE_zsysv_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zsysv_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsysv_Tile_Async( uplo, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -269,8 +270,8 @@ int MORSE_zsysv_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_csysv_Tile_Async * ******************************************************************************/ -int MORSE_zsysv_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zsysv_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zsytrf.c b/compute/zsytrf.c index 04cfd22f2772aac5354c75e8b7e873bf72e7df36..11d74b77576fe3474cfd29ef9396a1788abd88c9 100644 --- a/compute/zsytrf.c +++ b/compute/zsytrf.c @@ -78,8 +78,8 @@ * @sa MORSE_zpotrs * ******************************************************************************/ -int MORSE_zsytrf(MORSE_enum uplo, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_zsytrf( MORSE_enum uplo, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -120,7 +120,7 @@ int MORSE_zsytrf(MORSE_enum uplo, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -133,13 +133,13 @@ int MORSE_zsytrf(MORSE_enum uplo, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -196,14 +196,15 @@ int MORSE_zsytrf_Tile( MORSE_enum uplo, MORSE_desc_t *A ) morse_fatal_error("MORSE_zsytrf_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsytrf_Tile_Async( uplo, A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -234,8 +235,8 @@ int MORSE_zsytrf_Tile( MORSE_enum uplo, MORSE_desc_t *A ) * @sa MORSE_csytrf_Tile_Async * ******************************************************************************/ -int MORSE_zsytrf_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zsytrf_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zsytrs.c b/compute/zsytrs.c index cc95b1c19170bce3d97babffcc79e7601a485d55..7223993d95d2b857c8eebdb664f1c318fecf9333 100644 --- a/compute/zsytrs.c +++ b/compute/zsytrs.c @@ -79,9 +79,9 @@ * @sa MORSE_csytrs * ******************************************************************************/ -int MORSE_zsytrs(MORSE_enum uplo, int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_zsytrs( MORSE_enum uplo, int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -131,7 +131,7 @@ int MORSE_zsytrs(MORSE_enum uplo, int N, int NRHS, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -148,14 +148,14 @@ int MORSE_zsytrs(MORSE_enum uplo, int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -209,15 +209,16 @@ int MORSE_zsytrs_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) morse_fatal_error("MORSE_zsytrs_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zsytrs_Tile_Async( uplo, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -248,8 +249,8 @@ int MORSE_zsytrs_Tile( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B ) * @sa MORSE_csytrs_Tile_Async * ******************************************************************************/ -int MORSE_zsytrs_Tile_Async(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zsytrs_Tile_Async( MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/ztile.c b/compute/ztile.c index 94fff08a0db71a9745965b5261539c75f645e242..868a565e4178542648efe0234f848cd32a651df6 100644 --- a/compute/ztile.c +++ b/compute/ztile.c @@ -61,7 +61,7 @@ * @sa MORSE_sLapack_to_Tile * ******************************************************************************/ -int MORSE_zLapack_to_Tile( MORSE_Complex64_t *Af77, int LDA, MORSE_desc_t *A ) +int MORSE_zLapack_to_Tile( MORSE_Complex64_t *Af77, int LDA, MORSE_desc_t *A ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -86,7 +86,7 @@ int MORSE_zLapack_to_Tile( MORSE_Complex64_t *Af77, int LDA, MORSE_desc_t *A ) morse_getaddr_cm, morse_getblkldd_cm, NULL ); /* Start the computation */ - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); morse_pzlacpy( MorseUpperLower, B, A, sequence, &request ); @@ -99,7 +99,7 @@ int MORSE_zLapack_to_Tile( MORSE_Complex64_t *Af77, int LDA, MORSE_desc_t *A ) MORSE_Desc_Destroy( &B ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -138,7 +138,7 @@ int MORSE_zLapack_to_Tile( MORSE_Complex64_t *Af77, int LDA, MORSE_desc_t *A ) * @sa MORSE_sTile_to_Lapack * ******************************************************************************/ -int MORSE_zTile_to_Lapack( MORSE_desc_t *A, MORSE_Complex64_t *Af77, int LDA ) +int MORSE_zTile_to_Lapack( MORSE_desc_t *A, MORSE_Complex64_t *Af77, int LDA ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -163,7 +163,7 @@ int MORSE_zTile_to_Lapack( MORSE_desc_t *A, MORSE_Complex64_t *Af77, int LDA ) morse_getaddr_cm, morse_getblkldd_cm, NULL ); /* Start the computation */ - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); morse_pzlacpy( MorseUpperLower, A, B, sequence, &request ); @@ -175,6 +175,6 @@ int MORSE_zTile_to_Lapack( MORSE_desc_t *A, MORSE_Complex64_t *Af77, int LDA ) MORSE_Desc_Destroy( &B ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/ztpgqrt.c b/compute/ztpgqrt.c index 8bb2753958e9fdcb2095a7d2290dc3cef290cccc..e209f1438b3971421b148ebd53bdcbf357a6d3b2 100644 --- a/compute/ztpgqrt.c +++ b/compute/ztpgqrt.c @@ -208,7 +208,7 @@ int MORSE_ztpgqrt( int M, int N, int K, int L, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descV1l, &descV1t, MorseDescInput, MorseUpperLower, @@ -233,7 +233,7 @@ int MORSE_ztpgqrt( int M, int N, int K, int L, morse_ztile2lap( morse, &descQ2l, &descQ2t, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descV1l, &descV1t ); @@ -242,7 +242,7 @@ int MORSE_ztpgqrt( int M, int N, int K, int L, morse_ztile2lap_cleanup( morse, &descQ2l, &descQ2t ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -289,15 +289,20 @@ int MORSE_ztpgqrt_Tile( int L, morse_fatal_error("MORSE_ztpgqrt_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztpgqrt_Tile_Async( L, V1, T1, V2, T2, Q1, Q2, sequence, &request ); + + MORSE_Desc_Flush( V1, sequence ); + MORSE_Desc_Flush( T1, sequence ); + MORSE_Desc_Flush( V2, sequence ); + MORSE_Desc_Flush( T2, sequence ); MORSE_Desc_Flush( Q1, sequence ); MORSE_Desc_Flush( Q2, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/ztpqrt.c b/compute/ztpqrt.c index ea1db7270c048d22ef08cd402680c4e6d9725071..9f5f60bbffdeac37cbad20448128fae9a053466d 100644 --- a/compute/ztpqrt.c +++ b/compute/ztpqrt.c @@ -184,7 +184,7 @@ int MORSE_ztpqrt( int M, int N, int L, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, MorseUpper, @@ -201,14 +201,14 @@ int MORSE_ztpqrt( int M, int N, int L, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -251,7 +251,7 @@ int MORSE_ztpqrt( int M, int N, int L, * @sa MORSE_zgeqrs_Tile * ******************************************************************************/ -int MORSE_ztpqrt_Tile( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T ) +int MORSE_ztpqrt_Tile( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -263,15 +263,17 @@ int MORSE_ztpqrt_Tile( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T morse_fatal_error("MORSE_ztpqrt_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztpqrt_Tile_Async( L, A, B, T, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); + MORSE_Desc_Flush( T, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/ztradd.c b/compute/ztradd.c index ef66d816d55d1e9090ea7e05244f743b83ed5f34..bd88723f3df62ac702ab43c91ee3899661725b9c 100644 --- a/compute/ztradd.c +++ b/compute/ztradd.c @@ -97,9 +97,9 @@ * @sa MORSE_stradd * ******************************************************************************/ -int MORSE_ztradd(MORSE_enum uplo, MORSE_enum trans, int M, int N, - MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t beta, MORSE_Complex64_t *B, int LDB) +int MORSE_ztradd( MORSE_enum uplo, MORSE_enum trans, int M, int N, + MORSE_Complex64_t alpha, MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t beta, MORSE_Complex64_t *B, int LDB ) { int NB; int Am, An; @@ -162,7 +162,7 @@ int MORSE_ztradd(MORSE_enum uplo, MORSE_enum trans, int M, int N, /* Set MT & NT & KT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -179,14 +179,14 @@ int MORSE_ztradd(MORSE_enum uplo, MORSE_enum trans, int M, int N, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -244,9 +244,9 @@ int MORSE_ztradd(MORSE_enum uplo, MORSE_enum trans, int M, int N, * @sa MORSE_stradd_Tile * ******************************************************************************/ -int MORSE_ztradd_Tile(MORSE_enum uplo, MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *B) +int MORSE_ztradd_Tile( MORSE_enum uplo, MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -258,15 +258,16 @@ int MORSE_ztradd_Tile(MORSE_enum uplo, MORSE_enum trans, morse_fatal_error("MORSE_ztradd_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztradd_Tile_Async( uplo, trans, alpha, A, beta, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -299,10 +300,10 @@ int MORSE_ztradd_Tile(MORSE_enum uplo, MORSE_enum trans, * @sa MORSE_stradd_Tile_Async * ******************************************************************************/ -int MORSE_ztradd_Tile_Async(MORSE_enum uplo, MORSE_enum trans, - MORSE_Complex64_t alpha, MORSE_desc_t *A, - MORSE_Complex64_t beta, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_ztradd_Tile_Async( MORSE_enum uplo, MORSE_enum trans, + MORSE_Complex64_t alpha, MORSE_desc_t *A, + MORSE_Complex64_t beta, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; int M, N; diff --git a/compute/ztrmm.c b/compute/ztrmm.c index 099c9cf00bb0110a2c06099521736fc0791c8f54..3496308a72a917e781343b596c7836c4fe9a4d59 100644 --- a/compute/ztrmm.c +++ b/compute/ztrmm.c @@ -100,11 +100,11 @@ * @sa MORSE_strmm * ******************************************************************************/ -int MORSE_ztrmm(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - int N, int NRHS, MORSE_Complex64_t alpha, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_ztrmm( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + int N, int NRHS, MORSE_Complex64_t alpha, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB, NA; int status; @@ -176,7 +176,7 @@ int MORSE_ztrmm(MORSE_enum side, MORSE_enum uplo, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -193,14 +193,14 @@ int MORSE_ztrmm(MORSE_enum side, MORSE_enum uplo, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -266,9 +266,9 @@ int MORSE_ztrmm(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_strmm_Tile * ******************************************************************************/ -int MORSE_ztrmm_Tile(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B) +int MORSE_ztrmm_Tile( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -280,15 +280,16 @@ int MORSE_ztrmm_Tile(MORSE_enum side, MORSE_enum uplo, morse_fatal_error("MORSE_ztrmm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztrmm_Tile_Async(side, uplo, transA, diag, alpha, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -320,10 +321,10 @@ int MORSE_ztrmm_Tile(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_strmm_Tile_Async * ******************************************************************************/ -int MORSE_ztrmm_Tile_Async(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_ztrmm_Tile_Async( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/ztrsm.c b/compute/ztrsm.c index 1619d391b1659aa2ab6e17f3c267053b33f00560..053be9e3e2774c022818fbd7fdc6c5c783258285 100644 --- a/compute/ztrsm.c +++ b/compute/ztrsm.c @@ -101,11 +101,11 @@ * @sa MORSE_strsm * ******************************************************************************/ -int MORSE_ztrsm(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - int N, int NRHS, MORSE_Complex64_t alpha, - MORSE_Complex64_t *A, int LDA, - MORSE_Complex64_t *B, int LDB) +int MORSE_ztrsm( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + int N, int NRHS, MORSE_Complex64_t alpha, + MORSE_Complex64_t *A, int LDA, + MORSE_Complex64_t *B, int LDB ) { int NB, NA; int status; @@ -174,7 +174,7 @@ int MORSE_ztrsm(MORSE_enum side, MORSE_enum uplo, /* Set NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, uplo, @@ -191,14 +191,14 @@ int MORSE_ztrsm(MORSE_enum side, MORSE_enum uplo, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } /** @@ -264,9 +264,9 @@ int MORSE_ztrsm(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_strsm_Tile * ******************************************************************************/ -int MORSE_ztrsm_Tile(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B) +int MORSE_ztrsm_Tile( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -278,15 +278,16 @@ int MORSE_ztrsm_Tile(MORSE_enum side, MORSE_enum uplo, morse_fatal_error("MORSE_ztrsm_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztrsm_Tile_Async(side, uplo, transA, diag, alpha, A, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -318,10 +319,10 @@ int MORSE_ztrsm_Tile(MORSE_enum side, MORSE_enum uplo, * @sa MORSE_strsm_Tile_Async * ******************************************************************************/ -int MORSE_ztrsm_Tile_Async(MORSE_enum side, MORSE_enum uplo, - MORSE_enum transA, MORSE_enum diag, - MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_ztrsm_Tile_Async( MORSE_enum side, MORSE_enum uplo, + MORSE_enum transA, MORSE_enum diag, + MORSE_Complex64_t alpha, MORSE_desc_t *A, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/ztrsmpl.c b/compute/ztrsmpl.c index 50a37e418236a4b6cdb5c8284fdd3107160852af..d84670ad98dc9fcf87bcccfd63c73b3538380f2c 100644 --- a/compute/ztrsmpl.c +++ b/compute/ztrsmpl.c @@ -80,10 +80,10 @@ * @sa MORSE_zgetrf_incpiv * ******************************************************************************/ -int MORSE_ztrsmpl(int N, int NRHS, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descL, int *IPIV, - MORSE_Complex64_t *B, int LDB) +int MORSE_ztrsmpl( int N, int NRHS, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descL, int *IPIV, + MORSE_Complex64_t *B, int LDB ) { int NB; int status; @@ -129,7 +129,7 @@ int MORSE_ztrsmpl(int N, int NRHS, /* Set Mt, NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseLower, @@ -146,14 +146,14 @@ int MORSE_ztrsmpl(int N, int NRHS, morse_ztile2lap( morse, &descBl, &descBt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descBl, &descBt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -208,15 +208,17 @@ int MORSE_ztrsmpl_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_ morse_fatal_error("MORSE_ztrsmpl_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztrsmpl_Tile_Async( A, L, IPIV, B, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( L, sequence ); MORSE_Desc_Flush( B, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -250,8 +252,8 @@ int MORSE_ztrsmpl_Tile( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_ * @sa MORSE_zgetrf_incpiv_Tile_Async * ******************************************************************************/ -int MORSE_ztrsmpl_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_ztrsmpl_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_desc_t *B, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/ztrtri.c b/compute/ztrtri.c index d911cf16f1e93123b332e6a81a9835c496d6f1b7..e4d822d94fce6c2b348925b26faf29c0bed71cfa 100644 --- a/compute/ztrtri.c +++ b/compute/ztrtri.c @@ -82,8 +82,8 @@ * @sa MORSE_zpotri * ******************************************************************************/ -int MORSE_ztrtri(MORSE_enum uplo, MORSE_enum diag, int N, - MORSE_Complex64_t *A, int LDA) +int MORSE_ztrtri( MORSE_enum uplo, MORSE_enum diag, int N, + MORSE_Complex64_t *A, int LDA ) { int NB; int status; @@ -128,7 +128,7 @@ int MORSE_ztrtri(MORSE_enum uplo, MORSE_enum diag, int N, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInout, uplo, @@ -141,13 +141,13 @@ int MORSE_ztrtri(MORSE_enum uplo, MORSE_enum diag, int N, morse_ztile2lap( morse, &descAl, &descAt, MorseDescInout, uplo, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -214,14 +214,15 @@ int MORSE_ztrtri_Tile( MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A ) morse_fatal_error("MORSE_ztrtri_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_ztrtri_Tile_Async( uplo, diag, A, sequence, &request ); - MORSE_Desc_Flush( A, sequence ); - morse_sequence_wait(morse, sequence); + MORSE_Desc_Flush( A, sequence ); + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -255,8 +256,8 @@ int MORSE_ztrtri_Tile( MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A ) * @sa MORSE_zpotri_Tile_Async * ******************************************************************************/ -int MORSE_ztrtri_Tile_Async(MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_ztrtri_Tile_Async( MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; diff --git a/compute/zunglq.c b/compute/zunglq.c index 07546a9e4e5d95ba7f9e4b0255040f8a2f881c5c..a2307eb05f042f8158eb63d98251422142412f45 100644 --- a/compute/zunglq.c +++ b/compute/zunglq.c @@ -81,10 +81,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zunglq(int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *Q, int LDQ) +int MORSE_zunglq( int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *Q, int LDQ ) { int NB; int status; @@ -135,7 +135,7 @@ int MORSE_zunglq(int M, int N, int K, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpper, @@ -152,14 +152,14 @@ int MORSE_zunglq(int M, int N, int K, morse_ztile2lap( morse, &descQl, &descQt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descQl, &descQt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -210,15 +210,17 @@ int MORSE_zunglq_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q ) morse_fatal_error("MORSE_zunglq_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunglq_Tile_Async( A, T, Q, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( Q, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -250,8 +252,8 @@ int MORSE_zunglq_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q ) * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zunglq_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunglq_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zunglq_param.c b/compute/zunglq_param.c index 3f14ca1f38d386a1d24674070e688975d95784e2..af1096930c7b7663f2a6d952deaae1ad5834f1f6 100644 --- a/compute/zunglq_param.c +++ b/compute/zunglq_param.c @@ -78,10 +78,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zunglq_param(const libhqr_tree_t *qrtree, int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *Q, int LDQ) +int MORSE_zunglq_param( const libhqr_tree_t *qrtree, int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *Q, int LDQ ) { int NB; int status; @@ -132,7 +132,7 @@ int MORSE_zunglq_param(const libhqr_tree_t *qrtree, int M, int N, int K, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpper, @@ -149,14 +149,14 @@ int MORSE_zunglq_param(const libhqr_tree_t *qrtree, int M, int N, int K, morse_ztile2lap( morse, &descQl, &descQt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descQl, &descQt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -207,15 +207,18 @@ int MORSE_zunglq_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zunglq_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunglq_param_Tile_Async( qrtree, A, TS, TT, Q, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( Q, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -247,8 +250,8 @@ int MORSE_zunglq_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zunglq_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *Q, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunglq_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *Q, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zungqr.c b/compute/zungqr.c index ac15a493db86601fa3f477cdbeda6e6d723bff2b..eb4ecf58ecf238a6da607d04dd48c811959b78e6 100644 --- a/compute/zungqr.c +++ b/compute/zungqr.c @@ -81,10 +81,10 @@ * @sa MORSE_zgeqrf * ******************************************************************************/ -int MORSE_zungqr(int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *Q, int LDQ) +int MORSE_zungqr( int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *Q, int LDQ ) { int NB; int status; @@ -134,7 +134,7 @@ int MORSE_zungqr(int M, int N, int K, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseLower, @@ -151,14 +151,14 @@ int MORSE_zungqr(int M, int N, int K, morse_ztile2lap( morse, &descQl, &descQt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descQl, &descQt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -209,15 +209,17 @@ int MORSE_zungqr_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q ) morse_fatal_error("MORSE_zungqr_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zungqr_Tile_Async( A, T, Q, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( Q, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -249,8 +251,8 @@ int MORSE_zungqr_Tile( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q ) * @sa MORSE_zgeqrf_Tile_Async * ******************************************************************************/ -int MORSE_zungqr_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zungqr_Tile_Async( MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *Q, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zungqr_param.c b/compute/zungqr_param.c index 02b6079f997d32150f15ad821851301210344666..3607664f9371fbb07285811e2e60dec7281d6fb6 100644 --- a/compute/zungqr_param.c +++ b/compute/zungqr_param.c @@ -78,12 +78,12 @@ * @sa MORSE_zgeqrf * ******************************************************************************/ -int MORSE_zungqr_param(const libhqr_tree_t *qrtree, - int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, - MORSE_desc_t *descTT, - MORSE_Complex64_t *Q, int LDQ) +int MORSE_zungqr_param( const libhqr_tree_t *qrtree, + int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, + MORSE_desc_t *descTT, + MORSE_Complex64_t *Q, int LDQ ) { int NB; int status; @@ -133,7 +133,7 @@ int MORSE_zungqr_param(const libhqr_tree_t *qrtree, /* Set NT */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseLower, @@ -150,14 +150,14 @@ int MORSE_zungqr_param(const libhqr_tree_t *qrtree, morse_ztile2lap( morse, &descQl, &descQt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descQl, &descQt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -208,15 +208,18 @@ int MORSE_zungqr_param_Tile( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE morse_fatal_error("MORSE_zungqr_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zungqr_param_Tile_Async( qrtree, A, TS, TT, Q, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( Q, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } diff --git a/compute/zunmlq.c b/compute/zunmlq.c index 3cf4090f238233000a0c102cb031d49f332eb0ec..318906b5e9d0118d7f8b08c4720758d19a8ad0b1 100644 --- a/compute/zunmlq.c +++ b/compute/zunmlq.c @@ -105,10 +105,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zunmlq(MORSE_enum side, MORSE_enum trans, int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *C, int LDC) +int MORSE_zunmlq( MORSE_enum side, MORSE_enum trans, int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *C, int LDC ) { int NB, An; int status; @@ -172,7 +172,7 @@ int MORSE_zunmlq(MORSE_enum side, MORSE_enum trans, int M, int N, int K, /* Set MT, NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpper, @@ -189,14 +189,14 @@ int MORSE_zunmlq(MORSE_enum side, MORSE_enum trans, int M, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -249,8 +249,8 @@ int MORSE_zunmlq(MORSE_enum side, MORSE_enum trans, int M, int N, int K, * @sa MORSE_zgelqf_Tile * ******************************************************************************/ -int MORSE_zunmlq_Tile(MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C) +int MORSE_zunmlq_Tile( MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -262,15 +262,17 @@ int MORSE_zunmlq_Tile(MORSE_enum side, MORSE_enum trans, morse_fatal_error("MORSE_zunmlq_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunmlq_Tile_Async(side, trans, A, T, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -302,9 +304,9 @@ int MORSE_zunmlq_Tile(MORSE_enum side, MORSE_enum trans, * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zunmlq_Tile_Async(MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunmlq_Tile_Async( MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zunmlq_param.c b/compute/zunmlq_param.c index e856971648921fd53dca5e4c49a6181b8c4f694f..a272452ef9093ee3e591211faea19b070e4b5a6b 100644 --- a/compute/zunmlq_param.c +++ b/compute/zunmlq_param.c @@ -104,10 +104,10 @@ * @sa MORSE_zgelqf * ******************************************************************************/ -int MORSE_zunmlq_param(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *C, int LDC) +int MORSE_zunmlq_param( const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *C, int LDC ) { int NB, An; int status; @@ -171,7 +171,7 @@ int MORSE_zunmlq_param(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum /* Set MT, NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseUpper, @@ -188,14 +188,14 @@ int MORSE_zunmlq_param(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -248,8 +248,8 @@ int MORSE_zunmlq_param(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum * @sa MORSE_zgelqf_Tile * ******************************************************************************/ -int MORSE_zunmlq_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C) +int MORSE_zunmlq_param_Tile( const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -261,15 +261,18 @@ int MORSE_zunmlq_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_ morse_fatal_error("MORSE_zunmlq_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunmlq_param_Tile_Async( qrtree, side, trans, A, TS, TT, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -301,9 +304,9 @@ int MORSE_zunmlq_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_ * @sa MORSE_zgelqf_Tile_Async * ******************************************************************************/ -int MORSE_zunmlq_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunmlq_param_Tile_Async( const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zunmqr.c b/compute/zunmqr.c index d292d049bd5f71d0959f86105c63ff79cad70dfe..c859ad44ed1d8992fa960c836a338bacb18beb5d 100644 --- a/compute/zunmqr.c +++ b/compute/zunmqr.c @@ -107,10 +107,10 @@ * @sa MORSE_zgeqrf * ******************************************************************************/ -int MORSE_zunmqr(MORSE_enum side, MORSE_enum trans, int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descT, - MORSE_Complex64_t *C, int LDC) +int MORSE_zunmqr( MORSE_enum side, MORSE_enum trans, int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descT, + MORSE_Complex64_t *C, int LDC ) { int NB, Am; int status; @@ -175,7 +175,7 @@ int MORSE_zunmqr(MORSE_enum side, MORSE_enum trans, int M, int N, int K, /* Set MT, NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseLower, @@ -192,14 +192,14 @@ int MORSE_zunmqr(MORSE_enum side, MORSE_enum trans, int M, int N, int K, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -253,8 +253,8 @@ int MORSE_zunmqr(MORSE_enum side, MORSE_enum trans, int M, int N, int K, * @sa MORSE_zgeqrf_Tile * ******************************************************************************/ -int MORSE_zunmqr_Tile(MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C) +int MORSE_zunmqr_Tile( MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -266,15 +266,17 @@ int MORSE_zunmqr_Tile(MORSE_enum side, MORSE_enum trans, morse_fatal_error("MORSE_zunmqr_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunmqr_Tile_Async(side, trans, A, T, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( T, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -306,9 +308,9 @@ int MORSE_zunmqr_Tile(MORSE_enum side, MORSE_enum trans, * @sa MORSE_zgeqrf_Tile_Async * ******************************************************************************/ -int MORSE_zunmqr_Tile_Async(MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunmqr_Tile_Async( MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL; diff --git a/compute/zunmqr_param.c b/compute/zunmqr_param.c index f0ccc22d9da72ecde50ea454df13d7d064b336d2..7c1f2fb1571d79cfb0e876a4bf508f139edbb99f 100644 --- a/compute/zunmqr_param.c +++ b/compute/zunmqr_param.c @@ -107,11 +107,11 @@ * @sa MORSE_zgeqrf * ******************************************************************************/ -int MORSE_zunmqr_param(const libhqr_tree_t *qrtree, - MORSE_enum side, MORSE_enum trans, int M, int N, int K, - MORSE_Complex64_t *A, int LDA, - MORSE_desc_t *descTS, MORSE_desc_t *descTT, - MORSE_Complex64_t *C, int LDC) +int MORSE_zunmqr_param( const libhqr_tree_t *qrtree, + MORSE_enum side, MORSE_enum trans, int M, int N, int K, + MORSE_Complex64_t *A, int LDA, + MORSE_desc_t *descTS, MORSE_desc_t *descTT, + MORSE_Complex64_t *C, int LDC ) { int NB, Am; int status; @@ -176,7 +176,7 @@ int MORSE_zunmqr_param(const libhqr_tree_t *qrtree, /* Set MT, NT & NTRHS */ NB = MORSE_NB; - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); /* Submit the matrix conversion */ morse_zlap2tile( morse, &descAl, &descAt, MorseDescInput, MorseLower, @@ -193,14 +193,14 @@ int MORSE_zunmqr_param(const libhqr_tree_t *qrtree, morse_ztile2lap( morse, &descCl, &descCt, MorseDescInout, MorseUpperLower, sequence, &request ); - morse_sequence_wait(morse, sequence); + morse_sequence_wait( morse, sequence ); /* Cleanup the temporary data */ morse_ztile2lap_cleanup( morse, &descAl, &descAt ); morse_ztile2lap_cleanup( morse, &descCl, &descCt ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -254,8 +254,8 @@ int MORSE_zunmqr_param(const libhqr_tree_t *qrtree, * @sa MORSE_zgeqrf_Tile * ******************************************************************************/ -int MORSE_zunmqr_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C) +int MORSE_zunmqr_param_Tile( const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C ) { MORSE_context_t *morse; MORSE_sequence_t *sequence = NULL; @@ -267,15 +267,18 @@ int MORSE_zunmqr_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_ morse_fatal_error("MORSE_zunmqr_param_Tile", "MORSE not initialized"); return MORSE_ERR_NOT_INITIALIZED; } - morse_sequence_create(morse, &sequence); + morse_sequence_create( morse, &sequence ); + MORSE_zunmqr_param_Tile_Async( qrtree, side, trans, A, TS, TT, C, sequence, &request ); + MORSE_Desc_Flush( A, sequence ); + MORSE_Desc_Flush( TS, sequence ); + MORSE_Desc_Flush( TT, sequence ); MORSE_Desc_Flush( C, sequence ); - morse_sequence_wait(morse, sequence); - + morse_sequence_wait( morse, sequence ); status = sequence->status; - morse_sequence_destroy(morse, sequence); + morse_sequence_destroy( morse, sequence ); return status; } @@ -307,10 +310,10 @@ int MORSE_zunmqr_param_Tile(const libhqr_tree_t *qrtree, MORSE_enum side, MORSE_ * @sa MORSE_zgeqrf_Tile_Async * ******************************************************************************/ -int MORSE_zunmqr_param_Tile_Async(const libhqr_tree_t *qrtree, - MORSE_enum side, MORSE_enum trans, - MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C, - MORSE_sequence_t *sequence, MORSE_request_t *request) +int MORSE_zunmqr_param_Tile_Async( const libhqr_tree_t *qrtree, + MORSE_enum side, MORSE_enum trans, + MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_desc_t *C, + MORSE_sequence_t *sequence, MORSE_request_t *request ) { MORSE_context_t *morse; MORSE_desc_t D, *Dptr = NULL;