femwind/hexa_plot.m

   1 function hexa_plot(X1,X2,K)
   2 % hexa_plot(X,K)
   3 % plot matrix K on nodes X
   4 % in
   5 %   X1  size (3,m), node coordinates
   6 %   X2  size (3,n), node coordinates
   7 %   K  size (m,n), matrix
   8 h = ishold;
   9 [~,m]=size(X1);
  10 [~,n]=size(X2);
  11 if any([m,n]~=size(K))
  12     error('wrong sizes')
  13 end
  14 ml=8;
  15 tol=1e-6;
  16 symm = false;
  17 if m==n,
  18     symm = norm(K-K','fro')<tol;
  19 end
  20 for j=1:n
  21     % plot3(X(1,j),X(2,j),X(3,j),'o','Linewidth',ml); hold on
  22     for i=1:m
  23         for k=1:3,xij{k}=[X1(k,i),X2(k,j)];end
  24         if symm,
  25             s = K(i,j)/sqrt(K(i,i)*K(j,j));
  26         else
  27             s = K(i,j)/norm([K(:,j);K(i,:)'],inf);
  28         end
  29         if s>tol,
  30             plot3(xij{:},'--k','Linewidth',s*ml); hold on
  31         elseif s<-tol
  32             plot3(xij{:},'--r','Linewidth',-s*ml); hold on
  33         end
  34     end
  35 end
  36 if ~h, hold off, end
  37 grid on
  38 end