Mercurial > repos > public > sbplib

comparison +time/+expint/Magnus_mp.m @ 512:4ef2d2a493f1 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
add H_xyz to hypsyst to compute H-norm
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 26 Jun 2017 19:23:19 +0200
parents 566304a5733e
children bc39bb984d88
comparison
equal deleted inserted replaced
511:57f3493f851b 512:4ef2d2a493f1
3 % time derivatives. 3 % time derivatives.
4 function v = Magnus_mp(v,D, t , k) 4 function v = Magnus_mp(v,D, t , k)
5 5
6 if isa(D,'function_handle') 6 if isa(D,'function_handle')
7 % v = krylov(k*D(t +k/2*t),v); 7 % v = krylov(k*D(t +k/2*t),v);
8 v = expm(k*D(t +k/2))*v; 8 % v = expm(k*D(t +k/2))*v;
9 toler = 10^(-5);
10 expm_Arnoldi(-D,v,k,toler,100)
9 else 11 else
10 %v = krylov(k*D,v); 12 %v = krylov(k*D,v);
11 v = expm(k*D)*v; 13 % v = expm(k*D)*v;
12 end 14 end
13 15
14 end 16 end