Mercurial > repos > public > sbplib

annotate +time/+expint/Magnus_mp.m @ 509:566304a5733e feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
changed the magnus midpoint rule
author Ylva Rydin <ylva.rydin@telia.com>
date Fri, 02 Jun 2017 17:43:59 +0200
parents 49bc2929dfb1
children 4ef2d2a493f1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
505
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
1 % Takes one time step of size k using the magnus midpoinr
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2 % starting from v_0 and where the function F(v,t) gives the
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
3 % time derivatives.
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4 function v = Magnus_mp(v,D, t , k)
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
5
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
6 if isa(D,'function_handle')
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 % v = krylov(k*D(t +k/2*t),v);
509
566304a5733e changed the magnus midpoint rule
Ylva Rydin <ylva.rydin@telia.com>
parents: 505
diff changeset
8 v = expm(k*D(t +k/2))*v;
505
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 else
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10 %v = krylov(k*D,v);
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11 v = expm(k*D)*v;
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
12 end
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
13
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
14 end