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

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
line wrap: on
line diff
--- a/+scheme/Beam.m	Fri Jul 01 11:38:53 2016 +0200
+++ b/+scheme/Beam.m	Tue Jul 12 16:32:36 2016 +0200
@@ -20,7 +20,8 @@
     methods
         function obj = Beam(grid, order, alpha, opsGen)
             default_arg('alpha', 1);
-            default_arg('opsGen', @sbp.Higher);
+            % default_arg('opsGen', @sbp.Higher);
+            default_arg('opsGen', @sbp.HigherCompatibleVariable); % Supposed to be better
 
             if ~isa(grid, 'grid.Cartesian') || grid.D() ~= 1
                 error('Grid must be 1d cartesian');
@@ -107,8 +108,10 @@
             gamm_v = neighbour_scheme.gamm;
             delt_v = neighbour_scheme.delt;
 
-            % tuning = 2;
-            tuning = 1.1;
+            tuning = 2;
+            % tuning = 1.1;
+            % tuning = 0.5;
+            % tuning = 0.49998;
 
             alpha_u = obj.alpha;
             alpha_v = neighbour_scheme.alpha;