untitled3.m

五点差分型多重网格方法：各种插值算子的比较）

u_c_out=[1 2 3 4 5 6 7 8 9];
nx_c=3;
ny_c=3;
IX_c = 1:nx_c; IY_c = 1:ny_c;
%IX_f = 2:(nx_f+1); IY_f = 2:(ny_f+1);
coarse(IX_c,IY_c) = reshape(u_c_out,nx_c,ny_c);
disp(coarse);