#include <stdio.h>
#include "mpi.h"

#define SIZE 1000
char filename[256]="/mnt/fs3/devezep/VIEW_TEST";
int i, j, myid, dsp[3], lng[3], buffer[SIZE];
MPI_Status status;
MPI_File fh;
MPI_Datatype filetype;

int main(int argc, char **argv) {

 MPI_Init(&argc, &argv);
 MPI_Errhandler_set(MPI_COMM_WORLD, MPI_ERRORS_ARE_FATAL);
 MPI_Comm_rank(MPI_COMM_WORLD, &myid);
 
 if (!myid) {
    MPI_File_open(MPI_COMM_SELF, filename, MPI_MODE_CREATE | MPI_MODE_RDWR , MPI_INFO_NULL, &fh);
    MPI_File_set_errhandler(fh, MPI_ERRORS_ARE_FATAL);

    for (i=0; i<SIZE; i++) buffer[i] = i;
    MPI_File_write(fh, buffer, SIZE, MPI_INTEGER, &status);
    MPI_File_close(&fh);

// loop for test :
    for (i=2 ; i<4 ; i++) {
       lng[0]= 0;
       dsp[0]= 0;
       lng[1]= 1;
       dsp[1]= 1;
       lng[2]= 0;
       dsp[2]= i;

       MPI_Type_indexed(3, lng, dsp, MPI_INTEGER, &filetype);
       MPI_Type_commit(&filetype);
       MPI_File_open(MPI_COMM_SELF, filename, MPI_MODE_RDONLY , MPI_INFO_NULL, &fh);
       MPI_File_set_view(fh, 0, MPI_INTEGER, filetype,"native", MPI_INFO_NULL);
       MPI_File_read(fh, buffer, 5, MPI_INTEGER, &status);

       printf("data read with view for i=%d => lng[0]=%d dsp[0]=%d lng[1]=%d dsp[1]=%d lng[2]=%d dsp[2]=%d :",
             i, lng[0], dsp[0], lng[1], dsp[1], lng[2], dsp[2]);
       for (j=0 ; j<5 ; j++) printf("%x ", buffer[j]);
       if (buffer[1] != (i+2)) printf("\n =======>  test KO for i=%d : buffer[1]=%d instead of %d \n", i, buffer[1], i+2);
       else printf("\n =======> test OK for i=%d\n", i);
       MPI_Type_free(&filetype);
       MPI_File_close(&fh);
    }
 }
 MPI_Barrier(MPI_COMM_WORLD);
 MPI_Finalize();
}