From 2d95a8aa87b0133e84511ddb81543861f403b7f8 Mon Sep 17 00:00:00 2001 From: Edgar Costa Date: Wed, 10 Jul 2024 11:19:23 -0400 Subject: [PATCH] a more reasonable sorting of complex numbers --- src/acb/vec_sort_pretty.c | 16 +++------------- 1 file changed, 3 insertions(+), 13 deletions(-) diff --git a/src/acb/vec_sort_pretty.c b/src/acb/vec_sort_pretty.c index 0660eb1d22..1f50339891 100644 --- a/src/acb/vec_sort_pretty.c +++ b/src/acb/vec_sort_pretty.c @@ -25,24 +25,14 @@ int acb_cmp_pretty(const acb_t a, const acb_t b) arb_t t, u, v; int res; arb_init(t); - arb_init(u); - arb_init(v); - arb_abs(u, acb_imagref(a)); - arb_abs(v, acb_imagref(b)); - arb_sub(t, u, v, MAG_BITS); + arb_sub(t, acb_realref(a), acb_realref(b), MAG_BITS); res = 0; if (arb_contains_zero(t)) { - arb_sub(t, acb_realref(a), acb_realref(b), MAG_BITS); - res = arb_is_positive(t) ? 1 : -1; - } - else - { - res = arb_is_positive(t) ? 1 : -1; + arb_sub(t, acb_imagref(a), acb_imagref(b), MAG_BITS); } + res = arb_is_positive(t) ? 1 : -1; arb_clear(t); - arb_clear(u); - arb_clear(v); return res; }