p_part_gen_detailed_data.m

function detailed_data = p_part_gen_detailed_data(model,model_data)
%function detailed_data = p_part_gen_detailed_data(model,model_data)
%
% function computing the detailed data for a p_part_model. I.e. a
% list of detailed_data for all partitions plus some management information
% as p_part_mode 'uniform' or 'adaptive' can be chosen

% Bernard Haasdonk 12.2.2010
% Markus Dihlmann 17.2.2010

detailed_data = [];
% field RB is generated

detailed_data.pgrid = [];
%detailed_data.parts_ranges = {};
detailed_data.parts_detailed_data = {};

% initialize pgrid
if isfield(model,'p_part_ranges')
    params.range = model.p_part_ranges;
else
    params.range = model.mu_ranges;
end
params.numintervals = model.p_part_numintervals;
pgrid = cubegrid(params);

switch model.p_part_generation_mode
 case 'uniform'
  nelements = get(pgrid,'nelements');
  for p = 1:nelements
    base_model = model.base_model;
    ranges = get_ranges_of_element(pgrid, p);
    base_model.mu_ranges = ranges;
    detailed_data.parts_detailed_data{p} =  ...
        gen_detailed_data(base_model,model_data); 
    %detailed_data.parts_ranges{p} = ranges;
  end;
  
 case 'adaptive'
     p=1; %p indicates for which part of pgrid the detailed_data is generated.
  disp('Generating detailed_datas in adaptive p-partitions.');
  disp(['Maximum number of Basis per partition is : ',num2str(model.base_model.RB_stop_Nmax)]);
  while(p<=get(pgrid,'nelements'))
      base_model = model.base_model;
      if isfield(model,'p_part_early_refinement')
          base_model.p_part_early_refinement=model.p_part_early_refinement;
      end
      ranges = get_ranges_of_element(pgrid, p);
      base_model.mu_ranges = ranges;
      levels=get(pgrid,'level');
      disp(['Generating detailed_data in partition ',num2str(p),' with refinement level ',num2str(levels(p))]);
      
      %if not yet on last refinement level...
      if(levels(p)<base_model.p_part_max_refinement_level)
          %gen detailed data for partition
          detailed_data.parts_detailed_data{p}=gen_detailed_data(base_model, model_data);
      
            %if generation of detailed data stopped due to maximum Number of
            %Basis-vectors reached, then refine the actual partition
            if(((detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nmax)||detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nlimit_estimation)&& ~(detailed_data.parts_detailed_data{p}.RB_info.stopped_on_epsilon))
                %refine method must get a leaf_element_id, so conversion:
                lid=gid2lid(pgrid,p);
                pgrid=refine(pgrid,lid);
                if(detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nlimit_estimation)
                    disp(['Due to an estimation of needed basis vectors of ',num2str(detailed_data.parts_detailed_data{p}.RB_info.N_limit_est_av),'...']);
                end
                disp(['Refined partition ',num2str(p),' into leafelements of level ',num2str(levels(p)+1)]);
                
                %clean up to save capacity
                is_leaf = get(pgrid,'isleaf');
                if ~is_leaf(p)
                    fields = fieldnames(detailed_data.parts_detailed_data{p});
                    for f=1:length(fields)
                        if~(strcmp(fields{f},'RB_info'))
                            detailed_data.parts_detailed_data{p} = rmfield(detailed_data.parts_detailed_data{p},fields(f));
                        end
                    end
                end
            
            
           
            end
      else
         %last possible refinement level reached
         
         base_model.p_part_early_refinement=0;
         detailed_data.parts_detailed_data{p}=gen_detailed_data(base_model, model_data);
         
      
      end
      p=p+1;
  end   
  
 case 'adaptive_epsilon_fixed'
  
  p=1; %p indicates for which part of pgrid the detailed_data is generated.
  disp('Generating detailed_datas in adaptive p-partitions.');
  disp(['Maximum number of Basis per partition is : ',num2str(model.base_model.RB_stop_Nmax)]);
  while(p<=get(pgrid,'nelements'))
      base_model = model.base_model;
      if isfield(model,'p_part_early_refinement')
          base_model.p_part_early_refinement=model.p_part_early_refinement;
      end
      ranges = get_ranges_of_element(pgrid, p);
      base_model.mu_ranges = ranges;
      levels=get(pgrid,'level');
      disp(['Generating detailed_data in partition ',num2str(p),' with refinement level ',num2str(levels(p))]);
      
      %if not yet on last refinement level...
      if(levels(p)<base_model.p_part_max_refinement_level)
          %gen detailed data for partition
          detailed_data.parts_detailed_data{p}=gen_detailed_data(base_model, model_data);
      
            %if generation of detailed data stopped due to maximum Number of
            %Basis-vectors reached, then refine the actual partition
            if(((detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nmax)||detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nlimit_estimation)&& ~(detailed_data.parts_detailed_data{p}.RB_info.stopped_on_epsilon))
                %refine method must get a leaf_element_id, so conversion:
                lid=gid2lid(pgrid,p);
                pgrid=refine(pgrid,lid);
                if(detailed_data.parts_detailed_data{p}.RB_info.stopped_on_Nlimit_estimation)
                    disp(['Due to an estimation of needed basis vectors of ',num2str(detailed_data.parts_detailed_data{p}.RB_info.N_limit_est_av),'...']);
                end
                disp(['Refined partition ',num2str(p),' into leafelements of level ',num2str(levels(p)+1)]);
                
                %clean up to save capacity
                is_leaf = get(pgrid,'isleaf');
                if ~is_leaf(p)
                    fields = fieldnames(detailed_data.parts_detailed_data{p});
                    for f=1:length(fields)
                        if~(strcmp(fields{f},'RB_info'))
                            detailed_data.parts_detailed_data{p} = rmfield(detailed_data.parts_detailed_data{p},fields(f));
                        end
                    end
                end
            
            
           
            end
      else
         %last possible refinement level reached
         %so we can not prescribe a maximum number of Basisvectors to reach
         %a certain level in error tolerance, so it is set to an
         %arbitrarily high number
         base_model.RB_stop_Nmax = 500;
         base_model.p_part_early_refinement=0;
         detailed_data.parts_detailed_data{p}=gen_detailed_data(base_model, model_data);
         
      
      end
      p=p+1;
  end
      
 otherwise
  error('p_part_generation_mode not defined');
