Reset previous e_ops method for mesolve

Author: albertomercurio

src/time_evolution/callback_helpers/callback_helpers.jl

@@ -34,7 +34,7 @@ _get_e_ops_data(e_ops, ::Type{SaveFuncSESolve}) = get_data.(e_ops)
 _get_e_ops_data(e_ops, ::Type{SaveFuncMESolve}) = [_generate_mesolve_e_op(op) for op in e_ops] # Broadcasting generates type instabilities on Julia v1.10
 _get_e_ops_data(e_ops, ::Type{SaveFuncSSESolve}) = get_data.(e_ops)
 
-_generate_mesolve_e_op(op) = mat2vec(get_data(op))
+_generate_mesolve_e_op(op) = mat2vec(adjoint(get_data(op)))
 
 #=
 This function add the normalization callback to the kwargs. It is needed to stabilize the integration when using the ssesolve method.

src/time_evolution/callback_helpers/mesolve_callback_helpers.jl

@@ -21,7 +21,7 @@ function _save_func_mesolve(integrator, e_ops, progr, iter, expvals)
     # to reshape u to make it a matrix, but we reshape the e_ops once.
 
     ρ = integrator.u
-    _expect = op -> dot(ρ, op)
+    _expect = op -> dot(op, ρ)
     @. expvals[:, iter[]] = _expect(e_ops)
     iter[] += 1
 

src/time_evolution/smesolve.jl

@@ -103,7 +103,7 @@ function smesolveProblem(
     D_l = map(sc_ops_evo_data) do op
         # TODO: Implement the three-argument dot function for SciMLOperators.jl
         # Currently, we are assuming a time-independent MatrixOperator
-        op_vec = mat2vec(op.A)
+        op_vec = mat2vec(adjoint(op.A))
         return _spre(op, Id) + _spost(op', Id) + _smesolve_ScalarOperator(op_vec) * IdentityOperator(prod(dims)^2)
     end
     D = DiffusionOperator(D_l)