Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
AGULLO Emmanuel
Chameleon
Commits
35fe8c53
Commit
35fe8c53
authored
Jan 22, 2018
by
Mathieu Faverge
Browse files
Indent and whitespace cleanup
parent
4fd458c7
Changes
58
Hide whitespace changes
Inline
Side-by-side
compute/zgeadd.c
View file @
35fe8c53
...
...
@@ -164,10 +164,10 @@ int MORSE_zgeadd(MORSE_enum trans, int M, int N,
trans
,
alpha
,
&
descA
,
beta
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
compute/zgelqf.c
View file @
35fe8c53
...
...
@@ -76,8 +76,8 @@
*
******************************************************************************/
int
MORSE_zgelqf
(
int
M
,
int
N
,
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
)
{
int
NB
;
int
status
;
...
...
@@ -130,9 +130,9 @@ int MORSE_zgelqf(int M, int N,
MORSE_zgelqf_Tile_Async
(
&
descA
,
descT
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -230,7 +230,7 @@ int MORSE_zgelqf_Tile(MORSE_desc_t *A, MORSE_desc_t *T)
*
******************************************************************************/
int
MORSE_zgelqf_Tile_Async
(
MORSE_desc_t
*
A
,
MORSE_desc_t
*
T
,
MORSE_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
MORSE_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_context_t
*
morse
;
MORSE_desc_t
D
,
*
Dptr
=
NULL
;
...
...
@@ -269,10 +269,10 @@ int MORSE_zgelqf_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
#if defined(CHAMELEON_COPY_DIAG)
{
int
m
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
mb
;
...
...
compute/zgelqf_param.c
View file @
35fe8c53
...
...
@@ -127,9 +127,9 @@ int MORSE_zgelqf_param(const libhqr_tree_t *qrtree, int M, int N,
MORSE_zgelqf_param_Tile_Async
(
qrtree
,
&
descA
,
descTS
,
descTT
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -270,10 +270,10 @@ int MORSE_zgelqf_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
#if defined(CHAMELEON_COPY_DIAG)
{
int
m
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
mb
;
...
...
compute/zgelqs.c
View file @
35fe8c53
...
...
@@ -82,9 +82,9 @@
*
******************************************************************************/
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
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
,
MORSE_Complex64_t
*
B
,
int
LDB
)
{
int
NB
;
int
status
;
...
...
@@ -147,8 +147,8 @@ int MORSE_zgelqs(int M, int N, int NRHS,
MORSE_zgelqs_Tile_Async
(
&
descA
,
descT
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
...
...
@@ -251,7 +251,7 @@ int MORSE_zgelqs_Tile(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B)
*
******************************************************************************/
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_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_desc_t
*
subB
;
MORSE_desc_t
*
subA
;
...
...
@@ -296,11 +296,11 @@ int MORSE_zgelqs_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
/* subB = morse_desc_submatrix(B, A->m, 0, A->n-A->m, B->n);
morse_pzlaset( MorseUpperLower, 0., 0., subB, sequence, request );
free(subB); */
...
...
compute/zgelqs_param.c
View file @
35fe8c53
...
...
@@ -149,8 +149,8 @@ int MORSE_zgelqs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS,
MORSE_zgelqs_param_Tile_Async
(
qrtree
,
&
descA
,
descTS
,
descTT
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
...
...
@@ -305,14 +305,14 @@ int MORSE_zgelqs_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
/* subB = morse_desc_submatrix(B, A->m, 0, A->n-A->m, B->n);
morse_pzlaset( MorseUpperLower, 0., 0., subB, sequence, request);
free(subB); */
morse_pzlaset( MorseUpperLower, 0., 0., subB, sequence, request);
free(subB); */
subB
=
morse_desc_submatrix
(
B
,
0
,
0
,
A
->
m
,
B
->
n
);
subA
=
morse_desc_submatrix
(
A
,
0
,
0
,
A
->
m
,
A
->
m
);
...
...
compute/zgels.c
View file @
35fe8c53
...
...
@@ -108,9 +108,9 @@
*
******************************************************************************/
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
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
,
MORSE_Complex64_t
*
B
,
int
LDB
)
{
int
i
,
j
;
int
NB
;
...
...
@@ -189,11 +189,11 @@ int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS,
MORSE_zgels_Tile_Async
(
MorseNoTrans
,
&
descA
,
descT
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -255,7 +255,7 @@ int MORSE_zgels(MORSE_enum trans, int M, int N, int NRHS,
*
******************************************************************************/
int
MORSE_zgels_Tile
(
MORSE_enum
trans
,
MORSE_desc_t
*
A
,
MORSE_desc_t
*
T
,
MORSE_desc_t
*
B
)
MORSE_desc_t
*
T
,
MORSE_desc_t
*
B
)
{
MORSE_context_t
*
morse
;
MORSE_sequence_t
*
sequence
=
NULL
;
...
...
@@ -358,12 +358,12 @@ int MORSE_zgels_Tile_Async(MORSE_enum trans, MORSE_desc_t *A,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_NOT_SUPPORTED
);
}
/* Quick return - currently NOT equivalent to LAPACK's:
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
for (i = 0; i < chameleon_max(M, N); i++)
for (j = 0; j < NRHS; j++)
B[j*LDB+i] = 0.0;
return MORSE_SUCCESS;
}
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
for (i = 0; i < chameleon_max(M, N); i++)
for (j = 0; j < NRHS; j++)
B[j*LDB+i] = 0.0;
return MORSE_SUCCESS;
}
*/
if
(
A
->
m
>=
A
->
n
)
{
...
...
@@ -392,8 +392,8 @@ int MORSE_zgels_Tile_Async(MORSE_enum trans, MORSE_desc_t *A,
}
else
{
/* subB = morse_desc_submatrix(B, A->m, 0, A->n-A->m, B->n);
morse_pzlaset( MorseUpperLower, 0., 0., subB, sequence, request);
free(subB); */
morse_pzlaset( MorseUpperLower, 0., 0., subB, sequence, request);
free(subB); */
#if defined(CHAMELEON_COPY_DIAG)
{
int
m
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
mb
;
...
...
compute/zgemm.c
View file @
35fe8c53
...
...
@@ -131,9 +131,9 @@
*
******************************************************************************/
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
)
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
;
...
...
@@ -346,9 +346,9 @@ int MORSE_zgemm_Tile(MORSE_enum transA, MORSE_enum transB,
*
******************************************************************************/
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_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
;
...
...
compute/zgeqrf.c
View file @
35fe8c53
...
...
@@ -75,8 +75,8 @@
*
******************************************************************************/
int
MORSE_zgeqrf
(
int
M
,
int
N
,
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
)
{
int
NB
;
int
status
;
...
...
@@ -129,9 +129,9 @@ int MORSE_zgeqrf(int M, int N,
MORSE_zgeqrf_Tile_Async
(
&
descA
,
descT
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -229,7 +229,7 @@ int MORSE_zgeqrf_Tile(MORSE_desc_t *A, MORSE_desc_t *T)
*
******************************************************************************/
int
MORSE_zgeqrf_Tile_Async
(
MORSE_desc_t
*
A
,
MORSE_desc_t
*
T
,
MORSE_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
MORSE_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_context_t
*
morse
;
MORSE_desc_t
D
,
*
Dptr
=
NULL
;
...
...
@@ -268,10 +268,10 @@ int MORSE_zgeqrf_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
#if defined(CHAMELEON_COPY_DIAG)
{
int
n
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
nb
;
...
...
compute/zgeqrf_param.c
View file @
35fe8c53
...
...
@@ -129,12 +129,12 @@ int MORSE_zgeqrf_param(const libhqr_tree_t *qrtree, int M, int N,
A
,
NB
,
NB
,
LDA
,
N
,
M
,
N
,
sequence
,
&
request
);
/* Call the tile interface */
MORSE_zgeqrf_param_Tile_Async
(
qrtree
,
&
descA
,
descTS
,
descTT
,
sequence
,
&
request
);
MORSE_zgeqrf_param_Tile_Async
(
qrtree
,
&
descA
,
descTS
,
descTT
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -236,7 +236,7 @@ int MORSE_zgeqrf_param_Tile(const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_
*
******************************************************************************/
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_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_context_t
*
morse
;
MORSE_desc_t
D
,
*
Dptr
=
NULL
;
...
...
@@ -279,10 +279,10 @@ int MORSE_zgeqrf_param_Tile_Async(const libhqr_tree_t *qrtree, MORSE_desc_t *A,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
#if defined(CHAMELEON_COPY_DIAG)
{
int
n
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
nb
;
...
...
compute/zgeqrs.c
View file @
35fe8c53
...
...
@@ -82,9 +82,9 @@
*
*******************************************************************************/
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
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descT
,
MORSE_Complex64_t
*
B
,
int
LDB
)
{
int
NB
;
int
status
;
...
...
@@ -147,12 +147,12 @@ int MORSE_zgeqrs(int M, int N, int NRHS,
MORSE_zgeqrs_Tile_Async
(
&
descA
,
descT
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
return
status
;
...
...
@@ -251,7 +251,7 @@ int MORSE_zgeqrs_Tile(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B)
*
******************************************************************************/
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_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_desc_t
*
subA
;
MORSE_desc_t
*
subB
;
...
...
@@ -296,11 +296,11 @@ int MORSE_zgeqrs_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *B,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
/*
if (chameleon_min(M, chameleon_min(N, NRHS)) == 0) {
return MORSE_SUCCESS;
}
*/
#if defined(CHAMELEON_COPY_DIAG)
{
int
n
=
chameleon_min
(
A
->
mt
,
A
->
nt
)
*
A
->
nb
;
...
...
compute/zgeqrs_param.c
View file @
35fe8c53
...
...
@@ -142,11 +142,11 @@ int MORSE_zgeqrs_param(const libhqr_tree_t *qrtree, int M, int N, int NRHS,
MORSE_zgeqrs_param_Tile_Async
(
qrtree
,
&
descA
,
descTS
,
descTT
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
compute/zgesv_incpiv.c
View file @
35fe8c53
...
...
@@ -87,9 +87,9 @@
*
******************************************************************************/
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
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descL
,
int
*
IPIV
,
MORSE_Complex64_t
*
B
,
int
LDB
)
{
int
NB
;
int
status
;
...
...
@@ -146,12 +146,12 @@ int MORSE_zgesv_incpiv(int N, int NRHS,
MORSE_zgesv_incpiv_Tile_Async
(
&
descA
,
descL
,
IPIV
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
return
status
;
...
...
@@ -218,7 +218,7 @@ int MORSE_zgesv_incpiv_Tile(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_d
RUNTIME_desc_flush
(
A
,
sequence
);
RUNTIME_desc_flush
(
B
,
sequence
);
morse_sequence_wait
(
morse
,
sequence
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
return
status
;
...
...
@@ -255,7 +255,7 @@ int MORSE_zgesv_incpiv_Tile(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, MORSE_d
*
******************************************************************************/
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_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_context_t
*
morse
;
...
...
@@ -297,10 +297,10 @@ int MORSE_zgesv_incpiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV, M
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(N, NRHS) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(N, NRHS) == 0)
return MORSE_SUCCESS;
*/
morse_pzgetrf_incpiv
(
A
,
L
,
IPIV
,
sequence
,
request
);
...
...
compute/zgesv_nopiv.c
View file @
35fe8c53
...
...
@@ -145,12 +145,12 @@ int MORSE_zgesv_nopiv(int N, int NRHS,
MORSE_zgesv_nopiv_Tile_Async
(
&
descA
,
&
descB
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_zooptile2lap
(
descB
,
B
,
NB
,
NB
,
LDB
,
NRHS
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_desc_mat_free
(
&
descB
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
return
status
;
...
...
@@ -210,7 +210,7 @@ int MORSE_zgesv_nopiv_Tile(MORSE_desc_t *A, MORSE_desc_t *B)
RUNTIME_desc_flush
(
A
,
sequence
);
RUNTIME_desc_flush
(
B
,
sequence
);
morse_sequence_wait
(
morse
,
sequence
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
return
status
;
...
...
@@ -285,10 +285,10 @@ int MORSE_zgesv_nopiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *B,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(N, NRHS) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(N, NRHS) == 0)
return MORSE_SUCCESS;
*/
morse_pzgetrf_nopiv
(
A
,
sequence
,
request
);
...
...
compute/zgesvd.c
View file @
35fe8c53
...
...
@@ -148,12 +148,12 @@
*
******************************************************************************/
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
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
;
...
...
compute/zgetrf_incpiv.c
View file @
35fe8c53
...
...
@@ -79,8 +79,8 @@
*
******************************************************************************/
int
MORSE_zgetrf_incpiv
(
int
M
,
int
N
,
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descL
,
int
*
IPIV
)
MORSE_Complex64_t
*
A
,
int
LDA
,
MORSE_desc_t
*
descL
,
int
*
IPIV
)
{
int
NB
;
int
status
;
...
...
@@ -131,9 +131,9 @@ int MORSE_zgetrf_incpiv(int M, int N,
MORSE_zgetrf_incpiv_Tile_Async
(
&
descA
,
descL
,
IPIV
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);
status
=
sequence
->
status
;
morse_sequence_destroy
(
morse
,
sequence
);
...
...
@@ -234,7 +234,7 @@ int MORSE_zgetrf_incpiv_Tile(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV)
*
******************************************************************************/
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_sequence_t
*
sequence
,
MORSE_request_t
*
request
)
{
MORSE_context_t
*
morse
;
...
...
@@ -272,10 +272,10 @@ int MORSE_zgetrf_incpiv_Tile_Async(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV,
return
morse_request_fail
(
sequence
,
request
,
MORSE_ERR_ILLEGAL_VALUE
);
}
/* Quick return */
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
/*
if (chameleon_min(M, N) == 0)
return MORSE_SUCCESS;
*/
morse_pzgetrf_incpiv
(
A
,
L
,
IPIV
,
sequence
,
request
);
...
...
compute/zgetrf_nopiv.c
View file @
35fe8c53
...
...
@@ -124,9 +124,9 @@ int MORSE_zgetrf_nopiv(int M, int N,
MORSE_zgetrf_nopiv_Tile_Async
(
&
descA
,
sequence
,
&
request
);
/* Submit the matrix conversion */
morse_zooptile2lap
(
descA
,
A
,
NB
,
NB
,
LDA
,
N
,
sequence
,
&
request
);
morse_sequence_wait
(
morse
,
sequence
);
morse_desc_mat_free
(
&
descA
);