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
|
#include "fplll-interface.hpp"
#include "mutablemat.hpp"
#ifdef HAVE_FPLLL
#include <stddef.h>
#include <fplll.h>
#endif
bool fp_LLL(MutableMatrix *M, MutableMatrix *U, int strategy)
{
#ifndef HAVE_FPLLL
ERROR("fplll is not available (configure M2 with fplll!)");
return 0;
#else
assert(U == NULL);
double delta = 0.99;
double eta = 0.51;
LLLMethod method = LM_WRAPPER;
FloatType floatType = FT_DEFAULT;
int precision = 0;
int flags = LLL_DEFAULT;
int ncols = static_cast<int>(M->n_rows());
int nrows = static_cast<int>(M->n_cols());
ZZ_mat<mpz_t> mat(nrows, ncols);
for (int i = 0; i < nrows; i++)
for (int j = 0; j < ncols; j++)
{
ring_elem a;
if (M->get_entry(j, i, a))
{
mpz_set(mat[i][j].getData(), a.get_mpz());
}
}
int result =
lllReduction(mat, delta, eta, method, floatType, precision, flags);
switch (result)
{
case RED_SUCCESS:
break;
case RED_BABAI_FAILURE:
ERROR("Error in fpLLL");
return 0;
case RED_LLL_FAILURE:
ERROR("infinite loop in LLL");
return 0;
default:
ERROR("unknown error in fpLLL");
return 0;
}
/* Put this back into M */
mpz_t a;
mpz_init(a);
for (int j = 0; j < ncols; j++)
for (int i = 0; i < nrows; i++)
{
mpz_set(a, mat[i][j].getData());
ring_elem b = globalZZ->from_int(a);
M->set_entry(j, i, b);
}
mpz_clear(a);
return true;
#endif
}
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e "
// indent-tabs-mode: nil
// End:
|
