Function mod_sqrt 

pub fn mod_sqrt(
    x: &BigUint,
    modulus: &BigUint,
    non_qr: &BigUint,
) -> Option<BigUint>

Find the square root of x modulo modulus with non_qr a quadratic nonresidue of the field.