Fixed problem in the new liouvillian_generalized

qutip · Oct 7, 2023 · 2f06909 · 2f06909
1 parent feb2547
commit 2f06909
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/time_evolution/time_evolution.jl b/src/time_evolution/time_evolution.jl
@@ -735,16 +735,16 @@ function liouvillian_generalized(H::QuantumObject{<:AbstractArray, OperatorQuant
         end
 
         # Filter in the Hilbert space
-        σ = σ_filter == nothing ? 20 * sum(abs2, X_op) / length(X_op) : σ_filter
-        F1 = QuantumObject(gaussian.(Ω, 0, 2*κc), dims=dims)
+        σ = σ_filter == nothing ? 10 * sum(abs2, X_op) / length(X_op) : σ_filter
+        F1 = QuantumObject(gaussian.(Ω, 0, σ), dims=dims)
         F1 = dense_to_sparse(F1, tol)
 
         # Filter in the Liouville space
         M1 = ones(N_trunc, N_trunc)
         Ω1 = kron(Ω, M1)
         Ω2 = kron(M1, Ω)
         Ωdiff = Ω1 .- Ω2
-        F2 = QuantumObject(gaussian.(Ωdiff, 0, 2*κc), SuperOperatorQuantumObject, dims)
+        F2 = QuantumObject(gaussian.(Ωdiff, 0, σ), SuperOperatorQuantumObject, dims)
         F2 = dense_to_sparse(F2, tol)
 
         # Ohmic reservoir