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:
- with GMP kernel
- with GMP kernel using --mt=pthread
- with native kernel, including micro-assembler code
- with native kernel, without micro-assembler
- with native kernel, without micro-assembler, disabling GCC extensions
(DISABLE_INLINE)
- with GMP kernel, emulating an
x86_32 architecture at
Configure time via setarch
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 |
 |
|
|
|
|
Line data Source code
1 : #line 2 "../src/kernel/gmp/gcd.c"
2 : /* Copyright (C) 2000-2003 The PARI group.
3 :
4 : This file is part of the PARI/GP package.
5 :
6 : PARI/GP is free software; you can redistribute it and/or modify it under the
7 : terms of the GNU General Public License as published by the Free Software
8 : Foundation; either version 2 of the License, or (at your option) any later
9 : version. It is distributed in the hope that it will be useful, but WITHOUT
10 : ANY WARRANTY WHATSOEVER.
11 :
12 : Check the License for details. You should have received a copy of it, along
13 : with the package; see the file 'COPYING'. If not, write to the Free Software
14 : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
15 :
16 : /* assume y > x > 0. return y mod x */
17 : static ulong
18 36682646 : resiu(GEN y, ulong x)
19 : {
20 36682646 : return mpn_mod_1(LIMBS(y), NLIMBS(y), x);
21 : }
22 :
23 : GEN
24 830781619 : gcdii(GEN a, GEN b)
25 : {
26 : long v, w;
27 : pari_sp av;
28 : GEN t;
29 :
30 830781619 : switch (abscmpii(a,b))
31 : {
32 237252117 : case 0: return absi(a);
33 294511201 : case -1: swap(a,b);
34 : }
35 593477807 : if (!signe(b)) return absi(a);
36 : /* here |a|>|b|>0. Try single precision first */
37 336931613 : if (lgefint(a)==3)
38 277807906 : return igcduu((ulong)a[2], (ulong)b[2]);
39 59123707 : if (lgefint(b)==3)
40 : {
41 36682699 : ulong u = resiu(a,(ulong)b[2]);
42 36683210 : if (!u) return absi(b);
43 22101561 : return igcduu((ulong)b[2], u);
44 : }
45 : /* larger than gcd: "set_avma(av)" gerepile (erasing t) is valid */
46 22441008 : av = avma; (void)new_chunk(lgefint(b)+1); /* HACK */
47 22463685 : t = remii(a,b);
48 22463965 : if (!signe(t)) { set_avma(av); return absi(b); }
49 :
50 15497432 : a = b; b = t;
51 15497432 : v = vali(a); a = shifti(a,-v); setabssign(a);
52 15496996 : w = vali(b); b = shifti(b,-w); setabssign(b);
53 15496924 : if (w < v) v = w;
54 15496924 : switch(abscmpii(a,b))
55 : {
56 311909 : case 0: set_avma(av); a=shifti(a,v); return a;
57 3293189 : case -1: swap(a,b);
58 : }
59 15184966 : if (is_pm1(b)) { set_avma(av); return int2n(v); }
60 : {
61 : /* general case */
62 : /*This serve two purposes: 1) mpn_gcd destroy its input and need an extra
63 : * limb 2) this allows us to use icopy instead of gerepile later. NOTE: we
64 : * must put u before d else the final icopy could fail.
65 : */
66 15050537 : GEN res= cgeti(lgefint(a)+1);
67 15050334 : GEN ca = icopy_ef(a,lgefint(a)+1);
68 15049741 : GEN cb = icopy_ef(b,lgefint(b)+1);
69 15050065 : long l = mpn_gcd(LIMBS(res), LIMBS(ca), NLIMBS(ca), LIMBS(cb), NLIMBS(cb));
70 15051493 : res[1] = evalsigne(1)|evallgefint(l+2);
71 15051493 : set_avma(av);
72 15051398 : return shifti(res,v);
73 : }
74 : }
|