On Thursday, December 14, 2023 at 12:16:47 PM GMT+2, Bill Allombert <bill.allombert@math.u-bordeaux.fr> wrote:
On Wed, Dec 13, 2023 at 10:34:43PM -0800, Thomas D. Dean wrote:
> What purpose does F serve in qfsolve([matdiagonal([1,1,1,-n]),F])?
>
> ? threesquares_fact(n,F)
> %66 = [15, 0, 3]
> ? abs(qfsolve(matdiagonal([1,1,1,-n]))[1..3])~
> %67 = [15, 0, 3]
F is the factorization of discriminant of matdiagonal([1,1,1,-n]) which is n.
This is to avoid factoring n several time.
Cheers,
Bill