3 %
function computing the l2-error between the two fv-functions or function
7 % U1: matrix of column vectors representing finite volume function dof
8 % sequence of first function
9 % U2: matrix of column vectors representing finite volume function dof
10 % sequence of second function
11 % W: mass matrix for L^2 scalar product such that
12 % `u^t W v = \langle u, v \rangle_{L^2}`
15 % l2_error: a single value or sequence of values corresponding to the
16 % column-differences between U1 and U2 correct Omega-integrals are
19 % Bernard Haasdonk 20.7.2006
21 l2_error = sqrt(sum(W * ((U1-U2).^2)));