| Paul van Wamelen on Tue, 23 May 2000 17:29:10 -0500 (CDT) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| qfminim bug 1 |
Dear Pari developers,
In qfminim with the flag equal to 2 (that is for dealing with
non-integral coefficients), if the bound given is less than one it
gets changed to 1:
GP/PARI CALCULATOR Version 2.0.19 (beta)
UltraSparc (MicroSparc kernel) 32-bit version
(readline v2.2 enabled, extended help available)
Copyright (C) 1989-1999 by
C. Batut, K. Belabas, D. Bernardi, H. Cohen and M. Olivier.
Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.
realprecision = 28 significant digits
seriesprecision = 16 significant terms
format = g0.28
parisize = 4000000, primelimit = 500000
? Qmat = [0.1121838950156262993198861903, -0.002294825170315502480266502542, -0.0003341361628152951755992730843, -0.01013912749586145381073809834, 0.02116023427426229437120984486, -0.04256500745397313843968688094; -0.002294825170315502480266502542, 0.0101
5200067110896037957565105, -0.01062142221979604038700374101, 0.02221761484223348391636273181, -0.04692817364343353260667839410, 0.09682920303695533048890485886; -0.0003341361628152951755992730843, -0.01062142221979604038700374101, 0.0232236177136219324940
5380743, -0.04916043776688522464771472577, 0.1039418432338874002700813232, -0.2172034509288058753060698920; -0.01013912749586145381073809834, 0.02221761484223348391636273181, -0.04916043776688522464771472577, 0.1063794149634527727718385934, -0.22753531414
34766490036291424, 0.4811386293014932004986704852; 0.02116023427426229437120984486, -0.04692817364343353260667839410, 0.1039418432338874002700813232, -0.2275353141434766490036291424, 0.4923171763996275063527134313, -1.053129742820247211799751080; -0.04256
500745397313843968688094, 0.09682920303695533048890485886, -0.2172034509288058753060698920, 0.4811386293014932004986704852, -1.053129742820247211799751080, 2.278649858405744971446804119];
? qfminim(Qmat*1000,1,0,2)
%2 = [474, 0.9981077444117395979158268274, [;]]
? qfminim(Qmat*1000,0.5,0,2)
%3 = [474, 0.9981077444117395979158268274, [;]]
? qfminim(Qmat*1000,0.00001,0,2)
%4 = [474, 0.9981077444117395979158268274, [;]]
?
The fix is to change line 2861 in basemath/bibli1.c (in function
fincke_pohst) from
res = smallvectors(gram,gceil(bound? bound: gcoeff(gram,i,i)),
to
res = smallvectors(gram,bound? bound: gcoeff(gram,i,i),
(I'm also sending a bigger patch fixing other things also)
Happy computing!
Paul van Wamelen