eop_master_gui.m

function eop_master_gui(varargin)
%function eop_master_GUI(model,detailed_data,[sim_data.U],[reduced_data], [callbackfigure], [cbhandle]) 
%
% Builds a GUI for the european_option_pricing_model, which has 4 sliders
% for the 4 differnet parameters, another slider for the time and a last
% slider for N, the number of used reduced basis vectors. In the top right
% the reduced simulation to the current settings is plottet and in the top
% left the error ||u-u_N|| in a discrete L2-Norm is plottet, dependant on
% the number N of used reduced basis vectors. With the button "Knopf" in
% the center a detailed simulation is done and plottet to the current
% setting and theerror ||u-u_N|| for the current settings is updated.


%% construct the GUI
if nargin < 6 % no callback -> new figure
% create new figure g and extract the input and write the important parts in ud.
% ud will be passed on for callbacks etc.
  g = figure;
  ud = [];
  list = varargin;
  if nargin < 3
    list{3} = [];
  end
  
  if nargin < 4
    list{4} = [];
  end;

  if nargin >=2
    ud.model = list{1};
    ud.detailed_data = list{2};
  else
    error('see description for the right input arguments!');
  end
% write size of RB in ud.model
  if ~isfield(ud.model,'N')
    ud.model.N = ceil(ud.model.get_rb_size(ud.model,ud.detailed_data)/2);
  end
% write reduced data in ud
  if isempty(list{4})
    disp('Preparing offline-quantities...')
    ud.reduced_data = gen_reduced_data(ud.model, ud.detailed_data);
  else 
    ud.reduced_data = list{4};
  end
% if not in the input do a detailed_simulation for the current mu and save
% it. If in the input the detailed_simulation should have been done for the
% current mu.
  if isempty(list{3})
    detailed_simulation_data = detailed_simulation(ud.model, ud.detailed_data);
    ud.detailed_simulation_data = detailed_simulation_data.U;
  else 
    ud.detailed_simulation_data = list{3};
  end
  ud.first_mu = ud.model.get_mu(ud.model); % save this first mu, since the detailed_simulation was done with it (later used for the error)
  
  % setting GUI size parameters
  di = 0.05;
  textwidth= 0.15;
  textheight = 0.05;
  yscale_uicontrols = 1;
   % define some stuff so that the slider positioning works
   for i=1:length(ud.model.mu_names)
    Positioner1 = ceil(i/2);
    if i/2 ~= floor(i/2)
      Positioner2 = 0;
    elseif i/2 == floor(i/2)
      Positioner2 = 0.5;
    end
    % various sliders for the parameters - use current mus as first slider value/position
    mu = ud.model.get_mu(ud.model);
    val = mu(i);
    ra = ud.model.mu_ranges{i};
    norm_val = (val-ra(1))/(ra(2)-ra(1));
    txt(i) = uicontrol(g,'Style','text','Units','normalized',...
                       'Position',...
                       [di+Positioner2,(Positioner1*di+(Positioner1-1)*textheight)*yscale_uicontrols,...
                        textwidth,textheight*yscale_uicontrols],...
                       'String',[ud.model.mu_names{i},'=',num2str(val)],...
                       'Tag',['text_',ud.model.mu_names{i}]);
    slider(i) = uicontrol(g,'Style','slider','Units','normalized',...
                          'Position',...
                          [ 2*di+textwidth+Positioner2, ...
                          (Positioner1*di+(Positioner1-1)*textheight)*yscale_uicontrols, ...
                          (1-3*di-textwidth)/3,...
                          textheight*yscale_uicontrols],...
                          'Tag',['slider_',ud.model.mu_names{i}],...
                          'Min',0, ...
                          'Max',1, ...
                          'value',norm_val, ...                   
                          'Callback','eop_master_gui([],[],[],[],[],[],gcbf,gcbo)');
   end

  % N = number of RB slider
  i = ceil(length(ud.model.mu_names)/2)+1;
  txt(i) = uicontrol(g,'Style','text','Units','normalized',...
                     'Position', ...
                     [di,(i*di+(i-1)*textheight)*yscale_uicontrols,...
                      textwidth,textheight*yscale_uicontrols],...
                     'String',['RB-size N = ',num2str(ud.model.N)],...
                     'Tag','text_N');
  slider(i) = uicontrol(g,'Style','slider','Units','normalized',...
                        'Position',...
                        [2*di+textwidth, ...
                        (i*di + (i-1)*textheight)*yscale_uicontrols, ...
                        (1-3*di-textwidth)/3,textheight*yscale_uicontrols],...
                        'Tag','slider_N',...
                        'min',1,'max', ...
                        ud.model.get_rb_size(ud.model,ud.detailed_data),...
                        'value',ud.model.N, ...
                        'Callback','eop_master_gui([],[],[],[],[],[],gcbf,gcbo)',...
                        'sliderstep',...
                        [1/((ud.model.get_rb_size(ud.model,ud.detailed_data))-1),0.1]);
    %Time-slider
    txt(i) = uicontrol(g,'Style','text','Units','normalized',...
                       'Position', ...
                       [di+0.5,(i*di+(i-1)*textheight)*yscale_uicontrols,...
                        textwidth,textheight*yscale_uicontrols],...
                       'String','timestep 1','Tag','text_time');
    slider(i) = uicontrol(g,'Style','slider','Units','normalized',...
                          'Position',...
                          [2*di+textwidth+0.5, ...
                    (i*di + (i-1)*textheight)*yscale_uicontrols, ...
                    (1-3*di-textwidth)/3,textheight*yscale_uicontrols],...
                          'Tag','slider_time',...
                          'min',1,'max',ud.model.nt+1,'value',1, ...
                          'Callback','eop_master_gui([],[],[],[],[],[],gcbf,gcbo)',...
                          'sliderstep',[1/ud.model.nt,0.1]);

    % graphics axes
    axes('Units','normalized',...
            'Position',[di,((i+2)*di+i*textheight)*yscale_uicontrols,...
                    di+textwidth+(1-3*di-textwidth)/3,1-((i+3)*di + i*textheight)*yscale_uicontrols],...
            'Tag','axes_1');
    
    axes('Units','normalized',...
            'Position',[di+0.5,((i+2)*di+i*textheight)*yscale_uicontrols,...
                    di+textwidth+(1-3*di-textwidth)/3,1-((i+3)*di + i*textheight)*yscale_uicontrols],...
            'Tag','axes_2');
    Master_Button = uicontrol(g,'Style','pushbutton','Units','normalized','String','Knopf',...
                        'Position',[0.5-di/4,((i+2)*di+i*textheight)*yscale_uicontrols-6/4*di,di,di],...
                        'Tag','Master_button','Callback','eop_master_gui([],[],[],[],[],[],gcbf,gcbo)');

  set(g,'Color',[0.8147 0.9058 0.1270]); % fancy backgroundcolour
  set(g,'Userdata',ud); % take the structure ud into the Userdata in g, so it is enough to have g.
  replot(g,'new');  
