hp_simulation.m

 function sim_data = hp_simulation(hp_model)
 %function sim_data = hp_simulation_alt(tree,param)
 %
 % findes the correct subdomain and computes a reduced basis simutlation
 % need a hp_model with a generated reduced basis (function: hp_rb_generation);
 
 
 sim_data = navigation(hp_model.tree,get_mu(hp_model));
 
     function data = navigation(tree,param)
 test = iscell(tree);
        %Navigation
        if test ==1
            % Representanten bestimmen
            first_l=tree_first(tree{1});
            first_r=tree_first(tree{2});
            % Mit d-Funktion naviagtion durchführen
            
            
            if hp_model.error_distance_extionsion ==1
                
                distance_l = hp_model.distance_function(first_l.model,first_l.anchor,param,...
                    first_l.detailed_data,first_l.reduced_data);%norm(first_l-param,2);
                distance_r = hp_model.distance_function(first_r.model, first_r.anchor,param,...
                    first_r.detailed_data,first_r.reduced_data)  ;%norm(first_r-param,2);
                
            else
                
                distance_l = hp_model.distance_function(first_l.model,first_l.anchor,param);%norm(first_l-param,2);
                distance_r = hp_model.distance_function(first_r.model, first_r.anchor,param);%norm(first_r-param,2);
                
            end
            % Unterbaum finden und Rekursiv vortfahren
            if distance_l<distance_r
               data =  navigation(tree{1},param);
            else
               data =  navigation(tree{2},param);
            end;
        else
            % im passenden Gebiet angekommen - RB Simulation durchführen
            tree.model=set_mu(tree.model,param);
            data = rb_simulation(tree.model,tree.reduced_data);
            data = rb_reconstruction(tree.model,tree.detailed_data,data);
        end;
     end      
    %Hilfsfunktion für die Baumnavigation 
    function value= tree_first (tree)
        t=  iscell(tree);
        if t==1
            value=tree_first(tree{1});
        else
            value=tree;
        end;
        
    end

 end