byte_profile.c

Go to the documentation of this file.

/*
* File:    byte_profile.c
* Author:  Dan Terpstra
*          terpstra@cs.utk.edu
* Mods:    Maynard Johnson
*          maynardj@us.ibm.com
*/

/* This file profiles multiple events with byte level address resolution.
   It's patterned after code suggested by John Mellor-Crummey, Rob Fowler,
   and Nathan Tallent.
   It is intended to illustrate the use of Multiprofiling on a very tight
   block of code at byte level resolution of the instruction addresses.
*/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "papi.h"
#include "papi_test.h"
#include "prof_utils.h"

#include "do_loops.h"

#define PROFILE_ALL

static const PAPI_hw_info_t *hw_info;

static int num_events = 0;

#define N (1 << 23)
#define T (10)

double aa[N], bb[N];
double s = 0, s2 = 0;

static void
cleara( double a[N] )
{
    int i;

    for ( i = 0; i < N; i++ ) {
        a[i] = 0;
    }
}

static int
my_dummy( int i )
{
    return ( i + 1 );
}

static void
my_main( void )
{
    int i, j;

    for ( j = 0; j < T; j++ ) {
        for ( i = 0; i < N; i++ ) {
            bb[i] = 0;
        }
        cleara( aa );
        memset( aa, 0, sizeof ( aa ) );
        for ( i = 0; i < N; i++ ) {
            s += aa[i] * bb[i];
            s2 += aa[i] * aa[i] + bb[i] * bb[i];
        }
    }
}

static int
do_profile( caddr_t start, unsigned long plength, unsigned scale, int thresh,
        int bucket, unsigned int mask ) {

    int i, retval;
    unsigned long blength;
    int num_buckets,j=0;

    int num_bufs = num_events;
    int event = num_events;

    int events[MAX_TEST_EVENTS];
    char header[BUFSIZ];

    strncpy(header,"address\t\t",BUFSIZ);

    //= "address\t\t\tcyc\tins\tfp_ins\n";

    for(i=0;i<MAX_TEST_EVENTS;i++) {
      if (mask & test_events[i].mask) {
        events[j]=test_events[i].event;

        if (events[j]==PAPI_TOT_CYC) {
           strncat(header,"\tcyc",BUFSIZ-1);
        }
        if (events[j]==PAPI_TOT_INS) {
           strncat(header,"\tins",BUFSIZ-1);
        }
        if (events[j]==PAPI_FP_INS) {
           strncat(header,"\tfp_ins",BUFSIZ-1);
        }
        if (events[j]==PAPI_FP_OPS) {
           strncat(header,"\tfp_ops",BUFSIZ-1);
        }
        if (events[j]==PAPI_L2_TCM) {
           strncat(header,"\tl2_tcm",BUFSIZ-1);
        }

        j++;

      }
    }

    strncat(header,"\n",BUFSIZ-1);



    blength = prof_size( plength, scale, bucket, &num_buckets );
    prof_alloc( num_bufs, blength );

    if ( !TESTS_QUIET ) printf( "Overall event counts:\n" );

    for ( i = 0; i < num_events; i++ ) {
        if ( ( retval =
               PAPI_profil( profbuf[i], ( unsigned int ) blength, start, scale,
                            EventSet, events[i], thresh,
                            PAPI_PROFIL_POSIX | bucket ) ) != PAPI_OK ) {
               if (retval == PAPI_EINVAL) {
              test_warn( __FILE__, __LINE__, "Trying to profile with derived event", 1);
              num_events=i;
              break;
           }
                   else {
                printf("Failed with event %d %#x\n",i,events[i]);
            test_fail( __FILE__, __LINE__, "PAPI_profil", retval );
           }
        }
    }

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

    my_main(  );

    if ( ( retval = PAPI_stop( EventSet, values[0] ) ) != PAPI_OK )
        test_fail( __FILE__, __LINE__, "PAPI_stop", retval );

    if ( !TESTS_QUIET ) {
        printf( TAB1, "PAPI_TOT_CYC:", ( values[0] )[--event] );
        if ( strcmp( hw_info->model_string, "POWER6" ) != 0 ) {
            printf( TAB1, "PAPI_TOT_INS:", ( values[0] )[--event] );
        }
#if defined(__powerpc__)
        printf( TAB1, "PAPI_FP_INS", ( values[0] )[--event] );
#else
        if ( strcmp( hw_info->model_string, "Intel Pentium III" ) != 0 ) {
            printf( TAB1, "PAPI_FP_OPS:", ( values[0] )[--event] );
            printf( TAB1, "PAPI_L2_TCM:", ( values[0] )[--event] );
        }
#endif
    }

    for ( i = 0; i < num_events; i++ ) {
        if ( ( retval =
               PAPI_profil( profbuf[i], ( unsigned int ) blength, start, scale,
                            EventSet, events[i], 0,
                            PAPI_PROFIL_POSIX ) ) != PAPI_OK )
            test_fail( __FILE__, __LINE__, "PAPI_profil", retval );
    }

    if (!TESTS_QUIET) {
        prof_head( blength, bucket, num_buckets, header );
        prof_out( start, num_events, bucket, num_buckets, scale );
    }
    retval = prof_check( num_bufs, bucket, num_buckets );
    for ( i = 0; i < num_bufs; i++ ) {
        free( profbuf[i] );
    }
    return retval;
}



