Fix dynamical shifted fock

Author: albertomercurio

src/time_evolution/callback_helpers/callback_helpers.jl

@@ -32,8 +32,6 @@ function _generate_save_callback(e_ops, tlist, progress_bar, method)
     return PresetTimeCallback(tlist, _save_affect!, save_positions = (false, false))
 end
 
-_generate_mesolve_e_op(op) = mat2vec(adjoint(get_data(op)))
-
 ##
 
 # When e_ops is Nothing. Common for all solvers

src/time_evolution/callback_helpers/mesolve_callback_helpers.jl

@@ -39,3 +39,5 @@ function _mesolve_callbacks_new_e_ops!(integrator::AbstractODEIntegrator, e_ops)
         return nothing
     end
 end
+
+_generate_mesolve_e_op(op) = mat2vec(adjoint(get_data(op)))

src/time_evolution/time_evolution_dynamical.jl

@@ -522,7 +522,7 @@ function _DSF_mcsolve_Affect!(integrator)
     # e_ops0 = params.e_ops
     # c_ops0 = params.c_ops
 
-    e_ops0 = _get_e_ops(integrator, SaveFuncMCSEsolve)
+    e_ops0 = _get_e_ops(integrator, SaveFuncMCSolve)
     c_ops0, c_ops0_herm = _mcsolve_get_c_ops(integrator)
 
     copyto!(ψt, integrator.u)