Mercurial > repos > public > sbplib
view +sbp/Upwind.m @ 99:ecf77a50d4fe feature/arclen-param
The arclength parameter function was slow because it called fzero. Now it is constructed once and for all with splines. Much better performance.
author | Martin Almquist <martin.almquist@it.uu.se> |
---|---|
date | Sun, 06 Dec 2015 15:41:42 +0100 |
parents | 4f5a65f49035 |
children | c56437d097de |
line wrap: on
line source
classdef Upwind < sbp.OpSet properties norms % Struct containing norm matrices such as H,Q, M boundary % Struct contanging vectors for boundry point approximations derivatives % Struct containging differentiation operators borrowing % Struct with borrowing limits for different norm matrices m % Number of grid points. h % Step size end methods function obj = Upwind(m,h,order) if order == 3 [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 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h); elseif order == 8 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h); else error('Invalid operator order %d.',order); end obj.h = h; obj.m = m; obj.norms.H = H; obj.norms.HI = HI; obj.boundary.e_1 = e_1; obj.boundary.e_m = e_m; obj.derivatives.Dp = Dp; obj.derivatives.Dm = Dm; end end methods (Static) function lambda = smallestGrid(obj) error('Not implmented') end end end