pub fn phase_estimation(theta: f64, precision_bits: u8) -> (f64, f64)
Estimate phase using quantum phase estimation