Mercurial > repos > public > sbplib

changeset 239:ba56e0d621f2 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Beam: Added ability to pick alphaII and alphaIII to control tau and sigma.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 28 Jul 2016 12:50:43 +0200
parents 5f75d7094983
children 46256fffa329
files +scheme/Beam.m
diffstat 1 files changed, 10 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Beam.m	Thu Jul 28 11:45:02 2016 +0200
+++ b/+scheme/Beam.m	Thu Jul 28 12:50:43 2016 +0200
@@ -19,9 +19,11 @@
     end
 
     methods
-        function obj = Beam(grid, order, alpha, opsGen, interface_tuning)
+        function obj = Beam(grid, order, alpha, opsGen, interface_tuning, alphaII, alphaIII)
             default_arg('alpha', -1);
             default_arg('interface_tuning', 1.1);
+            default_arg('alphaII',  [])
+            default_arg('alphaIII', [])
 
             % default_arg('opsGen', @sbp.Higher);
             default_arg('opsGen', @sbp.HigherCompatibleVariable); % Supposed to be better
@@ -55,8 +57,13 @@
 
             obj.D = alpha*D4;
 
-            obj.gamm = h*ops.borrowing.N.S2/2;
-            obj.delt = h^3*ops.borrowing.N.S3/2;
+            if isempty(alphaII) && isempty(alphaIII)
+                alphaII = ops.borrowing.N.S2/2;
+                alphaIII = ops.borrowing.N.S3/2;
+            end
+
+            obj.gamm = h*alphaII;
+            obj.delt = h^3*alphaIII;
             obj.interface_tuning = interface_tuning;
         end