Mercurial > repos > public > sbplib
diff +sbp/Upwind.m @ 52:4f5a65f49035
Fixed the upwind operators.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Wed, 11 Nov 2015 15:43:15 -0800 |
parents | ebd25e5a481a |
children | c56437d097de |
line wrap: on
line diff
--- a/+sbp/Upwind.m Wed Nov 11 15:15:26 2015 -0800 +++ b/+sbp/Upwind.m Wed Nov 11 15:43:15 2015 -0800 @@ -1,4 +1,4 @@ -classdef Ordinary < sbp.OpSet +classdef Upwind < sbp.OpSet properties norms % Struct containing norm matrices such as H,Q, M boundary % Struct contanging vectors for boundry point approximations @@ -9,10 +9,10 @@ end methods - function obj = Ordinary(m,h,order) + function obj = Upwind(m,h,order) if order == 3 - [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3(m,h); + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3_3bp(m,h); elseif order == 4 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h); elseif order == 6 @@ -28,7 +28,6 @@ obj.norms.H = H; obj.norms.HI = HI; - obj.norms.Q = Q; obj.boundary.e_1 = e_1; obj.boundary.e_m = e_m;