Function square_montgomery_form 

pub(crate) const fn square_montgomery_form<const LIMBS: usize>(
    a: &Uint<LIMBS>,
    modulus: &Odd<Uint<LIMBS>>,
    mod_neg_inv: Limb,
) -> Uint<LIMBS>

Computes the Montgomery squaring of a modulo modulus where a is in Montgomery form.