quicwind/saddlepoint/plot_conditions.m

   1 function plot_conditions(n,v,C,D,desc)
   2
   3 figure, h = axes;
   4 plot(C*v,'Color','k'),
   5 hold on
   6 nc1 = (n(1)-1)*n(2)*n(3);
   7 line([nc1 nc1],get(h,'YLim'),'Color',[1 0 0],'LineWidth',2)
   8 nc2 = nc1 + n(1)*(n(2)-1)*n(3);
   9 line([nc2 nc2],get(h,'YLim'),'Color',[0 1 0],'LineWidth',2)
  10 nc3 = nc2 + n(1)*n(2)*(n(3)-1);
  11 line([nc3 nc3],get(h,'YLim'),'Color',[0 0 1],'LineWidth',2)
  12 hold off
  13 title(sprintf('Continuity conditions %s',desc))
  14
  15 figure,
  16 plot(D*v)
  17 title(sprintf('Divergence-free condition %s',desc))
  18
  19 end