xref: /phasta/phSolver/common/test/phIOwriteFields.cc (revision 2df2e51ff945c62919bcc7c5aa1185879485db34)

#include <mpi.h>
#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <cassert>
#include "phiostats.h"
#include "phIO.h"
#include "phstream.h" //for makeRStream and makeGRStream
#include "syncio.h"
#include "posixio.h"
#include "streamio.h"

int main(int argc, char* argv[]) {
  MPI_Init(&argc,&argv);
  int rank;
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  int size;
  MPI_Comm_size(MPI_COMM_WORLD, &size);
  if( argc != 2 ) {
    fprintf(stderr, "Usage: %s <numSyncFiles>\n",argv[0]);
    MPI_Finalize();
    return 1;
  }

  const char* phrase = "number of fishes";
  const char* type = "double";
  const char* iotype = "binary";
  double blockArray[4] = {2.0,4.0,8.0,16.0};
  int blockEntries = 4;
  int one = 1;
  int zero = 0;
  int numFish = 0;
  double fishWeight = 1.23;
  int nfiles = 1;
  int ppf = size/nfiles;
  const char* filename[3] = {"water-dat.", "water.dat.", "water."};
  rstream rs = makeRStream();
     phio_fp file[3];
  const char* modes[3]={"syncio", "posixio", "streamio"};
  fprintf(stderr,"nfiles %d\n", nfiles);
  fprintf(stderr,"ppf %d\n", ppf);
  syncio_setup_write(nfiles, one, ppf, &(file[0]));
  posixio_setup(&(file[1]), 'w');
  streamio_setup_r(&(file[2]), rs, 'w');//Check this _r? work
  fprintf(stderr, "%s\n" ,"Outside loop 1.0");
  for(int i=0; i<3; i++) {
    fprintf(stderr, "%s\n" ,"Within the i loop");
    if(!rank) fprintf(stderr, "%s\n", modes[i]);
    fprintf(stderr, "%s\n" ,"Before phastaio");
    phastaio_initStats();
    fprintf(stderr, "Opening files with %s\n", filename[i]);
    phio_openfile(filename[i], file[i]);
    char str [50];
    int n;
   for (int j = 0; j < 2 ; j++) {
      fprintf(stderr,"%s\n", "Inside loop");
      n = sprintf(str, " Number of times %d ", j);
      assert(n);
      fprintf(stderr,"str \'%s\'\n", str);
      fprintf(stderr,"Printing the int zero %d\n", zero);
      fprintf(stderr,"Printing the int one %d\n", one);
      fprintf(stderr,"blockentries %d\n",blockEntries);
      fprintf(stderr,"Printing the const char type %s\n", type);
      fprintf(stderr,"Printing the const char iotype %s\n", iotype);
      fprintf(stderr,"Calling writeheader\n");
      phio_writeheader(file[i], str, &blockEntries, &one, &blockEntries, "integer", iotype);
      fprintf(stderr,"Done calling writeheader\n");
      fprintf(stderr,"Calling writedatablock\n");
      phio_writedatablock(file[i], str, blockArray, &blockEntries, type, iotype);
      fprintf(stderr,"Done Calling writedatablock\n");
    }
    phio_closefile(file[i]);
    phastaio_printStats();
  }
  syncio_setup_read(nfiles, &(file[0]));
  posixio_setup(&(file[1]), 'r');
  streamio_setup_r(&(file[2]), rs, 'r');
  for(int i=0; i<3; i++) {
    if(!rank) fprintf(stderr, "%s\n", modes[i]);
    phastaio_initStats();
    phio_openfile(filename[i], file[i]);
    //Str was added
   // const char* str = "Number of times "+ nfiles ;
   // for (int j = 0; j < nfiles ; j++) {
    phio_readheader(file[i], phrase, &numFish, &one, type, iotype);
    //  phio_readheader(file[i], str, &numFish, &one, type, iotype);
     //Changing argument from file[i] to file[j]
    assert(!numFish);
    phio_readdatablock(file[i], phrase, &fishWeight, &numFish, type, iotype);
    //  phio_readdatablock(file[i], str, &fishWeight, &numFish, type, iotype);
    assert(fishWeight == 1.23);
    phio_closefile(file[i]);
    phastaio_printStats();
  }
  MPI_Finalize();
  return 0;
}