Mercurial > repos > public > sbplib
diff +time/Rungekutta.m @ 984:0585a2ee7ee7 feature/timesteppers
Inline the rk.rungekutta_4 function.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Tue, 08 Jan 2019 12:19:33 +0100 |
parents | 34b3d092a4d0 |
children | a32856fc2ad2 |
line wrap: on
line diff
--- a/+time/Rungekutta.m Tue Jan 08 12:07:57 2019 +0100 +++ b/+time/Rungekutta.m Tue Jan 08 12:19:33 2019 +0100 @@ -30,13 +30,7 @@ obj.coeffs = struct('s',s,'a',a,'b',b,'c',c); if isempty(discreteData) - % TODO: method "rk4" is also implemented in the butcher tableau, but the rungekutta_4.m implementation - % might be slightly more efficient. Need to do some profiling before deciding whether or not to keep it. - if (method == "rk4") - obj.scheme = @(v,t,n) time.rk.rungekutta_4(v ,t, dt, F); - else - obj.scheme = @(v,t,n) time.rk.rungekutta(v, t, dt, F, obj.coeffs); - end + obj.scheme = @(v,t,n) time.rk.rungekutta(v, t, dt, F, obj.coeffs); else obj.scheme = @(v,t,n) time.rk.rungekuttaDiscreteData(v, t, dt, F, obj.coeffs, discreteData, n); end