def METHOD_NAME():
    qubits = cirq.LineQubit.range(3)
    circuit = cirq.Circuit(cirq.Z(qubits[0]) ** 0.123, cirq.X(qubits[1]), cirq.X(qubits[2]))
    no_noise_simulator = cirq.DensityMatrixSimulator()
    estimated_fidelity, _ = dfe.direct_fidelity_estimation(
        circuit, qubits, no_noise_simulator, n_measured_operators=64, samples_per_term=0
    )
    assert np.isclose(estimated_fidelity, 1.0, atol=0.01)