Translating HPF program with inherited mapping
FORALL (I = 1:50) D(I) = 1.0 * I
void init(float* d, DAD* d_dad) {
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;
d [s * (b.sub_bas + b.sub_stp * l)] = 1.0 * i;