fn test_sqrt_qfe_2() {
        let x = super::BN254TwistCurveFieldElement::new([
            BN254FieldElement::from_hex_unchecked(
                "3010c68cb50161b7d1d96bb71edfec9880171954e56871abf3d93cc94d745fa1",
            ),
            BN254FieldElement::from_hex_unchecked(
                "0476be093a6d2b4bbf907172049874af11e1b6267606e00804d3ff0037ec57fd",
            ),
        ]);

        let qfe_b = BN254TwistCurve::b();

        let y_square = x.pow(3_u64) + qfe_b;
        let y = super::sqrt_qfe(&y_square, 0).unwrap();

        let y_expected = super::BN254TwistCurveFieldElement::new([
            BN254FieldElement::from_hex_unchecked(
                "01b33461f39d9e887dbb100f170a2345dde3c07e256d1dfa2b657ba5cd030427",
            ),
            BN254FieldElement::from_hex_unchecked(
                "14c059d74e5b6c4ec14ae5864ebe23a71781d86c29fb8fb6cce94f70d3de7a21",
            ),
        ]);

        let value_y = y.value();
        let value_y_expected = y_expected.value();

        assert_eq!(value_y[0].clone(), value_y_expected[0].clone());
        assert_eq!(value_y[1].clone(), value_y_expected[1].clone());
    }
We took the q1 point of the test two_pairs_of_points_match_2 from pairing.rs
crates/math/src/elliptic_curve/short_weierstrass/curves/bn_254/sqrt.rs