Mercurial > repos > public > sbplib
view +scheme/errorSbp.m @ 1119:8984b12feba6 feature/poroelastic
Add diffOp (currently scheme) for Divergence, not sure if this belongs in +scheme.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Mon, 06 May 2019 21:23:21 -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