Add tests, changelog and format code

Author: albertomercurio

CHANGELOG.md

@@ -9,6 +9,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 - Fix Dynamical Fock Dimension states saving due to wrong saving of dimensions. ([#375])
 - Support a list of observables for `expect`. ([#374], [#376])
+- Add checks for `tlist` in time evolution solvers. The checks are to ensure that `tlist` is not empty, the elements are in increasing order, and the elements are unique. ([#378])
 
 ## [v0.25.0]
 Release date: 2025-01-20
@@ -89,3 +90,4 @@ Release date: 2024-11-13
 [#374]: https://github.com/qutip/QuantumToolbox.jl/issues/374
 [#375]: https://github.com/qutip/QuantumToolbox.jl/issues/375
 [#376]: https://github.com/qutip/QuantumToolbox.jl/issues/376
+[#378]: https://github.com/qutip/QuantumToolbox.jl/issues/378

src/time_evolution/time_evolution.jl

@@ -213,7 +213,6 @@ struct DiscreteLindbladJumpCallback <: LindbladJumpCallbackType end
 ContinuousLindbladJumpCallback(; interp_points::Int = 10) = ContinuousLindbladJumpCallback(interp_points)
 
 function _check_tlist(tlist, T::Type)
-    
     tlist = convert(Vector{T}, tlist) # Convert it to support GPUs and avoid type instabilities for OrdinaryDiffEq.jl
 
     # Check if the list of times is not empty

test/core-test/time_evolution.jl

@@ -63,6 +63,13 @@
               "abstol = $(sol2.abstol)\n" *
               "reltol = $(sol2.reltol)\n"
 
+        tlist1 = Float64[]
+        tlist2 = [0, 0.2, 0.1]
+        tlist3 = [0, 0.1, 0.1, 0.2]
+        @test_throws ArgumentError sesolve(H, ψ0, tlist1, progress_bar = Val(false))
+        @test_throws ArgumentError sesolve(H, ψ0, tlist2, progress_bar = Val(false))
+        @test_throws ArgumentError sesolve(H, ψ0, tlist3, progress_bar = Val(false))
+
         @testset "Memory Allocations" begin
             allocs_tot = @allocations sesolve(H, ψ0, tlist, e_ops = e_ops, progress_bar = Val(false)) # Warm-up
             allocs_tot = @allocations sesolve(H, ψ0, tlist, e_ops = e_ops, progress_bar = Val(false))
@@ -190,6 +197,19 @@
               "abstol = $(sol_sse.abstol)\n" *
               "reltol = $(sol_sse.reltol)\n"
 
+        tlist1 = Float64[]
+        tlist2 = [0, 0.2, 0.1]
+        tlist3 = [0, 0.1, 0.1, 0.2]
+        @test_throws ArgumentError mesolve(H, ψ0, tlist1, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError mesolve(H, ψ0, tlist2, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError mesolve(H, ψ0, tlist3, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError mcsolve(H, ψ0, tlist1, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError mcsolve(H, ψ0, tlist2, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError mcsolve(H, ψ0, tlist3, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError ssesolve(H, ψ0, tlist1, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError ssesolve(H, ψ0, tlist2, c_ops, progress_bar = Val(false))
+        @test_throws ArgumentError ssesolve(H, ψ0, tlist3, c_ops, progress_bar = Val(false))
+
         # Time-Dependent Hamiltonian
         # ssesolve is slow to be run on CI. It is not removed from the test because it may be useful for testing in more powerful machines.