femwind/plot_mat.m

   1 function plot_mat(X1,X2,M,j)
   2 % plot_mat(X1,X2,M,i1,i2)
   3 % in:
   4 %   X1, X2  coordinates, cell size 3
   5 %   M       matrix
   6 %   j       index triple of the column of M to display
   7 % example: plot_mat(X,XC,P,[5,5,1])
   8
   9 m=size(X1{1});
  10 mm=prod(m);
  11 n=size(X2{2});
  12 nn=prod(n);
  13 if any(size(M)~=[mm,nn])
  14     error('wrong matrix M size')
  15 end
  16 if ~isvector(j) | length(j) ~= 3,
  17     error('j must be vector length 3')
  18 end
  19 XX2=zeros(3,1);
  20 for k=1:3,
  21     XX2(k,1)=[X2{k}(j(1),j(2),j(3))];
  22 end
  23 jx = sub2ind(n,j(1),j(2),j(3));
  24 ii=find(M(:,jx));
  25 mm=length(ii);
  26 XX1=zeros(3,mm);
  27 MM=zeros(mm,1);
  28 for i=1:mm
  29     ix = ii(i);
  30     [i1,i2,i3]=ind2sub(m,ix);
  31     for k=1:3,
  32         XX1(k,i)=[X1{k}(i1,i2,i3)];
  33     end
  34     MM(i,1)=M(ix,jx);
  35 end
  36 hexa_plot(XX1,XX2,MM)