Commit a2c4f527 authored by PRUVOST Florent's avatar PRUVOST Florent

memset T to 0 for parfb on gpu

parent bb558c25
......@@ -271,6 +271,7 @@ int testing_zgels(int argc, char **argv)
}
MORSE_Alloc_Workspace_zgels(M, N, &T);
memset(T->mat, 0, (T->llm*T->lln)*sizeof(MORSE_Complex64_t));
eps = BLAS_dfpinfo( blas_eps );
/*----------------------------------------------------------
......
......@@ -29,7 +29,7 @@ RunTest(int *iparam, double *dparam, morse_time_t *t_)
{
MORSE_desc_t *T;
PASTE_CODE_IPARAM_LOCALS( iparam );
if ( M != N ) {
fprintf(stderr, "This timing works only with M == N\n");
return -1;
......@@ -46,6 +46,7 @@ RunTest(int *iparam, double *dparam, morse_time_t *t_)
MORSE_zplrnt( M, NRHS, x, LDB, 5673 );
MORSE_Alloc_Workspace_zgels(M, N, &T);
memset(T->mat, 0, (T->llm*T->lln)*sizeof(MorseComplexDouble));
/* Save A and b */
if (check) {
......
......@@ -46,6 +46,7 @@ RunTest(int *iparam, double *dparam, morse_time_t *t_)
/* Allocate Workspace */
MORSE_Alloc_Workspace_zgels_Tile(M, N, &descT);
memset(descT->mat, 0, (descT->llm*descT->lln)*sizeof(MorseComplexDouble));
/* Save A and B for check */
if (check == 1){
......
......@@ -43,6 +43,7 @@ RunTest(int *iparam, double *dparam, morse_time_t *t_)
/* Allocate Workspace */
MORSE_Alloc_Workspace_zgels(M, N, &T);
memset(T->mat, 0, (T->llm*T->lln)*sizeof(MorseComplexDouble));
/* Save AT in lapack layout for check */
PASTE_CODE_ALLOCATE_COPY( Acpy, check, MORSE_Complex64_t, A, LDA, N );
......@@ -50,14 +51,14 @@ RunTest(int *iparam, double *dparam, morse_time_t *t_)
START_TIMING();
MORSE_zgeqrf( M, N, A, LDA, T );
STOP_TIMING();
/* Check the solution */
if ( check )
{
PASTE_CODE_ALLOCATE_MATRIX( X, 1, MORSE_Complex64_t, LDB, NRHS );
MORSE_zplrnt( N, NRHS, X, LDB, 5673 );
PASTE_CODE_ALLOCATE_COPY( B, 1, MORSE_Complex64_t, X, LDB, NRHS );
MORSE_zgeqrs(M, N, NRHS, A, LDA, T, X, LDB);
dparam[IPARAM_RES] = z_check_solution(M, N, NRHS, Acpy, LDA, B, X, LDB,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment