diff --git a/testing/CMakeLists.txt b/testing/CMakeLists.txt
index da4470ac393513e11fa2757611f018de3b3db017..df0573a94ebccdb71f57e7b24c607de3e60b5d21 100644
--- a/testing/CMakeLists.txt
+++ b/testing/CMakeLists.txt
@@ -96,6 +96,7 @@ set(ZSRC
# testing_zgeqrs_hqr.c
# testing_zgelqs_hqr.c
testing_zgels_hqr.c
+ testing_zgepdf_qdwh.c
testing_zgepdf_qr.c
testing_zplrnk.c
)
diff --git a/testing/CTestLists.cmake b/testing/CTestLists.cmake
index 06c2a819333690a760d2089b3620ff2360a3d3f4..007d8013938501d3e90f079eef6bb60d4740d3e2 100644
--- a/testing/CTestLists.cmake
+++ b/testing/CTestLists.cmake
@@ -65,7 +65,7 @@ if (NOT CHAMELEON_SIMULATION)
gels
gels_hqr )
set( TESTS ${TESTS}
- genm2 gepdf_qr )
+ genm2 gepdf_qr gepdf_qdwh )
foreach(cat ${TEST_CATEGORIES})
foreach(gpus ${N_GPUS})
diff --git a/testing/input/gepdf_qdwh.in b/testing/input/gepdf_qdwh.in
new file mode 100644
index 0000000000000000000000000000000000000000..6fd3322cf90c4dbddab6f253d95fe7e2eb89ec41
--- /dev/null
+++ b/testing/input/gepdf_qdwh.in
@@ -0,0 +1,17 @@
+# You can enumerate each parameter's values as an explicit list separated by commas or by a range start:end[:step]
+# Not given parameters will receive default values
+
+# GEPDF_QDWH
+
+# nb: Tile size
+# ib: Inner tile size
+# m: Number of rows of the A matrix
+# n: Number of columns of the A matrix
+# lda: Leading dimension of the A matrix
+
+op = gepdf_qdwh
+nb = 8
+ib = 3
+m = 8, 32, 64
+n = 8, 16, 32
+lda = 79
diff --git a/testing/testing_zgepdf_qdwh.c b/testing/testing_zgepdf_qdwh.c
new file mode 100644
index 0000000000000000000000000000000000000000..78b5b04c57409296cd90a4d43bf0e168658a1789
--- /dev/null
+++ b/testing/testing_zgepdf_qdwh.c
@@ -0,0 +1,133 @@
+/**
+ *
+ * @file testing_zgepdf_qdwh.c
+ *
+ * @copyright 2020-2020 Bordeaux INP, CNRS (LaBRI UMR 5800), Inria,
+ * Univ. Bordeaux. All rights reserved.
+ * @copyright 2020-2020 KAUST. All rights reserved.
+ *
+ ***
+ *
+ * @brief Chameleon zgepdf_qdwh testing
+ *
+ * @version 1.0.0
+ * @author Mathieu Faverge
+ * @author Hatem Ltaief
+ * @date 2020-10-12
+ * @precisions normal z -> c d s
+ *
+ */
+#include <chameleon.h>
+#if !defined(CHAMELEON_SIMULATION)
+#include <coreblas/lapacke.h>
+#include <coreblas/cblas.h>
+#include <coreblas.h>
+#endif
+#include "testings.h"
+#include "testing_zcheck.h"
+#include <chameleon/flops.h>
+
+int
+testing_zgepdf_qdwh( run_arg_list_t *args, int check )
+{
+ int hres = 0;
+ CHAM_desc_t *descA, *descA0, *descH;
+ gepdf_info_t info;
+
+ /* Reads arguments */
+ intptr_t mtxfmt = parameters_getvalue_int( "mtxfmt" );
+ int nb = run_arg_get_int( args, "nb", 320 );
+ int ib = run_arg_get_int( args, "ib", 48 );
+ int P = parameters_getvalue_int( "P" );
+ int N = run_arg_get_int( args, "N", 1000 );
+ int M = run_arg_get_int( args, "M", N );
+ int LDA = run_arg_get_int( args, "LDA", M );
+ int LDB = run_arg_get_int( args, "LDB", N );
+ int seedA = run_arg_get_int( args, "seedA", random() );
+ int Q = parameters_compute_q( P );
+ double cond = run_arg_get_double( args, "cond", 1.e16 );
+ int mode = run_arg_get_int( args, "mode", 4 );
+ int runtime;
+ cham_fixdbl_t t, gflops;
+
+ CHAMELEON_Set( CHAMELEON_TILE_SIZE, nb );
+ CHAMELEON_Set( CHAMELEON_INNER_BLOCK_SIZE, ib );
+
+ CHAMELEON_Get( CHAMELEON_RUNTIME, &runtime );
+ if ( runtime == RUNTIME_SCHED_PARSEC ) {
+ if ( CHAMELEON_Comm_rank() == 0 ) {
+ fprintf( stderr, "SKIPPED: The QDWH polar decomposition is not supported with PaRSEC\n" );
+ }
+ return -1;
+ }
+
+ if ( N > M ) {
+ if ( CHAMELEON_Comm_rank() == 0 ) {
+ fprintf( stderr, "SKIPPED: The QDWH polar decomposition is performed only when M >= N\n" );
+ }
+ return -1;
+ }
+
+ /* Creates the matrix */
+ CHAMELEON_Desc_Create(
+ &descA, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDA, N, 0, 0, M, N, P, Q );
+ CHAMELEON_Desc_Create(
+ &descH, (void*)(-mtxfmt), ChamComplexDouble, nb, nb, nb * nb, LDB, N, 0, 0, N, N, P, Q );
+
+ /* Fills the matrix with random values */
+ hres = CHAMELEON_zlatms_Tile(
+ ChamDistUniform, seedA, ChamNonsymPosv, NULL, mode, cond, 1., descA );
+ if ( hres != 0 ) {
+ return hres;
+ }
+
+ if ( check ) {
+ descA0 = CHAMELEON_Desc_Copy( descA, CHAMELEON_MAT_ALLOC_GLOBAL );
+ CHAMELEON_zlacpy_Tile( ChamUpperLower, descA, descA0 );
+ }
+
+ /* Calculates the norm */
+ START_TIMING( t );
+ hres = CHAMELEON_zgepdf_qdwh_Tile( descA, descH, &info );
+ STOP_TIMING( t );
+ gflops = info.flops * 1.e-9 / t;
+ run_arg_add_fixdbl( args, "time", t );
+ run_arg_add_fixdbl( args, "gflops", ( hres == CHAMELEON_SUCCESS ) ? gflops : -1. );
+
+ /* Checks the solution */
+ hres = 0;
+ if ( check ) {
+ hres += check_zxxpd( args, descA0, descA, descH );
+ hres += check_zortho( args, descA );
+
+ CHAMELEON_Desc_Destroy( &descA0 );
+ }
+
+ CHAMELEON_Desc_Destroy( &descA );
+ CHAMELEON_Desc_Destroy( &descH );
+
+ return hres;
+}
+
+testing_t test_zgepdf_qdwh;
+const char *zgepdf_qdwh_params[] = { "mtxfmt", "nb", "ib", "m", "n", "lda", "ldb", "seedA", "cond", "mode", NULL };
+const char *zgepdf_qdwh_output[] = { NULL };
+const char *zgepdf_qdwh_outchk[] = { "||A||", "||A-fact(A)||", "||I-QQ'||", "RETURN", NULL };
+
+/**
+ * @brief Testing registration function
+ */
+void testing_zgepdf_qdwh_init( void ) __attribute__( ( constructor ) );
+void
+testing_zgepdf_qdwh_init( void )
+{
+ test_zgepdf_qdwh.name = "zgepdf_qdwh";
+ test_zgepdf_qdwh.helper = "Polar decomposition factorization with QDWH algorithm";
+ test_zgepdf_qdwh.params = zgepdf_qdwh_params;
+ test_zgepdf_qdwh.output = zgepdf_qdwh_output;
+ test_zgepdf_qdwh.outchk = zgepdf_qdwh_outchk;
+ test_zgepdf_qdwh.fptr = testing_zgepdf_qdwh;
+ test_zgepdf_qdwh.next = NULL;
+
+ testing_register( &test_zgepdf_qdwh );
+}