Mercurial > repos > public > sbplib

changeset 236:8592358a3935 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Beam: Added tuning as a parameter to contructor.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 15 Jul 2016 11:22:12 +0200
parents d6344b549081
children 0fc1de5cd85b
files +scheme/Beam.m
diffstat 1 files changed, 8 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Beam.m	Thu Jul 14 10:41:32 2016 +0200
+++ b/+scheme/Beam.m	Fri Jul 15 11:22:12 2016 +0200
@@ -15,11 +15,14 @@
         d3_l, d3_r
         gamm
         delt
+        interface_tuning
     end
 
     methods
-        function obj = Beam(grid, order, alpha, opsGen)
-            default_arg('alpha', 1);
+        function obj = Beam(grid, order, alpha, opsGen, interface_tuning)
+            default_arg('alpha', -1);
+            default_arg('interface_tuning', 1.1);
+
             % default_arg('opsGen', @sbp.Higher);
             default_arg('opsGen', @sbp.HigherCompatibleVariable); % Supposed to be better
 
@@ -54,6 +57,7 @@
 
             obj.gamm = h*ops.borrowing.N.S2/2;
             obj.delt = h^3*ops.borrowing.N.S3/2;
+            obj.interface_tuning = interface_tuning;
         end
 
 
@@ -108,11 +112,8 @@
             gamm_v = neighbour_scheme.gamm;
             delt_v = neighbour_scheme.delt;
 
-            % tuning = 2;
-            tuning = 1.1;
-            % tuning = 0.5;
-            % tuning = 0.49998;
-            % tuning = 0.3;
+            tuning = obj.interface_tuning;
+
 
             alpha_u = obj.alpha;
             alpha_v = neighbour_scheme.alpha;