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
7289434f
Commit
7289434f
authored
Mar 02, 2017
by
THIBAULT Samuel
Browse files
fix build
parent
07d767f2
Changes
21
Hide whitespace changes
Inline
Side-by-side
compute/pzgelqf.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzgelqfrh.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzgeqrf.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzgeqrfrh.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
for
(
M
=
k
;
M
<
A
->
mt
;
M
+=
BS
)
{
...
...
compute/pzgetrf_incpiv.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzgetrf_nopiv.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzhetrd_he2hb.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
@@ -292,7 +292,7 @@ void morse_pzhetrd_he2hb(MORSE_enum uplo,
}
else
{
for
(
k
=
0
;
k
<
A
->
nt
-
1
;
k
++
){
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzpotrf.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -130,7 +130,7 @@ void morse_pzpotrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
else
{
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzpotrimm.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -277,7 +277,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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzsytrf.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -112,7 +112,7 @@ void morse_pzsytrf(MORSE_enum uplo, MORSE_desc_t *A,
*/
else
{
for
(
k
=
0
;
k
<
A
->
nt
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pztpgqrt.c
View file @
7289434f
...
...
@@ -101,7 +101,7 @@ void morse_pztpgqrt( int L,
#endif
for
(
k
=
V1
->
nt
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pztpqrt.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pztrtri.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -106,7 +106,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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzunglq.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzunglqrh.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzungqr.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
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
;
...
...
compute/pzungqrrh.c
View file @
7289434f
...
...
@@ -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
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
A
->
nt
-
1
?
A
->
n
-
k
*
A
->
nb
:
A
->
nb
;
lastRD
=
0
;
...
...
compute/pzunmlq.c
View file @
7289434f
...
...
@@ -105,7 +105,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseNoTrans
*/
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -157,7 +157,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseConjTrans
*/
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -211,7 +211,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseNoTrans
*/
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -263,7 +263,7 @@ void morse_pzunmlq(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseConjTrans
*/
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
mb
:
A
->
mb
;
...
...
compute/pzunmlqrh.c
View file @
7289434f
...
...
@@ -105,7 +105,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseNoTrans
*/
for
(
k
=
0
;
k
<
K
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -180,7 +180,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseConjTrans
*/
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -261,7 +261,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseNoTrans
*/
for
(
k
=
K
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
@@ -338,7 +338,7 @@ void morse_pzunmlqrh(MORSE_enum side, MORSE_enum trans,
* MorseRight / MorseConjTrans
*/
for
(
k
=
0
;
k
<
K
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
A
->
mt
-
1
?
A
->
m
-
k
*
A
->
mb
:
A
->
mb
;
ldak
=
BLKLDD
(
A
,
k
);
...
...
compute/pzunmqr.c
View file @
7289434f
...
...
@@ -105,7 +105,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
* MorseLeft / MorseConjTrans
*/
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -158,7 +158,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
*/
else
{
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkm
=
k
==
B
->
mt
-
1
?
B
->
m
-
k
*
B
->
mb
:
B
->
mb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -213,7 +213,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
else
{
if
(
trans
==
MorseConjTrans
)
{
for
(
k
=
minMT
-
1
;
k
>=
0
;
k
--
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
@@ -267,7 +267,7 @@ void morse_pzunmqr(MORSE_enum side, MORSE_enum trans,
*/
else
{
for
(
k
=
0
;
k
<
minMT
;
k
++
)
{
RUNTIME_set_iteration
(
k
);
RUNTIME_set_iteration
(
morse
,
k
);
tempkn
=
k
==
B
->
nt
-
1
?
B
->
n
-
k
*
B
->
nb
:
B
->
nb
;
tempkmin
=
k
==
minMT
-
1
?
minM
-
k
*
A
->
nb
:
A
->
nb
;
...
...
Prev
1
2
Next
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment