Bug numbers refer to the BTS at http://pari.math.u-bordeaux.fr/Bugs/

 [last column crossreferences current development release 2.18.1]
Done for version 2.17.4 (released 29/06/2026):
  Fixed
BA  1- call(printsep, [" "]) -> SEGV                                      [F44]
    2- forsubset: handle trivial case of negative cardinalities [#2652]   [F47]
BA  3- ellchangecurve(E).roots could be wrong [#2660]                     [F48]
    4- charpoly(Mod(x+y,y^2+x+1),y) -> invalid t_POL [#2668]              [F50]
    5- writebin(file) only saved polynomials variables in file            [F53]
BA  6- lfun(1/3+I*10^19) -> overflow in expo()                            [F54]
BA  7- fflog(a, b, o) was very slow for small o but large fields [#2676]  [F55]
LGr 8- psi(1+O(x^5)) -> wrong result [#2677]                              [F56]
    9- zetahurwitz(large s, |x| > 1) is inaccurate (correct absolute accuracy
        but zeta(s,x) is small and relative error is large) [#2673]       [F57]
BA  10- ellcard(ellinit([ffgen([53,36]),1])) -> wrong result              [F58]
    11- polcyclofactors(((x^6+x^3+1)*(3*x^3+1))^2) -> oo loop [#2680]     [F59]


 [last column crossreferences current development release 2.18.1]
Done for version 2.17.3 (released 24/11/2025):
  Fixed
    1- sumnumrat & friends: use full precision estimates to avoid failing with
        domain error when "close" to boundary                             [F17]
    2- zeta([3,2,1]) -> oo-loop                                           [F18]
BA  3- parforeach([1..1000],j,vector(10^5,i,i)) -> stack overflow         [F19]
    4- mfpetersson() incorrect when one form not cuspidal                 [F20]
    5- lift_to_frac: use of cgiv could corrupt PARI stack                 [F22]
    6- sumeulerrat / prodeulerrat fail in cases where they shouldn't [#2620]
       e.g., sumeulerrat(1/(p^2-5),1.1), prodeulerrat(1+1/(p^2-p-23),0.6) [F23]
    7- incgam(0,10^-50) -> GEN corruption                                 [F25]
BA  8- 12*harmonic(1) -> 12/1                                             [F26]
    9- hyperellratpoints(P, d) could miss points when P(-d) = 0           [F27]
   10- obscure bnrclassfield bug [internally creating an incorrect ideal
       factorization whose prime ideals were not sorted]                  [F28]
BA 11- qfeval(Qfb(2,3,4),[1,2;3,4]).disc was wrong [#2630]                [F29]
   12- bnfsunit: could fail to detect non-S-units                         [F30]
BA 13- lfunan: overflow error with large bad primes.                      [F32]
   14- mftonew: fix problem with mfbd forms [#2619]                       [F34]
   15- lfunthetainit: lost accuracy for large shifts.                     [F35]
HC 16- hypergeom([8,8.5],7,-0.4) was wrong [#2628]                        [F36]
HC 17- lfunmfspec(L) divided by zero if period vanishes                   [F37]
BA 18- fforder(,factor(1)) -> crash [#2618]                               [F38]
   19- mateigen must not call qfjacobi on non-real matrices [#2636]       [F39]
BA 20- ^C\n^C in parfor could crash GP                                    [F40]
   21- solve() iteration could incorrectly abort and return 0.0 [#2638]   [F41]
   22- zeta(1/2+10^22*I) -> precision too low                             [F42]
BA 23- lfuninit(lfunshift(1,3,1)) was incorrect [#2631]                   [F43]

 [last column crossreferences current development release 2.18.1]
Done for version 2.17.2 (released 01/03/2025):
  Fixed
BA   1- ellrootno over number fields could be slow [#2595]                 [F2]
BA   2- matsolvemod could fail [#2604]                                     [F3]
BA   3- nfcompositum(nfinit(y-2),x^2-y,x^3-x-1) -> error [#2605]           [F4]
     4- Mod(a,A) == Mod(b,B) was checking a == b && A == B instead         [F6]
        of the documented Mod(a,A) - Mod(b,B) == 0 [F6]
     5- normlp(, oo): support components equal to oo or -oo                [F8]
     6- sqrtn(27+O(3^4),3) -> 1 + O(3) (wrong valuation in Qp_sqrtn_ram)  [F10]
     7- vectorsmall(5, i, 2^63-1).normfu -> SEGV                          [F11]
AP   8- alginit(nf,[a,b]) would reject a in basis form                    [F13]
BA   9- for(i=1,10^6,bernpol(5,1)) -> "no more variables" [#2611]         [F14]
BA  10- Flx_mod_Xn1(T,n), deg T >= 2*n -> wrong result                    [F15]
    11- forprime(p=1099505335069-100000,1099505335069+10000,) -> crash    [F16]

 [last column crossreferences current development release 2.18.0]
Done for version 2.17.1 (released 24/12/2024):
  Fixed
BA   1- incorrect GP proto for ellpadicfrobenius, mftocoset, qfcvp [#2575] [F1]
AP   2- fix idealismaximal for ideal p*ZK [Sage issue 38766]               [F2]
AP   3- segfault in algsubalg when input is linearly dependent             [F3]
AP   4- fix alglatinter trivial intersection (misuse of matkermod(,1))     [F4]
BA   5- ellrank could try many more curves than requested by effort [#2577][F5]
     6- wrong result in issquare(Mod(a, composite)) [#2578]                [F6]
     7- gener_FpXQ_local(t_POL of degre 1) -> oo loop [#2581]              [F7]
     8- GC error when squaring integer 2x2 matrices [#2582]                [F8]
AP   9- algebras: fix splitting field for CSA, algtomatrix / algcharpoly
        were sometimes wrong by an automorphism of the center              [F9]
BA  10- paridecl.h could conflict with complex.h [#2583]                  [F10]
    11- forprime(p=524288,1048576,1) -> crash [#2584]
        forprime(p=1048607,1048617,1) -> oo loop                          [F11]
    12- printf("%10.6f", tiny rational) -> undefined result [#2576]       [F12]
AP  13- alginit(nf,[a,b]) was only allowing a,b integral over the power basis
        instead of arbitrary algebraic integers.                          [F13]
HC  14- zeta not continous around 1/2+3000*I [#2585]                      [F14]
    15- missing type checks for change of variables in ellchangecurve,    [F15]
        ellchangepoint, ellchangepointinv
BA  16- control-C could lead to segfault later [#2496]                    [F16]
    17- [libpari] missing normalization in FpX_mulu (when p is not prime) [F18]
HC  18- mfsymboleval could return wrong sign [#2586]                      [F19]
    19- missing type checks in factorback [#2590]                         [F20]
BA  20- charpoly(Mod(x,x^20+1/2)) -> crash [#2591]                        [F21]
    21- psi(x close to a nonpositive integer) -> error                    [F22]
BA  22- polresultant(t_POL with RFRAC coeffs) could be incorrect [#2588]  [F24]
BA  23- qflll(non-maximal rank) -> potential oo-loop [#2587]              [F25]

  Added
BA   1- [loongarch64] add loongarch64 assembly kernel

Done for version 2.17.0 (released 30/09/2024):
  Fixed
     1- rare SEGV in sumnumrat [#2566]
     2- gp -p N did not check for N < 2^32
     3- nflist([6,8],12167) -> SEGV [#2567]
     4- nflist([6,1],1229312), nflist([6,10],36192256) -> duplicate fields
AP   5- variable leak in algsplit
     6- mfslashexpansion: possible extra rational factor [#2572]
     7- nfeltsign: possible wrong result

  Added
BA   1- [libpari] Flm_sqr

  Changed
BA   1- qflll(,1), qfllgram(,1): disable FLATTER algorithm.
     2- tunings of forprime over "short" intervals (use sieve more frequently)
