idamnjanovic@51: /**************************************************************************
idamnjanovic@51:  *
idamnjanovic@51:  * File name: ompprof.c
idamnjanovic@51:  *
idamnjanovic@51:  * Ron Rubinstein
idamnjanovic@51:  * Computer Science Department
idamnjanovic@51:  * Technion, Haifa 32000 Israel
idamnjanovic@51:  * ronrubin@cs
idamnjanovic@51:  *
idamnjanovic@51:  * Last Updated: 11.4.2009
idamnjanovic@51:  *
idamnjanovic@51:  *************************************************************************/
idamnjanovic@51: 
idamnjanovic@51: 
idamnjanovic@51: #include "ompprof.h"
idamnjanovic@51: 
idamnjanovic@51: 
idamnjanovic@51: /* initialize profiling information */
idamnjanovic@51: 
idamnjanovic@51: void initprofdata(profdata *pd)
idamnjanovic@51: {
idamnjanovic@51:   pd->DtX_time = 0;
idamnjanovic@51:   pd->XtX_time = 0;
idamnjanovic@51:   pd->DtR_time = 0;
idamnjanovic@51:   pd->maxabs_time = 0;
idamnjanovic@51:   pd->DtD_time = 0;
idamnjanovic@51:   pd->Lchol_time = 0;
idamnjanovic@51:   pd->compcoef_time = 0;
idamnjanovic@51:   pd->update_DtR_time = 0;
idamnjanovic@51:   pd->update_resnorm_time = 0;
idamnjanovic@51:   pd->compres_time = 0;
idamnjanovic@51:   pd->indexsort_time = 0;
idamnjanovic@51:   
idamnjanovic@51:   pd->DtX_time_counted = 0;
idamnjanovic@51:   pd->XtX_time_counted = 0;
idamnjanovic@51:   pd->DtR_time_counted = 0;
idamnjanovic@51:   pd->DtD_time_counted = 0;
idamnjanovic@51:   pd->update_DtR_time_counted = 0;
idamnjanovic@51:   pd->resnorm_time_counted = 0;
idamnjanovic@51:   pd->compres_time_counted = 0;
idamnjanovic@51:   pd->indexsort_time_counted = 0;
idamnjanovic@51:   
idamnjanovic@51:   pd->prevtime = clock();
idamnjanovic@51: }
idamnjanovic@51: 
idamnjanovic@51: 
idamnjanovic@51: /* add elapsed time to profiling data according to specified computation */
idamnjanovic@51: 
idamnjanovic@51: void addproftime(profdata *pd, int comptype)
idamnjanovic@51: {
idamnjanovic@51:   switch(comptype) {
idamnjanovic@51:     case DtX_TIME:            pd->DtX_time            += clock()-pd->prevtime; pd->DtX_time_counted = 1; break;
idamnjanovic@51:     case XtX_TIME:            pd->XtX_time            += clock()-pd->prevtime; pd->XtX_time_counted = 1; break;
idamnjanovic@51:     case DtR_TIME:            pd->DtR_time            += clock()-pd->prevtime; pd->DtR_time_counted = 1; break;
idamnjanovic@51:     case DtD_TIME:            pd->DtD_time            += clock()-pd->prevtime; pd->DtD_time_counted = 1; break;
idamnjanovic@51:     case COMPRES_TIME:        pd->compres_time        += clock()-pd->prevtime; pd->compres_time_counted = 1; break;
idamnjanovic@51:     case UPDATE_DtR_TIME:     pd->update_DtR_time     += clock()-pd->prevtime; pd->update_DtR_time_counted = 1; break;
idamnjanovic@51:     case UPDATE_RESNORM_TIME: pd->update_resnorm_time += clock()-pd->prevtime; pd->resnorm_time_counted = 1; break;
idamnjanovic@51:     case INDEXSORT_TIME:      pd->indexsort_time      += clock()-pd->prevtime; pd->indexsort_time_counted = 1; break;
idamnjanovic@51:     case MAXABS_TIME:         pd->maxabs_time         += clock()-pd->prevtime; break;
idamnjanovic@51:     case LCHOL_TIME:          pd->Lchol_time          += clock()-pd->prevtime; break;
idamnjanovic@51:     case COMPCOEF_TIME:       pd->compcoef_time       += clock()-pd->prevtime; break;
idamnjanovic@51:   }
idamnjanovic@51:   pd->prevtime = clock();
idamnjanovic@51: }
idamnjanovic@51: 
idamnjanovic@51: 
idamnjanovic@51: /* print profiling info */
idamnjanovic@51: 
idamnjanovic@51: void printprofinfo(profdata *pd, int erroromp, int batchomp, int signum)
idamnjanovic@51: {
idamnjanovic@51:   clock_t tottime;
idamnjanovic@51:   
idamnjanovic@51:   tottime = pd->DtX_time + pd->XtX_time + pd->DtR_time + pd->DtD_time + pd->compres_time + pd->maxabs_time + 
idamnjanovic@51:             pd->Lchol_time + pd->compcoef_time + pd->update_DtR_time + pd->update_resnorm_time + pd->indexsort_time;
idamnjanovic@51:   
idamnjanovic@51:   mexPrintf("\n\n*****  Profiling information for %s  *****\n\n", erroromp? "OMP2" : "OMP");
idamnjanovic@51:   
idamnjanovic@51:   mexPrintf("OMP mode: %s\n\n", batchomp? "Batch-OMP" : "OMP-Cholesky");
idamnjanovic@51:   
idamnjanovic@51:   mexPrintf("Total signals processed: %d\n\n", signum);
idamnjanovic@51:   
idamnjanovic@51:   if (pd->DtX_time_counted) {
idamnjanovic@51:     mexPrintf("Compute DtX time:      %7.3lf seconds\n", pd->DtX_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   if (pd->XtX_time_counted) {
idamnjanovic@51:     mexPrintf("Compute XtX time:      %7.3lf seconds\n", pd->XtX_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   mexPrintf("Max abs time:          %7.3lf seconds\n", pd->maxabs_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   if (pd->DtD_time_counted) {
idamnjanovic@51:     mexPrintf("Compute DtD time:      %7.3lf seconds\n", pd->DtD_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   mexPrintf("Lchol update time:     %7.3lf seconds\n", pd->Lchol_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   mexPrintf("Compute coef time:     %7.3lf seconds\n", pd->compcoef_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   if (pd->compres_time_counted) {
idamnjanovic@51:     mexPrintf("Compute R time:        %7.3lf seconds\n", pd->compres_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   if (pd->DtR_time_counted) {
idamnjanovic@51:     mexPrintf("Compute DtR time:      %7.3lf seconds\n", pd->DtR_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   if (pd->update_DtR_time_counted) {
idamnjanovic@51:     mexPrintf("Update DtR time:       %7.3lf seconds\n", pd->update_DtR_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   if (pd->resnorm_time_counted) {
idamnjanovic@51:     mexPrintf("Update resnorm time:   %7.3lf seconds\n", pd->update_resnorm_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   if (pd->indexsort_time_counted) {
idamnjanovic@51:     mexPrintf("Index sort time:       %7.3lf seconds\n", pd->indexsort_time/(double)CLOCKS_PER_SEC);
idamnjanovic@51:   }
idamnjanovic@51:   mexPrintf("---------------------------------------\n");
idamnjanovic@51:   mexPrintf("Total time:            %7.3lf seconds\n\n", tottime/(double)CLOCKS_PER_SEC);
idamnjanovic@51: }
idamnjanovic@51: