elastic_membrane.m

function res = elastic_membrane(step,params)
%function res = elastic_membrane(params)
%
% function elastic membrane 

% B. Haasdonk 12.3.2012

res = [];

if nargin<1
  step = 4; % demo_rb_gui
end;

model = elastic_membrane_model;
model_data = gen_model_data(model);

switch step
  
 case 1 % detailed simulation with unconstrained fem solver 
  
  %  model = set_mu(model,[1,20]);
  model = set_mu(model,[1,1]);
  %  model = set_mu(model,[1,0]);
  model = set_mu(model,[20,0]);
  model = set_mu(model,[20,1]);
  model = set_mu(model,[0.3,0]);
  model.operators = @fem_operators;
  model.detailed_simulation = @lin_stat_detailed_simulation;
  model.compute_output_functional = 0;
  sim_data = filecache_function(@detailed_simulation,model,model_data);
  sim_data.U = sim_data.uh.dofs(:);
  sim_data.L = zeros(size(sim_data.U,1),1);
  plot_sim_data(model,model_data,sim_data,[]);
  view(-17,34);
 
 case 2 % detailed simulation with constrained fem solver 
  
  model = set_mu(model,[20,1]);
  model.dummy = 2;
%  model = set_mu(model,[0.3,0]);
  sim_data = filecache_function(@detailed_simulation,model,model_data);
  plot_sim_data(model,model_data,sim_data,[]);
  view(-17,34);
  
 case 3  % generate greedy reduced basis

  % rb-estimator with W subbasis for dual reduced basis:
  model.RB_generation_mode = 'greedy-rb-estimator';
  model.supremizer_enrichment = 1; % on/off
  model.enable_error_estimator = 1;
  model.RB_stop_Nmax = 20; % should suffice here!
  model.RB_numintervals = [16,16];
  model.dual_lin_independence_mode = 1; % use W  subbasis 
  
  disp('--------------------------------------------------');
  disp(model.RB_generation_mode);
  detailed_data = filecache_function(@gen_detailed_data,model,model_data);
  save('elastic_membrane_detailed_data',...
       'model','detailed_data');
  max_err_sequence7 = detailed_data.RB_info.max_err_sequence;
  plot(max_err_sequence7);
  set(gca,'Yscale','log');
  title('greedy rb-estimator, with supr.-enrichment, lin-indep=1');
  disp('press enter to continue');
  plot_detailed_data(model,detailed_data,[]);
  keyboard;
  
 case 4 % demo_rb_gui
  load('elastic_membrane_detailed_data') % from step 3
  model.enable_error_estimator = 0;
  plot_params = [];
  plot_params.plot_title = 'elastic membrane'; 
  plot_params.no_lines = 0;
  demo_rb_gui(model,detailed_data,[],plot_params);
  
%  demo_detailed_gui(model,model_data);
  view(-17,34);  

 case 5 % generate animation of demo_rb_gui

  load('elastic_membrane_detailed_data') % from step 3
  model.enable_error_estimator = 0;
  plot_params = [];
  plot_params.plot_title = 'elastic membrane'; 
  plot_params.no_lines = 0;
%  demo_rb_gui(model,detailed_data,[],plot_params);
  
  %  demo_detailed_gui(model,model_data);  
  %  loop_demo_rb_gui(gcf,'mu_1',0.1);
  
  % 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
    
  text_sequence = {''};
  
  outputfn = 'rb_elastic_membrane3.avi';
  mu_sequence = [0.3*ones(1,41); 0:0.025:1];
  
  reduced_data = gen_reduced_data(model,detailed_data);
  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);
  plot_params.show_colorbar = 1;
  
  for i = 1:length(mu_sequence)
    
    model = set_mu(model,mu_sequence(:,i));
    sim_data = rb_simulation(model,reduced_data);
    sim_data = rb_reconstruction(model,detailed_data,sim_data);
    
    cla;
    plot_sim_data(model,model_data,sim_data,plot_params);
    view(-17,34);  
    %    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);
    % -15,-78
  end;
  mov = close(mov);
    
 otherwise
  disp('step unknown')
end;