int
main( int argc, char **argv )
{
    long length;
    int mask;
    int retval;
    const PAPI_exe_info_t *prginfo;
    caddr_t start, end;
    int quiet;

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

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

    if ( ( prginfo = PAPI_get_executable_info(  ) ) == NULL ) {
        test_fail( __FILE__, __LINE__, "PAPI_get_executable_info", 1 );
    }

    hw_info = PAPI_get_hardware_info(  );
        if ( hw_info == NULL ) {
      test_fail( __FILE__, __LINE__, "PAPI_get_hardware_info", 2 );
    }

        mask = MASK_TOT_CYC | MASK_TOT_INS | MASK_FP_OPS | MASK_L2_TCM;

#if defined(__powerpc__)
    if ( strcmp( hw_info->model_string, "POWER6" ) == 0 )
        mask = MASK_TOT_CYC | MASK_FP_INS;
    else
        mask = MASK_TOT_CYC | MASK_TOT_INS | MASK_FP_INS;
#endif

#if defined(ITANIUM2)
    mask = MASK_TOT_CYC | MASK_FP_OPS | MASK_L2_TCM | MASK_L1_DCM;
#endif
    EventSet = add_test_events( &num_events, &mask, 0 );
    if (num_events==0) {
        if (!quiet) printf("Trouble adding events\n");
        test_skip(__FILE__,__LINE__,"add_test_events",2);
    }
    values = allocate_test_space( 1, num_events );

/* profile the cleara and my_main address space */
    start = ( caddr_t ) cleara;
    end = ( caddr_t ) my_dummy;

/* Itanium and PowerPC64 processors return function descriptors instead
 * of function addresses. You must dereference the descriptor to get the address.
*/
#if defined(ITANIUM1) || defined(ITANIUM2) \
    || (defined(__powerpc64__) && (_CALL_ELF != 2))
    start = ( caddr_t ) ( ( ( struct fdesc * ) start )->ip );
    end = ( caddr_t ) ( ( ( struct fdesc * ) end )->ip );
        /* PPC64 Big Endian is ELF version 1 which uses function descriptors.
         *  PPC64 Little Endian is ELF version 2 which does not use
         * function descriptors
         */
#endif

    /* call dummy so it doesn't get optimized away */
    retval = my_dummy( 1 );

    length = end - start;
    if ( length < 0 )
        test_fail( __FILE__, __LINE__, "Profile length < 0!", ( int ) length );

    if (!quiet) {
        prof_print_address( "Test case byte_profile: "
                "Multi-event profiling at byte resolution.\n",
                prginfo );
        prof_print_prof_info( start, end, THRESHOLD, event_name );
    }

    retval = do_profile( start, ( unsigned ) length,
                FULL_SCALE * 2, THRESHOLD,
                PAPI_PROFIL_BUCKET_32, mask );

    remove_test_events( &EventSet, mask );

    if (retval == 0) {
        test_fail( __FILE__, __LINE__, "No information in buffers", 1 );
    }

    test_pass( __FILE__ );

    return 0;
}