flops.c

Go to the documentation of this file.

/*
 * A simple example for the use of PAPI, the number of flops you should
 * get is about INDEX^3  on machines that consider add and multiply one flop
 * such as SGI, and 2*(INDEX^3) that don't consider it 1 flop such as INTEL
 * -Kevin London
 */

#include "papi_test.h"

#define INDEX 1000

char format_string[] =
    { "Real_time: %f Proc_time: %f Total flpins: %lld MFLOPS: %f\n" };
extern int TESTS_QUIET;                /* Declared in test_utils.c */


extern void dummy( void * );
float matrixa[INDEX][INDEX], matrixb[INDEX][INDEX], mresult[INDEX][INDEX];
int
main( int argc, char **argv )
{
    float real_time, proc_time, mflops;
    long long flpins;
    int retval;
    int i, j, k, fip = 0;

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

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

    if ( PAPI_query_event( PAPI_FP_INS ) == PAPI_OK )
        fip = 1;
    else if ( PAPI_query_event( PAPI_FP_OPS ) == PAPI_OK )
        fip = 2;
    else {
        if ( !TESTS_QUIET )
            printf
                ( "PAPI_FP_INS and PAPI_FP_OPS are not defined for this platform.\n" );
    }

    PAPI_shutdown(  );

    if ( fip > 0 ) {
        /* Initialize the Matrix arrays */
        for ( i = 0; i < INDEX; i++ ) {
          for ( j = 0; j < INDEX; j++) {
            mresult[j][i] = 0.0;
            matrixa[j][i] = matrixb[j][i] = ( float ) rand(  ) * ( float ) 1.1;
          }
        }

        /* Setup PAPI library and begin collecting data from the counters */
        if ( fip == 1 ) {
            if ( ( retval =
                   PAPI_flips( &real_time, &proc_time, &flpins,
                               &mflops ) ) < PAPI_OK )
                test_fail( __FILE__, __LINE__, "PAPI_flips", retval );
        } else {
            if ( ( retval =
                   PAPI_flops( &real_time, &proc_time, &flpins,
                               &mflops ) ) < PAPI_OK )
                test_fail( __FILE__, __LINE__, "PAPI_flops", retval );
        }

        /* Matrix-Matrix multiply */
        for ( i = 0; i < INDEX; i++ )
            for ( j = 0; j < INDEX; j++ )
                for ( k = 0; k < INDEX; k++ )
                    mresult[i][j] =
                        mresult[i][j] + matrixa[i][k] * matrixb[k][j];

        /* Collect the data into the variables passed in */
        if ( fip == 1 ) {
            if ( ( retval =
                   PAPI_flips( &real_time, &proc_time, &flpins,
                               &mflops ) ) < PAPI_OK )
                test_fail( __FILE__, __LINE__, "PAPI_flips", retval );
        } else {
            if ( ( retval =
                   PAPI_flops( &real_time, &proc_time, &flpins,
                               &mflops ) ) < PAPI_OK )
                test_fail( __FILE__, __LINE__, "PAPI_flops", retval );
        }
        dummy( ( void * ) mresult );

        if ( !TESTS_QUIET ) {
            if ( fip == 1 ) {
                printf( "Real_time: %f Proc_time: %f Total flpins: ", real_time,
                        proc_time );
            } else {
                printf( "Real_time: %f Proc_time: %f Total flpops: ", real_time,
                        proc_time );
            }
            printf( LLDFMT, flpins );
            printf( " MFLOPS: %f\n", mflops );
        }
    }
    test_pass( __FILE__, NULL, 0 );
    exit( 1 );
}