1 #include <mpi.h> 2 #include <stdio.h> 3 #include <stdlib.h> 4 #include <unistd.h> 5 #include <cassert> 6 #include "phIO.h" 7 #include "phstream.h" //for makeRStream and makeGRStream 8 #include "syncio.h" 9 #include "posixio.h" 10 #include "streamio.h" 11 12 int main(int argc, char* argv[]) { 13 MPI_Init(&argc,&argv); 14 int rank; 15 MPI_Comm_rank(MPI_COMM_WORLD, &rank); 16 int size; 17 MPI_Comm_size(MPI_COMM_WORLD, &size); 18 if( argc != 2 ) { 19 fprintf(stderr, "Usage: %s <numSyncFiles>\n",argv[0]); 20 MPI_Finalize(); 21 return 1; 22 } 23 const char* phrase = "number of fishes"; 24 const char* type = "double"; 25 const char* iotype = "binary"; 26 int zero = 0; 27 int one = 1; 28 int numFish = 0; 29 double fishWeight = 1.23; 30 int nfiles = atoi(argv[1]); 31 int ppf = size/nfiles; 32 const char* filename[2] = {"water-dat.", "water.dat."}; 33 rstream rs = makeRStream(); 34 grstream grs = makeGRStream(); 35 phio_fp file[3]; 36 syncio_setup_write(nfiles, one, ppf, &(file[0])); 37 posixio_setup(&(file[1]), 'w'); 38 streamio_setup_write(&(file[2]), rs); 39 for(int i=0; i<3; i++) { 40 phio_openfile(filename[i], file[i]); 41 phio_writeheader(file[i], phrase, &zero, &one, &zero, type, iotype); 42 phio_writedatablock(file[i], phrase, &fishWeight, &zero, type, iotype); 43 phio_closefile(file[i]); 44 } 45 syncio_setup_read(nfiles, &(file[0])); 46 posixio_setup(&(file[1]), 'r'); 47 streamio_read_r(&(file[2]), rs); 48 for(int i=0; i<3; i++) { 49 phio_openfile(filename[i], file[i]); 50 phio_readheader(file[i], phrase, &numFish, &one, type, iotype); 51 assert(!numFish); 52 phio_readdatablock(file[i], phrase, &fishWeight, &numFish, type, iotype); 53 assert(fishWeight == 1.23); 54 phio_closefile(file[i]); 55 } 56 MPI_Finalize(); 57 return 0; 58 } 59