quadr_nonlin_script.m

function quadr_nonlin_script(step)
%function quadr_nonlin_script(step)
%
% step 1: detailed simulation and plot
% step 10: demo_rb_gui 
%
% B. Haasdonk 8.12.2011

if nargin == 0
  step = 1;
end;

switch step
 case 1 
  model = quadratic_poisson_model;
  model_data = gen_model_data(model);
  model = set_mu(model,[0.1,1]);
  sim_data = detailed_simulation(model, model_data);
  p = plot_sim_data(model, model_data, sim_data,[]);

 case 2 % demo_rb_gui of detailed_simulation
  % setting rb-pointers to dummy-functions, such that
  % rb simulation does simply do a full simulation:
  model = quadratic_poisson_model;
  model.gen_detailed_data = @(model,model_data) model_data;
  model.gen_reduced_data  = @(model,detailed_data) detailed_data;
  model.rb_simulation  = model.detailed_simulation;
  model.rb_reconstruction = ...
     @(model,detailed_data,rb_sim_data) ...
    rb_sim_data;
  % end of setting pointers.
  model_data = gen_model_data(model);
  detailed_data = gen_detailed_data(model, model_data);
  demo_rb_gui(model, detailed_data);
   
 case 3 % simple reduced basis generation and plot
  model = quadratic_poisson_model;
  model_data = gen_model_data(model);
  detailed_data = gen_detailed_data(model, model_data);
  plot(model_data.grid.X(2:end-1),detailed_data.RB)
  title(['Reduced basis with ',...
         num2str(size(detailed_data.RB,2)),...
         ' basis vectors']);
  
 case 4 % reduced simulation and error for single parameter  
  
  model = quadratic_poisson_model;
  model_data = gen_model_data(model);
  detailed_data = gen_detailed_data(model, model_data);
  reduced_data = gen_reduced_data(model, detailed_data);
  mu =[1,1];
  model = set_mu(model,mu);
  rb_sim_data = rb_simulation(model, reduced_data);
  rb_sim_data = rb_reconstruction(model, detailed_data, rb_sim_data)
  sim_data = detailed_simulation(model, model_data);
  plot_params = [];
  err = [];
  err.u = sim_data.u - rb_sim_data.u;
  plot_params.title = ['error u-uN for mu=',num2str(mu)];
  plot_sim_data(model,detailed_data,err,plot_params);
  figure;
  plot_params.title = ['uN and u for mu=',num2str(mu)];
  plot_params.legend ={'u','u_N'};
  joint = [];
  joint.u = [rb_sim_data.u, sim_data.u];  
  plot_sim_data(model,detailed_data,joint,plot_params);
  
 case 5 % rb_demo_gui 
  model = quadratic_poisson_model;
  model_data = gen_model_data(model);
  detailed_data = gen_detailed_data(model, model_data);
  demo_rb_gui(model, detailed_data);

 otherwise
  error('step unknown');
end;