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
67
68
69
70
71
72
73
74
75
76
77
78

/*
    Copyright (C) 2024, 2026 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 <string.h>
#include <stdint.h>
#include "mpn_extras.h"
#include "radix.h"

FLINT_TLS_PREFIX radix_t default_decimal_ctx;
FLINT_TLS_PREFIX int default_decimal_ctx_initialized = 0;

static void
default_decimal_ctx_cleanup(void)
{
    if (default_decimal_ctx_initialized)
    {
        default_decimal_ctx_initialized = 0;
        radix_clear(default_decimal_ctx);
    }
}

radix_struct * get_default_decimal_ctx(void);

radix_struct * get_default_decimal_ctx(void)
{
    if (!default_decimal_ctx_initialized)
    {
        radix_init(default_decimal_ctx, 10, 0);
        flint_register_cleanup_function(default_decimal_ctx_cleanup);
        default_decimal_ctx_initialized = 1;
    }

    return default_decimal_ctx;
}



char * flint_mpn_get_str(char * res, int base, mp_srcptr x, mp_size_t xn, int negative)
{
    radix_struct * radix;
    nn_ptr t;
    slong alloc;
    slong decimal_limbs;
    TMP_INIT;

#if FLINT_HAVE_FFT_SMALL
    /* Todo: improve the radix code so that it beats GMP in this range too. */
    if ((xn > 40 && xn < 1400) || base != 10)
#endif
    {
        mpz_t tmp;
        while (xn > 0 && x[xn - 1] == 0)
            xn--;
        tmp->_mp_d = (mp_ptr) x;
        tmp->_mp_alloc = tmp->_mp_size = negative ? -xn : xn;
        if (res == NULL)
            res = flint_malloc(mpz_sizeinbase(tmp, base) + 2);
        return mpz_get_str(res, base, tmp);
    }

    TMP_START;
    radix = get_default_decimal_ctx();
    alloc = radix_set_mpn_need_alloc(xn, radix);
    t = TMP_ALLOC(alloc * sizeof(ulong));
    decimal_limbs = radix_set_mpn(t, x, xn, radix);
    res = radix_get_str_decimal(res, t, decimal_limbs, negative, radix);
    TMP_END;
    return res;
}