Closed amitjamadagni closed 9 years ago
I think it would be nice to give an example of how we think propagators should be used. For instance,
psi0 = QuArray([1.,0.])
ep = QuPropagator(sigmax, psi0, 0.:0.1:2pi, QuEuler())
for (t, psi) in ep
println("$t $(abs2(psi[1])) $(abs2(psi[2]))")
end
cnp = QuPropagator(sigmax, psi0, 0.:0.1:2pi, QuCrankNicolson())
for (t, psi) in cnp
println("$t $(abs2(psi[1])) $(abs2(psi[2]))")
end
kp = QuPropagator(sigmax, psi0, 0.:0.1:2pi, QuKrylov())
for (t, psi) in kp
println("$t $(abs2(psi[1])) $(abs2(psi[2]))")
end
This is not a very impressive example, but IMO it shows where we are going ...
@acroy a review would be helpful.
Have you already explained how to use QuDynamics in the last blog post? Maybe you can briefly mention that one has to clone QuBase and QuDynamics (basically copy from README.md).
Blog updates for the first two weeks.