in pq-crypto/sike_r3/sikep434r3_fpx.c [217:359]
static void fpinv_chain_mont(felm_t a)
{
unsigned int i, j;
felm_t t[31], tt;
/* Precomputed table */
fpsqr_mont(a, tt);
fpmul_mont(a, tt, t[0]);
for (i = 0; i <= 29; i++) {
fpmul_mont(t[i], tt, t[i + 1]);
}
fpcopy(a, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[5], tt, tt);
for (i = 0; i < 10; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[14], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[3], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[23], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[13], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[24], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[7], tt, tt);
for (i = 0; i < 8; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[12], tt, tt);
for (i = 0; i < 8; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[30], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[1], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[30], tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[21], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[2], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[19], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[1], tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[24], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[26], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[16], tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[10], tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[6], tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[0], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[20], tt, tt);
for (i = 0; i < 8; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[9], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[25], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[30], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[26], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(a, tt, tt);
for (i = 0; i < 7; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[28], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[6], tt, tt);
for (i = 0; i < 6; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[10], tt, tt);
for (i = 0; i < 9; i++) {
fpsqr_mont(tt, tt);
}
fpmul_mont(t[22], tt, tt);
for (j = 0; j < 35; j++) {
for (i = 0; i < 6; i++) fpsqr_mont(tt, tt);
fpmul_mont(t[30], tt, tt);
}
fpcopy(tt, a);
}