]> pd.if.org Git - zpackage/blob - crypto/ref10/fe_invert.c
add test vectors for zpm-sign
[zpackage] / 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 }