Translating simple HPF program to C++
!HPF$ DISTRIBUTE A(BLOCK) ONTO P
FORALL (I = 1:50) A(I) = 1.0 * I
BlockRange x(50, p.dim(0));
float* a = new float [x.volume()];
x.block(&b, p.dim(0).crd());
for (int l = 0; l < b.count; l++) {
const int i = b.glb_bas + b.glb_stp * l + 1;
a [b.sub_bas + b.sub_stp * l] = 1.0 * i;