aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorJ08nY2025-03-21 12:00:05 +0100
committerJ08nY2025-03-21 12:00:05 +0100
commit4bf51fc47b31cf03fd725fb90cc3ae53d178d172 (patch)
treeac0ee926d7eb59b4bd466db77eb6c4a02746a0b4 /src
parent0e91e62d6b9bc1fb7cbdea5b09eb0794d3870f25 (diff)
downloadecgen-4bf51fc47b31cf03fd725fb90cc3ae53d178d172.tar.gz
ecgen-4bf51fc47b31cf03fd725fb90cc3ae53d178d172.tar.zst
ecgen-4bf51fc47b31cf03fd725fb90cc3ae53d178d172.zip
Fix addprimes call in cm_any.
Diffstat (limited to 'src')
-rw-r--r--src/cm/cm.c9
-rw-r--r--src/gen/gens.c1
2 files changed, 8 insertions, 2 deletions
diff --git a/src/cm/cm.c b/src/cm/cm.c
index 79e0e0e..04b988a 100644
--- a/src/cm/cm.c
+++ b/src/cm/cm.c
@@ -181,10 +181,15 @@ static int cm_init(exhaustive_t *setup) {
for (long j = 1; j <= len; ++j) {
GEN factor = gel(factors, j);
if (isprime(factor)) {
+ debug_log("Adding prime factor %Pi", factor);
addprimes(gtovec(factor));
} else {
- GEN factored = Z_factor(order);
- addprimes(gel(factored, 1));
+ GEN factored = Z_factor(factor);
+ GEN primes = gel(factored, 1);
+ for (long k = 1; k <= glength(primes); ++k) {
+ debug_log("Adding prime factor %Pi", gel(primes, k));
+ addprimes(gel(primes, k));
+ }
}
order = mulii(order, factor);
}
diff --git a/src/gen/gens.c b/src/gen/gens.c
index a32f8b8..329149a 100644
--- a/src/gen/gens.c
+++ b/src/gen/gens.c
@@ -38,6 +38,7 @@ GENERATOR(gens_gen_any) {
GENERATOR(gens_gen_one) {
pari_sp ltop = avma;
GEN group = ellff_get_group(curve->curve);
+ debug_log("Group structure %Ps.", group);
long len = glength(group);
if (len == 2) {
avma = ltop;