/* file hl_rates.c

 * This test exercises the four PAPI High Level rate calls:

 *    PAPI_flops, PAPI_flips, PAPI_ipc, and PAPI_epc

 * flops and flips report cumulative real and process time since the first call,

 * and either floating point operations or instructions since the first call.

 * Also reported is incremental flop or flip rate since the last call.

 *

 * PAPI_ipc reports the same cumulative information, substituting

 * total instructions for flops or flips, and also reports

 * instructions per (process) cycle as a measure of execution efficiency.

 *

 * PAPI_epc is new in PAPI 5.2. It reports the same information as PAPI_IPC,

 * but for an arbitrary event instead of total cycles. It also reports

 * incremental core and (where available) reference cycles to allow the

 * computation of effective clock rates in the presence of clock scaling

 * like speed step or turbo-boost.

 *

 * This test computes a 1000 x 1000 matrix multiply for orders of indexing for

 * each of the four rate calls. It also accepts a command line parameter

 * for the event to be measured for PAPI_epc. If not provided, PAPI_TOT_INS

 * is measured.

 */

#include <stdio.h>

#include <stdlib.h>

#include "papi.h"

#include "papi_test.h"

#include "testcode.h"

int

main( int argc, char **argv )

{

	int retval, event = 0;

	float rtime, ptime, mflips, mflops, ipc, epc;

	long long flpins, flpops, ins, ref, core, evt;

	double mflips_classic,mflips_swapped;

	double mflops_classic,mflops_swapped;

	double ipc_classic,ipc_swapped;

	double epc_classic,epc_swapped;

	int quiet;

	/* Set TESTS_QUIET variable */

	quiet=tests_quiet( argc, argv );

	/* Initialize the test matrix */

	flops_float_init_matrix();

	/************************/

	/* FLIPS		*/

	/************************/

	if (!quiet) {

		printf( "\n----------------------------------\n" );

		printf( "PAPI_flips\n");

	}

	/* Run flips at start */

	retval=PAPI_flips(&rtime, &ptime, &flpins, &mflips);

	if (retval!=PAPI_OK) {

		if (!quiet) PAPI_perror( "PAPI_flips" );

		if (retval==PAPI_ENOEVNT) {

			test_skip(__FILE__,__LINE__,"Could not add event",0);

		}

	}

	if (!quiet) {

		printf( "\nStart\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Instructions: %lld\n", flpins);

		printf( "MFLIPS           %f\n", mflips);

	}

	/* Be sure we are all zero at beginning */

	if ((rtime!=0) || (ptime!=0) || (flpins!=0) || (mflips!=0)) {

		test_fail(__FILE__,__LINE__,"Not initialized to zero",0);

	}

	// Flips classic

	flops_float_matrix_matrix_multiply();

	if ( PAPI_flips(&rtime, &ptime, &flpins, &mflips)  != PAPI_OK )

		PAPI_perror( "PAPI_flips" );

	if (!quiet) {

		printf( "\nClassic\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Instructions: %lld\n", flpins);

		printf( "MFLIPS           %f\n", mflips);

	}

	mflips_classic=mflips;

	// Flips swapped

	flops_float_swapped_matrix_matrix_multiply();

	if ( PAPI_flips(&rtime, &ptime, &flpins, &mflips)  != PAPI_OK )

		PAPI_perror( "PAPI_flips" );

	if (!quiet) {

		printf( "\nSwapped\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Instructions: %lld\n", flpins);

		printf( "MFLIPS           %f\n", mflips);

	}

	mflips_swapped=mflips;

	// turn off flips

	if ( PAPI_stop_counters(NULL, 0)  != PAPI_OK ) {

		PAPI_perror( "PAPI_stop_counters" );

	}

	/************************/

	/* FLOPS		*/

	/************************/

	if (!quiet) {

		printf( "\n----------------------------------\n" );

		printf( "PAPI_flops\n");

	}

	// Start flops

	if ( PAPI_flops(&rtime, &ptime, &flpops, &mflops)  != PAPI_OK ) {

		PAPI_perror( "PAPI_flops" );

	}

	if (!quiet) {

		printf( "\nStart\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Operations:   %lld\n", flpops);

		printf( "MFLOPS           %f\n", mflops);

	}

	/* Be sure we are all zero at beginning */

	if ((rtime!=0) || (ptime!=0) || (flpops!=0) || (mflops!=0)) {

		test_fail(__FILE__,__LINE__,"Not initialized to zero",0);

	}

	// Classic flops

	flops_float_matrix_matrix_multiply();

	if ( PAPI_flops(&rtime, &ptime, &flpops, &mflops)  != PAPI_OK )

		PAPI_perror( "PAPI_flops" );

	if (!quiet) {

		printf( "\nClassic\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Operations:   %lld\n", flpops);

		printf( "MFLOPS           %f\n", mflops);

	}

	mflops_classic=mflops;

	// Swapped flops

	flops_float_swapped_matrix_matrix_multiply();

	if ( PAPI_flops(&rtime, &ptime, &flpops, &mflops)  != PAPI_OK )

		PAPI_perror( "PAPI_flops" );

	if (!quiet) {

		printf( "\nSwapped\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "FP Operations:   %lld\n", flpops);

		printf( "MFLOPS           %f\n", mflops);

	}

	mflops_swapped=mflops;

	// turn off flops

	if ( PAPI_stop_counters(NULL, 0)  != PAPI_OK ) {

		PAPI_perror( "PAPI_stop_counters" );

	}

	/************************/

	/* IPC  		*/

	/************************/

	if (!quiet) {

		printf( "\n----------------------------------\n" );

		printf( "PAPI_ipc\n");

	}

	// Start ipc

	if ( PAPI_ipc(&rtime, &ptime, &ins, &ipc)  != PAPI_OK )

		PAPI_perror( "PAPI_ipc" );

	if (!quiet) {

		printf( "\nStart\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Instructions:    %lld\n", ins);

		printf( "IPC              %f\n", ipc);

	}

	/* Be sure we are all zero at beginning */

	if ((rtime!=0) || (ptime!=0) || (ins!=0) || (ipc!=0)) {

		test_fail(__FILE__,__LINE__,"Not initialized to zero",0);

	}

	// Classic ipc

	flops_float_matrix_matrix_multiply();

	if ( PAPI_ipc(&rtime, &ptime, &ins, &ipc)  != PAPI_OK )

		PAPI_perror( "PAPI_ipc" );

	if (!quiet) {

		printf( "\nClassic\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Instructions:    %lld\n", ins);

		printf( "IPC              %f\n", ipc);

	}

	ipc_classic=ipc;

	// Swapped ipc

	flops_float_swapped_matrix_matrix_multiply();

	if ( PAPI_ipc(&rtime, &ptime, &ins, &ipc)  != PAPI_OK )

		PAPI_perror( "PAPI_ipc" );

	if (!quiet) {

		printf( "\nSwapped\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Instructions:    %lld\n", ins);

		printf( "IPC              %f\n", ipc);

	}

	ipc_swapped=ipc;

	// turn off ipc

	if ( PAPI_stop_counters(NULL, 0)  != PAPI_OK ) {

		PAPI_perror( "PAPI_stop_counters" );

	}

	/************************/

	/* EPC  		*/

	/************************/

	if (!quiet) {

		printf( "\n----------------------------------\n" );

		printf( "PAPI_epc\n");

	}

	/* This unfortunately conflicts a bit with the TESTS_QUIET */

	/* command line paramater nonsense.			   */

	if ( argc >= 2) {

		retval = PAPI_event_name_to_code( argv[1], &event );

		if (retval != PAPI_OK) {

		 	if (!quiet) printf("Can't find %s; Using PAPI_TOT_INS\n", argv[1]);

		 	event = PAPI_TOT_INS;

		} else {

		 	if (!quiet) printf("Using event %s\n", argv[1]);

		}

	}

	// Start epc

	if ( PAPI_epc(event, &rtime, &ptime, &ref, &core, &evt, &epc)  != PAPI_OK )

		PAPI_perror( "PAPI_epc" );

	if (!quiet) {

		printf( "\nStart\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Ref Cycles:      %lld\n", ref);

		printf( "Core Cycles:     %lld\n", core);

		printf( "Events:          %lld\n", evt);

		printf( "EPC:             %f\n", epc);

	}

	/* Be sure we are all zero at beginning */

	if ((rtime!=0) || (ptime!=0) || (ref!=0) || (core!=0)

			|| (evt!=0) || (epc!=0)) {

		test_fail(__FILE__,__LINE__,"Not initialized to zero",0);

	}

	// Classic epc

	flops_float_matrix_matrix_multiply();

	if ( PAPI_epc(event, &rtime, &ptime, &ref, &core, &evt, &epc)  != PAPI_OK )

		PAPI_perror( "PAPI_epc" );

	if (!quiet) {

		printf( "\nClassic\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Ref Cycles:      %lld\n", ref);

		printf( "Core Cycles:     %lld\n", core);

		printf( "Events:          %lld\n", evt);

		printf( "EPC:             %f\n", epc);

	}

	epc_classic=epc;

	// Swapped epc

	flops_float_swapped_matrix_matrix_multiply();

	if ( PAPI_epc(event, &rtime, &ptime, &ref, &core, &evt, &epc)  != PAPI_OK ) {

		PAPI_perror( "PAPI_epc" );

	}

	if (!quiet) {

		printf( "\nSwapped\n");

		printf( "real time:       %f\n", rtime);

		printf( "process time:    %f\n", ptime);

		printf( "Ref Cycles:      %lld\n", ref);

		printf( "Core Cycles:     %lld\n", core);

		printf( "Events:          %lld\n", evt);

		printf( "EPC:             %f\n", epc);

	}

	epc_swapped=epc;

	// turn off epc

	if ( PAPI_stop_counters(NULL, 0)  != PAPI_OK ) {

		PAPI_perror( "PAPI_stop_counters" );

	}

	if (!quiet) {

		printf( "\n----------------------------------\n" );

	}

	/* Validate */

	if (mflips_swapped

		test_fail(__FILE__,__LINE__,"FLIPS should be better when swapped",0);

	}

	if (mflops_swapped

		test_fail(__FILE__,__LINE__,"FLOPS should be better when swapped",0);

	}

	if (ipc_swapped

		test_fail(__FILE__,__LINE__,"IPC should be better when swapped",0);

	}

	if (epc_swapped

		test_fail(__FILE__,__LINE__,"EPC should be better when swapped",0);

	}

	test_pass( __FILE__ );

	return 0;

}