Macaulay2 Engine
Loading...
Searching...
No Matches

◆ ntl_LLL()

bool ntl_LLL ( MutableMatrix * M,
MutableMatrix * U,
long numer,
long denom,
int strategy )

Definition at line 69 of file ntl-interface.cpp.

74{
75 int nrows = static_cast<int>(M->n_rows());
76 int ncols = static_cast<int>(M->n_cols());
77
78 NTL::ZZ d;
79 // Note that the LLL routines all return the rank, but we ignore this return
80 // value.
81 double delta = static_cast<double>(numer) / static_cast<double>(denom);
82
83 if (M2_gbTrace >= 10) printf("LLL: using strategy %d\n", strategy);
84 NTL::mat_ZZ *A = mutableMatrix_to_NTL_mat_ZZ(M);
85 NTL::mat_ZZ *V = (U ? mutableMatrix_to_NTL_mat_ZZ(U) : nullptr);
86
87 switch (strategy)
88 {
89 case 2:
90 if (!V)
91 LLL(d, *A, numer, denom);
92 else
93 LLL(d, *A, *V, numer, denom);
94 break;
95
96 case useNTL + GS + useLLL + FP:
97 if (!V)
98 LLL_FP(*A, delta);
99 else
100 LLL_FP(*A, *V, delta);
101 break;
102 case useNTL + GS + useLLL + QP:
103 case useNTL + GS + useLLL + QP1:
104 if (!V)
105 LLL_QP(*A, delta);
106 else
107 LLL_QP(*A, *V, delta);
108 break;
109 case useNTL + GS + useLLL + XD:
110 if (!V)
111 LLL_XD(*A, delta);
112 else
113 LLL_XD(*A, *V, delta);
114 break;
115 case useNTL + GS + useLLL + useRR:
116 if (!V)
117 LLL_RR(*A, delta);
118 else
119 LLL_RR(*A, *V, delta);
120 break;
121
122 case useNTL + GS + useBKZ + FP:
123 if (!V)
124 BKZ_FP(*A, delta);
125 else
126 BKZ_FP(*A, *V, delta);
127 break;
128 case useNTL + GS + useBKZ + QP:
129 if (!V)
130 BKZ_QP(*A, delta);
131 else
132 BKZ_QP(*A, *V, delta);
133 break;
134 case useNTL + GS + useBKZ + QP1:
135 if (!V)
136 BKZ_QP1(*A, delta);
137 else
138 BKZ_QP1(*A, *V, delta);
139 break;
140 case useNTL + GS + useBKZ + XD:
141 if (!V)
142 BKZ_XD(*A, delta);
143 else
144 BKZ_XD(*A, *V, delta);
145 break;
146 case useNTL + GS + useBKZ + useRR:
147 if (!V)
148 BKZ_RR(*A, delta);
149 else
150 BKZ_RR(*A, *V, delta);
151 break;
152
153 case useNTL + Givens + useLLL + FP:
154 if (!V)
155 G_LLL_FP(*A, delta);
156 else
157 G_LLL_FP(*A, *V, delta);
158 break;
159 case useNTL + Givens + useLLL + QP:
160 case useNTL + Givens + useLLL + QP1:
161 if (!V)
162 G_LLL_QP(*A, delta);
163 else
164 G_LLL_QP(*A, *V, delta);
165 break;
166 case useNTL + Givens + useLLL + XD:
167 if (!V)
168 G_LLL_XD(*A, delta);
169 else
170 G_LLL_XD(*A, *V, delta);
171 break;
172 case useNTL + Givens + useLLL + useRR:
173 if (!V)
174 G_LLL_RR(*A, delta);
175 else
176 G_LLL_RR(*A, *V, delta);
177 break;
178
179 case useNTL + Givens + useBKZ + FP:
180 if (!V)
181 G_BKZ_FP(*A, delta);
182 else
183 G_BKZ_FP(*A, *V, delta);
184 break;
185 case useNTL + Givens + useBKZ + QP:
186 if (!V)
187 G_BKZ_QP(*A, delta);
188 else
189 G_BKZ_QP(*A, *V, delta);
190 break;
191 case useNTL + Givens + useBKZ + QP1:
192 if (!V)
193 G_BKZ_QP1(*A, delta);
194 else
195 G_BKZ_QP1(*A, *V, delta);
196 break;
197 case useNTL + Givens + useBKZ + XD:
198 if (!V)
199 G_BKZ_XD(*A, delta);
200 else
201 G_BKZ_XD(*A, *V, delta);
202 break;
203 case useNTL + Givens + useBKZ + useRR:
204 if (!V)
205 G_BKZ_RR(*A, delta);
206 else
207 G_BKZ_RR(*A, *V, delta);
208 break;
209 default:
210 delete A;
211 if (V) delete V;
212 ERROR("Strategy option to LLL not understood");
213 return false;
214 }
215
216 /* Put this back into M */
217 mpz_t a;
218 mpz_init(a);
219
220 for (int j = 0; j < ncols; j++)
221 for (int i = 0; i < nrows; i++)
222 {
223 mat_ZZ_get_entry(A, j, i, a);
224 ring_elem b = globalZZ->from_int(a);
225 M->set_entry(i, j, b);
226 }
227
228 if (U)
229 {
230 for (int j = 0; j < ncols; j++)
231 for (int i = 0; i < ncols; i++)
232 {
233 mat_ZZ_get_entry(V, j, i, a);
234 ring_elem b = globalZZ->from_int(a);
235 U->set_entry(i, j, b);
236 }
237 }
238 mpz_clear(a);
239 delete A;
240 if (V) delete V;
241 return true;
242}
MutableMatrix::n_rows
virtual size_t n_rows() const =0
MutableMatrix::n_cols
virtual size_t n_cols() const =0
MutableMatrix::set_entry
virtual bool set_entry(size_t r, size_t c, const ring_elem a)=0
globalZZ
RingZZ * globalZZ
Definition relem.cpp:13
ERROR
const int ERROR
Definition m2-mem.cpp:55
M2_gbTrace
int M2_gbTrace
Definition m2-types.cpp:52
useLLL
static const int useLLL
Definition ntl-interface.cpp:61
GS
static const int GS
Definition ntl-interface.cpp:59
Givens
static const int Givens
Definition ntl-interface.cpp:60
useNTL
static const int useNTL
Definition ntl-interface.cpp:58
useRR
static const int useRR
Definition ntl-interface.cpp:67
FP
static const int FP
Definition ntl-interface.cpp:63
QP
static const int QP
Definition ntl-interface.cpp:65
useBKZ
static const int useBKZ
Definition ntl-interface.cpp:62
mutableMatrix_to_NTL_mat_ZZ
NTL::mat_ZZ * mutableMatrix_to_NTL_mat_ZZ(const MutableMatrix *M)
Definition ntl-interface.cpp:6
XD
static const int XD
Definition ntl-interface.cpp:66
QP1
static const int QP1
Definition ntl-interface.cpp:64
mat_ZZ_get_entry
void mat_ZZ_get_entry(const NTL::mat_ZZ *A, long i, long j, mpz_t result)
Definition ntl-internal.cpp:61
ring_elem
Definition ringelem.hpp:85

References ERROR, FP, Givens, globalZZ, GS, M2_gbTrace, mat_ZZ_get_entry(), mutableMatrix_to_NTL_mat_ZZ(), MutableMatrix::n_cols(), MutableMatrix::n_rows(), QP, QP1, MutableMatrix::set_entry(), useBKZ, useLLL, useNTL, useRR, and XD.

Referenced by rawLLL().