#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "papi.h"
#include "papi_test.h"
#ifdef SETMAX
#define MAX SETMAX
#else
#define MAX 10000
#endif
#define TIMES 1000
#define PAPI_MAX_EVENTS 2
long long PAPI_values1[PAPI_MAX_EVENTS];
long long PAPI_values2[PAPI_MAX_EVENTS];
long long PAPI_values3[PAPI_MAX_EVENTS];
static int EventSet = PAPI_NULL;
void funcX( double a[MAX], double b[MAX], int n)
{
int i, k;
for ( k = 0; k < TIMES; k++ )
for ( i = 0; i < n; i++ )
a[i] = a[i] * b[i] + 1.;
}
void funcA( double a[MAX], double b[MAX], int n)
{
int i, k;
double t[MAX];
for ( k = 0; k < TIMES; k++ )
for ( i = 0; i < n; i++ ) {
t[i] = b[n - i];
b[i] = a[n - i];
a[i] = t[i];
}
}
int
main( int argc, char **argv )
{
int i, retval;
double a[MAX], b[MAX];
int quiet;
/* Set TESTS_QUIET variable */
quiet=tests_quiet( argc, argv );
for ( i = 0; i < MAX; i++ ) {
a[i] = 0.0;
b[i] = 0.;
}
for ( i = 0; i < PAPI_MAX_EVENTS; i++ )
PAPI_values1[i] = PAPI_values2[i] = 0;
retval = PAPI_library_init( PAPI_VER_CURRENT );
if ( retval != PAPI_VER_CURRENT )
test_fail( __FILE__, __LINE__, "PAPI_library_init", retval );
#ifdef MULTIPLEX
if ( !quiet ) {
printf( "Activating PAPI Multiplex\n" );
}
retval = PAPI_multiplex_init( );
if ( retval != PAPI_OK ) {
test_fail( __FILE__, __LINE__, "PAPI multiplex init fail\n",
retval );
}
#endif
retval = PAPI_create_eventset( &EventSet );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI set event fail\n", retval );
#ifdef MULTIPLEX
/* In Component PAPI, EventSets must be assigned a component index
before you can fiddle with their internals.
0 is always the cpu component */
retval = PAPI_assign_eventset_component( EventSet, 0 );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI_assign_eventset_component",
retval );
retval = PAPI_set_multiplex( EventSet );
if (retval == PAPI_ENOSUPP) {
test_skip( __FILE__, __LINE__, "Multiplex not supported", 1 );
}
else if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI_set_multiplex fails \n", retval );
#endif
retval = PAPI_add_event( EventSet, PAPI_FP_INS );
if ( retval < PAPI_OK ) {
retval = PAPI_add_event( EventSet, PAPI_TOT_INS );
if ( retval < PAPI_OK ) {
if (!quiet) printf("Trouble adding events\n");
test_skip( __FILE__, __LINE__,
"PAPI add PAPI_FP_INS or PAPI_TOT_INS fail\n", retval );
} else if ( !quiet ) {
printf( "PAPI_TOT_INS\n" );
}
} else if ( !quiet ) {
printf( "PAPI_FP_INS\n" );
}
retval = PAPI_add_event( EventSet, PAPI_TOT_CYC );
if ( retval < PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI add PAPI_TOT_CYC fail\n",
retval );
if ( !quiet ) {
printf( "PAPI_TOT_CYC\n" );
}
retval = PAPI_start( EventSet );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI start fail\n", retval );
funcX( a, b, MAX );
retval = PAPI_read( EventSet, PAPI_values1 );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI read fail \n", retval );
funcX( a, b, MAX );
retval = PAPI_read( EventSet, PAPI_values2 );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI read fail \n", retval );
#ifdef RESET
retval = PAPI_reset( EventSet );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI read fail \n", retval );
#endif
funcA( a, b, MAX );
retval = PAPI_stop( EventSet, PAPI_values3 );
if ( retval != PAPI_OK )
test_fail( __FILE__, __LINE__, "PAPI read fail \n", retval );
if ( !quiet ) {
printf( "values1 is:\n" );
for ( i = 0; i < PAPI_MAX_EVENTS; i++ )
printf( LLDFMT15, PAPI_values1[i] );
printf( "\nvalues2 is:\n" );
for ( i = 0; i < PAPI_MAX_EVENTS; i++ )
printf( LLDFMT15, PAPI_values2[i] );
printf( "\nvalues3 is:\n" );
for ( i = 0; i < PAPI_MAX_EVENTS; i++ )
printf( LLDFMT15, PAPI_values3[i] );
#ifndef RESET
printf( "\nPAPI value (2-1) is : \n" );
for ( i = 0; i < PAPI_MAX_EVENTS; i++ )
printf( LLDFMT15, PAPI_values2[i] - PAPI_values1[i] );
printf( "\nPAPI value (3-2) is : \n" );
for ( i = 0; i < PAPI_MAX_EVENTS; i++ ) {
long long diff;
diff = PAPI_values3[i] - PAPI_values2[i];
printf( LLDFMT15, diff);
if (diff<0) {
test_fail( __FILE__, __LINE__, "Multiplexed counter decreased", 1 );
}
}
#endif
printf( "\n\nVerification:\n" );
printf( "From start to first PAPI_read %d fp operations are made.\n",
2 * MAX * TIMES );
printf( "Between 1st and 2nd PAPI_read %d fp operations are made.\n",
2 * MAX * TIMES );
printf( "Between 2nd and 3rd PAPI_read %d fp operations are made.\n",
0 );
printf( "\n" );
}
test_pass( __FILE__ );
return 0;
}