out1[2] = x1010 out1[3] = x1011 out1[4] = x1012 out1[5] = x1013 out1[6] = x1014 out1[7] = x1015 out1[8] = x1016 } // p521Square squares a field element in the Montgomery domain. // // Preconditions: // // 0 ≤ eval arg1 < m // // Postconditions: // // eval (from_montgomery out1) mod m = (eval (from_montgomery arg1) * eval (from_montgomery arg1)) mod m // 0 ≤ eval out1 < m func p521Square(out1 *p521MontgomeryDomain, arg1 *p521MontgomeryDomain) { x1 := arg1[1] x2 := arg1[2] x3 := arg1[3] x4 := arg1[4]...