idealadd() details

[Jeroen Demeyer <jdemeyer@cage.ugent.be>]

Hello list,

I have a question about idealadd().  Either I'm misunderstanding 
something, or the implementation is wrong, or the documentation is wrong.

Consider the following:
gp> idealadd(nfinit(t^2+1), [1,0;0,2], [1,0;0,2])
%20 =
[1 0]

[0 1]

Here I expected to get [1,0;0,2] as a result, since the documentation 
states it just adds Z-modules.  Of course, as ideals, the answer is correct.


Cheers,
Jeroen