Karim Belabas on Wed, 28 Jan 2004 14:43:00 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: lindep |
* <Denis.Simon@math.unicaen.fr> [2004-01-28 11:14]: > maybe my question should be asked in pari-users, > but I can't understand the behaviour of lindep: > ? lindep([0.2527546621, -0.89865858,1.1733994]) > %1 = 66984134.27215961355 > > Here, I was expecting 0 or a vector ??? > > In the big help, it is written that I can use > lindep with a matrix if flag <0, but I wasn't able > to do it ! > > ? lindep(matrix(2,2,i,j,1)*1.23456789,-1) > *** incorrect type in lindep. > ? lindep(matrix(2,2,i,j,1)*1.23456789,-2) > %2 = [1.234567890000000000000000000, -1.234567890000000000000000000]~ > > Here, I was expecting a relation over Z !! > > ? lindep(matrix(2,2,i,j,1)*1.23456789,-3) > *** incorrect type in pslq. > ? lindep(matrix(2,2,i,j,1)*1.23456789,-4) > *** incorrect type in lindep2. > ? lindep(matrix(2,2,i,j,1)*1.23456789,-10) > *** incorrect type in lindep2. The help was completely outdated. I have fixed the long and short help in CVS. Note: the PSLQ implementation is not currently in a satisfactory shape. The original contributor (Henri Cohen) will not maintain that code. I have fixed some problems and rewrote part of it, but I won't invest much time in it either. (I don't use it at all and am not that interested in the associated maths.) The LLL option is currently, more flexible, reliable and much faster than PSLQ. This will probably change when (if) somebody improves the PSLQ code... I have added the "standard" example (2^(1/6)+3^(1/5), from Borwein & Lisonek's old paper) to algdep doc. Karim. -- Karim Belabas Tel: (+33) (0)1 69 15 57 48 Dep. de Mathematiques, Bat. 425 Fax: (+33) (0)1 69 15 60 19 Universite Paris-Sud http://www.math.u-psud.fr/~belabas/ F-91405 Orsay (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]