for(Idim=1; Idim <= Nlocal; Idim++) {
vecout[Idim] = 0.5*(vec[Idim-1] + vec[Idim+1]);
}
行列をベクトルに掛ける。
ここで vec[0] と vec[Nlocal+1] は隣の CPU から受信しなくてはならない。同時に vec[1] と vec[Nlocal] を隣の CPU へ送信しなくてはならない。