/*
* 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 <stdio.h>
#include <stdlib.h>
#include "papi.h"
#include "papi_test.h"
#include "testcode.h"
#include "display_error.h"
int
main( int argc, char **argv )
{
float real_time, proc_time, mflops;
long long flpins;
int retval;
int fip = 0;
int quiet=0;
long long expected;
double double_result,error;
/* Set TESTS_QUIET variable */
quiet=tests_quiet( argc, argv );
/* Initialize PAPI library */
retval = PAPI_library_init( PAPI_VER_CURRENT );
if ( retval != PAPI_VER_CURRENT ) {
test_fail( __FILE__, __LINE__, "PAPI_library_init", retval );
}
/* Try to use one of the FP events */
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 ( !quiet ) printf( "PAPI_FP_INS and PAPI_FP_OPS are not defined for this platform.\n" );
test_skip(__FILE__,__LINE__,"No FP events available",1);
}
/* Shutdown? */
/* I guess because it would interfere with the high-level interface? */
PAPI_shutdown( );
/* Initialize the Matrix arrays */
expected=flops_float_init_matrix();
/* Setup PAPI library and begin collecting data from the counters */
if ( fip == 1 ) {
retval = PAPI_flips( &real_time, &proc_time, &flpins, &mflops );
if (retval!=PAPI_OK) {
test_fail( __FILE__, __LINE__, "PAPI_flips", retval );
}
}
else {
retval = PAPI_flops( &real_time, &proc_time, &flpins, &mflops );
if (retval!=PAPI_OK) {
test_fail( __FILE__, __LINE__, "PAPI_flops", retval );
}
}
/* Matrix-Matrix multiply */
double_result=flops_float_matrix_matrix_multiply();
/* Collect the data into the variables passed in */
if ( fip == 1 ) {
retval = PAPI_flips( &real_time, &proc_time, &flpins, &mflops );
if (retval!=PAPI_OK) {
test_fail( __FILE__, __LINE__, "PAPI_flips", retval );
}
} else {
retval = PAPI_flops( &real_time, &proc_time, &flpins, &mflops );
if (retval!=PAPI_OK) {
test_fail( __FILE__, __LINE__, "PAPI_flops", retval );
}
}
if (!quiet) printf("result=%lf\n",double_result);
if ( !quiet ) {
printf( "Real_time: %f Proc_time: %f MFLOPS: %f\n",
real_time, proc_time, mflops );
if ( fip == 1 ) {
printf( "Total flpins: ");
} else {
printf( "Total flpops: ");
}
printf( "%lld\n\n", flpins );
}
error=display_error(flpins,flpins,flpins,expected,quiet);
if ((error > 1.0) || (error<-1.0)) {
if (!quiet) printf("Instruction count off by more than 1%%\n");
test_fail( __FILE__, __LINE__, "Validation failed", 1 );
}
test_pass( __FILE__ );
return 0;
}