t_part_initial_cond_POD.m

function U0=t_part_initial_cond_POD(model, detailed_data)
%function U0=t_part_initial_cond_POD(model, detailed_data)
%
% function calculates the final time step solutions for all parameters used for generating the 
% forgoing partition basis. 
% Then a POD over all these solutions is done and the
% model.basis_vector_overlap first modes are taken a initial basis for the
% new partition.

% Markus Dihlmann 25.02.2011

disp('calculating POD of initial data training set')

U0=[];

t_part_ind = model.t_part_for_simulation - 1;
if t_part_ind == 0
    U0 = RB_init_data_basis(model, detailed_data);
else

  mu_set = detailed_data.t_part_detailed_data{t_part_ind}.RB_info.mu_sequence;

  mus_done = [];

  for i = 1:size(mu_set,2)
    if ~isnan(mu_set(1,i))
      if ~ispartof(mu_set(:,i), mus_done)
        %detailed_simulation of mu
        mus_done = [mus_done, mu_set(:,i)];
        model = set_mu(model, mu_set(:,i));
        sim_data = detailed_simulation(model, detailed_data);

        U_last = sim_data.U(:,end);

        U0 = [U0, U_last];
      end
    end
  end
  U_init = U0(:,1);
  U0 = PCA_fixspace(U0(:,2:end), U_init, detailed_data.W, model.basis_vector_overlap);
  U0 = [U_init, U0];

end

end


function ret = ispartof(mu, M)

  ret =0;
  num_el = size(M,1);
  len=size(M,2);
  l=1;
  while (ret==0)&&(l<=len)

    v=(M(:,l)==mu);
    if sum(v)==num_el
      ret =1  ;    
    end
    l=l+1;
  end
end