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
main(int argc,char * argv[])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