conv_flux_affine_decomposed.m

Go to the documentation of this file.

function flux = conv_flux_affine_decomposed(glob, U, params)
%function flux = conv_flux_affine_decomposed(glob, U, params)
% function evaluating a function in the list of global coordinates
% specified in the columns of 'glob'.
%
% The result is a matrix of flux vector results as columns of flux.
%
% The function is given by function pointer to components and coefficients
% required fields of params:
%    conv_flux_coefficients_ptr: Function must
%                  allow call 'flux = f(params)' Return is a vector of length
%                  'Q' with (parameter dependent) coefficients in flux.
%    conv_flux_components_ptr: Functions must
%                  allow call 'flux = f(glob,U,params)' Return of components is
%                  a cell array of matrices of the same size as 'glob' with the
%                  point values of the flux.
%
% Linear combination of components by coefficients then yields the
% complete evaluation.

% Bernard Haasdonk 27.8.2009

% determine affine_decomposition_mode as integer  
% glob column check
if params.debug
  if ~isempty(glob) && size(glob,1) < size(glob,2)
    warning('coordinates in variable glob are given row-wise, but expected them to be column-wise');
    if params.debug > 2
      keyboard;
    end
  end
end
decomp_mode = params.decomp_mode;

if decomp_mode == 2
  flux = params.conv_flux_coefficients_ptr(params);
elseif decomp_mode == 1; % components
  flux = params.conv_flux_components_ptr(glob,params);
else % decomp_mode = 0, complete
  flux_coefficients = params.conv_flux_coefficients_ptr(params);
  flux_components = params.conv_flux_components_ptr(glob,params); 
  flux = lincomb_sequence(flux_components,flux_coefficients);     
end;