rb_init_values_default.m

Go to the documentation of this file.

function a0 = rb_init_values_default(model,detailed_data)
%function a0 = rb_init_values_default(models,detailed_data)
% function computing the reduced basis initial values. If the
% decomposition mode is 'coefficients', the detailed_data are
% superfluous, can (and should for H-independence) be empty.
%

% Bernard Haasdonk 23.7.2006

% determine affine_decomposition_mode as integer  
%model.decomp_mode = decomp_mode;
decomp_mode = model.decomp_mode;

if decomp_mode == 0 % complete: simple projection on RB
  Nmax = size(detailed_data.RB,2);
  u0 = model.init_values_algorithm(model,detailed_data);
  A = detailed_data.W;
  a0 = u0' * A * detailed_data.RB;
  %a0 = feval(params.inner_product_name, ...
  %      detailed_data.RB,u0,detailed_data.grid,params);
  %  equivalent for fv: a0 = (RB.*repmat(grid.A(:),1,Nmax))' * u0(:);
elseif decomp_mode == 1
  Nmax = size(detailed_data.RB,2);
  u0 = model.init_values_algorithm(model,detailed_data);
  Q_u0 = length(u0);
  a0 = cell(Q_u0,1);
  a0(:) = {zeros(Nmax,1)}; 
  %    a0     : initial data projected on RB set: starting coefficients
  A = detailed_data.W;
  for q=1:Q_u0
    a0{q} = u0{q}' * A *  detailed_data.RB;
  end;
  if model.debug && length(u0) > 1
    if max(abs(detailed_data.RB * a0{1}'-u0{1}))>0.001 
      U0 = [u0{1},u0{2}]; 
      U0_RB = [detailed_data.RB*a0{1}',detailed_data.RB*a0{2}']; 
      plot_params = [];
      plot_params.plot = model.plot;
      plot_params.title = 'U0 Components';
      plot_sequence(U0,detailed_data.grid,plot_params);
      plot_params.title = 'U0_RB';
      plot_sequence(U0_RB,detailed_data.grid,plot_params);
      error('rb initial data badly approximating initial data!!');
    end;
  end;
%  disp('halt in rb_init_values')
%  keyboard;
else % decomp_mode== 2 -> coefficients simply transfered from u0
  u0 = model.init_values_algorithm(model, []);
  a0 = u0;
end;