Re: Component [5] of prime ideal structure

[Karim Belabas <Karim.Belabas@math.u-bordeaux.fr>]

* Jeroen Demeyer [2014-02-11 14:53]:
> On 2014-02-11 14:07, Bill Allombert wrote:
> >Try idealprimedec(nfinit(x^2+1), 5)[1][5][,1]
> 
> This doesn't work for
> 
> gp> idealprimedec(nfinit(x^2+1), 3)[1][5]
> %1 = 1
> 
> One could of course easily wrap this in a (member) function.

Here's a possible implementation that caters for all versions

  pr.b = my(b=pr[5]);if(type(b)=="t_MAT",b[,1],b);

Cheers,

    K.B.
--
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~kbelabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`