Time dependence in QuDynamics

time_dependent_2.jl

julia> require("/home/amit/Downloads/SemVII/QuDynamics.jl/src/QuDynamics.jl")
WARNING: `require` is deprecated, use `using` or `import` instead
 in depwarn at deprecated.jl:73
 [inlined code] from deprecated.jl:694
 in require at no file:0
while loading no file, in expression starting on line 0

julia> using QuBase

julia> using QuDynamics

julia> immutable QuSchrodingerEqTD <: QuEquation{1}
           hamiltonian_const::QuBase.QuArray
           QuSchrodingerEqTD(hamiltonian_const) = new(hamiltonian_const)
       end

julia> function fn(H,t)
           expm(H*t)
       end
fn (generic function with 1 method)

julia> function QuDynamics.operator(QSTD::QuSchrodingerEqTD, t)
           fn(QSTD.hamiltonian_const,t)
       end
operator (generic function with 9 methods)

julia> qsetd = QuSchrodingerEqTD(sigmax + 2*sigmay + 3*sigmaz)
QuSchrodingerEqTD(2x2 QuMatrix in QuBase.FiniteBasis{QuBase.Orthonormal}:
...coefficients: Array{Complex{Float64},2}
Complex{Float64}[3.0 + 0.0im 1.0 - 2.0im
                 1.0 + 2.0im -3.0 + 0.0im])

julia> sys = QuPropagator(qsetd, statevec(1, FiniteBasis(2)), 0.:0.1:2*π,  QuExpmV())
QuDynamics.QuStateEvolution{QuDynamics.QuExpmV,QuBase.QuArray{QuBase.FiniteBasis{QuBase.Orthonormal},Float64,1,Array{Float64,1}},QuSchrodingerEqTD}(QuSchrodingerEqTD(2x2 QuMatrix in QuBase.FiniteBasis{QuBase.Orthonormal}:
...coefficients: Array{Complex{Float64},2}
Complex{Float64}[3.0 + 0.0im 1.0 - 2.0im
                 1.0 + 2.0im -3.0 + 0.0im]),2-element QuVector in QuBase.FiniteBasis{QuBase.Orthonormal}:
...coefficients: Array{Float64,1}
[1.0,0.0],0.0:0.1:6.2,QuDynamics.QuExpmV(Dict{Symbol,Any}()))

julia> for (t,psi) in sys
           println(t,psi)
       end