descr_default.m

Go to the documentation of this file.

function descr = descr_default
  % function model = descr_default
  % default discretization description for problems of type "NonlinEvol"
  %
  % generated fields of descr:
  %   init_values_algorithm            : fv_init_values()
  %   filecache_ignore_fields_in_model : filecache
  %   inner_product_matrix_algorithm   : fv_inner_product_matrix()
  %   inner_product                    : fv_inner_product()
  %   get_inner_product_matrix         : extracts inner product matrix from model_data.
  %   geometry_transformation          : 'none'
  %

% init values projection algorithm
descr.init_values        = SeparableFunctionDefault(@fv_init_values);

descr.filecache_ignore_fields_in_model = {'N','Nmax','M',...
  'filecache_ignore_fields_in_model',...
  'filecache_ignore_fields_in_detailed_data',...
  'M_by_N_ratio'};

descr.inner_product_matrix_algorithm = @fv_inner_product_matrix;
descr.inner_product                  = @fv_inner_product;
descr.get_inner_product_matrix       = @(model_data) model_data.W;

descr.geometry_transformation = 'none';

descr.decomp_mode    = 0;
descr.t              = 0;
descr.tstep          = 0;
descr.newton_steps   = 100;
descr.newton_epsilon = 1e-11;

descr.mean_ptr        = @(X) harmmean(max(X, 0), 2);
descr.mean_deriv1_ptr = @(X) 2./( (1./X(:,1) + 1./(X(:,2))).^2 .* X(:,1).^2 );
descr.mean_deriv2_ptr = @(X) 2./( (1./X(:,1) + 1./(X(:,2))).^2 .* X(:,2).^2 );

descr.L_saturation = Fv.TwoPhase.SaturationSpace;
descr.L_velocity   = Fv.TwoPhase.VelocitySpace;
descr.L_divergence = Fv.TwoPhase.DivergenceSpace;
descr.L_pressure   = Fv.TwoPhase.PressureMean;