Make changelog, add example on docs and remove spell check errors

Author: albertomercurio

CHANGELOG.md

@@ -14,6 +14,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 - Add more generic solver for `steadystate_floquet` to allow more linear solvers. ([#396])
 - Fix time evolution output when using `saveat` keyword argument. ([#398])
 - Improve ensemble generation of `ssesolve` and change parameters handling on stochastic processes. ([#403])
+- Introduce measurement on `ssesolve` and `smesolve`. ([#404])
 
 ## [v0.26.0]
 Release date: 2025-02-09
@@ -132,3 +133,4 @@ Release date: 2024-11-13
 [#398]: https://github.com/qutip/QuantumToolbox.jl/issues/398
 [#402]: https://github.com/qutip/QuantumToolbox.jl/issues/402
 [#403]: https://github.com/qutip/QuantumToolbox.jl/issues/403
+[#404]: https://github.com/qutip/QuantumToolbox.jl/issues/404

docs/src/users_guide/time_evolution/stochastic.md

@@ -105,12 +105,16 @@ sse_sol = ssesolve(
     sc_ops,
     e_ops = [x],
     ntraj = ntraj,
+    store_measurement = Val(true),
 )
 
+measurement_avg = sum(sse_sol.measurement, dims=2) / size(sse_sol.measurement, 2)
+measurement_avg = dropdims(measurement_avg, dims=2)
+
 # plot by CairoMakie.jl
 fig = Figure(size = (500, 350))
 ax = Axis(fig[1, 1], xlabel = "Time")
-#lines!(ax, tlist, real(sse_sol.xxxxxx), label = L"J_x", color = :red, linestyle = :solid) TODO: add this in the future
+lines!(ax, tlist[1:end-1], real(measurement_avg[1,:]), label = L"J_x", color = :red, linestyle = :solid)
 lines!(ax, tlist, real(sse_sol.expect[1,:]),  label = L"\langle x \rangle", color = :black, linestyle = :solid)
 
 axislegend(ax, position = :rt)
@@ -134,12 +138,16 @@ sme_sol = smesolve(
     sc_ops,
     e_ops = [x],
     ntraj = ntraj,
+    store_measurement = Val(true),
 )
 
+measurement_avg = sum(sme_sol.measurement, dims=2) / size(sme_sol.measurement, 2)
+measurement_avg = dropdims(measurement_avg, dims=2)
+
 # plot by CairoMakie.jl
 fig = Figure(size = (500, 350))
 ax = Axis(fig[1, 1], xlabel = "Time")
-#lines!(ax, tlist, real(sme_sol.xxxxxx), label = L"J_x", color = :red, linestyle = :solid) TODO: add this in the future
+lines!(ax, tlist[1:end-1], real(measurement_avg[1,:]), label = L"J_x", color = :red, linestyle = :solid)
 lines!(ax, tlist, real(sme_sol.expect[1,:]),  label = L"\langle x \rangle", color = :black, linestyle = :solid)
 
 axislegend(ax, position = :rt)

src/time_evolution/callback_helpers/smesolve_callback_helpers.jl

@@ -5,15 +5,15 @@ Helper functions for the smesolve callbacks. Almost equal to the mesolve case, b
 struct SaveFuncSMESolve{
     SM,
     TE,
-    TME,
+    TE2,
     PT<:Union{Nothing,ProgressBar},
     IT,
     TEXPV<:Union{Nothing,AbstractMatrix},
     TMEXPV<:Union{Nothing,AbstractMatrix},
 } <: AbstractSaveFunc
     store_measurement::Val{SM}
     e_ops::TE
-    m_ops::TME
+    m_ops::TE2
     progr::PT
     iter::IT
     expvals::TEXPV

src/time_evolution/callback_helpers/ssesolve_callback_helpers.jl

@@ -5,15 +5,15 @@ Helper functions for the ssesolve callbacks. Almost equal to the sesolve case, b
 struct SaveFuncSSESolve{
     SM,
     TE,
-    TME,
+    TE2,
     PT<:Union{Nothing,ProgressBar},
     IT,
     TEXPV<:Union{Nothing,AbstractMatrix},
     TMEXPV<:Union{Nothing,AbstractMatrix},
 } <: AbstractSaveFunc
     store_measurement::Val{SM}
     e_ops::TE
-    m_ops::TME
+    m_ops::TE2
     progr::PT
     iter::IT
     expvals::TEXPV