89 typedef uint128_t limb;  typedef
90 typedef limb felem[NLIMBS];
91 typedef limb longfelem[NLIMBS * 2];
215 #define two105m41m9 (((limb)1) << 105) - (((limb)1) << 41) - (((limb)1) << 9)
216 #define two105 (((limb)1) << 105)
217 #define two105m41p9 (((limb)1) << 105) - (((limb)1) << 41) + (((limb)1) << 9)
251 #define two107m43m11 (((limb)1) << 107) - (((limb)1) << 43) - (((limb)1) << 11)
252 #define two107 (((limb)1) << 107)
253 #define two107m43p11 (((limb)1) << 107) - (((limb)1) << 43) + (((limb)1) << 11)
283   static const limb two70m8p6 =  in longfelem_diff()
284       (((limb)1) << 70) - (((limb)1) << 8) + (((limb)1) << 6);  in longfelem_diff()
285   static const limb two70p40 = (((limb)1) << 70) + (((limb)1) << 40);  in longfelem_diff()
286   static const limb two70 = (((limb)1) << 70);  in longfelem_diff()
287   static const limb two70m40m38p6 = (((limb)1) << 70) - (((limb)1) << 40) -  in longfelem_diff()
288                                     (((limb)1) << 38) + (((limb)1) << 6);  in longfelem_diff()
289   static const limb two70m6 = (((limb)1) << 70) - (((limb)1) << 6);  in longfelem_diff()
312 #define two64m0 (((limb)1) << 64) - 1
313 #define two110p32m0 (((limb)1) << 110) + (((limb)1) << 32) - 1
314 #define two64m46 (((limb)1) << 64) - (((limb)1) << 46)
315 #define two64m32 (((limb)1) << 64) - (((limb)1) << 32)
347   tmp[3] += ((limb)a) << 32;  in felem_shrink()
355   tmp[3] += ((limb)a) << 32;  in felem_shrink()
361   tmp[1] -= (((limb)b) << 32);  in felem_shrink()
420   limb a;  in smallfelem_square()
509   limb a;  in smallfelem_mul()
635 #define two100m36m4 (((limb)1) << 100) - (((limb)1) << 36) - (((limb)1) << 4)
636 #define two100 (((limb)1) << 100)
637 #define two100m36p4 (((limb)1) << 100) - (((limb)1) << 36) + (((limb)1) << 4)
826 static limb smallfelem_is_zero(const smallfelem small) {  in smallfelem_is_zero()
827   limb result;  in smallfelem_is_zero()
854   result |= ((limb)is_zero) << 64;  in smallfelem_is_zero()
859   return (int)(smallfelem_is_zero(small) & ((limb)1));  in smallfelem_is_zero_int()
1078 static void copy_conditional(felem out, const felem in, limb mask) {  in copy_conditional()
1081     const limb tmp = mask & (in[i] ^ out[i]);  in copy_conditional()
1087 static void copy_small_conditional(felem out, const smallfelem in, limb mask) {  in copy_small_conditional()
1091     out[i] = ((limb)(in[i] & mask64)) | (out[i] & ~mask);  in copy_small_conditional()
1112   limb x_equal, y_equal, z1_is_zero, z2_is_zero;  in point_add()
1581         copy_small_conditional(ftmp, tmp[1], (((limb)sign) - 1));  in batch_mul()