view +scheme/errorSbp.m @ 256:87504db14c6d operator_remake

Bug fix in D2.m, added 12th order.
author Martin Almquist <martin.almquist@it.uu.se>
date Wed, 07 Sep 2016 16:41:31 +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