Mercurial > repos > public > sbplib

diff +time/+rk/ButcherTableau.m @ 994:2f89959fb9f0 feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implement method to get gain from butcher tableu
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 09 Jan 2019 12:14:30 +0100
parents 1066bb31bc95
children
line wrap: on
line diff
--- a/+time/+rk/ButcherTableau.m	Wed Jan 09 11:14:16 2019 +0100
+++ b/+time/+rk/ButcherTableau.m	Wed Jan 09 12:14:30 2019 +0100
@@ -27,6 +27,18 @@
             b = all(all(triu(obj.a)==0));
         end
 
+        function g = testEquationGain(obj, z)
+            default_arg('z', sym('z'));
+            s = obj.nStages();
+
+            b = sym(obj.b);
+            A = sym(obj.a);
+            one = sym(ones(s,1));
+            I = sym(eye(s));
+
+            g = abs(1 + z*b*inv(I-z*A)*one);
+        end
+
         % TBD: Add functions for checking accuracy, stability?
     end