vi_plot_sim_data.m

function p = vi_plot_sim_data(model,model_data,sim_data,plot_params)
%function p = vi_plot_sim_data(model,model_data,sim_data,plot_params)
%
% function plotting the solution of a Parametrized Variational
% Inequality.

% I.Maier 30.05.2011

if nargin < 4
  plot_params = [];
end;
if ~isfield(plot_params,'plot_title')
  plot_params.plot_title = '';
end;
if ~isfield(plot_params,'no_lines')
  plot_params.no_lines = '0';
end;

%keyboard;

if ~isfield(plot_params,'solution_component')
  plot_params.solution_component = 'u';
end;

if plot_params.solution_component == 'u'
  if isequal(model.gridtype, 'onedgrid')
    g = model.obstacle(model_data.grid.X,model);
    %subplot(1,2,1); 
    p1 = plot(model_data.grid.X(2:end-1),sim_data.U,'--k');
    hold on
    p2 = plot(model_data.grid.X,g,'c','Linewidth',4);
    title(plot_params.plot_title);
    hold off
    %subplot(1,2,2);p3 = plot(model_data.X,sim_data.L,'k');
    %p = [p1(:),p2(:),p3(:)];
    p = [p1(:),p2(:)];
  else % triagrid
    
    if ~isfield(plot_params,'subsampling_level');
      plot_params.subsampling_level = 0;
    end;
    if ~isfield(plot_params,'title');
      plot_params.title = '';
    end;
    
    uh = femdiscfunc([],model_data.df_info);
    uh.dofs = sim_data.U;
    plot_params.height_factor = 1;
    plot_params.no_lines = 1;
    p = plot(uh,plot_params);
    hold on;
    
    obstacle = femdiscfunc([],model_data.df_info);
    obstacle = fem_interpol_global(model.obstacle,obstacle,model);
    plot_params.plot_flat_colored = 1;
    plot_params.no_lines = 0;
    plot_params.flat_color = [0.75,0.75,0.75];    
    p = plot(obstacle,plot_params);
    title(plot_params.title);
    
  end;
else % plot lagrange multipliers
  if isequal(model.gridtype, 'onedgrid')
    p = plot(model_data.grid.X(2:end-1),sim_data.L,'--k');
  else % triagrid
    disp('plot_sim_data: please insert')
    keyboard;
  end    
  title(plot_params.plot_title);
end;