Bill Allombert on Tue, 11 Feb 2014 14:07:10 +0100

[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

 Re: Component  of prime ideal structure

```On Tue, Feb 11, 2014 at 11:33:48AM +0100, Jeroen Demeyer wrote:
> Consider
>
> gp> idealprimedec(nfinit(x^2+1), 5)
> %4 =
> [2 -1]
>
> [1  2]
>
> What is the meaning of this? The documentation ??idealprimedec says
> that it is supposed to be an algebraic integer, but I don't
> understand

Try idealprimedec(nfinit(x^2+1), 5)[,1]

> it. This changed from PARI-2.5 to PARI-2.6 and it breaks Denis
> Simon's GP scripts.

For reference, this is how Denis ell.gp script is using it:
valap = idealval(nf,a,p);
zlog = ideallog(nf,a*(nfbasistoalg(nf,p)/p.p)^valap,zinit)
...
f = idealval(nf,a,p);
a = nfalgtobasis(nf,a);
aaa = nfeltpow(nf,nfeltdiv(nf,a,p/p.p),f);
...
cont = idealval(nf,polcoeff(pol,0),p);
for( i = 1, deg,
if( cont, cont = min(cont,idealval(nf,polcoeff(pol,i),p))));
if( cont, pi = nfbasistoalg(nf,p/p.p));
if( cont > 1, pol *= pi^(2*(cont\2)));
...

Cheers,
Bill.

```