femwind/check_P.m

   1 function check_P(P,X,XC)
   2 n = size(X{1});
   3 nn=prod(n);
   4 nc=size(XC{1});
   5 YC = XC{3};
   6 for i=nc(3):-1:1
   7     YC(:,:,i)=YC(:,:,i)-YC(:,:,1);
   8 end
   9 for i=1:3;
  10     xmax(i)=max(XC{i}(:));
  11     xmin(i)=min(XC{i}(:));
  12 end
  13 f = @(x1,x2,y)(y-xmin(3))/(xmax(3)-xmin(3));
  14 FC = f(XC{1},XC{2},YC);
  15 F = reshape(P*FC(:),n);
  16 p.no_error_title=1;
  17 p.slice=0.5;
  18 p.err_slice_fig=1;
  19 plot_error_slice(F,[],X,'interpolation',p)
  20 p.err_slice_fig=2;
  21 plot_error_slice(FC,[],XC,'coarse',p)
  22 end