Add GEMM case for HCHAMELEON

include/main.h

@@ -233,6 +233,7 @@ enum algo {
   _gemvHLIBPRO,
   _solveHLIBPRO,
   _gemvHCHAMELEON,
+  _gemmHCHAMELEON,
   _solveHCHAMELEON
 } ;
 

src/main.c

@@ -43,6 +43,10 @@ int main(int argc, char **argv) {
     testedAlgo = _solveHLIBPRO ;
     printf("Testing : solve HLIBPRO.\n") ;
   }
+  if (MpfArgHasName(&argc, argv, 1, "-gemmhchameleon") > 0) {
+    testedAlgo = _gemmHCHAMELEON ;
+    printf("Testing : gemm H-CHAMELEON.\n") ;
+  }
   if (MpfArgHasName(&argc, argv, 1, "-gemvhchameleon") > 0) {
     testedAlgo = _gemvHCHAMELEON ;
     printf("Testing : gemv H-CHAMELEON.\n") ;
@@ -66,6 +70,7 @@ int main(int argc, char **argv) {
     case _solveHLIBPRO:
       break;
     case _gemvHCHAMELEON:
+    case _gemmHCHAMELEON:
     case _solveHCHAMELEON: // Add "--chameleon" in the command line if it is not yet there
       if (!MpfArgHasName(&argc, argv, 0, "--chameleon")) argv[argc++]="--chameleon";
       break;
@@ -96,6 +101,7 @@ int main(int argc, char **argv) {
       ierr=initHLIBPRO(&argc, &argv) ; CHKERRQ(ierr) ;
       break;
     case _gemvHCHAMELEON:
+    case _gemmHCHAMELEON:
     case _solveHCHAMELEON:
       mpf_hmat_settings.engine = mpf_hmat_seq;
       ierr = init_hmat_interface(); CHKERRQ(ierr);
@@ -301,6 +307,7 @@ int main(int argc, char **argv) {
       ierr = testHLIBPRO(); CHKERRQ(ierr);
       break;
     case _gemvHCHAMELEON:
+    case _gemmHCHAMELEON:
     case _solveHCHAMELEON:
       ierr = testHCHAMELEON(); CHKERRQ(ierr);
       break;

src/testHCHAMELEON.c

@@ -110,6 +110,38 @@ int testHCHAMELEON(void) {
   /* Switch on the various H-CHAMELEON tests */
   /* ===================================================================================== */
   switch(testedAlgo) {
+    case _gemmHCHAMELEON:
+	{
+	    HCHAM_desc_t hdescB;
+	    HCHAMELEON_generate_matrix( chameleonType, NB, dims, &hdescB, hi );
+	    HCHAM_desc_t hdescC;
+	    HCHAMELEON_generate_matrix( chameleonType, NB, dims, &hdescC, hi );
+
+	    /* ===================================================================================== */
+	    /* Appelle le produit mat.vec : A.rhs -> solHMAT */
+	    /* ===================================================================================== */
+	    Mpf_printf(MPI_COMM_WORLD, "\n**** Computing HMATxHMAT product...\n") ;
+
+	    temps_initial = getTime ();
+
+	    /* Compute the mat.vec product A.rhs -> solCHAM */
+	    if ( generate_traces ) {
+		CHAMELEON_Enable( CHAMELEON_PROFILING_MODE );
+	    }
+	    CHAMELEON_gemm_Tile( hdescA.super, hdescB.super, hdescC.super );
+	    if ( generate_traces ) {
+		CHAMELEON_Disable( CHAMELEON_PROFILING_MODE );
+	    }
+
+	    temps_final = getTime ();
+	    temps_cpu = temps_final - temps_initial;
+	    Mpf_printf( MPI_COMM_WORLD,"<PERFTESTS> TpsCpuGEMM = %f \n", temps_cpu );
+
+	    HCHAMELEON_destroy_matrix( &hdescB );
+	    HCHAMELEON_destroy_matrix( &hdescC );
+	}
+	break;
+
     case _gemvHCHAMELEON:
 
       /* ===================================================================================== */