From 4546846f0682f8dd7e327578b9f38234537d0af7 Mon Sep 17 00:00:00 2001
From: Kurt Stockinger <kstockinger@lbl.gov>
Date: Tue, 24 Apr 2007 22:48:53 +0000
Subject: [PATCH] Delete all files are benchmark is done.

---
 test/Bench.c | 265 +++++++++++++++++++++++++++------------------------
 1 file changed, 138 insertions(+), 127 deletions(-)

diff --git a/test/Bench.c b/test/Bench.c
index 446d7ab..b0df5a4 100644
--- a/test/Bench.c
+++ b/test/Bench.c
@@ -22,13 +22,14 @@
 
 /*
 
-  bench <nParticles>
+bench <nParticles>
 
 
 */
 
 
 int main(int argc,char *argv[]){
+  printf("Start benchmarks...\n");
   MPI_Info info;
   int nprocs,rank;
   int trial;
@@ -48,6 +49,7 @@ int main(int argc,char *argv[]){
   H5PartFile *f;
 #endif
   char newfilename[128];
+  char lastfilename[128];
   FILE *fd;
   MPI_File file;
   MPI_Offset foffset;
@@ -74,7 +76,7 @@ int main(int argc,char *argv[]){
 
   
   /* printf("about to call create subarray with nparticles=%u localnp=%u offset=%u\n",
-	nparticles,localnp,offset); */
+     nparticles,localnp,offset); */
   MPI_Type_create_subarray(1, /* rank */
 			   &nparticles, /* size of the global array */
 			   &localnp, /* size of my local chunk */
@@ -86,148 +88,157 @@ int main(int argc,char *argv[]){
   MPI_Info_create(&info);
 
   if(rank==0) printf("Nprocs=%u Particles=%u*6attribs*sizeof(double) Particles/proc=%u Nsteps=%u Ntrials=%u\n",
-	nprocs,nparticles,localnp,NSTEPS,NTRIALS);
+		     nprocs,nparticles,localnp,NSTEPS,NTRIALS);
 	
 
   for(trial=0;trial<NTRIALS;trial++){
-  if(rank==0) printf("---------------------- Trial %u of %u ---------------------\n",trial+1,NTRIALS);
+    if(rank==0) printf("---------------------- Trial %u of %u ---------------------\n",trial+1,NTRIALS);
 
 
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  sprintf(filename,"%s.%u.mpio.dat",FILENAME,nprocs);
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    sprintf(filename,"%s.%u.mpio.dat",FILENAME,nprocs);
+    sprintf(lastfilename,"%s.%u.mpio.dat",FILENAME,nprocs);
 
-  if(rank==0) unlink(filename);
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    if(rank==0) unlink(filename);
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
 
-  MPI_File_open(MPI_COMM_WORLD,filename,
-	MPI_MODE_CREATE | MPI_MODE_RDWR,
-	info,&file);
+    MPI_File_open(MPI_COMM_WORLD,filename,
+		  MPI_MODE_CREATE | MPI_MODE_RDWR,
+		  info,&file);
 	
-  MPI_File_set_view(file,0,MPI_DOUBLE,chunktype,"native",info);
-  /* now a barrier to get the start timers roughly synced*/
-  MPI_Barrier(MPI_COMM_WORLD);
-  curtime = starttime = MPI_Wtime();
-  endtime = starttime+5.0*60.0; /* end in 5 minutes */
-  MPI_Bcast(&endtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-  /* must touch the entire array after each write */
-  /* ensures cache-invalidation */
-  foffset=0;
-  i=0;
-  curtime=starttime;
-  for(i=0;i<NSTEPS;i++){
-    int n;
-    MPI_Status status;
-    for(j=0;j<6;j++){
-      /* touch data */
-      for(n=0;n<localnp;n++)
-	(data[j])[n]=(double)rank;
-      /* write to that file */
-      /*  MPI_File_set_view(file,foffset,MPI_DOUBLE,chunktype,"native",info);*/
-      MPI_File_write_at_all(file,
-			    foffset,
-			    data[j],
-			    localnp,
-			    MPI_DOUBLE,&status);
-      foffset+=nparticles/nprocs;
+    MPI_File_set_view(file,0,MPI_DOUBLE,chunktype,"native",info);
+    /* now a barrier to get the start timers roughly synced*/
+    MPI_Barrier(MPI_COMM_WORLD);
+    curtime = starttime = MPI_Wtime();
+    endtime = starttime+5.0*60.0; /* end in 5 minutes */
+    MPI_Bcast(&endtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
+    /* must touch the entire array after each write */
+    /* ensures cache-invalidation */
+    foffset=0;
+    i=0;
+    curtime=starttime;
+    for(i=0;i<NSTEPS;i++){
+      int n;
+      MPI_Status status;
+      for(j=0;j<6;j++){
+	/* touch data */
+	for(n=0;n<localnp;n++)
+	  (data[j])[n]=(double)rank;
+	/* write to that file */
+	/*  MPI_File_set_view(file,foffset,MPI_DOUBLE,chunktype,"native",info);*/
+	MPI_File_write_at_all(file,
+			      foffset,
+			      data[j],
+			      localnp,
+			      MPI_DOUBLE,&status);
+	foffset+=nparticles/nprocs;
+      }
+      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
+      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
+    }
+    MPI_File_close(&file);
+    MPI_Barrier(MPI_COMM_WORLD);
+    endtime=MPI_Wtime();
+    sprintf(filename,"%s.%u.h5.dat",FILENAME,nprocs);
+    /* foffset*=nprocs; if we want total megabytes written */
+    if(rank==0){
+      puts("*");
+      unlink(filename);
+      puts("======================================================");
+      printf("Raw MPI-IO Total Duration %lf seconds, iterations=%u %lf Megabytes written per processor Nprocs= %u \n",
+	     (endtime-starttime),i,((double)foffset)/(1024.0*1024.0),nprocs);
+      printf("Raw MPI-IO Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
+	     (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
+	     (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
+      puts("======================================================");
     }
-    curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-    MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-  }
-  MPI_File_close(&file);
-  MPI_Barrier(MPI_COMM_WORLD);
-  endtime=MPI_Wtime();
-  sprintf(filename,"%s.%u.h5.dat",FILENAME,nprocs);
-  /* foffset*=nprocs; if we want total megabytes written */
-  if(rank==0){
-    puts("*");
-    unlink(filename);
-    puts("======================================================");
-    printf("Raw MPI-IO Total Duration %lf seconds, iterations=%u %lf Megabytes written per processor Nprocs= %u \n",
-	   (endtime-starttime),i,((double)foffset)/(1024.0*1024.0),nprocs);
-    printf("Raw MPI-IO Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	   (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
-	   (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-    puts("======================================================");
-  }
 
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  /* OK, now we do this using H5Part */
-  sprintf(newfilename,"testio%u.%u.dat",rank,nprocs);
-  unlink(newfilename);
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  fd = fopen(newfilename,"w");
-  /* start the timer */
-  starttime=endtime=MPI_Wtime();
-  for(i=0;i<NSTEPS;i++){
-    for(j=0;j<6;j++){
-      /* touch data */
-      for(n=0;n<localnp;n++)
-	(data[j])[n]=(double)rank;
-      fwrite(data[j],sizeof(double),localnp,fd);
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    /* OK, now we do this using H5Part */
+    sprintf(newfilename,"testio%u.%u.dat",rank,nprocs);
+    unlink(newfilename);
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    fd = fopen(newfilename,"w");
+    /* start the timer */
+    starttime=endtime=MPI_Wtime();
+    for(i=0;i<NSTEPS;i++){
+      for(j=0;j<6;j++){
+	/* touch data */
+	for(n=0;n<localnp;n++)
+	  (data[j])[n]=(double)rank;
+	fwrite(data[j],sizeof(double),localnp,fd);
+      }
+      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
+      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
+    }
+    fclose(fd);
+    MPI_Barrier(MPI_COMM_WORLD);
+    endtime=MPI_Wtime();
+    if(rank==0) puts("*");
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    unlink(newfilename);
+    MPI_Barrier(MPI_COMM_WORLD);
+    if(rank==0){
+      puts("======================================================");
+      printf("Raw 1-file-per-proc Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
+	     (endtime-starttime),NSTEPS,((double)foffset)/(1024.0*1024.0),nprocs);
+      printf("Raw 1-file-per-proc Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
+	     (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
+	     (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
+      puts("======================================================");
     }
-    curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-    MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-  }
-  fclose(fd);
-  MPI_Barrier(MPI_COMM_WORLD);
-  endtime=MPI_Wtime();
-  if(rank==0) puts("*");
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  unlink(newfilename);
-  MPI_Barrier(MPI_COMM_WORLD);
-  if(rank==0){
-    puts("======================================================");
-    printf("Raw 1-file-per-proc Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
-	   (endtime-starttime),NSTEPS,((double)foffset)/(1024.0*1024.0),nprocs);
-    printf("Raw 1-file-per-proc Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	   (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
-	   (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-    puts("======================================================");
-  }
 
 #ifndef DISABLE_H5PART
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  /* OK, now we do this using H5Part */
-  f = H5PartOpenFileParallel(filename,H5PART_WRITE,MPI_COMM_WORLD);
-  MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-  /* start the timer */
-  starttime=endtime=MPI_Wtime();
-  H5PartSetNumParticles(f,localnp);
-  for(i=0;i<NSTEPS;i++){
-    for(j=0;j<6;j++){
-      /* touch data */
-      for(n=0;n<localnp;n++)
-	(data[j])[n]=(double)rank;
-    }
-    H5PartSetStep(f,i);
-    H5PartWriteDataFloat64(f,"x",x);
-    H5PartWriteDataFloat64(f,"y",y);
-    H5PartWriteDataFloat64(f,"z",z);
-    H5PartWriteDataFloat64(f,"px",px);
-    H5PartWriteDataFloat64(f,"py",py);
-    H5PartWriteDataFloat64(f,"pz",pz);
+    //printf("H5Part benchmark...");
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    /* OK, now we do this using H5Part */
+    f = H5PartOpenFileParallel(filename,H5PART_WRITE,MPI_COMM_WORLD);
+    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
+    /* start the timer */
+    starttime=endtime=MPI_Wtime();
+    H5PartSetNumParticles(f,localnp);
+    for(i=0;i<NSTEPS;i++){
+      for(j=0;j<6;j++){
+	/* touch data */
+	for(n=0;n<localnp;n++)
+	  (data[j])[n]=(double)rank;
+      }
+      H5PartSetStep(f,i);
+      H5PartWriteDataFloat64(f,"x",x);
+      H5PartWriteDataFloat64(f,"y",y);
+      H5PartWriteDataFloat64(f,"z",z);
+      H5PartWriteDataFloat64(f,"px",px);
+      H5PartWriteDataFloat64(f,"py",py);
+      H5PartWriteDataFloat64(f,"pz",pz);
 
-    curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-    MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-  }
-  H5PartCloseFile(f);
-  MPI_Barrier(MPI_COMM_WORLD);
-  endtime=MPI_Wtime();
-  if(rank==0){
-    puts("*");
-    unlink(filename);
-    puts("======================================================");
-    printf("H5Part Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
-	   (endtime-starttime),NSTEPS,((double)foffset)/(1024.0*1024.0),nprocs);
-    printf("H5Part Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	   (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
-	   (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-    puts("======================================================");
-  }
-  MPI_Barrier(MPI_COMM_WORLD);
+      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
+      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
+    }
+    H5PartCloseFile(f);
+    MPI_Barrier(MPI_COMM_WORLD);
+    endtime=MPI_Wtime();
+    if(rank==0){
+      puts("*");
+      unlink(filename);
+      puts("======================================================");
+      printf("H5Part Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
+	     (endtime-starttime),NSTEPS,((double)foffset)/(1024.0*1024.0),nprocs);
+      printf("H5Part Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
+	     (double)(nprocs*localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),
+	     (double)(localnp*sizeof(double))*((double)NSTEPS)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
+      puts("======================================================");
+    }
+    MPI_Barrier(MPI_COMM_WORLD);
 #endif
   } /* trials */
+
+  if (rank==0) {
+    printf("lastfilename: %s \n", lastfilename);
+    unlink(lastfilename);
+  }
+
   MPI_Finalize();
 
+
   return 0;
 }