1.13.10/doc/gen__burgers__movie_8m_source.html

% small script generating a movie of simulation sequences

%

% the problem is, that powerpoint cannot properly play these clips

% within a slide. So No compression is used to generate the avi.

% Then, for powerpoint-includion, use Programme/avi2mpeg to generate a

% mpeg movie


% Bernard Haasdonk 25.5.2008


load burgers_fv_detailed_data


movie = 5

text_sequence = {''};


switch movie

  case 1

   outputfn = 'rb_burgers1.avi';

   mu_sequence = {[1,0]};

  case 2

   outputfn = 'rb_burgers2.avi';

   mu_sequence = {[1,-pi/4]};

  case 3

   outputfn = 'rb_burgers3.avi';

   mu_sequence = {[-1,0]};

  case 4

   outputfn = 'rb_burgers4.avi';

   mu_sequence = {[-1,-pi/4]};

  case 5

   outputfn = 'simtech.avi';

   mu_sequence = {[1,0,0]};

   load(fullfile(rbmatlabresult,'simtech_detailed'));

   detailed_data.QM = [];

 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 = [-1 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


for i = 1:length(mu_sequence)


  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);


  for j = 1:size(Uappr,2)

    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 = getframe(gca);

    pause(0.01);

    mov = addframe(mov,F);

  end;

end;


mov = close(mov);


%| \docupdate