1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
/*
* ecgen, tool for generating Elliptic curve domain parameters
* Copyright (C) 2017 J08nY
*/
#include "order.h"
#include "io/input.h"
GEN order_factors(curve_t *curve, const config_t *cfg) {
if (cfg->prime) {
return gtovec(curve->order);
} else {
GEN factors = Z_factor(curve->order);
return gel(factors, 1);
}
}
GEN order_groups(curve_t *curve, const config_t *cfg, GEN factors) {
long nprimes = glength(factors);
if (cfg->prime) {
return gtovec(curve->order);
} else {
GEN amount = int2n(nprimes);
GEN groups = gtovec0(gen_0, itos(amount) - 1);
for (size_t count = 1; count < (size_t)(1 << nprimes); ++count) {
GEN result = gen_1;
for (long bit = 0; bit < nprimes; ++bit) {
size_t mask = (size_t)(1 << bit);
if (count & mask) {
result = mulii(result, gel(factors, bit + 1));
}
}
gel(groups, count) = result;
}
// TODO: sort this, as it is not necessarily sorted, in fact most likely
// not
return groups;
}
}
GENERATOR(order_gen_input) {
pari_sp ltop = avma;
GEN ord = input_int("order", cfg->bits);
if (gequalm1(ord)) {
avma = ltop;
return -4;
} else {
curve->order = ord;
obj_insert_shallow(curve->curve, 1, ord);
return 1;
}
}
GENERATOR(order_gen_any) {
GEN ord = ellff_get_card(curve->curve);
if (isclone(ord)) {
curve->order = gcopy(ord);
} else {
curve->order = ord;
}
return 1;
}
GENERATOR(order_gen_sea) {
pari_sp ltop = avma;
GEN order = ellsea(curve->curve, 0);
if (gequal0(order)) {
avma = ltop;
return -4;
} else {
curve->order = order;
obj_insert_shallow(curve->curve, 1, order);
return 1;
}
}
GENERATOR(order_gen_smallfact) {
if (!args) {
fprintf(stderr, "No args to an arged function. order_gen_smallfact\n");
return INT_MIN;
}
pari_ulong smallfact = *(pari_ulong *)args->args;
pari_sp ltop = avma;
GEN fact = mpfact(smallfact);
if (lgefint(fact) > 3) {
smallfact = 0;
} else {
smallfact = itou(fact);
}
GEN order = ellsea(curve->curve, smallfact);
if (gequal0(order) || gequal1(gcdii(order, fact))) {
avma = ltop;
return -4;
} else {
curve->order = order;
obj_insert_shallow(curve->curve, 1, curve->order);
return 1;
}
}
GENERATOR(order_gen_prime) {
pari_sp ltop = avma;
GEN order = ellsea(curve->curve, 1);
if (gequal0(order) || !(isprime(order))) {
avma = ltop;
return -4;
} else {
curve->order = order;
obj_insert_shallow(curve->curve, 1, curve->order);
return 1;
}
}
|
