level1.c

shpf 1.9一个并行编译器

//       PROGRAM MAIN
//       REAL X(2:99), Y(100)
// !HPF$ PROCESSORS P(4)
// !HPF$ TEMPLATE TX(400),TY(-200:199)
// !HPF$ DISTRIBUTE TX(BLOCK) ONTO P
// !HPF$ DISTRIBUTE TY(BLOCK) ONTO P
// !HPF$ ALIGN X(i) WITH TX(2*i+1)   
// !HPF$ ALIGN Y(i) WITH TY(3*i-150)
//
//       FORALL (i=2:99) X(i) = 0
//       FORALL (i=1:100) Y(i) = i
//
//       FORALL (i=2:99) X(i) = Y(i+1) + Y(i-1)
//
//       PRINT *, (X(i),i=2,99)       
//       END


#include "ad++.h"
#include "admacros.h"

void main(int argc, char* argv[]) {
  AdlibInit(argc, argv) ;

  Procs1 p(4) ;

  on(p) {

    Range itx(400, p.dim(0), BLK) ;
    Range ix = itx.subrng(98, 4, 2) ;
    Array1<float> x(ix) ;

    Range ity(400, p.dim(0), BLK) ;
    Range iy = ity.subrng(100, 56, 3) ;
    Array1<float> y(iy) ;

    where(ix)
      x(ix) = 0 ;
    erewh(ix) ;

    where(iy)
      y(iy) = iy + 1 ;
    erewh(iy) ;

    Array1<float> tmp1(ix), tmp2(ix) ;
    remap(tmp1, y.sect(iy.subrng(98, 2))) ;
    remap(tmp2, y.sect(iy.subrng(98, 0))) ;

    where(ix)
      x(ix) = tmp1(ix) + tmp2(ix) ;
    erewh(ix) ;

    aprintf("%12.2f\n", x) ;

  } no(p) ;

  AdlibFinalize() ;
}