1b184aa6dSCameron Smith #include <mpi.h> 2b184aa6dSCameron Smith #include <stdio.h> 3b184aa6dSCameron Smith #include <stdlib.h> 4b184aa6dSCameron Smith #include <unistd.h> 5b184aa6dSCameron Smith #include "phIO.h" 6b184aa6dSCameron Smith #include "syncio.h" 7b184aa6dSCameron Smith #include "posixio.h" 81a21b96dSCameron Smith #include <assert.h> 9b184aa6dSCameron Smith 10b184aa6dSCameron Smith int main(int argc, char* argv[]) { 11b184aa6dSCameron Smith MPI_Init(&argc,&argv); 12b184aa6dSCameron Smith int rank; 13b184aa6dSCameron Smith MPI_Comm_rank(MPI_COMM_WORLD, &rank); 14b184aa6dSCameron Smith const char* iotype = "binary"; 15b184aa6dSCameron Smith int seven = 7; 16b184aa6dSCameron Smith int headerData[7] = {0,0,0,0,0,0,0}; 17b184aa6dSCameron Smith int blocksRead = 0; 18b184aa6dSCameron Smith phio_fp file; 19b184aa6dSCameron Smith posixio_setup(&file, 'r'); 20b184aa6dSCameron Smith phio_openfile("geombc.dat.", file); 211a21b96dSCameron Smith for(int i=0;i<2;i++) { 22b184aa6dSCameron Smith phio_readheader(file, "connectivity interior", 23b184aa6dSCameron Smith headerData, &seven, "integer", iotype); 241a21b96dSCameron Smith assert(headerData[0] > 0 && headerData[3] > 0); 251a21b96dSCameron Smith int size = headerData[0]*headerData[3]; /* neltp*nshl */ 261a21b96dSCameron Smith int* vals = (int*) calloc(size,sizeof(int)); 271a21b96dSCameron Smith phio_readdatablock(file,"connectivity interior",vals,&size,"integer",iotype); 281a21b96dSCameron Smith free(vals); 29b184aa6dSCameron Smith blocksRead += (headerData[0] > 0); 301a21b96dSCameron Smith } 31b184aa6dSCameron Smith phio_closefile(file); 32b184aa6dSCameron Smith MPI_Finalize(); 33*310e2401SCameron Smith return !(blocksRead==2); 34b184aa6dSCameron Smith } 35