flint3-sys 3.5.1

Rust 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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96

/*
    Copyright (C) 2023 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"
#include "fft_small.h"
#include "crt_helpers.h"

#define TEST_ADD(n) \
        flint_mpn_copyi(a, d, n); \
        mpn_add_n(c, a, b, n); \
        multi_add_ ## n(a, b); \
        if (mpn_cmp(a, c, n) != 0) \
        { \
            flint_printf("FAIL: add_%i\n", n); \
            flint_abort(); \
        }

#define TEST_SUB(n) \
        flint_mpn_copyi(a, d, n); \
        mpn_sub_n(c, a, b, n); \
        multi_sub_ ## n(a, b); \
        if (mpn_cmp(a, c, n) != 0) \
        { \
            flint_printf("FAIL: sub_%i\n", n); \
            flint_abort(); \
        }

TEST_FUNCTION_START(flint_mpn_add_inplace_c, state)
{
    slong iter;

    for (iter = 0; iter < 1000 * flint_test_multiplier(); iter++)
    {
        ulong a[10], b[10], c[10];
        slong an, bn;
        unsigned char cf, c1, c2;

        bn = 1 + n_randint(state, 4);
        an = bn + n_randint(state, 4);

        flint_mpn_rrandom(a, state, an);
        flint_mpn_rrandom(b, state, bn);
        flint_mpn_copyi(c, a, an);
        cf = n_randint(state, 2);

        c1 = flint_mpn_add_inplace_c(a, an, b, bn, cf);

        c2 = mpn_add(c, c, an, b, bn);
        c2 += mpn_add_1(c, c, an, cf);

        if (c1 != c2 || mpn_cmp(a, c, an) != 0)
        {
            flint_printf("FAIL\n");
            flint_abort();
        }
    }

    for (iter = 0; iter < 1000 * flint_test_multiplier(); iter++)
    {
        ulong a[8], b[8], c[8], d[8];

        flint_mpn_rrandom(a, state, 8);
        flint_mpn_rrandom(b, state, 8);
        flint_mpn_copyi(c, a, 8);
        flint_mpn_copyi(d, a, 8);

        TEST_ADD(1)
        TEST_ADD(2)
        TEST_ADD(3)
        TEST_ADD(4)
        TEST_ADD(5)
        TEST_ADD(6)
        TEST_ADD(7)
        TEST_ADD(8)

        TEST_SUB(1)
        TEST_SUB(2)
        TEST_SUB(3)
        TEST_SUB(4)
        TEST_SUB(5)
        TEST_SUB(6)
        TEST_SUB(7)
        TEST_SUB(8)
    }

    TEST_FUNCTION_END(state);
}