Mercurial > repos > public > sbplib

diff +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
line wrap: on
line diff
--- a/+time/+expint/Magnus_mp.m	Thu Jun 08 10:33:36 2017 +0200
+++ b/+time/+expint/Magnus_mp.m	Mon Jun 26 19:23:19 2017 +0200
@@ -5,10 +5,12 @@
 
 if isa(D,'function_handle')
    % v = krylov(k*D(t +k/2*t),v);
-   v = expm(k*D(t +k/2))*v;
+%   v = expm(k*D(t +k/2))*v;
+      toler = 10^(-5);
+   expm_Arnoldi(-D,v,k,toler,100)
 else
    %v = krylov(k*D,v);
-   v = expm(k*D)*v;
+   % v = expm(k*D)*v;
 end
 
 end
\ No newline at end of file