Minor fixes

Author: albertomercurio

Project.toml

@@ -1,7 +1,7 @@
 name = "QuantumToolbox"
 uuid = "6c2fb7c5-b903-41d2-bc5e-5a7c320b9fab"
 authors = ["Alberto Mercurio", "Yi-Te Huang"]
-version = "0.33.1"
+version = "0.33.0"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
@@ -61,7 +61,7 @@ OrdinaryDiffEqTsit5 = "1"
 Pkg = "1"
 Random = "1"
 SciMLBase = "2"
-SciMLOperators = "1"
+SciMLOperators = "0.3, 0.4, 1"
 SparseArrays = "1"
 SpecialFunctions = "2"
 StaticArraysCore = "1"

src/qobj/quantum_object_evo.jl

@@ -512,12 +512,7 @@ function (A::QuantumObjectEvolution)(
         )
     end
 
-    # We put `nothing` in the place of `u` because the time-dependence doesn't depend on the state
-    ## TODO: should use SciMLOperators v1 API: A.data(ψout.data, ψin.data, nothing, p, t)
-    ## But somehow that doesn't work for some specific cases in our runtests.
-    ## Therefore, we use the most basic way to do it for now
-    update_coefficients!(A.data, nothing, p, t)
-    mul!(ψout.data, A.data, ψin.data)
+    A.data(ψout.data, ψin.data, nothing, p, t)
 
     return ψout
 end