From 3657f2303dc96818f62a7ed01a34fb7cc3143f56 Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@t002.ph.ed.ac.uk>
Date: Sun, 10 May 2015 15:18:04 +0100
Subject: [PATCH] ET ready benchmark with bytes counted assuming loop
 interchange

---
 benchmarks/Grid_memory_bandwidth.cc | 14 ++++++--------
 1 file changed, 6 insertions(+), 8 deletions(-)

diff --git a/benchmarks/Grid_memory_bandwidth.cc b/benchmarks/Grid_memory_bandwidth.cc
index 3254e95b..fb6db6d8 100644
--- a/benchmarks/Grid_memory_bandwidth.cc
+++ b/benchmarks/Grid_memory_bandwidth.cc
@@ -38,7 +38,6 @@ int main (int argc, char ** argv)
 
       double start=usecond();
       for(int i=0;i<Nloop;i++){
-	//	z=a*x+y;
 	//   inline void axpy(Lattice<vobj> &ret,double a,const Lattice<vobj> &lhs,const Lattice<vobj> &rhs){
 	axpy(z,a,x,y);
       }
@@ -53,7 +52,7 @@ int main (int argc, char ** argv)
   std::cout << "===================================================================================================="<<std::endl;
   std::cout << "= Benchmarking a*x + y bandwidth"<<std::endl;
   std::cout << "===================================================================================================="<<std::endl;
-  std::cout << "  L  "<<"\t\t"<<"bytes"<<"\t\t"<<"GB/s"<<"\t\tGB/s (eff)"<<std::endl;
+  std::cout << "  L  "<<"\t\t"<<"bytes"<<"\t\t"<<"GB/s"<<std::endl;
   std::cout << "----------------------------------------------------------"<<std::endl;
 
   
@@ -73,14 +72,13 @@ int main (int argc, char ** argv)
 
       double start=usecond();
       for(int i=0;i<Nloop;i++){
-	z=a*x+y;
+	z=a*x-y;
       }
       double stop=usecond();
       double time = (stop-start)/Nloop*1000;
-      
-      double effbytes=3*lat*lat*lat*lat*Nvec*sizeof(Real);
-      double bytes=5*lat*lat*lat*lat*Nvec*sizeof(Real);
-      std::cout<<std::setprecision(2) << lat<<"\t\t"<<bytes<<"\t\t"<<bytes/time<<"\t\t"<<effbytes/time<<std::endl;
+     
+      double bytes=3*lat*lat*lat*lat*Nvec*sizeof(Real);
+      std::cout<<std::setprecision(2) << lat<<"\t\t"<<bytes<<"\t\t"<<bytes/time<<std::endl;
 
     }
 
@@ -105,7 +103,7 @@ int main (int argc, char ** argv)
 
       double start=usecond();
       for(int i=0;i<Nloop;i++){
-	x=a*z;
+	z=a*x;
       }
       double stop=usecond();
       double time = (stop-start)/Nloop*1000;