global.openmploops.h

Go to the documentation of this file.

#define OPENMPNUMCHUNKS 100 // when large blocksize, if too small chunks, then overhead is large
#define MINCHUNKSIZE 10 // if too few chunks, then overhead is large (GODMARK: should worry about if <10 iterations? -- not for now)




#if(N1*N2*N3<200) // overhead becomes a problem for too few iterations in loops
#define OPENMPSCHEDULE(arg) static
#define OPENMPCHUNKSIZE(blocksize) (MAX(blocksize/numopenmpthreads,MINCHUNKSIZE))
#else
#define OPENMPSCHEDULE(arg) guided
#define OPENMPCHUNKSIZE(blocksize) (MAX(blocksize/OPENMPNUMCHUNKS,MINCHUNKSIZE))
#endif


#define OPENMPNOVARYSCHEDULE(arg) static

#define OPENMPFULLNOVARYSCHEDULE(arg) static

#define OPENMPVARYENDTIMESCHEDULE(arg) guided

//#define OPENMP3DLOOPBLOCK     for(blockijk=0;blockijk<blocksize;blockijk++)
#define OPENMP3DLOOPBLOCK     for(blockijk=1;blockijk<=blocksize;blockijk++)



#define OPENMP3DLOOPVARSDEFINE int nxsize, nxshift, nysize, nyshift, nzsize, nzshift, blocksize, blockijk


#define OPENMP3DLOOPSETUPSUPERGENCOMP(i,j,k,istart,istop,jstart,jstop,kstart,kstop,di,dj,dk) \
  {nxsize=(  ((istop+SHIFTX1UP) - (istart+SHIFTX1DN))*di +1 );          \
    nxshift=istart+SHIFTX1DN;                                           \
    nysize=( ((jstop+SHIFTX2UP) - (jstart+SHIFTX2DN))*dj +1 );          \
    nyshift=jstart+SHIFTX2DN;                                           \
    nzsize=( ((kstop+SHIFTX3UP) - (kstart+SHIFTX3DN))*dk +1 );          \
    nzshift=kstart+SHIFTX3DN;                                           \
    blocksize=nxsize*nysize*nzsize;}

#define OPENMP3DLOOPSETUPSUPERGEN(istart,istop,jstart,jstop,kstart,kstop,di,dj,dk) \
  {nxsize=(  (istop - istart)*di +1 );                                  \
    nxshift=istart;                                                     \
    nysize=( (jstop - jstart)*dj +1 );                                  \
    nyshift=jstart;                                                     \
    nzsize=( (kstop - kstart)*dk +1 );                                  \
    nzshift=kstart;                                                     \
    blocksize=nxsize*nysize*nzsize;}

#define OPENMP3DLOOPSETUP(istart,istop,jstart,jstop,kstart,kstop)       \
  {nxsize=((istop+SHIFTX1UP) - (istart+SHIFTX1DN) +1 );                 \
    nxshift=istart+SHIFTX1DN;                                           \
    nysize=((jstop+SHIFTX2UP) - (jstart+SHIFTX2DN) +1 );                \
    nyshift=jstart+SHIFTX2DN;                                           \
    nzsize=((kstop+SHIFTX3UP) - (kstart+SHIFTX3DN) +1 );                \
    nzshift=kstart+SHIFTX3DN;                                           \
    blocksize=nxsize*nysize*nzsize;}


#define OPENMP3DLOOPSETUPFULL OPENMP3DLOOPSETUP(-N1BND,N1-1+N1BND,-N2BND,N2-1+N2BND,-N3BND,N3-1+N3BND)

#define OPENMP3DLOOPSETUPFULLINOUT2 OPENMP3DLOOPSETUP(-N1BND+2*SHIFT1,N1-1+N1BND-2*SHIFT1,-N2BND+2*SHIFT2,N2-1+N2BND-2*SHIFT2,-N3BND+2*SHIFT3,N3-1+N3BND-2*SHIFT3)

