  
 +#ifndef NO_GEOPM
 +  uint64_t R_ID;
-+  geopm_prof_region( "PHASE_MATVEC", 0, &R_ID);
++  geopm_prof_region( "PHASE_MATVEC", GEOPM_REGION_HINT_COMPUTE, &R_ID);
 +  geopm_prof_enter( R_ID);
 +#endif
    #pragma omp parallel for

 
+#ifndef NO_GEOPM
+  uint64_t R_ID;
+  geopm_prof_region( "PHASE_WAXPBY", 0, &R_ID);
+  geopm_prof_enter( R_ID);
+#endif
   if(beta == 0.0) {

 
+#ifndef NO_GEOPM
+  uint64_t R_ID;
+  geopm_prof_region( "PHASE_DOT", GEOPM_REGION_HINT_MEMORY, &R_ID);
+  geopm_prof_enter( R_ID);
+#endif
   #pragma omp parallel for reduction(+:result)
