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 36334723 : resiu(GEN y, ulong x)
19 : {
20 36334723 : return mpn_mod_1(LIMBS(y), NLIMBS(y), x);
21 : }
22 :
23 : GEN
24 851985214 : gcdii(GEN a, GEN b)
25 : {
26 : long v, w;
27 : pari_sp av;
28 : GEN t;
29 :
30 851985214 : switch (abscmpii(a,b))
31 : {
32 235976686 : case 0: return absi(a);
33 300449321 : case -1: swap(a,b);
34 : }
35 615966823 : if (!signe(b)) return absi(a);
36 : /* here |a|>|b|>0. Try single precision first */
37 364353730 : if (lgefint(a)==3)
38 305427101 : return igcduu((ulong)a[2], (ulong)b[2]);
39 58926629 : if (lgefint(b)==3)
40 : {
41 36334793 : ulong u = resiu(a,(ulong)b[2]);
42 36334956 : if (!u) return absi(b);
43 22010713 : return igcduu((ulong)b[2], u);
44 : }
45 : /* larger than gcd: GC using set_avma(av) (erasing t) is valid */
46 22591836 : av = avma; (void)new_chunk(lgefint(b)+1); /* HACK */
47 22614212 : t = remii(a,b);
48 22614724 : if (!signe(t)) { set_avma(av); return absi(b); }
49 :
50 15605409 : a = b; b = t;
51 15605409 : v = vali(a); a = shifti(a,-v); setabssign(a);
52 15604866 : w = vali(b); b = shifti(b,-w); setabssign(b);
53 15604696 : if (w < v) v = w;
54 15604696 : switch(abscmpii(a,b))
55 : {
56 300655 : case 0: set_avma(av); a=shifti(a,v); return a;
57 3300466 : case -1: swap(a,b);
58 : }
59 15303972 : if (is_pm1(b)) { set_avma(av); return int2n(v); }
60 : {
61 : /* Serves two purposes: 1) mpn_gcd destroys its input and needs an extra limb
62 : * 2) allows us to use icopy for GC. */
63 15179698 : GEN res= cgeti(lgefint(a)+1);
64 15179473 : GEN ca = icopy_ef(a,lgefint(a)+1);
65 15178446 : GEN cb = icopy_ef(b,lgefint(b)+1);
66 15179275 : long l = mpn_gcd(LIMBS(res), LIMBS(ca), NLIMBS(ca), LIMBS(cb), NLIMBS(cb));
67 15180736 : res[1] = evalsigne(1)|evallgefint(l+2);
68 15180736 : set_avma(av); return shifti(res,v);
69 : }
70 : }
|