Code coverage tests

This page documents the degree to which the PARI/GP source code is tested by our public test suite, distributed with the source distribution in directory src/test/. This is measured by the gcov utility; we then process gcov output using the lcov frond-end.

We test a few variants depending on Configure flags on the pari.math.u-bordeaux.fr machine (x86_64 architecture), and agregate them in the final report:

The target is to exceed 90% coverage for all mathematical modules (given that branches depending on DEBUGLEVEL or DEBUGMEM are not covered). This script is run to produce the results below.

LCOV - code coverage report
Current view: top level - basemath - buch4.c (source / functions) Hit Total Coverage
Test: PARI/GP v2.16.2 lcov report (development 29115-f22e516b23) Lines: 140 146 95.9 %
Date: 2024-03-18 08:03:28 Functions: 8 8 100.0 %
Legend: Lines: hit not hit

          Line data    Source code
       1             : /* Copyright (C) 2000  The PARI group.
       2             : 
       3             : This file is part of the PARI/GP package.
       4             : 
       5             : PARI/GP is free software; you can redistribute it and/or modify it under the
       6             : terms of the GNU General Public License as published by the Free Software
       7             : Foundation; either version 2 of the License, or (at your option) any later
       8             : version. It is distributed in the hope that it will be useful, but WITHOUT
       9             : ANY WARRANTY WHATSOEVER.
      10             : 
      11             : Check the License for details. You should have received a copy of it, along
      12             : with the package; see the file 'COPYING'. If not, write to the Free Software
      13             : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
      14             : 
      15             : /*******************************************************************/
      16             : /*                        RNFISNORM                                */
      17             : /*     (Adapted from Denis Simon's original implementation)        */
      18             : /*******************************************************************/
      19             : #include "pari.h"
      20             : #include "paripriv.h"
      21             : 
      22             : static void
      23         609 : p_append(GEN p, hashtable *H, hashtable *H2)
      24             : {
      25         609 :   ulong h = H->hash(p);
      26         609 :   hashentry *e = hash_search2(H, (void*)p, h);
      27         609 :   if (!e)
      28             :   {
      29         539 :     hash_insert2(H, (void*)p, NULL, h);
      30         539 :     if (H2) hash_insert2(H2, (void*)p, NULL, h);
      31             :   }
      32         609 : }
      33             : 
      34             : /* N a t_INT */
      35             : static void
      36         196 : Zfa_append(GEN N, hashtable *H, hashtable *H2)
      37             : {
      38         196 :   if (!is_pm1(N))
      39             :   {
      40         126 :     GEN v = gel(absZ_factor(N),1);
      41         126 :     long i, l = lg(v);
      42         308 :     for (i=1; i<l; i++) p_append(gel(v,i), H, H2);
      43             :   }
      44         196 : }
      45             : /* N a t_INT or t_FRAC or ideal in HNF*/
      46             : static void
      47         140 : fa_append(GEN N, hashtable *H, hashtable *H2)
      48             : {
      49         140 :   switch(typ(N))
      50             :   {
      51         112 :     case t_INT:
      52         112 :       Zfa_append(N,H,H2);
      53         112 :       break;
      54           0 :     case t_FRAC:
      55           0 :       Zfa_append(gel(N,1),H,H2);
      56           0 :       Zfa_append(gel(N,2),H,H2);
      57           0 :       break;
      58          28 :     default: /*t_MAT*/
      59          28 :       Zfa_append(gcoeff(N,1,1),H,H2);
      60          28 :       break;
      61             :   }
      62         140 : }
      63             : 
      64             : /* apply lift(rnfeltup) to all coeffs, without rnf structure */
      65             : static GEN
      66           7 : nfX_eltup(GEN nf, GEN rnfeq, GEN x)
      67             : {
      68             :   long i, l;
      69           7 :   GEN y = cgetg_copy(x, &l), zknf = nf_nfzk(nf, rnfeq);
      70          35 :   for (i=2; i<l; i++) gel(y,i) = nfeltup(nf, gel(x,i), zknf);
      71           7 :   y[1] = x[1]; return y;
      72             : }
      73             : 
      74             : static hashtable *
      75         196 : hash_create_INT(ulong s)
      76         196 : { return hash_create(s, (ulong(*)(void*))&hash_GEN,
      77             :                         (int(*)(void*,void*))&equalii, 1); }
      78             : GEN
      79          56 : rnfisnorminit(GEN T, GEN R, int galois)
      80             : {
      81          56 :   pari_sp av = avma;
      82             :   long i, l, dR;
      83             :   GEN S, gen, cyc, bnf, nf, nfabs, rnfeq, bnfabs, k, polabs;
      84          56 :   GEN y = cgetg(9, t_VEC);
      85          56 :   hashtable *H = hash_create_INT(100UL);
      86             : 
      87          56 :   if (galois < 0 || galois > 2) pari_err_FLAG("rnfisnorminit");
      88          56 :   T = get_bnfpol(T, &bnf, &nf);
      89          56 :   if (!bnf) bnf = Buchall(nf? nf: T, nf_FORCE, DEFAULTPREC);
      90          56 :   if (!nf) nf = bnf_get_nf(bnf);
      91             : 
      92          56 :   R = get_bnfpol(R, &bnfabs, &nfabs);
      93          56 :   if (!gequal1(leading_coeff(R))) pari_err_IMPL("non monic relative equation");
      94          56 :   dR = degpol(R);
      95          56 :   if (dR <= 2) galois = 1;
      96             : 
      97          56 :   R = RgX_nffix("rnfisnorminit", T, R, 1);
      98          56 :   if (nf_get_degree(nf) == 1) /* over Q */
      99          35 :     rnfeq = mkvec5(R,gen_0,gen_0,T,R);
     100          21 :   else if (galois == 2) /* needs eltup+abstorel */
     101           7 :     rnfeq = nf_rnfeq(nf, R);
     102             :   else /* needs abstorel */
     103          14 :     rnfeq = nf_rnfeqsimple(nf, R);
     104          56 :   polabs = gel(rnfeq,1);
     105          56 :   k = gel(rnfeq,3);
     106          56 :   if (!bnfabs || !gequal0(k))
     107          28 :     bnfabs = Buchall(polabs, nf_FORCE, nf_get_prec(nf));
     108          56 :   if (!nfabs) nfabs = bnf_get_nf(bnfabs);
     109             : 
     110          56 :   if (galois == 2)
     111             :   {
     112          21 :     GEN P = polabs==R? leafcopy(R): nfX_eltup(nf, rnfeq, R);
     113          21 :     setvarn(P, fetch_var_higher());
     114          21 :     galois = !!nfroots_if_split(&nfabs, P);
     115          21 :     (void)delete_var();
     116             :   }
     117             : 
     118          56 :   cyc = bnf_get_cyc(bnfabs);
     119          56 :   gen = bnf_get_gen(bnfabs); l = lg(cyc);
     120          84 :   for(i=1; i<l; i++)
     121             :   {
     122          35 :     GEN g = gel(gen,i);
     123          35 :     if (ugcdiu(gel(cyc,i), dR) == 1) break;
     124          28 :     Zfa_append(gcoeff(g,1,1), H, NULL);
     125             :   }
     126          56 :   if (!galois)
     127             :   {
     128          21 :     GEN Ndiscrel = diviiexact(nf_get_disc(nfabs), powiu(nf_get_disc(nf), dR));
     129          21 :     Zfa_append(Ndiscrel, H, NULL);
     130             :   }
     131          56 :   S = hash_keys_GEN(H);
     132          56 :   gel(y,1) = bnf;
     133          56 :   gel(y,2) = bnfabs;
     134          56 :   gel(y,3) = R;
     135          56 :   gel(y,4) = rnfeq;
     136          56 :   gel(y,5) = S;
     137          56 :   gel(y,6) = nf_pV_to_prV(nf, S);
     138          56 :   gel(y,7) = nf_pV_to_prV(nfabs, S);
     139          56 :   gel(y,8) = stoi(galois); return gerepilecopy(av, y);
     140             : }
     141             : 
     142             : /* T as output by rnfisnorminit
     143             :  * if flag=0 assume extension is Galois (==> answer is unconditional)
     144             :  * if flag>0 add to S all primes dividing p <= flag
     145             :  * if flag<0 add to S all primes dividing abs(flag)
     146             : 
     147             :  * answer is a vector v = [a,b] such that
     148             :  * x = N(a)*b and x is a norm iff b = 1  [assuming S large enough] */
     149             : GEN
     150          70 : rnfisnorm(GEN T, GEN x, long flag)
     151             : {
     152          70 :   pari_sp av = avma;
     153             :   GEN bnf, rel, R, rnfeq, nfpol;
     154             :   GEN nf, aux, H, U, Y, M, A, bnfS, sunitrel, futu, S, S1, S2, Sx;
     155             :   long L, i, itu;
     156             :   hashtable *H0, *H2;
     157          70 :   if (typ(T) != t_VEC || lg(T) != 9)
     158           0 :     pari_err_TYPE("rnfisnorm [please apply rnfisnorminit()]", T);
     159          70 :   bnf = gel(T,1);
     160          70 :   rel = gel(T,2);
     161          70 :   bnf = checkbnf(bnf);
     162          70 :   rel = checkbnf(rel);
     163          70 :   nf = bnf_get_nf(bnf);
     164          70 :   x = nf_to_scalar_or_alg(nf,x);
     165          70 :   if (gequal0(x)) { set_avma(av); return mkvec2(gen_0, gen_1); }
     166          70 :   if (gequal1(x)) { set_avma(av); return mkvec2(gen_1, gen_1); }
     167          70 :   R = gel(T,3);
     168          70 :   rnfeq = gel(T,4);
     169          70 :   if (gequalm1(x) && odd(degpol(R)))
     170           0 :   { set_avma(av); return mkvec2(gen_m1, gen_1); }
     171             : 
     172             :   /* build set T of ideals involved in the solutions */
     173          70 :   nfpol = nf_get_pol(nf);
     174          70 :   S = gel(T,5);
     175          70 :   H0 = hash_create_INT(100UL);
     176          70 :   H2 = hash_create_INT(100UL);
     177          70 :   L = lg(S);
     178         147 :   for (i = 1; i < L; i++) p_append(gel(S,i),H0,NULL);
     179          70 :   S1 = gel(T,6);
     180          70 :   S2 = gel(T,7);
     181          70 :   if (flag > 0)
     182             :   {
     183             :     forprime_t T;
     184             :     ulong p;
     185          14 :     u_forprime_init(&T, 2, flag);
     186         364 :     while ((p = u_forprime_next(&T))) p_append(utoipos(p), H0,H2);
     187             :   }
     188          56 :   else if (flag < 0)
     189           7 :     Zfa_append(utoipos(-flag),H0,H2);
     190             :   /* overkill: prime ideals dividing x would be enough */
     191          70 :   A = idealnumden(nf, x);
     192          70 :   fa_append(gel(A,1), H0,H2);
     193          70 :   fa_append(gel(A,2), H0,H2);
     194          70 :   Sx = hash_keys_GEN(H2); L = lg(Sx);
     195          70 :   if (L > 1)
     196             :   { /* new primes */
     197          49 :     S1 = shallowconcat(S1, nf_pV_to_prV(nf, Sx));
     198          49 :     S2 = shallowconcat(S2, nf_pV_to_prV(rel, Sx));
     199             :   }
     200             : 
     201             :   /* computation on T-units */
     202          70 :   futu = shallowconcat(bnf_get_fu(rel), bnf_get_tuU(rel));
     203          70 :   bnfS = bnfsunit(bnf,S1,LOWDEFAULTPREC);
     204          70 :   sunitrel = shallowconcat(futu, gel(bnfsunit(rel,S2,LOWDEFAULTPREC), 1));
     205             : 
     206          70 :   A = lift_shallow(bnfissunit(bnf,bnfS,x));
     207          70 :   L = lg(sunitrel);
     208          70 :   itu = lg(nf_get_roots(nf))-1; /* index of torsion unit in bnfsunit(nf) output */
     209          70 :   M = cgetg(L+1,t_MAT);
     210        1449 :   for (i=1; i<L; i++)
     211             :   {
     212        1379 :     GEN u = eltabstorel(rnfeq, gel(sunitrel,i));
     213        1379 :     gel(sunitrel,i) = u;
     214        1379 :     u = bnfissunit(bnf,bnfS, gnorm(u));
     215        1379 :     if (lg(u) == 1) pari_err_BUG("rnfisnorm");
     216        1379 :     gel(u,itu) = lift_shallow(gel(u,itu)); /* lift root of 1 part */
     217        1379 :     gel(M,i) = u;
     218             :   }
     219          70 :   aux = zerocol(lg(A)-1); gel(aux,itu) = utoipos( bnf_get_tuN(rel) );
     220          70 :   gel(M,L) = aux;
     221          70 :   H = ZM_hnfall(M, &U, 2);
     222          70 :   Y = RgM_RgC_mul(U, inverseimage(H,A));
     223             :   /* Y: sols of MY = A over Q */
     224          70 :   setlg(Y, L);
     225          70 :   aux = factorback2(sunitrel, gfloor(Y));
     226          70 :   x = mkpolmod(x,nfpol);
     227          70 :   if (!gequal1(aux)) x = gdiv(x, gnorm(aux));
     228          70 :   x = lift_if_rational(x);
     229          70 :   if (typ(aux) == t_POLMOD && degpol(nfpol) == 1)
     230          28 :     gel(aux,2) = lift_if_rational(gel(aux,2));
     231          70 :   return gerepilecopy(av, mkvec2(aux, x));
     232             : }
     233             : 
     234             : GEN
     235          28 : bnfisnorm(GEN bnf, GEN x, long flag)
     236             : {
     237          28 :   pari_sp av = avma;
     238          28 :   GEN T = rnfisnorminit(pol_x(fetch_var()), bnf, flag == 0? 1: 2);
     239          28 :   GEN r = rnfisnorm(T, x, flag == 1? 0: flag);
     240          28 :   (void)delete_var();
     241          28 :   return gerepileupto(av,r);
     242             : }

Generated by: LCOV version 1.14