fft/poisson_fft2_test.m

   1 function poisson_test\r
   2 \r
   3 n=[111,155]; % gridpoints in each direction excluding boundaries\r
   4 h=rand(1,2); % meshstep\r
   5 \r
   6 F=ones(n);\r
   7 F=rand(n);\r
   8 \r
   9 U=poisson_fft2(F,h,-1);\r
  10 FF=mlap(U,h);\r
  11 R=F-FF;\r
  12 err=norm(R,inf)\r
  13 F2=poisson_fft2(U,h,1/2);\r
  14 F2=poisson_fft2(F2,h,1/2);\r
  15 R=F-F2;\r
  16 err=norm(R,inf)\r
  17 \r
  18 return