Pheryus

	for (i = 0; i < NGLOB; i += BS_NGLOB) {
actual_size = ((NGLOB - i) >= BS_NGLOB ? BS_NGLOB : (NGLOB - i));
#pragma omp task firstprivate(actual_size, i, bs_nglob, dim) depend(inout: displ[0:bs_nglob][0:dim], veloc[0:bs_nglob][0:dim], accel[0:bs_nglob][0:dim])
update_disp_vel (actual_size, (void *) &displ[i][0], (void *) &veloc[i][0], (void *) &accel[i][0]);
}