in src/integer_arith/scalar.rs [198:210]
fn pow_mod(base: &Scalar, b: &Scalar, q: &Scalar) -> Scalar {
let bits: Vec<bool> = b.get_bits();
let mut res = Self::one();
res = Self::modulus(&res, q);
let mut pow = Scalar::new(base.rep);
for bit in bits.iter() {
if *bit {
res = Self::mul_mod(&res, &pow, q);
}
pow = Self::mul_mod(&pow, &pow, q);
}
res
}