#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;
}
}