Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 407:e1db62d14835 feature/SBPInTimeGauss

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug with initial data and made gauss the default type
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 02 Feb 2017 19:30:16 +0100
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