Stepsize.m

From Spinach Documentation Wiki
Jump to: navigation, search

This function estimates the norm of the Liouvillian supplied and returns the time step that ensures monotonic convergence of the powder series in the exponential propagator. From the numerical accuracy point of view, this is the optimal time step for the evolution. Syntax:

    [timestep,nsteps]=stepsize(L,interval)

Parameters:

    L   -      Liouvillian operator or superoperator inder which the 
               evolution is taking place.
    interval - the total evolution time, must be specified if the op-
               timal number of steps is also wanted in the output.

Outputs:

    timestep - optimal time step
    nsteps   - number of steps covering the interval, if the interval
               was specified

Note: if a zero Liouvillian is supplied, the function would return the following values: timestep=0, nsteps=1