replay_instant.m

来自「解决麦克斯韦的matlab源文件」· M 代码 · 共 24 行

function replay_instant(mesh,t,times,solv,sole,step,filename)
% Rendering of nodal/edge element solution for
% transient Maxwell problem in a cavity
% at a particular instant t of time

N = length(times);
if (nargin < 6), filename = []; end;
if (nargin < 5), step = 1; end;
if (size(solv,2) ~= N), error('Wrong number of samples in solv'); end
if (size(sole,2) ~= N), error('Wrong number of samples in sole'); end

fprintf('Displaying %d data samples\n',N);

figure; clf; figno = gcf;
idx = find(times == t);
for i=idx 
  plotiterate(mesh,sole(:,i),solv(:,i),times(i),figno);
  if (~isempty(filename))
    epsfile = [filename , '_t' ,  num2str(times(i),3) , '.eps'];
    fprintf('Printing to file %s\n',epsfile);
    print('-depsc','-r600',epsfile);
  end
end