sion_fortran_mpi.c

Go to the documentation of this file.

/****************************************************************************
**  SIONLIB     http://www.fz-juelich.de/jsc/sionlib                       **
*****************************************************************************
**  Copyright (c) 2008-2016                                                **
**  Forschungszentrum Juelich, Juelich Supercomputing Centre               **
**                                                                         **
**  See the file COPYRIGHT in the package base directory for details       **
****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "mpi.h"

/* #define _FORTRANUNDERSCORE
   #define _FORTRANNOCAPS
   #define _FORTRANCAPS
   #define _FORTRANDOUBLEUNDERSCORE */
#include "sion.h"
#include "sion_debug.h"
#include "sion_fd.h"
#include "sion_filedesc.h"
#include "sion_mpi.h"
#include "sion_fortran_mpi.h"

void fsion_paropen_mpi_c(char *fname,
             char *file_mode,
             int  *numFiles,
             MPI_Fint * fgComm,
             MPI_Fint * flComm,
             sion_int64 *chunksize, 
             sion_int32 *fsblksize, 
             int *globalrank, 
             char *newfname,
             int *sid, 
             int  fname_len, 
             int  file_mode_len,
             int  newfname_len) {

  FILE     *fileptr;
  char     *fname_tmp,*fmode_tmp;
  MPI_Comm  cgComm,clComm;


#ifdef SION_DEBUG_BUG_OPENMPI
  {
    int grank;
    MPI_Comm newcomm;
    MPI_Comm_rank(MPI_COMM_WORLD,&grank);
    fprintf(stderr,"fsion_paropen_mpi: grank: %d\n", (int) grank);
    MPI_Comm_dup(MPI_COMM_WORLD,&newcomm);
    MPI_Comm_split(MPI_COMM_WORLD,0,grank,&newcomm);
  }
#endif

  fname_tmp = (char *) malloc((size_t) ((fname_len + 1) * sizeof(char)));
  fmode_tmp = (char *) malloc((size_t) ((file_mode_len + 1) * sizeof(char)));

  /* Copy the strings to the new buffers and pad with nulls */
  strncpy(fname_tmp, fname, fname_len);
  strncpy(fmode_tmp, file_mode, file_mode_len);

  fname_tmp[fname_len] = '\0';
  fmode_tmp[file_mode_len] = '\0';
/*   newfname_tmp[newfname_len] = '\0'; */
/*   strncpy(newfname_tmp, newfname, newfname_len); */

  cgComm = MPI_Comm_f2c(*fgComm);
  clComm = MPI_Comm_f2c(*flComm); 

#ifdef SION_DEBUG_UNUSED
  if (*globalrank == 0){
    if (cgComm == MPI_COMM_WORLD) fprintf(stderr,"COMM_WORLD selected\n");
    fprintf(stderr,"fsion_paropen_mpi: filename_len: %d\n", (int) fname_len);
    fprintf(stderr,"fsion_paropen_mpi: filename: %s\n",fname_tmp);
    /*
      fprintf(stderr,"fsion_paropen_mpi: newfname_len: %d\n",(int) newfname_len);
      fprintf(stderr,"fsion_paropen_mpi: newfname: %s\n",newfname_tmp);
    */
    fprintf(stderr,"fsion_paropen_mpi: filemode_len: %d\n", (int) file_mode_len);
    fprintf(stderr,"fsion_paropen_mpi: filemode: %s\n",fmode_tmp);
    
    fprintf(stderr,"fsion_paropen_mpi: chunksize: %lld\n",*chunksize);
    fprintf(stderr,"fsion_paropen_mpi: fsblksize: %d\n",*fsblksize);
    fprintf(stderr,"fsion_paropen_mpi: numFiles:  %d\n",*numFiles);
  }
#endif

  (*sid) = sion_paropen_mpi(fname_tmp, fmode_tmp, numFiles, cgComm, &clComm, chunksize, fsblksize, globalrank, &fileptr, NULL);

  /* int help=MPI_Comm_c2f(clComm); */
  /* printf("WF: flComm=%d (%x), help=%d\n",*flComm,flComm,help); */
/*   *flComm = MPI_Comm_c2f(clComm); WF */
/*   strncpy(newfname, newfname_tmp, newfname_len); */

  /* Free the used memory */
  free(fname_tmp);
  free(fmode_tmp);

}

void fsion_paropen_multi_mpi_c(char *fname,
                   char *file_mode,
                   int *numFiles,
                   MPI_Fint * fgComm,
                   MPI_Fint * flComm,
                   sion_int64 *chunksize,
                   sion_int32 *fsblksize,
                   int *globalrank, int *sid, char *newfname, int fname_len, int file_mode_len, int newfname_len)
{
  FILE     *fileptr;
  MPI_Comm  cgComm, clComm;
  
  /* Fortran strings are not NULL-terminated => insert \0 at the end */
  char     *fname_tmp, *fmode_tmp, *newfname_tmp;
  
  
  fname_tmp = (char *) malloc((size_t) ((fname_len + 1) * sizeof(char)));
  fmode_tmp = (char *) malloc((size_t) ((file_mode_len + 1) * sizeof(char)));
  newfname_tmp = (char *) malloc((size_t) ((newfname_len + 1) * sizeof(char)));
  /* Copy the strings to the new buffers and pad with nulls */
  strncpy(fname_tmp, fname, fname_len);
  strncpy(fmode_tmp, file_mode, file_mode_len);

  fname_tmp[fname_len] = '\0';
  fmode_tmp[file_mode_len] = '\0';
  newfname_tmp[newfname_len] = '\0';
  strncpy(newfname_tmp, newfname, newfname_len);

  cgComm = MPI_Comm_f2c(*fgComm);
  /*        clComm = MPI_Comm_f2c(*flComm); */

  /*
     #ifdef SION_DEBUG
     if (*globalrank == 0){
     if (cgComm == MPI_COMM_WORLD) fprintf(stderr,"COMM_WORLD selected\n");
     fprintf(stderr,"fsion_paropen_mpi: filename_len: %d\n", fname_len);
     fprintf(stderr,"fsion_paropen_mpi: filename: %s\n",fname_tmp);

     fprintf(stderr,"fsion_paropen_mpi: newfname_len: %d\n", newfname_len);
     fprintf(stderr,"fsion_paropen_mpi: newfname: %s\n",newfname_tmp);

     fprintf(stderr,"fsion_paropen_mpi: filemode_len: %d\n", file_mode_len);
     fprintf(stderr,"fsion_paropen_mpi: filemode: %s\n",fmode_tmp);

     fprintf(stderr,"fsion_paropen_mpi: chunksize: %lld\n",*chunksize);
     fprintf(stderr,"fsion_paropen_mpi: fsblksize: %d\n",*fsblksize);
     fprintf(stderr,"fsion_paropen_mpi: numFiles:  %d\n",*numFiles);
     }
     #endif
   */
  (*sid) = sion_paropen_multi_mpi(fname_tmp, fmode_tmp, numFiles, cgComm, &clComm, chunksize, fsblksize, globalrank, &fileptr, newfname_tmp);

  *flComm = MPI_Comm_c2f(clComm);
  strncpy(newfname, newfname_tmp, newfname_len);

  /* Free the used memory */
  free(fname_tmp);
  free(fmode_tmp);
  free(newfname_tmp);
}

void fsion_parclose_mpi_c(int *sid, int *ierr)
{
  (*ierr) = sion_parclose_mpi(*sid);
}