demo_adaptive_basisgen.m

Go to the documentation of this file.

% demo_adaptive_basisgen 
%
% small animation of the adaptive basis generation for the
% gdl-convection-diffusion problem. I.e. the parameter domain is
% plotted with the adaptive grid that is used in the basis
% generation procedure. The selected parameter is indicated with a
% point for each new basis vector and the colors indicate the
% number of basis functions selected for the parameter (as we have
% a time-dependent problem, many snapshots for identical parameter
% can possibly be end up in the basis: essentially all
% time-instants of the trajectory.)
%
% The nice aspect here is, that the physics of the problem is
% recovered: The beta-value (x-axis) represents a weight in the
% convex combination defining the boundary data. Hence the most
% important snapshots are correctly found for beta=0 and 1 (and
% almost reflection symmetric!)
% The y-axis represents the diffusivity: We also recover the
% optimal logarithmical distribution of snapshot points, which is
% known from theory. But here it is found completely automatically.

% Bernard Haasdonk 12.12.2007

help demo_adaptive_basisgen

fn = 'RB_adaptive_refined_2_2_r1_smax1000_theta5e-2.mat';
load(fn);

% change/add some settings, if datafile has been generated by earlier versions
%extend = 1;
extend = 0;
if extend
  params.flux_linear = 1;
  for i = 1:length(detailed_data.RB_info.MMesh_list)
    detailed_data.RB_info.MMesh_list{i} = ...
        extend_class_from_struct('cubegrid',detailed_data.RB_info.MMesh_list{i});
  end;
  save(fullfile(rbmatlabhome,'demos','datafiles',fn),'detailed_data','params');
end;
  
animate_basisgen(detailed_data,params)
%| \docupdate