Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 695:ba0d31ce4121 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed jacobian at right side of the penalty term in 2D, seems like it is working
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 19 Sep 2017 09:51:08 +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