lin_evol_get_output.m

Go to the documentation of this file.

function [y,data] = lin_evol_get_output(model, varargin)
%lin_evol_get_output_detailed(model, varargin)
%
%Function returns the output which should be optimized. 
%It performs a detailed/reduced simulation (considering varagin = reduced_data or detailed_data) for the mu set in model and returns the
%output (the average concentration in the lower right quadrant).
% 
% in the struct data additional data can be returned
%
% Markus Dihlmann 04.02.2010
data = [];

if isfield(varargin{1},'grid')
    model_data = varargin{1};
    if length(varargin)>1
        mu= varargin{2};
%         if (mu(1)<0)||(mu(1)>1)||(mu(2)<0)||(mu(2)>1)
%             keyboard;
%         end
        model = set_mu_to_optimize(model, mu);
    end
     %perform detailed simulation with actual parameter set
        model.compute_derivative_info = 0;
        sim_data = detailed_simulation(model,model_data);

    y = sim_data.y(end);
    %y=y(2);
    if(model.verbose>=8)
        disp('detailed simulation output result:')
        y
    end
    
else 
   reduced_data = varargin{1};
    %perform reduced simulation
    if length(varargin)>1
        mu= varargin{2};
%         if (mu(1)<0)||(mu(1)>1)||(mu(2)<0)||(mu(2)>1)
%             keyboard;
%         end
        model = set_mu_to_optimize(model, mu);
    end
    model.compute_derivative_info = 0;
    sim_data =rb_simulation(model, reduced_data);
    %if isempty(model_data)
    %    model_data = gen_model_data(model);
    %end
    %disp('rewrite lin_evol_get_output for faster evaluation');
  %  sim_data = model.rb_reconstruction(model, model_data, detailed_data, sim_data);
    y=sim_data.s(end);
    %y=sim_data.s(end);
    if(model.verbose>=8)
        disp('reduced simulation output result:')
        y
    end
    
    data.nRB = sim_data.nRB;
    
end