Mercurial > repos > public > sbplib
view +scheme/errorSbp.m @ 852:fbb8be3177c8 feature/burgers1d
Fix bug in SAT terms for upwind operators
- Set the appropriate boundary derivatives in respect to choice of upwind discretization of D2
- Distinguish between upwind discretizations e.g upwind+ vs upwind-
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Thu, 27 Sep 2018 09:30:21 +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