/*
 This program is for testing large stride data transfer

I put this program under the GPL.

Contact:Naohiko Shimizu, 
	School of Engineering, Tokai University.
	1117 Kitakaname, Kanagawa 259-12 Japan
	email:nshimizu@keyaki.cc.u-tokai.ac.jp
	TEL: +81-463-58-1211(ext.4084)
	FAX: +81-463-58-8320
	<URL:http://shimizu-lab.et.u-tokai.ac.jp/>
*/


#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <math.h>
#include <string.h>
#include <sys/time.h>
#include <sys/mman.h>
#include <sys/resource.h>

#ifndef DIMENSION
#define DIMENSION 1000
#endif

#define ITR 10*4000*4000/(DIMENSION*DIMENSION)

double a[DIMENSION][DIMENSION];
double b[DIMENSION][DIMENSION];

struct rusage rusage;
double dtime()
{
 double q;

 getrusage(RUSAGE_SELF,&rusage);

 q = (double)(rusage.ru_utime.tv_sec);
 q = q + (double)(rusage.ru_utime.tv_usec) * 1.0e-06;
	
 return q;
}

int main(int argc, char **argv)
{
 int i,j,k;
 int    dim=DIMENSION;
 double starttime,endtime,ss,ls;
 for(i=0; i<dim; i++) {
	for(j=0; j< dim; j++) {
		a[j][i] = b[i][j] = 0.0;
		  }
	 }

 starttime = dtime();
   for(k=0; k<ITR; k++)
	for(i=0; i<dim; i++) {
		for(j=0; j< dim; j++) {
			a[j][i] = b[i][j];
		  }
	 }
 endtime = dtime();
 ss=endtime-starttime;
// printf("store stride: %5.2fMB/S\n", ITR*dim*dim*sizeof(double)*1e-6/(endtime-starttime));
 starttime = dtime();
   for(k=0; k<ITR; k++)
	for(j=0; j< dim; j++) {
		for(i=0; i<dim; i++) {
			a[j][i] = b[i][j];
		  }
	 }
 endtime = dtime();
 ls=endtime-starttime;
// printf("load stride: %5.2fMB/S\n", ITR*dim*dim*sizeof(double)*1e-6/(endtime-starttime));
printf("%d %5.2f %5.2f\n", dim,
 ITR*dim*dim*sizeof(double)*1e-6/ss,
 ITR*dim*dim*sizeof(double)*1e-6/ls
);
 exit(0);
}