Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 680:cd1a76c38565 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove schemes elasticDilation and elasticShear because they are subsets of elastic.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 14:58:48 -0800
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