Function square_redc 

pub fn square_redc<const N: usize>(
    a: [u64; N],
    modulus: [u64; N],
    inv: u64,
) -> [u64; N]

⚠️ Computes a^2 * 2^(-BITS) mod modulus

⚠️ This function is not part of the stable API.

Requires that `inv` is the inverse of `-modulus[0]` modulo `2^64`. Requires that `a` is less than `modulus`.