Steady state

quantum/mesolve.py

@@ -284,7 +284,6 @@ def mesolve(H, rho0, tlist, c_ops=[], e_ops=[], args={}, options=None,
     #check if should use OPENMP
     check_use_openmp(options)
 
-
     res = None
 
     #

quantum/steadystate.py

@@ -960,7 +960,7 @@ def _pseudo_inverse_dense(L, rhoss, w=None, **pseudo_args):
         L = L
     else:
         L = 1.0j * w * spre(tr_mat) + L
-        # It's possible that there's an error here!
+    # It's possible that there's an error here because of the use of the spre function.
 
     if pseudo_args['method'] == 'direct':
         try:

quantum/superoperator.py

@@ -197,6 +197,8 @@ def spre(A):
 
 # In the following, the uncommented lines correspond to the QuTip original
 # assignations
+# The creation of a vector from a matrix in QuTip approach is ordering in
+# columns, while our approach is to order in rows.
 #------------------------------------------------------------------------------
 def operator_to_vector(op):
     """