Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 1282:a52033540dd9 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make forcingSetup return empty array instead of zero function if all boundary data functions are empty.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 15 Jun 2020 16:11:35 -0700
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