Mercurial > repos > public > sbplib

comparison +time/+expint/Magnus_4.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 71908bbce2e1
children bc39bb984d88
comparison
equal deleted inserted replaced
511:57f3493f851b 512:4ef2d2a493f1
9 c2 = 1/2 + sqrt(3)/6; 9 c2 = 1/2 + sqrt(3)/6;
10 10
11 A1 = D(t +c1*k); 11 A1 = D(t +c1*k);
12 A2 = D(t + c2*k); 12 A2 = D(t + c2*k);
13 Omega = k/2*(A1 + A2) + sqrt(3)*k^2/12*(A1*A2-A2*A1); 13 Omega = k/2*(A1 + A2) + sqrt(3)*k^2/12*(A1*A2-A2*A1);
14 v = expm(Omega)*v; 14 % v = expm(Omega)*v;
15 toler = 10^(-8);
16 v = time.expint.expm_Arnoldi(-Omega,v,k,toler,100);
15 else 17 else
16 %v = krylov(k*D,v); 18 %v = krylov(k*D,v);
17 v = expm(k*D)*v; 19 v = expm(k*D)*v;
18 end 20 end
19 21