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

/*
    Copyright (C) 2011, 2021 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 "nmod.h"
#include "nmod_vec.h"
#include "nmod_mat.h"

#if FLINT_USES_BLAS || (FLINT_BITS == 32) || !defined(__AVX2__)

/* Tuned for Zen 3 with BLAS */
static const slong lu_cutoff_tab[64] = { 64, 64, 244, 280, 296, 320, 332,
    792, 800, 800, 800, 400, 400, 400, 400, 404, 404, 400, 408, 416, 400,
    412, 424, 408, 484, 1352, 1032, 260, 68, 56, 56, 56, 148, 160, 148, 156,
    156, 156, 120, 168, 160, 124, 156, 124, 156, 148, 112, 156, 148, 136,
    156, 160, 116, 136, 168, 148, 156, 160, 120, 128, 68, 64, 104, 104 };

#else

/* Tuned for Zen 3 without BLAS */
static const slong lu_cutoff_tab[64] = { 64, 64, 212, 260, 280, 316, 344,
    792, 872, 904, 856, 1016, 1136, 1456, 1440, 1464, 1376, 1392, 1448, 1448,
    1360, 1392, 1400, 1392, 1448, 1416, 1032, 260, 68, 56, 56, 60, 168, 164,
    152, 152, 156, 148, 148, 148, 152, 164, 160, 148, 164, 148, 164, 160, 160,
    148, 156, 156, 148, 164, 148, 164, 148, 148, 148, 128, 68, 64, 96, 96 };

#endif

slong
nmod_mat_lu(slong * P, nmod_mat_t A, int rank_check)
{
    slong nrows = A->r, ncols = A->c, n;
    n = FLINT_MIN(nrows, ncols);

    if (n <= 3 || (NMOD_BITS(A->mod) > 28 && n <= 7))
    {
        return nmod_mat_lu_classical(P, A, rank_check);
    }
    else
    {
        if (n < lu_cutoff_tab[NMOD_BITS(A->mod) - 1])
            return nmod_mat_lu_classical_delayed(P, A, rank_check);
        else
            return nmod_mat_lu_recursive(P, A, rank_check);
    }
}