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