mendes-alt.c

Go to the documentation of this file.

#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;
 
}