| Phil Carmody on Thu, 27 Jul 2006 09:17:14 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| divisors/factor problem in 2.4.0 (development CHANGES-1.1570) |
In some contexts, divisors(3)=[1], and factor(3)=[;] alas.
Not had time to check other versions or anything.
Once you have a value that refuses to factor, it 'infects' all other values it
touches (so factor(round(broken*210.0))=[;] too).
Divisors() not including the constant factors I can live with, but being unable
to put them back manually is shall we say more of an inconvenience.
bash-3.1$ gp
GP/PARI CALCULATOR Version 2.4.0 (development CHANGES-1.1570)
i686 running linux (ix86/GMP-4.1.4 kernel) 32-bit version
compiled: Jul 22 2006, gcc-4.1.2 20060613 (prerelease) (Debian 4.1.1-5)
(readline v5.1 enabled, extended help enabled)
...
parisize = 4000000, primelimit = 500000
? pdiv
%60 = [1, x, x^2, x^3]
? poly=3*x^3
%61 = 3*x^3
? pdiv=divisors(poly)
%62 = [1, x, x^2, x^3]
? const=poly/pdiv[length(pdiv)]
%63 = 3
? cdiv=divisors(const)
%64 = [1]
? factor(const)
%65 = [;]
() ASCII ribbon campaign () Hopeless ribbon campaign
/\ against HTML mail /\ against gratuitous bloodshed
[stolen with permission from Daniel B. Cristofani]
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com