end

%% callback

if nargin >=6 % callback
% get the handles of the current figure and the current callback object 
  cbf = varargin{7};
  cbo = varargin{8};
  ud = get(cbf,'Userdata'); % extract data
  tag = get(cbo,'Tag'); % tag of the callback object
% updates the correct callbackobject and gives back a sourcestring for the
% replot function to work properly
  if isequal(tag,'slider_time')
    val = round(get(cbo,'value'));
    set(cbo,'value',val);
    ud.timestep = val;
    txt = findobj(cbf,'Tag','text_time');
    set(txt,'String',['timestep = ',num2str(val)]);
    sourcestr = 'time';
  elseif isequal(tag,'slider_N')
    val = round(get(cbo,'value'));
    set(cbo,'value',val);
    ud.model.N = val;
    txt = findobj(cbf,'Tag','text_N');
    set(txt,'String',['RB-size N = ',num2str(val)]);
    sourcestr = 'N';
  elseif isequal(tag,'Master_button')
    sourcestr = 'Master';
  else
    found_cbo = 0;
    for i=1:length(ud.model.mu_names)
      if isequal(tag,['slider_',ud.model.mu_names{i}])
        found_cbo = 1;
        % rescaling of the mu values since the slider is normed to [0,1]
        val = get(cbo,'value')*...
                  (ud.model.mu_ranges{i}(2)-...
                   ud.model.mu_ranges{i}(1)) + ... 
                   ud.model.mu_ranges{i}(1);
        mu = ud.model.get_mu(ud.model);
        mu(i) = val;
        ud.model = ud.model.set_mu(ud.model,mu);
        txt = findobj(cbf,'Tag',['text_',ud.model.mu_names{i}]);
        set(txt,'String',[ud.model.mu_names{i},' = ',num2str(val)]);
        sourcestr = 'mu';
      end
    end
    if ~found_cbo
      error('callback object unknown!!');
    end
  end
  set(cbf,'Userdata',ud); % update userdata
  replot(cbf, sourcestr);