#define OPENMP3DLOOPSETUPFULLINOUT2DIR1 OPENMP3DLOOPSETUP(-N1BND+2*SHIFT1,N1-1+N1BND-2*SHIFT1,-N2BND,N2-1+N2BND,-N3BND,N3-1+N3BND)

#define OPENMP3DLOOPSETUPFULLINOUT2DIR2 OPENMP3DLOOPSETUP(-N1BND,N1-1+N1BND,-N2BND+2*SHIFT2,N2-1+N2BND-2*SHIFT2,-N3BND,N3-1+N3BND)

#define OPENMP3DLOOPSETUPFULLINOUT2DIR3 OPENMP3DLOOPSETUP(-N1BND,N1-1+N1BND,-N2BND,N2-1+N2BND,-N3BND+2*SHIFT3,N3-1+N3BND-2*SHIFT3)


#define OPENMP3DLOOPSETUPFULLP1 OPENMP3DLOOPSETUP(-N1BND,N1-1+N1BND+SHIFT1,-N2BND,N2-1+N2BND+SHIFT2,-N3BND,N3-1+N3BND+SHIFT3)

#define OPENMP3DLOOPSETUPFULLP2 OPENMP3DLOOPSETUP(-N1BND-SHIFT1,N1-1+N1BND+SHIFT1*2,-N2BND-SHIFT2,N2-1+N2BND+SHIFT2*2,-N3BND-SHIFT3,N3-1+N3BND+SHIFT3*2)

#define OPENMP3DLOOPSETUPFULLP1EXCEPTX2 OPENMP3DLOOPSETUP(-N1BND,N1-1+N1BND+SHIFT1,-N2BND,N2-1+N2BND,-N3BND,N3-1+N3BND+SHIFT3)

#define OPENMP3DLOOPSETUPZLOOP OPENMP3DLOOPSETUP(0,N1-1,0,N2-1,0,N3-1)



#define OPENMPBCLOOPSETUP(istart,istop,jstart,jstop,kstart,kstop)       \
  {nxsize=((istop) - (istart) +1 );                                     \
    nxshift=istart;                                                     \
    nysize=((jstop) - (jstart) +1 );                                    \
    nyshift=jstart;                                                     \
    nzsize=((kstop) - (kstart) +1 );                                    \
    nzshift=kstart;                                                     \
    blocksize=nxsize*nysize*nzsize;}


#define OPENMPBCLOOPVARSDEFINELOOPX1DIR int nxsize, nxshift, nysize, nyshift, nzsize, nzshift, blocksize, blockijk, fooi

#define OPENMPBCLOOPVARSDEFINELOOPX2DIR int nxsize, nxshift, nysize, nyshift, nzsize, nzshift, blocksize, blockijk, fooj

#define OPENMPBCLOOPVARSDEFINELOOPX3DIR int nxsize, nxshift, nysize, nyshift, nzsize, nzshift, blocksize, blockijk, fook

#define OPENMPBCLOOPSETUPLOOPX1DIR OPENMPBCLOOPSETUP(0,0,innormalloop[2],outnormalloop[2],innormalloop[3],outnormalloop[3])

#define OPENMPBCLOOPSETUPLOOPX2DIR OPENMPBCLOOPSETUP(inboundloop[1],outboundloop[1],0,0,innormalloop[3],outnormalloop[3])

#define OPENMPBCLOOPSETUPLOOPX3DIR OPENMPBCLOOPSETUP(inboundloop[1],outboundloop[1],inboundloop[2],outboundloop[2],0,0)

#define OPENMPBCLOOPBLOCK2IJKLOOPX1DIR(j,k) OPENMP3DLOOPBLOCK2IJK(fooi,j,k);

#define OPENMPBCLOOPBLOCK2IJKLOOPX2DIR(i,k) OPENMP3DLOOPBLOCK2IJK(i,fooj,k);

#define OPENMPBCLOOPBLOCK2IJKLOOPX3DIR(i,j) OPENMP3DLOOPBLOCK2IJK(i,j,fook);

#define OPENMPBCLOOPBLOCK OPENMP3DLOOPBLOCK