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
62
63
64
65
66

/*
    Copyright (C) 2024 Albin Ahlbäck

    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_mullow_n, state)
{
    slong ix;
    int result;

    mp_ptr rp, rpf, xp, yp;

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

        /* Trigger full multiplication in mulhigh */
        if (n_randint(state, 1000) == 0)
            n = 1 + FLINT_MPN_MULHIGH_MUL_CUTOFF +  n_randint(state, 50);
        else if (n_randint(state, 100) == 0)
            n = 1 + n_randint(state, FLINT_MPN_MULHIGH_MUL_CUTOFF);
        else
            n = 1 + n_randint(state, 2 * FLINT_MPN_MULHIGH_MULDERS_CUTOFF);

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

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

        ret = flint_mpn_mullow_n(rp, xp, yp, n);
        flint_mpn_mul_n(rpf, xp, yp, n);

        result = (mpn_cmp(rp, rpf, n) == 0 && ret == rpf[n]);
        if (!result)
            TEST_FUNCTION_FAIL(
                    "ix = %wd\n"
                    "n = %wd\n"
                    "xp = %{ulong*}\n"
                    "yp = %{ulong*}\n"
                    "Exp ret: %{ulong}\n"
                    "Got ret: %{ulong}\n"
                    "Expected: %{ulong*}\n"
                    "Got:      %{ulong*}\n",
                    ix, n, xp, n, yp, n, rpf[n], ret, rpf, n, rp, n);

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

    TEST_FUNCTION_END(state);
}