comsol_RB_extension_greedy_lin_stat.m

function [RBext,par_dummy] = comsol_RB_extension_greedy_lin_stat(model,detailed_data)
%function [RBext,par_dummy] = RB_extension_simple_greedy_lin_stat_comsol(model,detailed_data)
%
% this function is needed for the greedy-algorithm for stationary functions.
% It returns an orthonormalized basis function
% 
% Oliver Zeeb, 18.09.2012

par_dummy=[];

%allgemeiner:
sim_data=detailed_simulation(model,detailed_data);
RBext = model.get_dofs_from_sim_data(sim_data);

% orthonormalize new basis vector
InProd_matrix = model.get_inner_product_matrix(detailed_data);
Uproj = RBext - ...
       detailed_data.RB * (detailed_data.RB' * InProd_matrix * ...
                           RBext); %projection
norm_Uproj = sqrt(max(Uproj'*InProd_matrix*Uproj,0)); %norm
if norm_Uproj > 1e-10 
      % if nonzero projection error -> end with new snapshot
      Uproj = Uproj/norm_Uproj;
else
      % due to scaling, no precise orthonormalization is obtained. So perform
      % explicit orthonormalization
      RBtmp = orthonormalize([detailed_data.RB,Uproj],InProd_matrix);
      if size(RBtmp,2)> size(detailed_data.RB,2)
        Uproj = RBtmp(:,end); 
      else
        error('accuracy problem, should not reach this point!!');
      end;
end;

%return orthonormalized basis vector
RBext=Uproj;