Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
AGULLO Emmanuel
Chameleon
Commits
a3481c8e
Commit
a3481c8e
authored
Mar 03, 2017
by
THIBAULT Samuel
Browse files
Rather use iteration_push/pop
parent
67246b24
Changes
27
Hide whitespace changes
Inline
Side-by-side
CMakeLists.txt
View file @
a3481c8e
...
...
@@ -706,10 +706,10 @@ if( CHAMELEON_SCHED_STARPU )
set
(
HAVE_STARPU_IDLE_PREFETCH 1
)
message
(
"--
${
Blue
}
Add definition HAVE_STARPU_IDLE_PREFETCH
${
ColourReset
}
"
)
endif
()
check_function_exists
(
starpu_
set_
iteration STARPU_
SET_
ITERATION_FOUND
)
if
(
STARPU_
SET_
ITERATION_FOUND
)
set
(
HAVE_STARPU_
SET_
ITERATION 1
)
message
(
"--
${
Blue
}
Add definition HAVE_STARPU_
SET_
ITERATION
${
ColourReset
}
"
)
check_function_exists
(
starpu_iteration
_push
STARPU_ITERATION_
PUSH_
FOUND
)
if
(
STARPU_ITERATION_
PUSH_
FOUND
)
set
(
HAVE_STARPU_ITERATION
_PUSH
1
)
message
(
"--
${
Blue
}
Add definition HAVE_STARPU_ITERATION
_PUSH
${
ColourReset
}
"
)
endif
()
check_function_exists
(
starpu_data_wont_use STARPU_DATA_WONT_USE_FOUND
)
if
(
STARPU_DATA_WONT_USE_FOUND
)
...
...
compute/pzgelqf.c
View file @
a3481c8e
...
...
@@ -115,7 +115,7 @@ void morse_pzgelqf(MORSE_desc_t *A, MORSE_desc_t *T,
#endif
for
(
k
=
0
;
k
<
minMNT
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -173,6 +173,8 @@ void morse_pzgelqf(MORSE_desc_t *A, MORSE_desc_t *T,
T
(
k
,
n
),
T
->
mb
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzgelqfrh.c
View file @
a3481c8e
...
...
@@ -119,7 +119,7 @@ void morse_pzgelqfrh(MORSE_desc_t *A, MORSE_desc_t *T, int BS,
#endif
for
(
k
=
0
;
k
<
chameleon_min
(
A
->
mt
,
A
->
nt
);
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -205,6 +205,8 @@ void morse_pzgelqfrh(MORSE_desc_t *A, MORSE_desc_t *T, int BS,
}
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzgeqrf.c
View file @
a3481c8e
...
...
@@ -110,7 +110,7 @@ void morse_pzgeqrf(MORSE_desc_t *A, MORSE_desc_t *T,
#endif
for
(
k
=
0
;
k
<
minMNT
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -167,6 +167,8 @@ void morse_pzgeqrf(MORSE_desc_t *A, MORSE_desc_t *T,
T
(
m
,
k
),
T
->
mb
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzgeqrfrh.c
View file @
a3481c8e
...
...
@@ -118,7 +118,7 @@ void morse_pzgeqrfrh(MORSE_desc_t *A, MORSE_desc_t *T, int BS,
K
=
chameleon_min
(
A
->
mt
,
A
->
nt
);
for
(
k
=
0
;
k
<
K
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
for
(
M
=
k
;
M
<
A
->
mt
;
M
+=
BS
)
{
...
...
@@ -204,6 +204,8 @@ void morse_pzgeqrfrh(MORSE_desc_t *A, MORSE_desc_t *T, int BS,
}
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzgetrf_incpiv.c
View file @
a3481c8e
...
...
@@ -78,7 +78,7 @@ void morse_pzgetrf_incpiv(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV,
morse_zdesc_alloc_diag
(
*
DIAG
,
A
->
mb
,
A
->
nb
,
chameleon_min
(
A
->
m
,
A
->
n
),
A
->
nb
,
0
,
0
,
chameleon_min
(
A
->
m
,
A
->
n
),
A
->
nb
,
A
->
p
,
A
->
q
);
for
(
k
=
0
;
k
<
minMNT
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -135,6 +135,8 @@ void morse_pzgetrf_incpiv(MORSE_desc_t *A, MORSE_desc_t *L, int *IPIV,
IPIV
(
m
,
k
));
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzgetrf_nopiv.c
View file @
a3481c8e
...
...
@@ -54,7 +54,7 @@ void morse_pzgetrf_nopiv(MORSE_desc_t *A,
ib
=
MORSE_IB
;
for
(
k
=
0
;
k
<
chameleon_min
(
A
->
mt
,
A
->
nt
);
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -100,6 +100,8 @@ void morse_pzgetrf_nopiv(MORSE_desc_t *A,
zone
,
A
(
m
,
n
),
ldam
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_finalize
(
&
options
,
morse
);
MORSE_TASK_dataflush_all
();
...
...
compute/pzhetrd_he2hb.c
View file @
a3481c8e
...
...
@@ -137,7 +137,7 @@ void morse_pzhetrd_he2hb(MORSE_enum uplo,
if
(
uplo
==
MorseLower
)
{
for
(
k
=
0
;
k
<
A
->
nt
-
1
;
k
++
){
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
+
1
==
A
->
mt
-
1
?
A
->
m
-
(
k
+
1
)
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -288,11 +288,13 @@ void morse_pzhetrd_he2hb(MORSE_enum uplo,
T
(
m
,
k
),
T
->
mb
);
options
.
priority
=
0
;
}
RUNTIME_iteration_pop
(
morse
);
}
}
else
{
for
(
k
=
0
;
k
<
A
->
nt
-
1
;
k
++
){
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
+
1
==
A
->
nt
-
1
?
A
->
n
-
(
k
+
1
)
*
A
->
nb
:
A
->
nb
;
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
...
...
@@ -440,6 +442,8 @@ void morse_pzhetrd_he2hb(MORSE_enum uplo,
T
(
k
,
n
),
T
->
mb
);
}
options
.
priority
=
0
;
RUNTIME_iteration_pop
(
morse
);
}
}
...
...
compute/pzpotrf.c
View file @
a3481c8e
...
...
@@ -71,7 +71,7 @@ void morse_pzpotrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
if
(
uplo
==
MorseLower
)
{
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -123,6 +123,7 @@ void morse_pzpotrf(MORSE_enum uplo, MORSE_desc_t *A,
}
MORSE_TASK_dataflush
(
&
options
,
A
(
n
,
k
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
/*
...
...
@@ -130,7 +131,7 @@ void morse_pzpotrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
else
{
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -181,6 +182,8 @@ void morse_pzpotrf(MORSE_enum uplo, MORSE_desc_t *A,
}
MORSE_TASK_dataflush
(
&
options
,
A
(
k
,
m
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
...
...
compute/pzpotrimm.c
View file @
a3481c8e
...
...
@@ -76,7 +76,7 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
* ZPOTRF
*/
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -121,12 +121,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
}
MORSE_TASK_dataflush
(
&
options
,
A
(
n
,
k
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZTRTRI
*/
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
A
->
nt
+
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -169,12 +171,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
uplo
,
MorseNonUnit
,
tempkn
,
A
->
mb
,
A
(
k
,
k
),
ldak
,
A
->
nb
*
k
);
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZLAUUM
*/
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
2
*
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
2
*
A
->
nt
+
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -212,12 +216,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
&
options
,
uplo
,
tempkm
,
A
->
mb
,
A
(
k
,
k
),
ldak
);
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZSYMM Right / Lower
*/
for
(
k
=
0
;
k
<
C
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
3
*
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
3
*
A
->
nt
+
k
);
tempkn
=
k
==
C
->
nt
-
1
?
C
->
n
-
k
*
C
->
nb
:
C
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -267,6 +273,8 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
for
(
n
=
0
;
n
<=
k
;
n
++
)
{
MORSE_TASK_dataflush
(
&
options
,
A
(
k
,
n
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
/*
...
...
@@ -277,7 +285,7 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
* ZPOTRF
*/
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -322,12 +330,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
}
MORSE_TASK_dataflush
(
&
options
,
A
(
k
,
m
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZTRTRI
*/
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
A
->
nt
+
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -370,12 +380,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
uplo
,
MorseNonUnit
,
tempkm
,
A
->
mb
,
A
(
k
,
k
),
ldak
,
A
->
mb
*
k
);
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZLAUUM
*/
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
2
*
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
2
*
A
->
nt
+
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -415,12 +427,14 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
&
options
,
uplo
,
tempkn
,
A
->
mb
,
A
(
k
,
k
),
ldak
);
RUNTIME_iteration_pop
(
morse
);
}
/*
* ZSYMM Right / Upper
*/
for
(
k
=
0
;
k
<
C
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
3
*
A
->
nt
+
k
);
RUNTIME_iteration
_push
(
morse
,
3
*
A
->
nt
+
k
);
tempkn
=
k
==
C
->
nt
-
1
?
C
->
n
-
k
*
C
->
nb
:
C
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -470,6 +484,8 @@ void morse_pzpotrimm(MORSE_enum uplo, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_de
for
(
m
=
0
;
m
<=
k
;
m
++
)
{
MORSE_TASK_dataflush
(
&
options
,
A
(
m
,
k
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
...
...
compute/pzsytrf.c
View file @
a3481c8e
...
...
@@ -60,7 +60,7 @@ void morse_pzsytrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
if
(
uplo
==
MorseLower
)
{
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -105,6 +105,8 @@ void morse_pzsytrf(MORSE_enum uplo, MORSE_desc_t *A,
}
MORSE_TASK_dataflush
(
&
options
,
A
(
n
,
k
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
/*
...
...
@@ -112,7 +114,7 @@ void morse_pzsytrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
else
{
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -157,6 +159,8 @@ void morse_pzsytrf(MORSE_enum uplo, MORSE_desc_t *A,
}
MORSE_TASK_dataflush
(
&
options
,
A
(
k
,
m
)
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
...
...
compute/pztpgqrt.c
View file @
a3481c8e
...
...
@@ -101,7 +101,7 @@ void morse_pztpgqrt( int L,
#endif
for
(
k
=
V1
->
nt
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
V1
->
mt
-
1
?
V1
->
m
-
k
*
V1
->
mb
:
V1
->
mb
;
tempkk
=
k
==
V1
->
nt
-
1
?
V1
->
n
-
k
*
V1
->
nb
:
V1
->
nb
;
...
...
@@ -171,6 +171,8 @@ void morse_pztpgqrt( int L,
T1
(
k
,
k
),
T1
->
mb
,
Q1
(
k
,
n
),
ldqk
);
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
...
...
compute/pztpqrt.c
View file @
a3481c8e
...
...
@@ -91,7 +91,7 @@ void morse_pztpqrt( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T,
RUNTIME_options_ws_alloc
(
&
options
,
ws_worker
,
ws_host
);
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -122,6 +122,8 @@ void morse_pztpqrt( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T,
}
maxmt
=
chameleon_min
(
B
->
mt
,
maxmt
+
1
);
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pztrtri.c
View file @
a3481c8e
...
...
@@ -56,7 +56,7 @@ void morse_pztrtri(MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A,
*/
if
(
uplo
==
MorseLower
)
{
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -99,6 +99,8 @@ void morse_pztrtri(MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A,
uplo
,
diag
,
tempkn
,
A
->
mb
,
A
(
k
,
k
),
ldak
,
A
->
nb
*
k
);
RUNTIME_iteration_pop
(
morse
);
}
}
/*
...
...
@@ -106,7 +108,7 @@ void morse_pztrtri(MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A,
*/
else
{
for
(
k
=
0
;
k
<
A
->
mt
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -149,6 +151,8 @@ void morse_pztrtri(MORSE_enum uplo, MORSE_enum diag, MORSE_desc_t *A,
uplo
,
diag
,
tempkm
,
A
->
mb
,
A
(
k
,
k
),
ldak
,
A
->
mb
*
k
);
RUNTIME_iteration_pop
(
morse
);
}
}
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzunglq.c
View file @
a3481c8e
...
...
@@ -98,7 +98,7 @@ void morse_pzunglq(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T,
#endif
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempAkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempAkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -145,6 +145,8 @@ void morse_pzunglq(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T,
T
(
k
,
k
),
T
->
mb
,
Q
(
m
,
k
),
ldqm
);
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzunglqrh.c
View file @
a3481c8e
...
...
@@ -99,7 +99,7 @@ void morse_pzunglqrh(MORSE_desc_t *A, MORSE_desc_t *Q,
K
=
chameleon_min
(
A
->
mt
,
A
->
nt
);
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -171,6 +171,8 @@ void morse_pzunglqrh(MORSE_desc_t *A, MORSE_desc_t *Q,
Q
(
m
,
N
),
ldqm
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzungqr.c
View file @
a3481c8e
...
...
@@ -98,7 +98,7 @@ void morse_pzungqr(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T,
#endif
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempAkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
tempAkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -146,6 +146,8 @@ void morse_pzungqr(MORSE_desc_t *A, MORSE_desc_t *Q, MORSE_desc_t *T,
T
(
k
,
k
),
T
->
mb
,
Q
(
k
,
n
),
ldqk
);
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzungqrrh.c
View file @
a3481c8e
...
...
@@ -101,7 +101,7 @@ void morse_pzungqrrh(MORSE_desc_t *A, MORSE_desc_t *Q,
K
=
chameleon_min
(
A
->
mt
,
A
->
nt
);
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
lastRD
=
0
;
...
...
@@ -176,6 +176,8 @@ void morse_pzungqrrh(MORSE_desc_t *A, MORSE_desc_t *Q,
Q
(
M
,
n
),
ldqM
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
RUNTIME_options_ws_free
(
&
options
);
RUNTIME_options_finalize
(
&
options
,
morse
);
...
...
compute/pzunmlq.c
View file @
a3481c8e
...
...
@@ -105,7 +105,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseNoTrans
*/
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -150,6 +150,8 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
T
(
k
,
m
),
T
->
mb
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
}
else
{
...
...
@@ -157,7 +159,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseConjTrans
*/
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -202,6 +204,8 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
T
(
k
,
k
),
T
->
mb
,
B
(
k
,
n
),
ldbk
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
}
...
...
@@ -211,7 +215,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseNoTrans
*/
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -256,6 +260,8 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
T
(
k
,
k
),
T
->
mb
,
B
(
m
,
k
),
ldbm
);
}
RUNTIME_iteration_pop
(
morse
);
}
}
else
{
...
...
@@ -263,7 +269,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseConjTrans
*/
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
mb
:
A
->
mb
;
...
...
@@ -308,6 +314,8 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
T
(
k
,
n
),
T
->
mb
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
}
}
...
...
compute/pzunmlqrh.c
View file @
a3481c8e
...
...
@@ -105,7 +105,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseNoTrans
*/
for
(
k
=
0
;
k
<
K
;
k
++
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -174,13 +174,15 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
}
}
}
RUNTIME_iteration_pop
(
morse
);
}
}
else
{
/*
* MorseLeft / MorseConjTrans
*/
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -252,6 +254,8 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
B
(
N
,
n
),
ldbN
);
}
}
RUNTIME_iteration_pop
(
morse
);
}
}
...
...
@@ -261,7 +265,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseNoTrans
*/
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_
set_
iteration
(
morse
,
k
);
RUNTIME_iteration
_push
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -332,13 +336,15 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
B
(
m
,
N
),
ldbm
);