Include dependency graph for libmsr_basic.c:

Data Structures
union	ll_dbl_union_t

Macros
#define	MAX_LIBMSR_EVENTS 64
	Tests basic functionality of libmsr component. More...

#define	MATRIX_SIZE 1024

Functions
void	run_test (int quiet)

int	main (int argc, char **argv)

Variables
static double	a [MATRIX_SIZE][MATRIX_SIZE]

static double	b [MATRIX_SIZE][MATRIX_SIZE]

static double	c [MATRIX_SIZE][MATRIX_SIZE]

Macro Definition Documentation

◆ MATRIX_SIZE

#define MATRIX_SIZE 1024

Definition at line 36 of file libmsr_basic.c.

◆ MAX_LIBMSR_EVENTS

#define MAX_LIBMSR_EVENTS 64

Author: Asim YarKhan; Vince Weaver (original rapl version)

Test case for libmsr component

Definition at line 22 of file libmsr_basic.c.

Function Documentation

◆ main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 76 of file libmsr_basic.c.

 {
     int retval,cid,libmsr_cid=-1,numcmp;
     int EventSet = PAPI_NULL;
     long long *values;
     int num_events=0;
     int code;
     char event_names[MAX_LIBMSR_EVENTS][PAPI_MAX_STR_LEN];
     char units[MAX_LIBMSR_EVENTS][PAPI_MIN_STR_LEN];
     int data_type[MAX_LIBMSR_EVENTS];
     int r,i, do_wrap = 0;
     const PAPI_component_info_t *cmpinfo = NULL;
     PAPI_event_info_t evinfo;
     long long before_time,after_time;
     double elapsed_time;
     ll_dbl_union_t tmp_ll_dbl;
     int repeat;
 
     /* Set TESTS_QUIET variable */
     tests_quiet( argc, argv );
     if ( argc > 1 && strstr( argv[1], "-w" ) )
         do_wrap = 1;
 
     /* PAPI Initialization */
     retval = PAPI_library_init( PAPI_VER_CURRENT );
     if ( retval != PAPI_VER_CURRENT ) 
         test_fail( __FILE__, __LINE__,"PAPI_library_init failed\n",retval );
     
     if ( !TESTS_QUIET ) 
         printf( "Trying all LIBMSR events\n" );
 
     numcmp = PAPI_num_components();
 
     for( cid=0; cid<numcmp; cid++ ) {
 
         if ( ( cmpinfo = PAPI_get_component_info( cid ) ) == NULL ) 
             test_fail( __FILE__, __LINE__,"PAPI_get_component_info failed\n", 0 );
         
         if ( strstr( cmpinfo->name,"libmsr" ) ) {
             libmsr_cid=cid;
             if ( !TESTS_QUIET ) printf( "Found libmsr component at cid %d\n",libmsr_cid );
             if ( cmpinfo->disabled ) {
                 if ( !TESTS_QUIET ) printf( "libmsr component disabled: %s\n",cmpinfo->disabled_reason );
                 test_skip( __FILE__,__LINE__, "libmsr component disabled", 0 );
             }
             break;
         }
     }
     
     /* Component not found */
     if ( cid==numcmp ) 
         test_skip( __FILE__,__LINE__,"No libmsr component found\n",0 );
 
     /* Create EventSet */
     retval = PAPI_create_eventset( &EventSet );
     if ( retval != PAPI_OK ) 
         test_fail( __FILE__, __LINE__, "PAPI_create_eventset()",retval );
 
     /* Add all events */
     num_events = 0;
     code = PAPI_NATIVE_MASK;
     r = PAPI_enum_cmp_event( &code, PAPI_ENUM_FIRST, libmsr_cid );
     while ( r == PAPI_OK ) {
         retval = PAPI_event_code_to_name( code, event_names[num_events] );
         if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__, "PAPI_event_code_to_name", retval );
         retval = PAPI_get_event_info( code,&evinfo );
         if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__,"Error getting event info\n",retval );
         strncpy( units[num_events], evinfo.units, sizeof( units[0] )-1 );
         // buffer must be null terminated to safely use strstr operation on it below
         units[num_events][sizeof( units[0] )-1] = '\0';
         data_type[num_events] = evinfo.data_type;
 
         retval = PAPI_add_event( EventSet, code );
         if ( retval != PAPI_OK ) 
             break; /* We've hit an event limit */
         num_events++;
 
         r = PAPI_enum_cmp_event( &code, PAPI_ENUM_EVENTS, libmsr_cid );
     }
 
     values=calloc( num_events,sizeof( long long ) );
     if ( values==NULL ) 
         test_fail( __FILE__, __LINE__,"No memory - calloc failed",retval );
     
     if ( !TESTS_QUIET ) printf( "Starting measurements...\n" );
 
     /* Start Counting */
     before_time=PAPI_get_real_nsec();
     retval = PAPI_start( EventSet );
     if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__, "PAPI_start()",retval );
 
     for ( repeat=0; repeat<3; repeat++ ) {
 
         /* Run test */
         run_test( TESTS_QUIET );
         
         /* Stop Counting */
         after_time=PAPI_get_real_nsec();
         retval = PAPI_read( EventSet, values );
         if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__, "PAPI_read()",retval );
         
         elapsed_time=( ( double )( after_time-before_time ) )/1.0e9;
         
         if ( !TESTS_QUIET ) {
             printf( "Stopping measurements, took %.3fs, gathering results...\n", elapsed_time );
             for( i=0; i<num_events; i++ ) {
                 tmp_ll_dbl.ll = values[i];
                 if ( data_type[i]==PAPI_DATATYPE_FP64 ) {
                     printf( "%-40s %12.6f %s\n", event_names[i], tmp_ll_dbl.dbl, units[i] );
                 } else if ( data_type[i]==PAPI_DATATYPE_UINT64 ) {
                     printf( "%-40s %12ld %s\n", event_names[i], tmp_ll_dbl.ll, units[i] );
                 }
             }            
             printf( "\n" );
         }
     }
 
     retval = PAPI_stop( EventSet, values );
     if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__, "PAPI_stop()",retval );
      
     retval = PAPI_cleanup_eventset( EventSet );
     if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__, "PAPI_cleanup_eventset()",retval );
 
     retval = PAPI_destroy_eventset( &EventSet );
     if ( retval != PAPI_OK ) test_fail( __FILE__, __LINE__,"PAPI_destroy_eventset()",retval );
 
     test_pass( __FILE__ );
 
     return 0;
 }

Here is the call graph for this function:

◆ run_test()

void run_test ( int quiet )

Definition at line 43 of file libmsr_basic.c.

 {
     double s;
     int i,j,k;
     if ( !quiet ) 
         printf( "Doing a naive %dx%d MMM...\n",MATRIX_SIZE,MATRIX_SIZE );
     for( i=0; i<MATRIX_SIZE; i++ ) {
         for( j=0; j<MATRIX_SIZE; j++ ) {
             a[i][j]=( double )i*( double )j;
             b[i][j]=( double )i/( double )( j+5 );
         }
     }
     for( j=0; j<MATRIX_SIZE; j++ ) {
         for( i=0; i<MATRIX_SIZE; i++ ) {
             s=0;
             for( k=0; k<MATRIX_SIZE; k++ ) {
                 s+=a[i][k]*b[k][j];
             }
             c[i][j] = s;
         }
     }
     s=0.0;
     for( i=0; i<MATRIX_SIZE; i++ ) {
         for( j=0; j<MATRIX_SIZE; j++ ) {
             s+=c[i][j];
         }
     }
     if ( !quiet ) printf( "Matrix multiply sum: s=%lf\n",s );
 }