Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 549:ae905a11e32c feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change away from matrix formulation in Rk4 second order form The new formulation can be up to 8 times faster
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Aug 2017 11:06:00 +0200
parents ed6a704b028d
children
line wrap: on
line source

function e = errorSbp(discr, v1, v2)
    % If v1 and v2 are more complex types, something like grid functions... Then we may use .getVectorFrom here!
    H = discr.H;
    err = v2 - v1;
    e = sqrt(err'*H*err);
end