fv_l2_error.m

Go to the documentation of this file.

function l2_error = fv_l2_error(U1,U2,W)
%function l2_error = fv_l2_error(U1,U2,W)
% function computing the l2-error between the two fv-functions or function
% sequences in U1,U2.
%
% Parameters:
%   U1: matrix of column vectors representing finite volume function dof
%       sequence of first function
%   U2: matrix of column vectors representing finite volume function dof
%       sequence of second function
%   W:  mass matrix for L^2 scalar product such that
%       `u^t W v = \langle u, v \rangle_{L^2}`
%
% Return values:
%   l2_error: a single value or sequence of values corresponding to the
%             column-differences between U1 and U2 correct Omega-integrals are
%             computed.

% Bernard Haasdonk 20.7.2006

l2_error = sqrt(sum(W * ((U1-U2).^2)));