papi_l2_dca.c

Go to the documentation of this file.

/* This code attempts to test the L2 Data Cache Acceesses   */
/* performance counter PAPI_L2_DCA              */
 
/* Notes:                           */
/*  Should this be equivelent to PAPI_L1_DCM?       */
/*      (on IVY it is)                  */
/*  On Haswell/Broadwell/Skylake this maps to :     */
/*      L2_RQSTS:ALL_DEMAND_REFERENCES          */
 
/*  Should this include *all* L2 accesses or just those */
/*      caused by the user?  Prefetch? MESI?        */
 
 
 
/* by Vince Weaver, vincent.weaver@maine.edu            */
 
#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
 
#include "papi.h"
#include "papi_test.h"
 
#include "cache_helper.h"
#include "display_error.h"
 
#include "testcode.h"
 
#define NUM_RUNS    100
 
int main(int argc, char **argv) {
 
    int i;
    int eventset=PAPI_NULL;
    int num_runs=NUM_RUNS;
    long long high,low,average,expected;
    long long count,total;
 
    int retval;
    int l1_size,l2_size,l1_linesize,l2_linesize,l2_entries;
    int arraysize;
    int quiet,errors=0;
 
    double error;
    double *array;
    double aSumm = 0.0;
 
    quiet=tests_quiet(argc,argv);
 
    if (!quiet) {
        printf("Testing the PAPI_L2_DCA event\n");
    }
 
    /* Init the PAPI library */
    retval = PAPI_library_init(PAPI_VER_CURRENT);
    if (retval != PAPI_VER_CURRENT) {
        test_fail(__FILE__,__LINE__,"PAPI_library_init",retval);
    }
 
    retval=PAPI_create_eventset(&eventset);
    if (retval!=PAPI_OK) {
        test_fail( __FILE__, __LINE__, "PAPI_create_eventset", retval );
    }
 
    retval=PAPI_add_named_event(eventset,"PAPI_L2_DCA");
    if (retval!=PAPI_OK) {
        test_skip( __FILE__, __LINE__, "adding PAPI_L2_DCA", retval );
    }
 
    l1_size=get_cachesize(L1D_CACHE);
    l1_linesize=get_linesize(L1D_CACHE);
    l2_size=get_cachesize(L2_CACHE);
    l2_linesize=get_linesize(L2_CACHE);
    l2_entries=get_entries(L2_CACHE);
 
    if ((l2_size==0) || (l2_linesize==0)) {
        if (!quiet) {
            printf("Unable to determine size of L2 cache!\n");
        }
        test_skip( __FILE__, __LINE__, "adding PAPI_L2_DCA", retval );
    }
 
    if (!quiet) {
        printf("\tDetected %dk L1 DCache, %dB linesize\n",
            l1_size/1024,l1_linesize);
        printf("\tDetected %dk L2 DCache, %dB linesize, %d entries\n",
            l2_size/1024,l2_linesize,l2_entries);
    }
 
    arraysize=l2_size/sizeof(double);
 
    if (!quiet) {
        printf("\tAllocating %zu bytes of memory (%d doubles)\n",
            arraysize*sizeof(double),arraysize);
    }
 
    array=calloc(arraysize,sizeof(double));
    if (array==NULL) {
        test_fail(__FILE__,__LINE__,"Can't allocate memory",0);
    }
 
    /******************/
    /* Testing Writes */
    /******************/
 
    if (!quiet) {
        printf("\nWrite Test: Initializing an array of %d doubles:\n",
            arraysize);
        }
 
    high=0; low=0; total=0;
 
    for(i=0;i<num_runs;i++) {
 
        PAPI_reset(eventset);
        PAPI_start(eventset);
 
        cache_write_test(array,arraysize);
 
        retval=PAPI_stop(eventset,&count);
 
        if (retval!=PAPI_OK) {
            test_fail( __FILE__, __LINE__,
                "reading PAPI_L2_DCA", retval );
        }
 
        if (count>high) high=count;
        if ((low==0) || (count<low)) low=count;
        total+=count;
    }
 
    average=(total/num_runs);
 
    expected=arraysize/(l2_linesize/sizeof(double));
 
    if (!quiet) {
        printf("\tShould be roughly "
            "arraysize/L2_linesize/double_size (%d/%d/%zu): "
            "%lld\n\n",
            arraysize,l2_linesize,sizeof(double),
            expected);
    }
 
    error=display_error(average,high,low,expected,quiet);
 
    if ((error > 1.0) || (error<-1.0)) {
        if (!quiet) printf("Instruction count off by more than 1%%\n");
        errors++;
    }
 
    if (!quiet) printf("\n");
 
    /******************/
    /* Testing Reads  */
    /******************/
 
    if (!quiet) {
        printf("\nRead Test: Summing an array of %d doubles:\n",
            arraysize);
        }
 
    high=0; low=0; total=0;
 
    for(i=0;i<num_runs;i++) {
 
        PAPI_reset(eventset);
        PAPI_start(eventset);
 
        aSumm+=cache_read_test(array,arraysize);
 
        retval=PAPI_stop(eventset,&count);
 
        if (retval!=PAPI_OK) {
            test_fail( __FILE__, __LINE__,
                "reading PAPI_L2_DCA", retval );
        }
 
        if (count>high) high=count;
        if ((low==0) || (count<low)) low=count;
        total+=count;
    }
 
    average=(total/num_runs);
 
    expected=arraysize/(l2_linesize/sizeof(double));
 
    if (!quiet) {
        printf("\tShould be roughly "
            "arraysize/L2_linesize/double_size (%d/%d/%zu): "
            "%lld\n\n",
            arraysize,l2_linesize,sizeof(double),
            expected);
    }
 
    error=display_error(average,high,low,expected,quiet);
 
    if ((error > 1.0) || (error<-1.0)) {
        if (!quiet) printf("Instruction count off by more than 1%%\n");
        errors++;
    }
 
    if (!quiet) {
        printf("\n");
    }
 
    /* Warn for now, as we get errors we can't easily */
    /* explain on haswell and more recent Intel chips */
    if (errors) {
        test_warn( __FILE__, __LINE__, "Error too high", 1 );
    }
 
    test_pass(__FILE__);
 
    return 0;
}