eop_fd_norm.m

Go to the documentation of this file.

function norm = eop_fd_norm(fd_function1, fd_function2, grid, ~)
%norm = eop_fd_norm(fd_function1, fd_function2, grid, ~)
%
% function that computes the discrete L2-norm of the error ||fd_function1 -
% fd_function2|| in every timestep, i.e. a sequence. Choose one of the
% functions to be zero in order to recieve the norm of a single function.

% fd_function1 and fd_function2 have to have the same size and they are
% assumed to be in matrix-form, i.e. one column represents the function in
% one timeslice.

% grid is for example model_data.grid

% Dominik Garmatter 10.04 2012

A = speye(grid.nvertices)*1/grid.nvertices; % mass matrix of the discrete L2-SP
C = fd_function1-fd_function2;
norm = zeros(1,size(C,2));

for i = 1:size(C,2)
    skalarproduct = C(:,i)'*A*C(:,i);
    norm(1,i) = sqrt(skalarproduct);
end