pd.if.org Git - zpackage/blob - crypto/ref10/fe_invert.c

   1 #include <stdint.h>
   2
   3 #include "fe.h"
   4
   5 void fe_invert(int32_t out[10], int32_t z[10]) {
   6         int32_t t0[10];
   7         int32_t t1[10];
   8         int32_t t2[10];
   9         int32_t t3[10];
  10         int i;
  11
  12         /* qhasm: z2 = z1^2^1 */
  13         fe_sq(t0,z); for (i = 1;i < 1;++i) fe_sq(t0,t0);
  14
  15         /* qhasm: z8 = z2^2^2 */
  16         fe_sq(t1,t0); for (i = 1;i < 2;++i) fe_sq(t1,t1);
  17
  18         /* qhasm: z9 = z1*z8 */
  19         fe_mul(t1,z,t1);
  20
  21         /* qhasm: z11 = z2*z9 */
  22         fe_mul(t0,t0,t1);
  23
  24         /* qhasm: z22 = z11^2^1 */
  25         fe_sq(t2,t0); for (i = 1;i < 1;++i) fe_sq(t2,t2);
  26
  27         /* qhasm: z_5_0 = z9*z22 */
  28         fe_mul(t1,t1,t2);
  29
  30         /* qhasm: z_10_5 = z_5_0^2^5 */
  31         fe_sq(t2,t1); for (i = 1;i < 5;++i) fe_sq(t2,t2);
  32
  33         /* qhasm: z_10_0 = z_10_5*z_5_0 */
  34         fe_mul(t1,t2,t1);
  35
  36         /* qhasm: z_20_10 = z_10_0^2^10 */
  37         fe_sq(t2,t1); for (i = 1;i < 10;++i) fe_sq(t2,t2);
  38
  39         /* qhasm: z_20_0 = z_20_10*z_10_0 */
  40         fe_mul(t2,t2,t1);
  41
  42         /* qhasm: z_40_20 = z_20_0^2^20 */
  43         fe_sq(t3,t2); for (i = 1;i < 20;++i) fe_sq(t3,t3);
  44
  45         /* qhasm: z_40_0 = z_40_20*z_20_0 */
  46         fe_mul(t2,t3,t2);
  47
  48         /* qhasm: z_50_10 = z_40_0^2^10 */
  49         fe_sq(t2,t2); for (i = 1;i < 10;++i) fe_sq(t2,t2);
  50
  51         /* qhasm: z_50_0 = z_50_10*z_10_0 */
  52         fe_mul(t1,t2,t1);
  53
  54         /* qhasm: z_100_50 = z_50_0^2^50 */
  55         fe_sq(t2,t1); for (i = 1;i < 50;++i) fe_sq(t2,t2);
  56
  57         /* qhasm: z_100_0 = z_100_50*z_50_0 */
  58         fe_mul(t2,t2,t1);
  59
  60         /* qhasm: z_200_100 = z_100_0^2^100 */
  61         fe_sq(t3,t2); for (i = 1;i < 100;++i) fe_sq(t3,t3);
  62
  63         /* qhasm: z_200_0 = z_200_100*z_100_0 */
  64         fe_mul(t2,t3,t2);
  65
  66         /* qhasm: z_250_50 = z_200_0^2^50 */
  67         fe_sq(t2,t2); for (i = 1;i < 50;++i) fe_sq(t2,t2);
  68
  69         /* qhasm: z_250_0 = z_250_50*z_50_0 */
  70         fe_mul(t1,t2,t1);
  71
  72         /* qhasm: z_255_5 = z_250_0^2^5 */
  73         fe_sq(t1,t1); for (i = 1;i < 5;++i) fe_sq(t1,t1);
  74
  75         /* qhasm: z_255_21 = z_255_5*z11 */
  76         fe_mul(out,t1,t0);
  77
  78         return;
  79 }