cpi.c

Go to the documentation of this file.

/* From Dave McNamara at PSRV. Thanks! */

/* If an event is countable but you've exhausted the counter resources
and you try to add an event, it seems subsequent PAPI_start and/or
PAPI_stop will causes a Seg. Violation.

   I got around this by calling PAPI to get the # of countable events,
then making sure that I didn't try to add more than these number of
events. I still have a problem if someone adds Level 2 cache misses
and then adds FLOPS 'cause I didn't count FLOPS as actually requiring
2 counters. */

#include "papi_test.h"
#include <mpi.h>
#include <stdio.h>
#include <math.h>

extern int TESTS_QUIET;                /* Declared in test_utils.c */
char *netevents[] =
    { "LO_RX_PACKETS", "LO_TX_PACKETS", "ETH0_RX_PACKETS", "ETH0_TX_PACKETS" };

double
f( double a )
{
    return ( 4.0 / ( 1.0 + a * a ) );
}

int
main( int argc, char **argv )
{
    int EventSet = PAPI_NULL, EventSet1 = PAPI_NULL;
    int evtcode;
    int retval, i, ins = 0;
    long long g1[2], g2[2];

    int done = 0, n, myid, numprocs;
    double PI25DT = 3.141592653589793238462643;
    double mypi, pi, h, sum, x;
    double startwtime = 0.0, endwtime;
    int namelen;
    char processor_name[MPI_MAX_PROCESSOR_NAME];

    tests_quiet( argc, argv );  /* Set TESTS_QUIET variable */

    if ( ( retval =
           PAPI_library_init( PAPI_VER_CURRENT ) ) != PAPI_VER_CURRENT )
        test_fail( __FILE__, __LINE__, "PAPI_library_init", retval );

    if ( ( retval = PAPI_create_eventset( &EventSet ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_create_eventset", retval );

    if ( ( retval = PAPI_create_eventset( &EventSet1 ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_create_eventset", retval );

    PAPI_event_name_to_code( netevents[2], &evtcode );
    if ( ( retval = PAPI_query_event( evtcode ) ) != PAPI_OK ) {
        if ( retval != PAPI_ECNFLCT )
            test_fail( __FILE__, __LINE__, "PAPI_aquery_event", retval );
    }
    if ( ( retval = PAPI_add_event( EventSet, evtcode ) ) != PAPI_OK ) {
        if ( retval != PAPI_ECNFLCT )
            test_fail( __FILE__, __LINE__, "PAPI_add_event", retval );
    }

    PAPI_event_name_to_code( netevents[3], &evtcode );
    if ( ( retval = PAPI_query_event( evtcode ) ) != PAPI_OK ) {
        if ( retval != PAPI_ECNFLCT )
            test_fail( __FILE__, __LINE__, "PAPI_aquery_event", retval );
    }
    if ( ( retval = PAPI_add_event( EventSet, evtcode ) ) != PAPI_OK ) {
        if ( retval != PAPI_ECNFLCT )
            test_fail( __FILE__, __LINE__, "PAPI_add_event", retval );
    }

    if ( ( retval = PAPI_query_event( PAPI_FP_INS ) ) != PAPI_OK ) {
        if ( ( retval = PAPI_query_event( PAPI_FP_OPS ) ) == PAPI_OK ) {
            ins = 2;
            if ( ( retval =
                   PAPI_add_event( EventSet1, PAPI_FP_OPS ) ) != PAPI_OK ) {
                if ( retval != PAPI_ECNFLCT )
                    test_fail( __FILE__, __LINE__, "PAPI_add_event", retval );
            }
        }
    } else {
        ins = 1;
        if ( ( retval = PAPI_add_event( EventSet1, PAPI_FP_INS ) ) != PAPI_OK ) {
            if ( retval != PAPI_ECNFLCT )
                test_fail( __FILE__, __LINE__, "PAPI_add_event", retval );
        }
    }

    if ( ( retval = PAPI_add_event( EventSet1, PAPI_TOT_CYC ) ) != PAPI_OK ) {
        if ( retval != PAPI_ECNFLCT )
            test_fail( __FILE__, __LINE__, "PAPI_add_event", retval );
    }

    MPI_Init( &argc, &argv );

    MPI_Comm_size( MPI_COMM_WORLD, &numprocs );
    MPI_Comm_rank( MPI_COMM_WORLD, &myid );
    MPI_Get_processor_name( processor_name, &namelen );

    fprintf( stdout, "Process %d of %d on %s\n",
             myid, numprocs, processor_name );
    fflush( stdout );

    if ( ( retval = PAPI_start( EventSet ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_start", retval );

    if ( ( retval = PAPI_start( EventSet1 ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_start", retval );

    n = 0;
    while ( !done ) {
        if ( myid == 0 ) {
            if ( n == 0 )
                n = 1000000;
            else
                n = 0;

            startwtime = MPI_Wtime(  );
        }
        MPI_Bcast( &n, 1, MPI_INT, 0, MPI_COMM_WORLD );
        if ( n == 0 )
            done = 1;
        else {
            h = 1.0 / ( double ) n;
            sum = 0.0;
            /* A slightly better approach starts from large i and works back */
            for ( i = myid + 1; i <= n; i += numprocs ) {
                x = h * ( ( double ) i - 0.5 );
                sum += f( x );
            }
            mypi = h * sum;

            MPI_Reduce( &mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD );

            if ( myid == 0 ) {
                printf( "pi is approximately %.16f, Error is %.16f\n",
                        pi, fabs( pi - PI25DT ) );
                endwtime = MPI_Wtime(  );
                printf( "wall clock time = %f\n", endwtime - startwtime );
                fflush( stdout );
            }
        }
    }

    if ( ( retval = PAPI_stop( EventSet1, g1 ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_stop", retval );

    if ( ( retval = PAPI_stop( EventSet, g2 ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_stop", retval );

    MPI_Finalize(  );


    printf( "ETH0_RX_BYTES: %lld   ETH0_TX_BYTES: %lld\n", g2[0], g2[1] );
    if ( ins == 0 ) {
        printf( "PAPI_TOT_CYC : %lld\n", g1[0] );
    } else if ( ins == 1 ) {
        printf( "PAPI_FP_INS  : %lld   PAPI_TOT_CYC : %lld\n", g1[0], g1[1] );
    } else if ( ins == 2 ) {
        printf( "PAPI_FP_OPS  : %lld   PAPI_TOT_CYC : %lld\n", g1[0], g1[1] );
    }
    test_pass( __FILE__, NULL, 0 );
    return 0;
}