Skip to content

Support time-dependent operators on dsf_mesolve #610

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
albertomercurio merged 2 commits into from
Dec 18, 2025
Merged

Support time-dependent operators on dsf_mesolve #610

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
3d5a040
Support time-dependent operators on dsf_mesolve
albertomercurio Dec 10, 2025
130aa78
Make changelog
albertomercurio Dec 10, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions CHANGELOG.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,6 +7,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0

## [Unreleased](https://github.com/qutip/QuantumToolbox.jl/tree/main)
- Add error message for bad input in state/operator generating functions ([#603])
- Support time-dependent operators on `dsf_mesolve` ([#610])

## [v0.39.1]
Release date: 2025-11-19
Expand Down Expand Up @@ -383,3 +384,4 @@ Release date: 2024-11-13
[#591]: https://github.com/qutip/QuantumToolbox.jl/issues/591
[#596]: https://github.com/qutip/QuantumToolbox.jl/issues/596
[#603]: https://github.com/qutip/QuantumToolbox.jl/issues/603
[#610]: https://github.com/qutip/QuantumToolbox.jl/issues/610
13 changes: 9 additions & 4 deletions src/time_evolution/time_evolution_dynamical.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -332,8 +332,10 @@ function _DSF_mesolve_Affect!(integrator)
op_l2 = op_list .+ αt_list
e_ops2 = e_ops(op_l2, dsf_params)
_mesolve_callbacks_new_e_ops!(integrator, [_generate_mesolve_e_op(op) for op in e_ops2])
# By doing this, we are assuming that the system is time-independent and f is a MatrixOperator
copyto!(integrator.f.f.A, liouvillian(H(op_l2, dsf_params), c_ops(op_l2, dsf_params)).data)

# By doing this, we are assuming that all the arguments of ODEFunction are the default ones
integrator.f =
ODEFunction{true,FullSpecialize}(_mesolve_make_L_QobjEvo(H(op_l2, dsf_params), c_ops(op_l2, dsf_params)).data)
return u_modified!(integrator, true)
end

Expand All @@ -360,9 +362,12 @@ function dsf_mesolveProblem(

αt_list = convert(Vector{T}, α0_l)
op_l_vec = map(op -> mat2vec(get_data(op)'), op_list)
# Create the Krylov subspace with kron(H₀.data, H₀.data) just for initialize
expv_cache = arnoldi(kron(H₀.data, H₀.data), mat2vec(ket2dm(ψ0).data), krylov_dim)

dsf_identity = Eye(prod(H₀.dimensions))

# Create the Krylov subspace just for initialize
expv_cache = arnoldi(kron(dsf_identity, dsf_identity), mat2vec(ket2dm(ψ0).data), krylov_dim)

dsf_displace_cache_left = sum(op -> ScalarOperator(one(T)) * MatrixOperator(kron(op.data, dsf_identity)), op_list)
dsf_displace_cache_left_dag =
sum(op -> ScalarOperator(one(T)) * MatrixOperator(kron(sparse(op.data'), dsf_identity)), op_list)
Expand Down