Re: foursquares-2.16.1.gp threesquares_fact

[Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>]

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