Function square_repeat_montgomery_form 

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

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

This method is variable time in n.