The function divstepPrecomp returns the precomputed value for Bernstein-Yang-inversion (in montgomery form).
Postconditions: eval (from_montgomery out1) = ⌊(m - 1) / 2⌋^(if ⌊log2 m⌋ + 1 < 46 then ⌊(49 * (⌊log2 m⌋ + 1) + 80) / 17⌋ else ⌊(49 * (⌊log2 m⌋ + 1) + 57) / 17⌋) 0 ≤ eval out1 < m
Output Bounds: out1: [[0x0 ~> 0xffffffffffffffff], [0x0 ~> 0xffffffffffffffff], [0x0 ~> 0xffffffffffffffff], [0x0 ~> 0xffffffffffffffff], [0x0 ~> 0xffffffffffffffff], [0x0 ~> 0xffffffffffffffff]]
pub fn divstepPrecomp(out1: *[6]u64) void
pub fn divstepPrecomp(out1: *[6]u64) void {
@setRuntimeSafety(mode == .debug);
out1[0] = 0xfff69400fff18fff;
out1[1] = 0x2b7feffffd3ff;
out1[2] = 0xfffedbfffffe97ff;
out1[3] = 0x2840000002fff;
out1[4] = 0x6040000050400;
out1[5] = 0xfffc480000038000;
}