diff --git a/testing/testing_zhemm.c b/testing/testing_zhemm.c
index 9fc0006429ccc736f71de8b0d27cb726a76d014d..b89b2bb3ce1d32f322547ec4a91f500e9fa3a088 100644
--- a/testing/testing_zhemm.c
+++ b/testing/testing_zhemm.c
@@ -156,8 +156,8 @@ testing_zhemm_std( run_arg_list_t *args, int check )
/* Create the matrices */
A = malloc( LDA*Am*sizeof(CHAMELEON_Complex64_t) );
- B = malloc( LDB*N*sizeof(CHAMELEON_Complex64_t) );
- C = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
+ B = malloc( LDB*N *sizeof(CHAMELEON_Complex64_t) );
+ C = malloc( LDC*N *sizeof(CHAMELEON_Complex64_t) );
/* Fills the matrix with random values */
CHAMELEON_zplghe( bump, uplo, N, A, LDA, seedA );
@@ -175,8 +175,7 @@ testing_zhemm_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplrnt( M, N, Cinit, LDC, seedC );
- // hres += check_zsymm( args, ChamHermitian, side, uplo, alpha, descA, descB,
- // beta, descCinit, descC );
+ hres += check_zsymm_std( args, ChamHermitian, side, uplo, M, N, alpha, A, LDA, B, LDB, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_zher2k.c b/testing/testing_zher2k.c
index 1507b35e719c4c394215aaee6c362046d37158fd..79509cd0809c5fbf7493c5e36edeb7685c53bb3a 100644
--- a/testing/testing_zher2k.c
+++ b/testing/testing_zher2k.c
@@ -161,11 +161,11 @@ testing_zher2k_std( run_arg_list_t *args, int check )
/* Create the matrices */
A = malloc( LDA*An*sizeof(CHAMELEON_Complex64_t) );
B = malloc( LDB*An*sizeof(CHAMELEON_Complex64_t) );
- C = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
+ C = malloc( LDC*N *sizeof(CHAMELEON_Complex64_t) );
/* Fill the matrix with random values */
CHAMELEON_zplrnt( Am, An, B, LDA, seedA );
- CHAMELEON_zplrnt( K, An, B, LDB, seedB );
+ CHAMELEON_zplrnt( K, An, B, LDB, seedB );
CHAMELEON_zplghe( bump, uplo, N, C, LDC, seedC );
/* Calculate the product */
@@ -179,8 +179,7 @@ testing_zher2k_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplghe( bump, uplo, N, Cinit, LDC, seedC );
- // hres += check_zsyrk( args, ChamHermitian, uplo, trans, alpha, descA, descB,
- // beta, descCinit, descC );
+ hres += check_zsyrk_std( args, ChamHermitian, uplo, trans, N, K, alpha, A, B, LDA, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_zherk.c b/testing/testing_zherk.c
index 3f6a02f32af8b7e1b70b66d2c5b4fd90ae194ef2..83ade364b25c4e19b90d1ddad0e11a257bc44f10 100644
--- a/testing/testing_zherk.c
+++ b/testing/testing_zherk.c
@@ -152,7 +152,7 @@ testing_zherk_std( run_arg_list_t *args, int check )
/* Creates the matrices */
A = malloc( LDA*An*sizeof(CHAMELEON_Complex64_t) );
- C = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
+ C = malloc( LDC*N *sizeof(CHAMELEON_Complex64_t) );
/* Fills the matrix with random values */
CHAMELEON_zplrnt( Am, An, A, LDA, seedA );
@@ -169,8 +169,7 @@ testing_zherk_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplghe( bump, uplo, N, Cinit, LDC, seedC );
- // hres += check_zsyrk( args, ChamHermitian, uplo, trans, alpha, descA, NULL,
- // beta, descCinit, descC );
+ hres += check_zsyrk_std( args, ChamHermitian, uplo, trans, N, K, alpha, A, NULL, LDA, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_zsymm.c b/testing/testing_zsymm.c
index a8d8bf0c0a9668f62cda2de9c27fabe660a59da5..f87066a874863646f478a512ae5e527d82db09bb 100644
--- a/testing/testing_zsymm.c
+++ b/testing/testing_zsymm.c
@@ -175,8 +175,7 @@ testing_zsymm_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplrnt( M, N, Cinit, LDC, seedC );
- // hres += check_zsymm( args, ChamSymmetric, side, uplo, alpha, descA, descB,
- // beta, descCinit, descC );
+ hres += check_zsymm_std( args, ChamSymmetric, side, uplo, M, N, alpha, A, LDA, B, LDB, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_zsyr2k.c b/testing/testing_zsyr2k.c
index 49f2047908ca69f0f748fa760cb04d7d04e821c2..1db847448215cc36c36887e438c7d3c2bed9ec40 100644
--- a/testing/testing_zsyr2k.c
+++ b/testing/testing_zsyr2k.c
@@ -161,7 +161,7 @@ testing_zsyr2k_std( run_arg_list_t *args, int check )
/* Create the matrices */
A = malloc( LDA*An*sizeof(CHAMELEON_Complex64_t) );
B = malloc( LDB*An*sizeof(CHAMELEON_Complex64_t) );
- C = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
+ C = malloc( LDC*N *sizeof(CHAMELEON_Complex64_t) );
/* Fill the matrix with random values */
CHAMELEON_zplrnt( Am, An, A, LDA, seedA );
@@ -179,8 +179,7 @@ testing_zsyr2k_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplgsy( bump, uplo, N, Cinit, LDC, seedC );
- // hres += check_zsyrk( args, ChamSymmetric, uplo, trans, alpha, descA, descB,
- // beta, descCinit, descC );
+ hres += check_zsyrk_std( args, ChamSymmetric, uplo, trans, N, K, alpha, A, B, LDA, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_zsyrk.c b/testing/testing_zsyrk.c
index 164d8ff5fed4aab73a3ddb7cba6c05a2cb5894e3..87be0bd64561a4cff298d4ce36b3b3fb1be252b8 100644
--- a/testing/testing_zsyrk.c
+++ b/testing/testing_zsyrk.c
@@ -151,7 +151,7 @@ testing_zsyrk_std( run_arg_list_t *args, int check )
/* Creates the matrices */
A = malloc( LDA*An*sizeof(CHAMELEON_Complex64_t) );
- C = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
+ C = malloc( LDC*N *sizeof(CHAMELEON_Complex64_t) );
/* Fills the matrix with random values */
CHAMELEON_zplrnt( Am, An, A, LDA, seedA );
@@ -168,8 +168,7 @@ testing_zsyrk_std( run_arg_list_t *args, int check )
Cinit = malloc( LDC*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplgsy( bump, uplo, N, Cinit, LDC, seedC );
- // hres += check_zsyrk( args, ChamSymmetric, uplo, trans, alpha, descA, NULL,
- // beta, descCinit, descC );
+ hres += check_zsyrk_std( args, ChamSymmetric, uplo, trans, N, K, alpha, A, NULL, LDA, beta, Cinit, C, LDC );
free( Cinit );
}
diff --git a/testing/testing_ztrmm.c b/testing/testing_ztrmm.c
index dced32a2e51699cbea23e26bbaae0455835279ea..9809931fc2d5bf586bcabb1e9ed549e1b40ea6dc 100644
--- a/testing/testing_ztrmm.c
+++ b/testing/testing_ztrmm.c
@@ -39,16 +39,15 @@ testing_ztrmm_desc( run_arg_list_t *args, int check )
cham_uplo_t uplo = run_arg_get_uplo( args, "uplo", ChamUpper );
cham_diag_t diag = run_arg_get_diag( args, "diag", ChamNonUnit );
int N = run_arg_get_int( args, "N", 1000 );
- int K = run_arg_get_int( args, "K", N );
+ int M = run_arg_get_int( args, "M", N );
int LDA = run_arg_get_int( args, "LDA", N );
- int LDB = run_arg_get_int( args, "LDB", N );
+ int LDB = run_arg_get_int( args, "LDB", M );
CHAMELEON_Complex64_t alpha = testing_zalea();
int seedA = run_arg_get_int( args, "seedA", random() );
int seedB = run_arg_get_int( args, "seedB", random() );
int Q = parameters_compute_q( P );
/* Descriptors */
- int Bm, Bn;
CHAM_desc_t *descA, *descB, *descBinit;
alpha = run_arg_get_complex64( args, "alpha", alpha );
@@ -57,19 +56,19 @@ testing_ztrmm_desc( run_arg_list_t *args, int check )
/* Calculates the dimensions according to the side */
if ( side == ChamLeft ) {
- Bm = N;
- Bn = K;
+ LDA = M;
+ CHAMELEON_Desc_Create(
+ &descA, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDA, M, 0, 0, M, M, P, Q );
}
else {
- Bm = K;
- Bn = N;
+ LDA = N;
+ CHAMELEON_Desc_Create(
+ &descA, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDA, N, 0, 0, N, N, P, Q );
}
/* Creates the matrices */
CHAMELEON_Desc_Create(
- &descA, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDA, N, 0, 0, N, N, P, Q );
- CHAMELEON_Desc_Create(
- &descB, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDB, Bn, 0, 0, Bm, Bn, P, Q );
+ &descB, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDB, N, 0, 0, M, N, P, Q );
/* Fills the matrix with random values */
CHAMELEON_zplrnt_Tile( descA, seedA );
@@ -87,12 +86,12 @@ testing_ztrmm_desc( run_arg_list_t *args, int check )
hres = CHAMELEON_ztrmm_Tile( side, uplo, trans, diag, alpha, descA, descB );
}
test_data.hres = hres;
- testing_stop( &test_data, flops_ztrmm( side, N, K ) );
+ testing_stop( &test_data, flops_ztrmm( side, N, M ) );
/* Checks the solution */
if ( check ) {
CHAMELEON_Desc_Create(
- &descBinit, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDB, Bn, 0, 0, Bm, Bn, P, Q );
+ &descBinit, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDB, N, 0, 0, M, N, P, Q );
CHAMELEON_zplrnt_Tile( descBinit, seedB );
hres += check_ztrmm( args, CHECK_TRMM, side, uplo, trans, diag,
@@ -120,15 +119,15 @@ testing_ztrmm_std( run_arg_list_t *args, int check )
cham_uplo_t uplo = run_arg_get_uplo( args, "uplo", ChamUpper );
cham_diag_t diag = run_arg_get_diag( args, "diag", ChamNonUnit );
int N = run_arg_get_int( args, "N", 1000 );
- int K = run_arg_get_int( args, "K", N );
+ int M = run_arg_get_int( args, "M", N );
int LDA = run_arg_get_int( args, "LDA", N );
- int LDB = run_arg_get_int( args, "LDB", N );
+ int LDB = run_arg_get_int( args, "LDB", M );
CHAMELEON_Complex64_t alpha = testing_zalea();
int seedA = run_arg_get_int( args, "seedA", random() );
int seedB = run_arg_get_int( args, "seedB", random() );
/* Descriptors */
- int Bm, Bn;
+ int An;
CHAMELEON_Complex64_t *A, *B, *Binit;
alpha = run_arg_get_complex64( args, "alpha", alpha );
@@ -137,35 +136,35 @@ testing_ztrmm_std( run_arg_list_t *args, int check )
/* Calculates the dimensions according to the side */
if ( side == ChamLeft ) {
- Bm = N;
- Bn = K;
+ An = M;
+ LDA = M;
+ A = malloc( LDA*N *sizeof(CHAMELEON_Complex64_t) );
}
else {
- Bm = K;
- Bn = N;
+ An = N;
+ LDA = N;
+ A = malloc( LDA*N *sizeof(CHAMELEON_Complex64_t) );
}
/* Creates the matrices */
- A = malloc( LDA*N*sizeof(CHAMELEON_Complex64_t) );
- B = malloc( LDB*Bn*sizeof(CHAMELEON_Complex64_t) );
+ B = malloc( LDB*N*sizeof(CHAMELEON_Complex64_t) );
/* Fills the matrix with random values */
- CHAMELEON_zplrnt( K, N, A, LDA, seedA );
- CHAMELEON_zplrnt( Bm, Bn, B, LDB, seedB );
+ CHAMELEON_zplrnt( An, An, A, LDA, seedA );
+ CHAMELEON_zplrnt( M, N, B, LDB, seedB );
/* Calculates the product */
testing_start( &test_data );
- hres = CHAMELEON_ztrmm( side, uplo, trans, diag, N, K, alpha, A, LDA, B, LDB );
+ hres = CHAMELEON_ztrmm( side, uplo, trans, diag, N, M, alpha, A, LDA, B, LDB );
test_data.hres = hres;
- testing_stop( &test_data, flops_ztrmm( side, N, K ) );
+ testing_stop( &test_data, flops_ztrmm( side, N, M ) );
/* Checks the solution */
if ( check ) {
- Binit = malloc( LDB*Bn*sizeof(CHAMELEON_Complex64_t) );
- CHAMELEON_zplrnt( Bm, Bn, Binit, LDB, seedB );
+ Binit = malloc( LDB*N*sizeof(CHAMELEON_Complex64_t) );
+ CHAMELEON_zplrnt( M, N, Binit, LDB, seedB );
- // hres += check_ztrmm( args, CHECK_TRMM, side, uplo, trans, diag,
- // alpha, descA, descB, descBinit );
+ hres += check_ztrmm_std( args, CHECK_TRMM, side, uplo, trans, diag, N, M, alpha, A, LDA, B, Binit, LDB );
free( Binit );
}
@@ -177,8 +176,8 @@ testing_ztrmm_std( run_arg_list_t *args, int check )
}
testing_t test_ztrmm;
-const char *ztrmm_params[] = { "mtxfmt", "nb", "trans", "side", "uplo", "diag", "n",
- "k", "lda", "ldb", "alpha", "seedA", "seedB", NULL };
+const char *ztrmm_params[] = { "mtxfmt", "nb", "trans", "side", "uplo", "diag", "m",
+ "n", "lda", "ldb", "alpha", "seedA", "seedB", NULL };
const char *ztrmm_output[] = { NULL };
const char *ztrmm_outchk[] = { "RETURN", NULL };
diff --git a/testing/testing_ztrsm.c b/testing/testing_ztrsm.c
index 851df92b88960469fbcf3515ecbdd73f2b2ee037..89cd4f73354a488301ab2d5b266f7726b1626996 100644
--- a/testing/testing_ztrsm.c
+++ b/testing/testing_ztrsm.c
@@ -128,7 +128,7 @@ testing_ztrsm_std( run_arg_list_t *args, int check )
/* Creates the matrices */
A = malloc( LDA*Ak*sizeof(CHAMELEON_Complex64_t) );
- B = malloc( LDB*N*sizeof(CHAMELEON_Complex64_t) );
+ B = malloc( LDB*N *sizeof(CHAMELEON_Complex64_t) );
/* Fills the matrix with random values */
/* We bump a little bit the diagonal to make it stable */
@@ -146,8 +146,7 @@ testing_ztrsm_std( run_arg_list_t *args, int check )
Binit = malloc( LDB*N*sizeof(CHAMELEON_Complex64_t) );
CHAMELEON_zplrnt( M, N, Binit, LDB, seedB );
- // hres += check_ztrmm( args, CHECK_TRSM, side, uplo, trans, diag,
- // alpha, descA, descB, descBinit );
+ hres += check_ztrmm_std( args, CHECK_TRSM, side, uplo, trans, diag, M, N, alpha, A, LDA, B, Binit, LDB );
free( Binit );
}