/* This program shows how to use PAPI_sprofil */
#include <stdlib.h>
#include <stdio.h>
#include <unistd.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <string.h>
#include "papi.h" /* This needs to be included every time you use PAPI */
#define NUM_FLOPS 20000000
#define NUM_ITERS 100000
#define THRESHOLD 100000
#define ERROR_RETURN(retval) { fprintf(stderr, "Error %d %s:line %d: \n", retval,__FILE__,__LINE__); exit(retval); }
#if (defined(linux) && defined(__ia64__)) || (defined(_AIX))
#define DO_FLOPS1 (caddr_t)(*(void **)do_flops1)
#define DO_FLOPS2 (caddr_t)(*(void **)do_flops2)
#else
#define DO_FLOPS1 (caddr_t)(do_flops1)
#define DO_FLOPS2 (caddr_t)(do_flops2)
#endif
void do_flops2(int);
volatile double t1 = 0.8, t2 = 0.9;
void do_flops1(int n)
{
int i;
double c = 22222.11;
for (i = 0; i < n; i++)
c -= t1 * t2;
}
void do_both(int n)
{
int i;
const int flops2 = NUM_FLOPS / n;
const int flops1 = NUM_FLOPS / n;
for (i = 0; i < n; i++)
{
do_flops1(flops1);
do_flops2(flops2);
}
}
int main(int argc, char **argv)
{
int i , PAPI_event;
int EventSet = PAPI_NULL;
unsigned short *profbuf;
unsigned short *profbuf2;
unsigned short *profbuf3;
unsigned long length;
caddr_t start, end;
long long values[2];
const PAPI_exe_info_t *prginfo = NULL;
PAPI_sprofil_t sprof[3];
int retval;
/* initializaion */
if ((retval = PAPI_library_init(PAPI_VER_CURRENT)) != PAPI_VER_CURRENT)
{
printf("Library initialization error! \n");
exit(1);
}
if ((prginfo = PAPI_get_executable_info()) == NULL)
ERROR_RETURN(1);
start = prginfo->address_info.text_start;
end = prginfo->address_info.text_end;
length = (end - start)/sizeof(unsigned short) * sizeof(unsigned short);
printf("start= %p end =%p \n", start, end);
profbuf = (unsigned short *) malloc(length);
if (profbuf == NULL)
ERROR_RETURN(PAPI_ESYS);
memset(profbuf, 0x00, length );
profbuf2 = (unsigned short *) malloc(length);
if (profbuf2 == NULL)
ERROR_RETURN(PAPI_ESYS);
memset(profbuf2, 0x00, length );
profbuf3 = (unsigned short *) malloc(1 * sizeof(unsigned short));
if (profbuf3 == NULL)
ERROR_RETURN(PAPI_ESYS);
memset(profbuf3, 0x00, 1 * sizeof(unsigned short));
/* First half */
sprof[0].pr_base = profbuf;
sprof[0].pr_size = length / 2;
sprof[0].pr_off = DO_FLOPS2;
fprintf(stderr, "do_flops is at %p %lx\n", &do_flops2, strtoul(sprof[0].pr_off,NULL,0));
sprof[0].pr_scale = 65536; /* constant needed by PAPI_sprofil */
/* Second half */
sprof[1].pr_base = profbuf2;
sprof[1].pr_size = length / 2;
sprof[1].pr_off = DO_FLOPS1;
fprintf(stderr, "do_flops1 is at %p %lx\n", &do_flops1, strtoul(sprof[1].pr_off,NULL,0));
sprof[1].pr_scale = 65536; /* constant needed by PAPI_sprofil */
/* Overflow bin */
sprof[2].pr_base = profbuf3;
sprof[2].pr_size = 1;
sprof[2].pr_off = 0;
sprof[2].pr_scale = 0x2; /* constant needed by PAPI_sprofil */
/* Creating the eventset */
if ( (retval = PAPI_create_eventset(&EventSet)) != PAPI_OK)
ERROR_RETURN(retval);
PAPI_event = PAPI_TOT_CYC;
/* Add Total Instructions Executed to our EventSet */
if ( (retval = PAPI_add_event(EventSet, PAPI_event)) != PAPI_OK)
ERROR_RETURN(retval);
/* Add Total Instructions Executed to our EventSet */
if ( (retval = PAPI_add_event(EventSet, PAPI_TOT_INS)) != PAPI_OK)
ERROR_RETURN(retval);
/* set profile flag */
if ((retval = PAPI_sprofil(sprof, 3, EventSet, PAPI_event, THRESHOLD,
PAPI_PROFIL_POSIX)) != PAPI_OK)
ERROR_RETURN(retval);
if ((retval = PAPI_start(EventSet)) != PAPI_OK)
ERROR_RETURN(retval);
do_both(NUM_ITERS);
if ((retval = PAPI_stop(EventSet, values)) != PAPI_OK)
ERROR_RETURN(retval);
/* to clear the profile flag before removing the events */
if ((retval = PAPI_sprofil(sprof, 3, EventSet, PAPI_event, 0,
PAPI_PROFIL_POSIX)) != PAPI_OK)
ERROR_RETURN(retval);
/* free the resources hold by PAPI */
PAPI_shutdown();
printf("Test case: PAPI_sprofil()\n");
printf("---------Buffer 1--------\n");
for (i = 0; i < length / 2; i++)
{
if (profbuf[i])
printf("%#lx\t%d\n", strtoul(DO_FLOPS2,NULL,0) + 2 * i, profbuf[i]);
}
printf("---------Buffer 2--------\n");
for (i = 0; i < length / 2; i++)
{
if (profbuf2[i])
printf("%#lx\t%d\n", strtoul(DO_FLOPS1,NULL,0) + 2 * i, profbuf2[i]);
}
printf("-------------------------\n");
printf("%u samples fell outside the regions.\n", *profbuf3);
exit(0);
}
/* Declare a and b to be volatile.
This is to try to keep the
compiler from optimizing the loop */
volatile double a = 0.5, b = 2.2;
void do_flops2(int n)
{
int i;
double c = 0.11;
for (i = 0; i < n; i++)
c += a * b;
}