[no ci] Check tlist properties

Author: albertomercurio

src/time_evolution/lr_mesolve.jl

@@ -409,7 +409,7 @@ function lr_mesolveProblem(
     c_ops = get_data.(c_ops)
     e_ops = get_data.(e_ops)
 
-    t_l = convert(Vector{_FType(H)}, tlist)
+    t_l = _check_tlist(tlist, _FType(H))
 
     # Initialization of Arrays
     expvals = Array{ComplexF64}(undef, length(e_ops), length(t_l))

src/time_evolution/mcsolve.jl

@@ -167,7 +167,7 @@ function mcsolveProblem(
     c_ops isa Nothing &&
         throw(ArgumentError("The list of collapse operators must be provided. Use sesolveProblem instead."))
 
-    tlist = convert(Vector{_FType(ψ0)}, tlist) # Convert it to support GPUs and avoid type instabilities for OrdinaryDiffEq.jl
+    tlist = _check_tlist(tlist, _FType(ψ0))
 
     H_eff_evo = _mcsolve_make_Heff_QobjEvo(H, c_ops)
 

src/time_evolution/mesolve.jl

@@ -68,7 +68,7 @@ function mesolveProblem(
     haskey(kwargs, :save_idxs) &&
         throw(ArgumentError("The keyword argument \"save_idxs\" is not supported in QuantumToolbox."))
 
-    tlist = convert(Vector{_FType(ψ0)}, tlist) # Convert it to support GPUs and avoid type instabilities for OrdinaryDiffEq.jl
+    tlist = _check_tlist(tlist, _FType(ψ0))
 
     L_evo = _mesolve_make_L_QobjEvo(H, c_ops)
     check_dimensions(L_evo, ψ0)

src/time_evolution/sesolve.jl

@@ -59,7 +59,7 @@ function sesolveProblem(
     haskey(kwargs, :save_idxs) &&
         throw(ArgumentError("The keyword argument \"save_idxs\" is not supported in QuantumToolbox."))
 
-    tlist = convert(Vector{_FType(ψ0)}, tlist) # Convert it to support GPUs and avoid type instabilities for OrdinaryDiffEq.jl
+    tlist = _check_tlist(tlist, _FType(ψ0))
 
     H_evo = _sesolve_make_U_QobjEvo(H) # Multiply by -i
     isoper(H_evo) || throw(ArgumentError("The Hamiltonian must be an Operator."))

src/time_evolution/ssesolve.jl

@@ -162,7 +162,7 @@ function ssesolveProblem(
     sc_ops isa Nothing &&
         throw(ArgumentError("The list of collapse operators must be provided. Use sesolveProblem instead."))
 
-    tlist = convert(Vector{Float64}, tlist) # Convert it into Float64 to avoid type instabilities for StochasticDiffEq.jl
+    tlist = _check_tlist(tlist, _FType(ψ0))
 
     H_eff_evo = _mcsolve_make_Heff_QobjEvo(H, sc_ops)
     isoper(H_eff_evo) || throw(ArgumentError("The Hamiltonian must be an Operator."))

src/time_evolution/time_evolution.jl

@@ -212,6 +212,20 @@ 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
+    isempty(tlist) && throw(ArgumentError("The list of times must not be empty."))
+    # Check if the list of times is sorted
+    !issorted(tlist) && throw(ArgumentError("The list of times must be sorted."))
+    # Check if the list of times is unique
+    length(tlist) != length(unique(tlist)) && throw(ArgumentError("The list of times must be unique."))
+
+    return tlist
+end
+
 #######################################
 
 function liouvillian_floquet(