end


function replot(g,tag)
ud = get(g,'Userdata'); % extract userdata

if ismember(tag,{'new','mu','N','time'})
    if strcmp(tag,'new') == 1
        ud.timestep = 1;
    end
    reduced_data = extract_reduced_data_subset(ud.model,ud.reduced_data); % reduced data of the desired size
    simulation_data = rb_simulation(ud.model,reduced_data);
    rb_sim_data=rb_reconstruction(ud.model,ud.detailed_data,simulation_data);
    ud.sim_data = rb_sim_data; % save the new rb_simulation
    % go to the correct axis
    ax2 = findobj(g,'Tag','axes_2');
    set(gcf,'CurrentAxes',ax2);
    % to plot the new data of the correct time step
    plot_vertex_data(ud.detailed_data.grid,ud.sim_data.U(:,ud.timestep),[]);
    set(gca,'Title',text('String','the reduced simulation')) % has to be set again :(
    if strcmp(tag,'N') == 1 || strcmp(tag,'new') == 1 % plot the error above N the actual size of the RB
        disp(['performing RB simulation with N = ',num2str(ud.model.N)]);
        if strcmp(tag,'new') == 1
            ud.model = ud.model.set_mu(ud.model,ud.first_mu);
            ud.current_mu = ud.first_mu;
        else
            ud.model = ud.model.set_mu(ud.model,ud.current_mu);
        end
        reduced_data = extract_reduced_data_subset(ud.model,ud.reduced_data);
        simulation_data = rb_simulation(ud.model,reduced_data);
        rb_sim_data=rb_reconstruction(ud.model,ud.detailed_data,simulation_data);
        ud.sim_data = rb_sim_data;
        % new error
        Fehler = eop_fd_norm(ud.detailed_simulation_data, ud.sim_data.U, ud.detailed_data.grid);
        % correct axis
        ax = findobj(g,'Tag','axes_1');
        set(gcf,'CurrentAxes',ax);
        fprintf('plotting...\n');
        % plot error
        plot(0:ud.model.deltaT:ud.model.T,Fehler','r o :');
        set(gca,'Tag','axes_1');
        set(gca,'Title',text('String','the error ||u-u_N|| above N'))
        set(get(gca,'XLabel'),'String','time t');
        set(get(gca,'YLabel'),'String','error');
        set(gca,'XTick',0:2*ud.model.deltaT:ud.model.T);
    end
    
elseif isequal(tag,'Master')
    % master button performs a detailed_simulation to the current mu and
    % updates the error
    detailed_simulation_data = detailed_simulation(ud.model, ud.detailed_data);
    ud.detailed_simulation_data = detailed_simulation_data.U;
    if ud.model.compute_output_functional == 1
        params.plot_output_title = ['the output functional for mu = ' mat2str(get_mu(ud.model))];
        params.plot_title = ['the detailed simulation for the current parameter mu = ', mat2str(round((get_mu(ud.model).*1000))/1000)];
    else
        params.plot_title = ['the detailed simulation for the current parameter mu = ', mat2str(round((get_mu(ud.model).*1000))/1000)];
    end
    % plots it in a seperate window
    lin_evol_plot_sim_data(ud.model,ud.detailed_data,detailed_simulation_data,params);
    ud.current_mu = get_mu(ud.model);
    reduced_data = extract_reduced_data_subset(ud.model,ud.reduced_data);
    simulation_data = rb_simulation(ud.model,reduced_data);
    rb_sim_data=rb_reconstruction(ud.model,ud.detailed_data,simulation_data);
    ud.sim_data = rb_sim_data;
    % error update in the GUI
    Fehler = eop_fd_norm(ud.detailed_simulation_data, ud.sim_data.U, ud.detailed_data.grid);
    ax = findobj(g,'Tag','axes_1');
    set(g,'CurrentAxes',ax);
    figure(g);
    fprintf('plotting...\n');
    plot(0:ud.model.deltaT:ud.model.T,Fehler','r o :');
    set(gca,'Tag','axes_1');
    set(gca,'Title',text('String','the error ||u-u_N|| above N'))
    set(get(gca,'XLabel'),'String','time t');
    set(get(gca,'YLabel'),'String','error');
    set(gca,'XTick',0:2*ud.model.deltaT:ud.model.T);
end
set(g,'Userdata',ud); % update user data