end;

detailed_data.pgrid = pgrid;

%Test of consistency: In all nodepoints of the p-grid the estimated error
%must be below etol..
% to be implemented...
disp('Consistency check...')
vertexes = get(pgrid,'vertex');
reduced_data=gen_reduced_data(model, detailed_data);
consistency_ok = 1;
i=1;
num_vertices = get(pgrid,'nvertices');
try
    while(consistency_ok && (i<=num_vertices))
        model = set_mu(model,vertexes(i,:));
        sim_data = rb_simulation(model,reduced_data);

        if(isfield(sim_data,'DeltaX'))
            sim_error = sim_data.DeltaX(end);
        elseif(isfield(sim_data,'Delta'))
            sim_error = sim_data.Delta(end);
        else
            sim_error = 0;
            error('It was not possible to identify error estimator during p-part consistency check')
        end

        if(model.base_model.RB_stop_epsilon < sim_error)
            consitency_ok = 0;
            disp('!!!!!!!!!!!!!!!!!!!!!!!!!!!');
            error('Consistency Error');

         end
        i=i+1;
    end
catch ME
    disp('Error in consistency check. Define Delta in p_part_gen_detailed_data.m !! Ignore if fixed p-part grid with Nmax...')
end




%detailed_data.number_parts = length(detailed_data.parts_detailed_data); 





%| \docupdate