flint-sys 0.9.0

Bindings to the FLINT C library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61

/*
    Copyright (C) 2024 Albin Ahlbäck
    Copyright (C) 2024 Fredrik Johansson

    This file is part of FLINT.

    FLINT is free software: you can redistribute it and/or modify it under
    the terms of the GNU Lesser General Public License (LGPL) as published
    by the Free Software Foundation; either version 3 of the License, or
    (at your option) any later version.  See <https://www.gnu.org/licenses/>.
*/

#include "test_helpers.h"
#include "mpn_extras.h"

TEST_FUNCTION_START(flint_mpn_mulhigh_n_recursive, state)
{
    slong ix;
    int result;
    mp_ptr rp, rpc, xp, yp;

    for (ix = 0; ix < 10000 * flint_test_multiplier(); ix++)
    {
        mp_size_t n;
        mp_limb_t ret0, ret1;

        n = 1 + n_randint(state, 100);

        xp = flint_malloc(sizeof(mp_limb_t) * n);
        yp = flint_malloc(sizeof(mp_limb_t) * n);
        rp = flint_malloc(sizeof(mp_limb_t) * n);
        rpc = flint_malloc(sizeof(mp_limb_t) * n);

        flint_mpn_rrandom(xp, state, n);
        flint_mpn_rrandom(yp, state, n);

        ret0 = _flint_mpn_mulhigh_n_naive(rpc, xp, yp, n);
        ret1 = _flint_mpn_mulhigh_n_recursive(rp, xp, yp, n);

        result = (mpn_cmp(rp, rpc, n) == 0 && ret0 == ret1);
        if (!result)
            TEST_FUNCTION_FAIL(
                    "Not coherent with generic version.\n"
                    "ix = %wd\n"
                    "n = %wd\n"
                    "xp = %{ulong*}\n"
                    "yp = %{ulong*}\n"
                    "Expected ret: %{ulong}\n"
                    "Got ret:      %{ulong}\n"
                    "Expected limbs: %{ulong*}\n"
                    "Got limbs:      %{ulong*}\n",
                    ix, n, xp, n, yp, n, ret0, ret1, rpc, n, rp, n);

        flint_free(xp);
        flint_free(yp);
        flint_free(rp);
        flint_free(rpc);
    }

    TEST_FUNCTION_END(state);
}