Make changelog and format

Author: albertomercurio

CHANGELOG.md

@@ -7,6 +7,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ## [Unreleased](https://github.com/qutip/QuantumToolbox.jl/tree/main)
 
+- Add support for `OperatorKet` state input for `mesolve` and `smesolve`. ([#423])
+
 ## [v0.28.0]
 Release date: 2025-02-22
 
@@ -171,3 +173,4 @@ Release date: 2024-11-13
 [#419]: https://github.com/qutip/QuantumToolbox.jl/issues/419
 [#420]: https://github.com/qutip/QuantumToolbox.jl/issues/420
 [#421]: https://github.com/qutip/QuantumToolbox.jl/issues/421
+[#423]: https://github.com/qutip/QuantumToolbox.jl/issues/423

src/time_evolution/mesolve.jl

@@ -79,7 +79,7 @@ function mesolveProblem(
     ρ0 = if isoperket(ψ0) # Convert it to dense vector with complex element type
         to_dense(_CType(T), copy(ψ0.data))
     else
-        to_dense(_CType(T), mat2vec(ket2dm(ψ0).data)) 
+        to_dense(_CType(T), mat2vec(ket2dm(ψ0).data))
     end
     L = L_evo.data
 
@@ -89,7 +89,7 @@ function mesolveProblem(
     tspan = (tlist[1], tlist[end])
     prob = ODEProblem{getVal(inplace),FullSpecialize}(L, ρ0, tspan, params; kwargs3...)
 
-    return TimeEvolutionProblem(prob, tlist, L_evo.dimensions, (isoperket=Val(isoperket(ψ0)),))
+    return TimeEvolutionProblem(prob, tlist, L_evo.dimensions, (isoperket = Val(isoperket(ψ0)),))
 end
 
 @doc raw"""

src/time_evolution/smesolve.jl

@@ -104,7 +104,7 @@ function smesolveProblem(
     ρ0 = if isoperket(ψ0) # Convert it to dense vector with complex element type
         to_dense(_CType(T), copy(ψ0.data))
     else
-        to_dense(_CType(T), mat2vec(ket2dm(ψ0).data)) 
+        to_dense(_CType(T), mat2vec(ket2dm(ψ0).data))
     end
 
     progr = ProgressBar(length(tlist), enable = getVal(progress_bar))
@@ -148,7 +148,7 @@ function smesolveProblem(
         kwargs3...,
     )
 
-    return TimeEvolutionProblem(prob, tlist, dims, (isoperket=Val(isoperket(ψ0)),))
+    return TimeEvolutionProblem(prob, tlist, dims, (isoperket = Val(isoperket(ψ0)),))
 end
 
 @doc raw"""

test/core-test/time_evolution.jl

@@ -159,9 +159,29 @@
         sol_sme3 = smesolve(H, ψ0, tlist, c_ops_sme2, sc_ops_sme2, e_ops = e_ops, progress_bar = Val(false))
 
         # For testing the `OperatorKet` input
-        sol_me4 = mesolve(H, operator_to_vector(ket2dm(ψ0)), tlist, c_ops, saveat=saveat, progress_bar = Val(false))
-        sol_sme4 = smesolve(H, ψ0, tlist, c_ops_sme, sc_ops_sme, saveat=saveat, ntraj=10, progress_bar = Val(false), rng = MersenneTwister(12))
-        sol_sme5 = smesolve(H, operator_to_vector(ket2dm(ψ0)), tlist, c_ops_sme, sc_ops_sme, saveat=saveat, ntraj=10, progress_bar = Val(false), rng = MersenneTwister(12))
+        sol_me4 = mesolve(H, operator_to_vector(ket2dm(ψ0)), tlist, c_ops, saveat = saveat, progress_bar = Val(false))
+        sol_sme4 = smesolve(
+            H,
+            ψ0,
+            tlist,
+            c_ops_sme,
+            sc_ops_sme,
+            saveat = saveat,
+            ntraj = 10,
+            progress_bar = Val(false),
+            rng = MersenneTwister(12),
+        )
+        sol_sme5 = smesolve(
+            H,
+            operator_to_vector(ket2dm(ψ0)),
+            tlist,
+            c_ops_sme,
+            sc_ops_sme,
+            saveat = saveat,
+            ntraj = 10,
+            progress_bar = Val(false),
+            rng = MersenneTwister(12),
+        )
 
         ρt_mc = [ket2dm.(normalize.(states)) for states in sol_mc_states.states]
         expect_mc_states = mapreduce(states -> expect.(Ref(e_ops[1]), states), hcat, ρt_mc)
@@ -208,7 +228,9 @@
         @test isnothing(sol_sme.measurement)
         @test size(sol_sse2.measurement) == (length(c_ops), 20, length(tlist) - 1)
         @test size(sol_sme2.measurement) == (length(sc_ops_sme), 20, length(tlist) - 1)
-        @test all([sol_sme4.states[j][i] ≈ vector_to_operator(sol_sme5.states[j][i]) for i in eachindex(saveat), j in 1:10])
+        @test all([
+            sol_sme4.states[j][i] ≈ vector_to_operator(sol_sme5.states[j][i]) for i in eachindex(saveat), j in 1:10
+        ])
 
         @test sol_me_string ==
               "Solution of time evolution\n" *