Mercurial > repos > public > sbplib
diff +time/Magnus4.m @ 510:71908bbce2e1 feature/quantumTriangles
tried to add 4th order magnus integrator which does not work
author | Ylva Rydin <ylva.rydin@telia.com> |
---|---|
date | Fri, 02 Jun 2017 17:46:01 +0200 |
parents | |
children | bc39bb984d88 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+time/Magnus4.m Fri Jun 02 17:46:01 2017 +0200 @@ -0,0 +1,43 @@ +classdef Magnus4 < time.Timestepper + properties + D + S + F + k + t + v + m + n + end + + + methods + function obj = Magnus4(D, k, t0, v0) + obj.D = D; + obj.k = k; + obj.t = t0; + obj.v = v0; + obj.m = length(v0); + obj.n = 0; + end + + function [v,t] = getV(obj) + v = obj.v; + t = obj.t; + end + + function obj = step(obj) + obj.v = time.expint.Magnus_4(obj.v,obj.D, obj.t, obj.k); + obj.t = obj.t + obj.k; + obj.n = obj.n + 1; + end + end + + + methods (Static) + function k = getTimeStep(lambda) + k = rk4.get_rk4_time_step(lambda); + end + end + +end \ No newline at end of file