Fix cut paste typo

qutip/solver/brmesolve.py

@@ -102,7 +102,7 @@ def brmesolve(H, psi0, tlist, a_ops=(), e_ops=(), c_ops=(),
             On `None` the states will be saved if no expectation operators are
             given.
         - | normalize_output : bool
-          | Normalize output state to hide ODE numerical errors. Only normalize 
+          | Normalize output state to hide ODE numerical errors. Only normalize
             the state if the initial state is already normalized.
         - | progress_bar : str {'text', 'enhanced', 'tqdm', ''}
           | How to present the solver progress.

qutip/solver/solver_base.py

@@ -99,7 +99,9 @@ def _prepare_state(self, state):
             norm = state.tr()
         else:
             norm = state.norm()
-        self._normalized = np.abs( - 1) <= settings.core["atol"]
+        # Use the settings atol instead of the solver one since the second
+        # refer to the ODE tolerance and some integrator do not use it.
+        self._normalized = np.abs(norm - 1) <= settings.core["atol"]
         if self.rhs.dims[1] == state.dims:
             return stack_columns(state.data)
         return state.data