#include <med.h>
#define MESGERR 1
#include "med_utils.h"
void generateFullIDatas(const int myrank, const int lastrank, const int sizeoftype,
                        const med_storage_mode profilemode, const med_size profilesize, const med_int * const profilearray,
                        const med_size start, const med_size stride, const med_size count, const med_size blocksize, const med_size lastblocksize,
                        const int nentities, const int nvaluesperentity, const int nconstituentpervalue,
                        med_float ** valuesarray ) {
  med_size _start=start-1,_blockstart = 0,_blocksize=blocksize,_allblocksize=0,_index=0;
  med_int  (*_profilearrayfunc)(int)=0;
  inline med_int _identity  (int i)  { return i; }
  inline med_int _withprofilearray(int i) { return (profilearray[i]-1); }
  int _blocknum=0,_i=0,_j=0,_k=0;
  if (profilesize) {
    if ( profilearray == NULL ) {MESSAGE("Error, profilesize > 0 && profilearray == 0"); }
    MESSAGE("Using a profile...");
    _profilearrayfunc = _withprofilearray;
  } else {
    _profilearrayfunc = _identity;
  }
  switch(profilemode) {
  case MED_GLOBAL_PFLMODE :
    
    
    *valuesarray = (med_float *) calloc(nentities*nvaluesperentity*nconstituentpervalue,sizeoftype);
    for (_blocknum=0; _blocknum< count; ++_blocknum) {
      _blockstart=_blocknum*stride;
      
      if ( (count > 1) && (_blocknum == (count-1) ) && (myrank == lastrank) ) _blocksize=lastblocksize;
      
      for (_i=0; _i<_blocksize; ++_i)
        for (_j=0; _j < nvaluesperentity; ++_j)
          for (_k=0; _k < nconstituentpervalue; ++_k) {
            _index = _profilearrayfunc(_start+_blockstart+_i)*nvaluesperentity*nconstituentpervalue
              +_j*nconstituentpervalue+_k;
            (*valuesarray)[_index]= (myrank+1)*1000+_blocknum*100+_i+0.1*_j+0.01*_k;
            
            
          }
    }
    break;
  case MED_COMPACT_PFLMODE :
    
    if ( (myrank == lastrank) ) _allblocksize=blocksize*count+lastblocksize; else _allblocksize = blocksize*count;
    *valuesarray = (med_float *) calloc(_allblocksize*nvaluesperentity*nconstituentpervalue,sizeoftype);
    _index = 0;
    for (_blocknum=0; _blocknum< count; ++_blocknum) {
      if ( (count > 1) && (_blocknum == (count-1) ) && (myrank == lastrank) ) _blocksize=lastblocksize;
      for (_i=0; _i<_blocksize; ++_i)
        for (_j=0; _j < nvaluesperentity; ++_j)
          for (_k=0; _k < nconstituentpervalue; ++_k) {
            (*valuesarray)[_index]= (myrank+1)*1000+_blocknum*100+_i+0.1*_j+0.01*_k;
            
            
            _index++;
          }
    }
    break;
  default:
    break;
  }
}
void generateNoIDatas(const int myrank, const int lastrank, const int sizeoftype,
                      const med_storage_mode storagemode, const med_size profilearraysize, const med_int * const profilearray,
                      const med_size start, const med_size stride, const med_size count, const med_size blocksize, const med_size lastblocksize,
                      const int nentities, const int nvaluesperentity, const int nconstituentpervalue,
                      med_float ** valuesarray ) {
  med_size _start=start-1,_blockstart = 0,_blocksize=blocksize,_allblocksize=0,_index=0,_dim=0;
  med_int  (*_profilearrayfunc)(int)=0;
  inline med_int _identity  (int i)  { return i; }
  inline med_int _withprofilearray(int i) { return (profilearray[i]-1); }
  int _blocknum=0,_i=0,_j=0,_k=0;
  if (profilearraysize) {
    MESSAGE("Using a profile...");
    if ( profilearray == NULL ) {MESSAGE("Error, profilearraysize > 0 && profilearray == 0"); }
    _profilearrayfunc = _withprofilearray;
  } else {
    _profilearrayfunc = _identity;
  }
  switch(storagemode) {
  case MED_GLOBAL_PFLMODE :
    
    *valuesarray = (med_float *) calloc(nentities*nvaluesperentity*nconstituentpervalue,sizeoftype);
    for (_dim=0; _dim< nconstituentpervalue; ++_dim) {
      _blocksize = blocksize;
      for (_blocknum=0; _blocknum< count; ++_blocknum) {
        _blockstart=_blocknum*stride;
        if ( (count > 1) && (_blocknum == (count-1) ) && (myrank == lastrank) ) _blocksize=lastblocksize;
        for (_i=0; _i<_blocksize; ++_i)
          for (_j=0; _j < nvaluesperentity; ++_j) {
            _index = ( _dim*nentities
                       +_profilearrayfunc(_start+_blockstart+_i) )
              *nvaluesperentity + _j;
            (*valuesarray)[_index]= (myrank+1)*1000+_blocknum*100+_i+0.1*_j+0.01*_dim;
            
            
          }
      }
    }
    break;
  case MED_COMPACT_PFLMODE :
    if ( (myrank == lastrank) ) _allblocksize=blocksize*count+lastblocksize; else _allblocksize = blocksize*count;
    *valuesarray = (med_float *) calloc(_allblocksize*nvaluesperentity*nconstituentpervalue,sizeoftype);
    _index = 0;
    for (_dim=0; _dim< nconstituentpervalue; ++_dim) {
      _blocksize = blocksize;
      for (_blocknum=0; _blocknum< count; ++_blocknum) {
        if ( (count > 1) && (_blocknum == (count-1) ) && (myrank == lastrank) ) _blocksize=lastblocksize;
        for (_i=0; _i<_blocksize; ++_i)
          for (_j=0; _j < nvaluesperentity; ++_j) {
            (*valuesarray)[_index]= (myrank+1)*1000+_blocknum*100+_i+0.1*_j+0.01*_dim;
            
            
            _index++;
          }
      }
    }
    break;
  default:
    break;
  }
}