gen_nonlin_symmetry_mpg.m

Go to the documentation of this file.

% small script generating a movie of simulation sequences
%
% the problem with avi capturing is, that powerpoint cannot
% properly play these clips. So using package mpgwrite here.
% the first frame is chosen as the snapshot at end time in order 
% to have a useful preview. As the initial dist will be identical
% for different parameters

% Bernard Haasdonk 25.5.2008

load nonlin_symmetry_detailed
params.no_lines = 1;

movie = 2
text_sequence = {''};

switch movie
  case 1
   outputfn = 'rb_nonlin_symmetry1.mpg';
   mu_sequence = {[1]};
  case 2
   outputfn = 'rb_nonlin_symmetry2.mpg';
   mu_sequence = {[2]};
%  case 3
%   outputfn = 'rb_burgers3.mpg';
%   mu_sequence = {[-1,0]};
%  case 4
%   outputfn = 'rb_burgers4.mpg';
%   mu_sequence = {[-1,-pi/4]};
 otherwise
  error('movie number not specified');
end;
%outputfn = 'rb_burgers1.avi';
%mu_sequence = {[1,0],[1,-pi/8],[1,-pi/4]};
%text_sequence = {'\phi = 0', '\phi = -\pi/8','\phi = -\pi/4'};

offline_data = rb_offline_prep(detailed_data,params);
params.N = size(detailed_data.RB,2);
params.M = size(detailed_data.QM{1},2);
%params.M = 45;
%params.vrot_angle = 0;
%params.c_init_lo = 1;
params.clim = [0 1];
%quality = 75;
%compression = 'None';
%compression = 'Cinepak';
%compression = 'Indeo5';
%compression = 'IV50';
%compression = 'Indeo5'; % Indeo3, Cinepak, MSVC, RLE, None
%compression='TSCC';
fps = 15;
%mov = avifile(outputfn,'Fps',fps,'Quality',quality,...
%             'Compression',compression);

%if movie == 5
%  params.clim = [0,1];
%  fps = 25;
%end;
params.show_colorbar = 1

% cinit, vrot

f = [];

for i = 1:length(mu_sequence)
%for i = 1:25
  
  params = params.set_mu(mu_sequence{i},params);
  reduced_data = rb_online_prep(offline_data,params);
  simulation_data = rb_simulation(reduced_data,params);
  Uappr = rb_reconstruction(detailed_data,simulation_data);

  cla;
  plot_element_data(detailed_data.grid,Uappr(:,end),params);
  text(0.1,0.9,text_sequence{i},'FontSize',30,'Color',[0,0,1])
  %    F = getframe(gcf);
  F(1) = getframe(gca);
   
  for j = 1:size(Uappr,2)
%  for j = 1:25
    cla;
    plot_element_data(detailed_data.grid,Uappr(:,j),params);
    text(0.1,0.9,text_sequence{i},'FontSize',30,'Color',[0,0,1])
%    F = getframe(gcf);
     F(j+1) = getframe(gca);
%    F = getframe(gca);
%    pause(0.01);
%    mov = addframe(mov,F);
  end;
end;

C = colormap;
%mov = close(mov);
mpgwrite(F,C,outputfn,...
    [1, 0, 1, 0, 10, 1, 1, 1]);
%keyboard;
disp('finished!');

%| \docupdate