Mentions légales du service

Skip to content
Snippets Groups Projects

Hotfix/sonarqube

Merged Hotfix/sonarqube
faverge/chameleon:hotfix/sonarqube into master
Merged Mathieu Faverge requested to merge faverge/chameleon:hotfix/sonarqube into master
Compare
and
7 files
+ 80
− 81
Compare changes
  • Side-by-side
  • Inline
Files
7
compute/pzgebrd_ge2gb.c
+ 65
− 65
@@ -20,86 +20,86 @@
@@ -20,86 +20,86 @@
*/
*/
#include "control/common.h"
#include "control/common.h"
void morse_pzgebrd_ge2gb(MORSE_desc_t A, MORSE_desc_t T, MORSE_desc_t D,
void morse_pzgebrd_ge2gb(MORSE_desc_t *A, MORSE_desc_t *T, MORSE_desc_t *D,
MORSE_sequence_t *sequence, MORSE_request_t *request)
MORSE_sequence_t *sequence, MORSE_request_t *request)
{
{
int k;
int k;
int tempkm, tempkn;
int tempkm, tempkn;
if (A.m >= A.n){
MORSE_desc_t *A1, *A2, *T1, *D1 = NULL;
for (k = 0; k < A.nt; 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;
morse_pzgeqrf(
if (A->m >= A->n){
morse_desc_submatrix(&A, k*A.mb, k*A.nb, A.m-k*A.mb, tempkn),
for (k = 0; k < A->nt; k++) {
morse_desc_submatrix(&T, k*T.mb, k*T.nb, T.m-k*T.mb, tempkn),
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
morse_desc_submatrix(&D, k*T.mb, k*T.nb, T.m-k*T.mb, tempkn),
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
sequence, request);
morse_pzunmqr(
A1 = morse_desc_submatrix(A, k*A->mb, k*A->nb, A->m-k*A->mb, tempkn);
MorseLeft,
A2 = morse_desc_submatrix(A, k*A->mb, (k+1)*A->nb, A->m-k*A->mb, A->n-(k+1)*A->nb);
MorseConjTrans,
T1 = morse_desc_submatrix(T, k*T->mb, k*T->nb, T->m-k*T->mb, T->nb );
morse_desc_submatrix(&A, k*A.mb, k*A.nb, A.m-k*A.mb, tempkn),
if ( D != NULL ) {
morse_desc_submatrix(&A, k*A.mb, (k+1)*A.nb, A.m-k*A.mb, A.n-(k+1)*A.nb),
D1 = morse_desc_submatrix(D, k*D->mb, k*D->nb, D->m-k*D->mb, tempkn);
morse_desc_submatrix(&T, k*T.mb, k*T.nb, T.m-k*T.mb, tempkn),
}
morse_desc_submatrix(&D, k*T.mb, k*T.nb, T.m-k*T.mb, tempkn),
sequence, request);
if (k+1 < A.nt){
morse_pzgeqrf( A1, T1, D1,
tempkn = k+1 == A.nt-1 ? A.n-(k+1)*A.nb : A.nb;
sequence, request);
morse_pzgelqf(
morse_pzunmqr( MorseLeft, MorseConjTrans,
morse_desc_submatrix(&A, k*A.mb, (k+1)*A.nb, tempkm, A.n-(k+1)*A.nb),
A1, A2, T1, D1,
morse_desc_submatrix(&T, k*T.mb, (k+1)*T.nb, T.mb, T.n-(k+1)*T.nb),
sequence, request);
morse_desc_submatrix(&D, k*T.mb, (k+1)*T.nb, T.mb, T.n-(k+1)*T.nb),
sequence, request);
morse_pzunmlq(
if (k+1 < A->nt){
MorseRight, MorseConjTrans,
tempkn = k+1 == A->nt-1 ? A->n-(k+1)*A->nb : A->nb;
morse_desc_submatrix(&A, k*A.mb, (k+1)*A.nb, tempkm, A.n-(k+1)*A.nb),
morse_desc_submatrix(&A, (k+1)*A.mb, (k+1)*A.nb, A.m-(k+1)*A.mb, A.n-(k+1)*A.nb),
A1 = morse_desc_submatrix(A, k*A->mb, (k+1)*A->nb, tempkm, A->n-(k+1)*A->nb);
morse_desc_submatrix(&T, k*T.mb, (k+1)*T.nb, T.mb, T.n-(k+1)*T.nb),
A2 = morse_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
morse_desc_submatrix(&D, k*T.mb, (k+1)*T.nb, T.mb, T.n-(k+1)*T.nb),
T1 = morse_desc_submatrix(T, k*T->mb, (k+1)*T->nb, T->mb, T->n-(k+1)*T->nb);
sequence, request);
if ( D != NULL ) {
}
D1 = morse_desc_submatrix(D, k*D->mb, (k+1)*D->nb, tempkm, D->n-(k+1)*D->nb);
}
}
 
 
morse_pzgelqf( A1, T1, D1,
 
sequence, request);
 
 
morse_pzunmlq( MorseRight, MorseConjTrans,
 
A1, A2, T1, D1,
 
sequence, request);
 
}
 
}
}
}
else{
else{
for (k = 0; k < A.mt; k++) {
for (k = 0; k < A->mt; k++) {
tempkm = k == A.mt-1 ? A.m-k*A.mb : A.mb;
tempkm = k == A->mt-1 ? A->m-k*A->mb : A->mb;
tempkn = k == A.nt-1 ? A.n-k*A.nb : A.nb;
tempkn = k == A->nt-1 ? A->n-k*A->nb : A->nb;
 
 
A1 = morse_desc_submatrix(A, k*A->mb, k*A->nb, tempkm, A->n-k*A->nb);
 
A2 = morse_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, A->n-k*A->nb);
 
T1 = morse_desc_submatrix(T, k*T->mb, k*T->nb, T->mb, T->n-k*T->nb);
 
if ( D != NULL ) {
 
D1 = morse_desc_submatrix(D, k*D->mb, k*D->nb, tempkm, D->n-k*D->nb);
 
}
 
morse_pzgelqf( A1, T1, D1,
 
sequence, request);
morse_pzgelqf(
morse_pzunmlq( MorseRight, MorseConjTrans,
morse_desc_submatrix(&A, k*A.mb, k*A.nb, tempkm, A.n-k*A.nb),
A1, A2, T1, D1,
morse_desc_submatrix(&T, k*T.mb, k*T.nb, T.mb, T.n-k*T.nb),
sequence, request);
morse_desc_submatrix(&D, k*T.mb, k*T.nb, T.mb, T.n-k*T.nb),
sequence, request);
morse_pzunmlq(
if (k+1 < A->mt){
MorseRight, MorseConjTrans,
tempkm = k+1 == A->mt-1 ? A->m-(k+1)*A->mb : A->mb;
morse_desc_submatrix(&A, k*A.mb, k*A.nb, tempkm, A.n-k*A.nb),
morse_desc_submatrix(&A, (k+1)*A.mb, k*A.nb, A.m-(k+1)*A.mb, A.n-k*A.nb),
morse_desc_submatrix(&T, k*T.mb, k*T.nb, T.mb, T.n-k*T.nb),
morse_desc_submatrix(&D, k*T.mb, k*T.nb, T.mb, T.n-k*T.nb),
sequence, request);
if (k+1 < A.mt){
A1 = morse_desc_submatrix(A, (k+1)*A->mb, k*A->nb, A->m-(k+1)*A->mb, tempkn);
tempkm = k+1 == A.mt-1 ? A.m-(k+1)*A.mb : A.mb;
A2 = morse_desc_submatrix(A, (k+1)*A->mb, (k+1)*A->nb, A->m-(k+1)*A->mb, A->n-(k+1)*A->nb);
 
T1 = morse_desc_submatrix(T, (k+1)*T->mb, k*T->nb, T->m-(k+1)*T->mb, T->nb );
 
if ( D != NULL ) {
 
D1 = morse_desc_submatrix(D, (k+1)*D->mb, k*D->nb, D->m-(k+1)*D->mb, tempkn);
 
}
morse_pzgeqrf(
morse_pzgeqrf( A1, T1, D1,
morse_desc_submatrix(&A, (k+1)*A.mb, k*A.nb, A.m-(k+1)*A.mb, tempkn),
sequence, request);
morse_desc_submatrix(&T, (k+1)*T.mb, k*T.nb, T.m-(k+1)*T.mb, tempkn),
morse_desc_submatrix(&D, (k+1)*T.mb, k*T.nb, T.m-(k+1)*T.mb, tempkn),
sequence, request);
morse_pzunmqr(
morse_pzunmqr( MorseLeft, MorseConjTrans,
MorseLeft, MorseConjTrans,
A1, A2, T1, D1,
morse_desc_submatrix(&A, (k+1)*A.mb, k*A.nb, A.m-(k+1)*A.mb, tempkn),
sequence, request);
morse_desc_submatrix(&A, (k+1)*A.mb, (k+1)*A.nb, A.m-(k+1)*A.mb, A.n-(k+1)*A.nb),
}
morse_desc_submatrix(&T, (k+1)*T.mb, k*T.nb, T.m-(k+1)*T.mb, tempkn),
}
morse_desc_submatrix(&D, (k+1)*T.mb, k*T.nb, T.m-(k+1)*T.mb, tempkn),
sequence, request);
}
}
}
}
}
}