Mercurial > repos > public > sbplib

comparison +scheme/Beam.m @ 232:443b14ea705b feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
scheme.Beam: changes to possibly better operators.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 12 Jul 2016 16:32:36 +0200
parents 05304d7792dc
children d6344b549081
comparison
equal deleted inserted replaced
231:7f6f04bfc007 232:443b14ea705b
18 end 18 end
19 19
20 methods 20 methods
21 function obj = Beam(grid, order, alpha, opsGen) 21 function obj = Beam(grid, order, alpha, opsGen)
22 default_arg('alpha', 1); 22 default_arg('alpha', 1);
23 default_arg('opsGen', @sbp.Higher); 23 % default_arg('opsGen', @sbp.Higher);
24 default_arg('opsGen', @sbp.HigherCompatibleVariable); % Supposed to be better
24 25
25 if ~isa(grid, 'grid.Cartesian') || grid.D() ~= 1 26 if ~isa(grid, 'grid.Cartesian') || grid.D() ~= 1
26 error('Grid must be 1d cartesian'); 27 error('Grid must be 1d cartesian');
27 end 28 end
28 29
105 delt_u = obj.delt; 106 delt_u = obj.delt;
106 107
107 gamm_v = neighbour_scheme.gamm; 108 gamm_v = neighbour_scheme.gamm;
108 delt_v = neighbour_scheme.delt; 109 delt_v = neighbour_scheme.delt;
109 110
110 % tuning = 2; 111 tuning = 2;
111 tuning = 1.1; 112 % tuning = 1.1;
113 % tuning = 0.5;
114 % tuning = 0.49998;
112 115
113 alpha_u = obj.alpha; 116 alpha_u = obj.alpha;
114 alpha_v = neighbour_scheme.alpha; 117 alpha_v = neighbour_scheme.alpha;
115 118
116 tau1 = ((alpha_u/2)/delt_u + (alpha_v/2)/delt_v)/2*tuning; 119 tau1 = ((alpha_u/2)/delt_u + (alpha_v/2)/delt_v)/2*tuning;