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