/* $Header$ */

#include <assert.h>
#include <stdlib.h>
#include <string.h>

#include "cctk.h"
#include "cctk_Parameters.h"

#include "Hyperslab.h"

#if 0
int Hyperslab_GetHyperslab (const cGH *GH,
			    int target_proc,
			    int vindex,
			    int vtimelvl,
                            int hdim,
                            const int global_startpoint [/*vdim*/],
                            const int directions [/*vdim*/],
                            const int lengths [/*hdim*/],
                            const int downsample_ [/*hdim*/],
                            void **hdata,
			    int hsize [/*hdim*/]);
#endif

int Hyperslab_FillHyperslab (const cGH *cctkGH,
			     int target_proc,
			     int vindex,
			     int vtimelevel,
			     int hdim,
			     const int global_startpoint [/*vdim*/],
			     const int directions [/*vdim*/],
			     const int lengths [/*hdim*/],
			     const int downsample [/*hdim*/],
			     int nhdata,
			     void *hdata,
			     int hsize [/*hdim*/])
{
  void *myhdata;
  int mynhdata;
  int vtype;
  int vtypesize;
  int d;
  int ierr;
  myhdata = 0;
  ierr = Hyperslab_GetHyperslab
    (cctkGH, target_proc, vindex, vtimelevel, hdim,
     global_startpoint, directions, lengths, downsample,
     &myhdata, hsize);
  if (myhdata) {
    mynhdata = 1;
    for (d=0; d<hdim; ++d) mynhdata *= hsize[d];
    if (mynhdata > nhdata) {
      CCTK_WARN (1, "Not enough space in hdata array");
      mynhdata = nhdata;
    }
    assert (hdata);
    vtype = CCTK_VarTypeI(vindex);
    assert (vtype>=0);
    vtypesize = CCTK_VarTypeSize(vtype);
    assert (vtypesize>0);
    memcpy (hdata, myhdata, mynhdata * vtypesize);
    free (myhdata);
  }
  return ierr;
}



void CCTK_FCALL CCTK_FNAME(Hyperslab_FillHyperslab)
     (int *ierr,
      const cGH *cctkGH,
      const int *target_proc,
      const int *vindex,
      const int *vtimelevel,
      const int *hdim,
      const int global_startpoint [/*vdim*/],
      const int directions [/*vdim*/],
      const int lengths [/*hdim*/],
      const int downsample [/*hdim*/],
      const int *nhdata,
      void *hdata,
      int hsize [/*hdim*/])
{
  *ierr = Hyperslab_FillHyperslab
    (cctkGH, *target_proc, *vindex, *vtimelevel, *hdim,
     global_startpoint, directions, lengths, downsample,
     *nhdata, hdata, hsize);
}