Mercurial > repos > public > sbplib

annotate +time/+expint/Magnus_mp.m @ 704:111fcbcff2e9 feature/optim

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
merg with featuew grids
author Ylva Rydin <ylva.rydin@telia.com>
date Fri, 03 Nov 2017 10:53:15 +0100
parents 8434063ed162
children
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.
513
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
4 function v = Magnus_mp(v,D, t , k,matrixexp,tol)
505
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')
513
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
7 Omega = D(t +k/2);
505
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8 else
513
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
9 Omega = D;
505
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10 end
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11
513
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
12 switch matrixexp
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
13 case 'expm'
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
14 v = expm(k*Omega)*v;
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
15 case 'Arnoldi'
697
8434063ed162 1 objective optimization finished
Ylva Rydin <ylva.rydin@telia.com>
parents: 513
diff changeset
16 v = time.expint.expm_Arnoldi(-Omega,v,k,tol,1000);
513
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
17 otherwise
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
18 error('No such matrix exponential evaluation')
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
19
bc39bb984d88 Added arnoldi krylov subspace approximation
Ylva Rydin <ylva.rydin@telia.com>
parents: 512
diff changeset
20 end
505
49bc2929dfb1 Added the midpont magnus exponential integrator
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
21 end