quicwind/mlap3z.m

   1 function y=mlap3z(x,h)
   2 % evaluate minus Laplacian y= -x_11 - x_22 - x_33 on rectangular uniform grid
   3 % assuming boundary values of x zero wrapped around
   4 % in:
   5 %   x   values on 3d grid
   6 n=size(x);
   7 xx=zeros(n+2);
   8 xx(2:n(1)+1,2:n(2)+1,2:n(3)+1)=x;
   9 mid1=2:n(1)+1;
  10 mid2=2:n(2)+1;
  11 mid3=2:n(3)+1;
  12 y=(2*xx(mid1,mid2,mid3)-xx(mid1-1,mid2,mid3)-xx(mid1+1,mid2,mid3))/(h(1)*h(1))...
  13  +(2*xx(mid1,mid2,mid3)-xx(mid1,mid2-1,mid3)-xx(mid1,mid2+1,mid3))/(h(2)*h(2))...
  14  +(2*xx(mid1,mid2,mid3)-xx(mid1,mid2,mid3-1)-xx(mid1,mid2,mid3+1))/(h(3)*h(3));