Mercurial > repos > public > sbplib

view +time/+expint/Magnus_mp.m @ 508:508b7493be94 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
moved all multiplications and divisions bu J and Ji to the siscretization to easier find the error
author Ylva Rydin <ylva.rydin@telia.com>
date Fri, 02 Jun 2017 11:33:04 +0200
parents 49bc2929dfb1
children 566304a5733e
line wrap: on
line source

% Takes one time step of size k using the magnus midpoinr
% starting from v_0 and where the function F(v,t) gives the
% time derivatives.
function v = Magnus_mp(v,D, t , k)

if isa(D,'function_handle')
   % v = krylov(k*D(t +k/2*t),v);
   v = expm(k*D(t +k/2*t))*v;
else
   %v = krylov(k*D,v);
   v = expm(k*D)*v;
end

end