def METHOD_NAME():
    qubits = cirq.LineQubit.range(1)
    circuit = cirq.Circuit(cirq.I(qubits[0]))
    no_noise_simulator = cirq.DensityMatrixSimulator()
    _, intermediate_result = dfe.direct_fidelity_estimation(
        circuit, qubits, no_noise_simulator, n_measured_operators=1, samples_per_term=0
    )
    # We only test a few fields to be sure that they are set properly. In
    # particular, some of them are random, and so we don't test them.
    assert str(intermediate_result.clifford_tableau) == "+ Z "
    np.testing.assert_equal(len(intermediate_result.pauli_traces), 1)
    assert np.isclose(intermediate_result.pauli_traces[0].rho_i, 1.0)
    assert np.isclose(intermediate_result.pauli_traces[0].Pr_i, 0.5)
    np.testing.assert_equal(len(intermediate_result.trial_results), 1)
    assert np.isclose(intermediate_result.trial_results[0].sigma_i, 1.0)
    assert np.isclose(intermediate_result.std_dev_estimate, 0.0)
    assert np.isclose(intermediate_result.std_dev_bound, 0.5)