Mercurial > repos > public > sbplib

view +scheme/errorSbp.m @ 678:06676c40e77f feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme for complete elastic operator. Traction BC working with MMS variable coefficient.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 11:06:15 -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