def METHOD_NAME(self):
    """
    Superoperator: Efficient numerics and reference return same result,
    acting on non-composite system
    """
    tol = 1e-10
    rho_list = list(map(rand_dm, [2, 3, 2, 3, 2]))
    rho_input = tensor(rho_list)
    superop = kraus_to_super(rand_kraus_map(3))
    analytic_result = rho_list
    analytic_result[1] = Qobj(
        vector_to_operator(
            superop @ operator_to_vector(analytic_result[1])
        ))
    analytic_result[3] = Qobj(
        vector_to_operator(
            superop @ operator_to_vector(analytic_result[3])
        ))
    analytic_result = tensor(analytic_result)
    naive_result = subsystem_apply(rho_input, superop,
                                   [False, True, False, True, False],
                                   reference=True)
    naive_diff = (analytic_result - naive_result).full()
    naive_diff_norm = norm(naive_diff)
    assert naive_diff_norm < tol
    efficient_result = subsystem_apply(rho_input, superop,
                                       [False, True, False, True, False])
    efficient_diff = (efficient_result - analytic_result).full()
    efficient_diff_norm = norm(efficient_diff)
    assert efficient_diff_norm < tol