Bug numbers refer to the BTS at http://pari.math.u-bordeaux.fr/Bugs/
Done for version 2.13.3 (released 25/10/2021):
[last column crossreferences current development release 2.14.0]
  Fixed
     1- zeta(-2 + O(x)) => division by 0                                  [F25]
     2- (-1)^(-1/2) -> I instead of -I                                    [F56]
     3- GC error in ZX_Uspensky(, flag = 0)                               [F58]
BA   4- ellisomat(ell/nf not integral) -> oo-loop or wrong result         [F59]
BA   5- (f(~v)=v[1]=1);(g(v) = f(~v)); my(V=[0]);g(V);V -> [1] instead
         of [0] due to missing copy on write [#2304]                      [F60]
     6- hypergeom([1/12,1/12],2/3,3) -> wrong result                      [F61]
     7- overflow in mulu_interval(,ULONG_MAX) [#2309]                     [F62]
     8- ellE(1) -> domain error [#2312]                                   [F63]
     9- log1p(-1 + 10^-10) -> oo time                                     [F64]
    10- bernvec entries were wrong from B_{2*4064} on [#2311]             [F65]
    11- plotrecthraw for multi-curves                                     [F66]
    12- RgXn_powu(,,1) wrong result [#2314]                               [F67]
    13- erfc incorrect on negative imaginary axis [#2315]                 [F68]
    14- mfgalpoistype bug for dihedral types (when conductor at oo is not
        stable under Galois) [#2293]                                      [F69]
BA  15- [32bit] forfactored(n=223092868, 223092871, ) -> SEGV [#2318]     [F70]
    16- Fp_issquare(0, p) would return 0                                  [F71]

Done for version 2.13.2 (released 24/06/2021):
[last column crossreferences current development release 2.14.0]
  Fixed
     1- rnfidealup(, non-integral ideal) -> division by 0                 [F22]
     2- memory leak in ellheight                                          [F23]
     3- memory leak in bnfinit                                            [F24]
     4- nfeltdiv(nfinit(x^2+1), 0, 1+x/2) -> SEGV [#2277]                 [F26]
     5- znchargauss(g,chi,0) -> oo loop                                   [F28]
     6- bnflog(bnr, ...) -> SEGV                                          [F29]
     7- elliptic functions: incorrect reduction of z/w1 mod [1,tau] =>
        catastrophic cancellation                                         [F30]
     8- powcx(x, t_COMPLEX with 0 imaginary part) -> crash                [F31]
     9- nfsubfields(t_POL in other var than x,,1) -> SEGV                 [F32]
    10- extra terms in power t_SER substitutions [#2281]                  [F33]
    11- k=bnfinit(quadpoly(-2923,a)); bnrclassfield(k,3) -> SEGV [#2280]  [F34]
BA  12- rnfdisc mishandled polynomials with denominators [#2284]          [F35]
BA  13- elltamagawa(ellinit([-1456/243,93184/19683])) -> wrong result     [F36]
BA  14- mfsearch would sometimes miss the trivial form                    [F37]
    15- typo in allh [affects ECPP tunings]                               [F38]
HC  16- mfisetaquo(mffrometaquo([1,12;193,12])) -> 0                      [F39]
    17- polroots(1.E-47*t^6+1.E-47*t^5+4.E-46*t^4) -> bug in type [#2286] [F40]
    18- ceil(quadgen(5)) -> 1 instead of 0 [#2288]                        [F41]
    19- floor(Mod(0,3)*quadgen(5)) -> crash                               [F42]
    20- wrong ispower(t_RFRAC) and issquare(t_RFRAC, &z) [#2289]          [F45]
    21- 1+O(x)==O(x)                                                      [F47]
BA  22- qfauto([[1,0;0,1],[0,0;1,0]]) -> SEGV                             [F49]
    23- bnfsunit(...)[4] (the S-regulator) did not correspond to its
        definition in the documentation [ was using log(p) instead of
        log(Norm P)) for P in S above p ]                                 [F50]
BA  24- rnfdisc could return corrupted result                             [F51]
BA  25- concat(List([1])[1..0]) ->SEGV (instead of error) [#2299]         [F53]
BA  26- List([1])[1..0]==List() -> 0 instead of 1                         [F54]

Done for version 2.13.1 (released 25/01/2021):
[last column crossreferences current development release 2.14.0]
  Fixed
     1- GC error in idealinv [#2258]                                      [F2]
     2- qfminim inefficiency (initial bound from non-LLL-reduced matrix)  [F3]
     3- mfshimura in weight 3/2 => oo loop                                [F4]
     4- matsolve([1,1,1/2; 1,0,1/2; 2,1,0], [-1;1;-1]) -> SEGV [#2261]    [F6]
     5- qfminim(Mat(1),,,2) -> precision error [#2262]                    [F7]
     6- subst(p, v, t_SER) very slow if p contains many variables [#2264] [F8]
     7- mfsymboleval for trivial path returned 0 instead of a vector of 0s
        when f has multiple embeddings [#2260]                            [F9]
     8- 2^(1 + O(3^4)) -> type error                                      [F10]
     9- Zn_quad_roots(8,0,-1) to compute roots of x^2-1 mod 8 -> [4,[1,3]]
        (1 or 3 mod 4) instead of the expected [2,[1]] (1 mod 2) [#2265]  [F11]
    10- tan(1+10^20*I) -> overflow (same for cotan and psi)               [F12]
BA  11- Mod(2,3)^1000000000000000000001 -> Mod(-1,3) [#2266]              [F13]
    12- subst(O(y) + x, x, (1 - y + O(y^2))*x + O(x^2)) -> SEGV [#2267]   [F14]
    13- (Mod(0,3) + x) + O(x^2) -> x + O(x^2)  [now Mod(1,3)*x + O(x^2)]  [F15]
    14- precision too low for units in mfgaloisprojrep                    [F16]
    15- missing GC at the end of RgXn_inv                                 [F17]
    16- (-1)^(1/3) -> -1 instead of exp(log(-1)/3)                        [F18]
    17- mfeisenstein(k, trivial characters mod N > 1) was incorrect       [F19]
BA  18- missing GC in qfsolve                                             [F20]
