fv_error.m

function errs = fv_error(U1,U2,W, params)
%function errs = fv_error(U1,U2,W, params);
%
% function computing the error between the two time-dependent
% fv-functions in the columns of U1,U2. Result is a single value or
% sequence of values representing either the l2 differences of the
% temporal sequences or the energy error.
% correct Omega-integrals are computed by respecting the cell-areas
% defined in grid
%
% required fields of model:
%   error_norm: 'l2l2', 'l1l2', 'l2' or 'energy'
%
% plus additional fields required by fv_l2_error, fv_energy_error.
  
% Bernard Haasdonk 30.8.2006
  
if isequal(params.error_norm,'l2')
  errs = fv_l2_error(U1,U2,W);
elseif isequal(params.error_norm, 'energy')
  errs = fv_energy_error(U1,U2,W,params);
elseif isequal(params.error_norm, 'l2l2')
  errs = fv_l2l2_error(U1,U2,W,params);
elseif isequal(params.error_norm, 'l1l2')
  errs = fv_l1l2_error(U1,U2,W,params);
else
  error('error_norm unknown!');
end;