diff options
Diffstat (limited to 'src/gen/gens.c')
| -rw-r--r-- | src/gen/gens.c | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/src/gen/gens.c b/src/gen/gens.c index f5f9fcb..d40de67 100644 --- a/src/gen/gens.c +++ b/src/gen/gens.c @@ -3,6 +3,8 @@ * Copyright (C) 2017 J08nY */ #include "gens.h" +#include <misc/types.h> +#include "exhaustive/arg.h" #include "point.h" static int gens_put(curve_t *curve, GEN generators, long len) { @@ -37,6 +39,44 @@ GENERATOR(gens_gen_one) { return gens_put(curve, generators, len); } +CHECK(gens_check_anomalous) { + if (cfg->field == FIELD_BINARY) return 1; + pari_sp ltop = avma; + for (size_t i = 0; i < curve->ngens; ++i) { + if (mpcmp(curve->field, curve->generators[i]->order) == 0) { + avma = ltop; + return -5; + } + } + return 1; +} + +CHECK(gens_check_embedding) { + HAS_ARG(args); + if (cfg->field == FIELD_BINARY) return 1; + pari_sp ltop = avma; + + const char *min_degree = args->args; + GEN mind = strtoi(min_degree); + + for (size_t i = 0; i < curve->ngens; ++i) { + GEN power = gen_0; + GEN pm; + do { + power = addii(power, gen_1); + GEN ppow = powii(curve->field, power); + pm = subii(ppow, gen_1); + } while (!dvdii(pm, curve->generators[i]->order)); + + if (mpcmp(power, mind) <= 0) { + avma = ltop; + return -5; + } + } + avma = ltop; + return 1; +} + UNROLL(gens_unroll) { if (curve->generators) { points_free_deep(&curve->generators, curve->ngens); |