Mercurial > repos > public > sbplib

diff +time/Ode45.m @ 1113:47e86b5270ad feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change name of property k to dt in time.Timestepper
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 10 Apr 2019 22:40:55 +0200
parents f5e14e5986b5
children
line wrap: on
line diff
--- a/+time/Ode45.m	Wed Apr 10 22:22:46 2019 +0200
+++ b/+time/Ode45.m	Wed Apr 10 22:40:55 2019 +0200
@@ -1,7 +1,7 @@
 classdef Ode45 < time.Timestepper
     properties
         F
-        k
+        dt
         t
         w
         m
@@ -15,7 +15,7 @@
 
 
     methods
-        function obj = Ode45(D, E, S, k, t0, v0, v0t)
+        function obj = Ode45(D, E, S, dt, t0, v0, v0t)
             obj.D = D;
             obj.E = E;
             obj.S = S;
@@ -31,7 +31,7 @@
                 obj.C = [zeros(obj.m,1), S];
             end
 
-            obj.k = k;
+            obj.dt = dt;
             obj.t = t0;
             obj.w = [v0; v0t];
 
@@ -49,7 +49,7 @@
         end
 
         function obj = step(obj)
-            [t,w] = ode45(@(t,w)(obj.F(w,t)),[obj.t obj.t+obj.k],obj.w);
+            [t,w] = ode45(@(t,w)(obj.F(w,t)),[obj.t obj.t+obj.dt],obj.w);
 
             obj.t = t(end);
             obj.w = w(end,:)';
@@ -59,8 +59,8 @@
 
 
     methods (Static)
-        function k = getTimeStep(lambda)
-            k = rk4.get_rk4_time_step(lambda);
+        function dt = getTimeStep(lambda)
+            dt = rk4.get_rk4_time_step(lambda);
         end
     end