blob: 95d44d7452f9c6b162bfca17dab7d45370162795 [file] [log] [blame]
caryclark@google.com639df892012-01-10 21:46:10 +00001#include "CubicIntersection_TestData.h"
2#include <limits>
3
4const Cubic pointDegenerates[] = {
5 {{0, 0}, {0, 0}, {0, 0}, {0, 0}},
6 {{1, 1}, {1, 1}, {1, 1}, {1, 1}},
7 {{1 + PointEpsilon - std::numeric_limits<double>::epsilon(), 1},
8 {1, 1 + PointEpsilon - std::numeric_limits<double>::epsilon()}, {1, 1}, {1, 1}},
9 {{1 + PointEpsilon/2 - std::numeric_limits<double>::epsilon(), 1},
10 {1 - (PointEpsilon/2 - std::numeric_limits<double>::epsilon()), 1}, {1, 1}, {1, 1}}
11};
12
13const size_t pointDegenerates_count = sizeof(pointDegenerates) / sizeof(pointDegenerates[0]);
14
15const Cubic notPointDegenerates[] = {
16 {{1 + PointEpsilon + std::numeric_limits<double>::epsilon(), 1}, {1, 1 + PointEpsilon}, {1, 1}, {1, 1}},
17 {{1 + PointEpsilon/2 + std::numeric_limits<double>::epsilon(), 1}, {1 - PointEpsilon/2, 1}, {1, 1}, {1, 1}}
18};
19
20const size_t notPointDegenerates_count = sizeof(notPointDegenerates) / sizeof(notPointDegenerates[0]);
21
22// from http://www.truetex.com/bezint.htm
23const Cubic tests[][2] = {
24 { // intersects in one place (data gives bezier clip fits
25 {{0, 45},
26 {6.0094158284751593, 51.610357411322688},
27 {12.741093228940867, 55.981703949474607},
28 {20.021417396476362, 58.652245509710262}},
29 {{2.2070737699246674, 52.703494107327209},
30 {31.591482272629477, 23.811002295222025},
31 {76.824588616426425, 44.049473790502674},
32 {119.25488947221436, 55.599248272955073}}
33 },
34 { // intersects in three places
35 {{0, 45}, {50, 100}, {150, 0}, {200, 55}},
36 {{0, 55}, {50, 0}, {150, 100}, {200, 45}}
37 },
38 { // intersects in one place, cross over is nearly parallel
39 {{0, 0}, {0, 100}, {200, 0}, {200, 100}},
40 {{0, 100}, {0, 0}, {200, 100}, {200, 0}}
41 },
42 { // intersects in two places
43 {{0, 0}, {0, 100}, {200, 100}, {200, 0}},
44 {{0, 100}, {0, 0}, {200, 0}, {200, 100}}
45 },
46 {
47 {{150, 100}, {150 + 0.1, 150}, {150, 200}, {150, 250}},
48 {{250, 150}, {200, 150 + 0.1}, {150, 150}, {100, 150}}
49 },
50 { // single intersection around 168,185
51 {{200, 100}, {150, 100}, {150, 150}, {200, 150}},
52 {{250, 150}, {250, 100}, {100, 100}, {100, 150}}
53 },
54 {
55 {{1.0, 1.5}, {15.5, 0.5}, {-8.0, 3.5}, {5.0, 1.5}},
56 {{4.0, 0.5}, {5.0, 15.0}, {2.0, -8.5}, {4.0, 4.5}}
57 },
58 {
59 {{664.00168, 0}, {726.11545, 124.22757}, {736.89069, 267.89743}, {694.0017, 400.0002}},
60 {{850.66843, 115.55563}, {728.515, 115.55563}, {725.21347, 275.15309}, {694.0017, 400.0002}}
61 },
62 {
63 {{1, 1}, {12.5, 6.5}, {-4, 6.5}, {7.5, 1}},
64 {{1, 6.5}, {12.5, 1}, {-4, 1}, {.5, 6}}
65 },
66 {
67 {{315.748, 312.84}, {312.644, 318.134}, {305.836, 319.909}, {300.542, 316.804}},
68 {{317.122, 309.05}, {316.112, 315.102}, {310.385, 319.19}, {304.332, 318.179}}
69 },
70 {
71 {{1046.604051, 172.937967}, {1046.604051, 178.9763059}, {1041.76745, 183.9279165}, {1035.703842, 184.0432409}},
72 {{1046.452235, 174.7640504}, {1045.544872, 180.1973817}, {1040.837966, 184.0469882}, {1035.505925, 184.0469882}}
73 },
74 {
75 {{125.79356, 199.57382}, {51.16556, 128.93575}, {87.494, 16.67848}, {167.29361, 16.67848}},
76 {{167.29361, 55.81876}, {100.36128, 55.81876}, {68.64099, 145.4755}, {125.7942, 199.57309}}
77 }
78};
79
80const size_t tests_count = sizeof(tests) / sizeof(tests[0]);
81
82Cubic hexTests[][2] = {
83 {
84 {{0}} // placeholder for hex converted below
85 }
86};
87
88const size_t hexTests_count = sizeof(hexTests) / sizeof(hexTests[0]);
89
90static const uint64_t testx[2][8] = {
91 {
92 0xf0d0d1ca63075a40LLU, 0x9408ce996a237740LLU, 0x6d5675460fbe5e40LLU, 0x6ef501e1b7487940LLU,
93 0x9a71d2f8143d6540LLU, 0x6bc18bbe02907a40LLU, 0x5b94d92093aa6b40LLU, 0x6ac18bbe02907a40LLU
94 },
95 {
96 0x92c56ed7b6145d40LLU, 0xede4f1255edb7740LLU, 0x1138c1101af75940LLU, 0x42e4f1255edb7740LLU,
97 0x408e51603ad95640LLU, 0x1e2e8fe9dd927740LLU, 0x1cb4777cd3a75440LLU, 0x212e1390de017740LLU
98 }
99};
100
101void convert_testx() {
102 const uint64_t* inPtr = testx[0];
103 double* outPtr = &hexTests[sizeof(tests) / sizeof(tests[0]) - 1][0][0].x;
104 for (unsigned index = 0; index < sizeof(testx) / sizeof(testx[0][0]); ++index) {
105 uint64_t input = *inPtr++;
106 unsigned char* output = (unsigned char*) outPtr++;
107 for (unsigned byte = 0; byte < sizeof(input); ++byte) {
108 output[byte] = input >> (7 - byte) * 8;
109 }
110 }
111}
112
113const Cubic lines[] = {
114 {{0, 0}, {0, 0}, {0, 0}, {1, 0}}, // 0: horizontal
115 {{0, 0}, {0, 0}, {1, 0}, {0, 0}},
116 {{0, 0}, {1, 0}, {0, 0}, {0, 0}},
117 {{1, 0}, {0, 0}, {0, 0}, {0, 0}},
118 {{1, 0}, {2, 0}, {3, 0}, {4, 0}},
119 {{0, 0}, {0, 0}, {0, 0}, {0, 1}}, // 5: vertical
120 {{0, 0}, {0, 0}, {0, 1}, {0, 0}},
121 {{0, 0}, {0, 1}, {0, 0}, {0, 0}},
122 {{0, 1}, {0, 0}, {0, 0}, {0, 0}},
123 {{0, 1}, {0, 2}, {0, 3}, {0, 4}},
124 {{0, 0}, {0, 0}, {0, 0}, {1, 1}}, // 10: 3 coincident
125 {{0, 0}, {0, 0}, {1, 1}, {0, 0}},
126 {{0, 0}, {1, 1}, {0, 0}, {0, 0}},
127 {{1, 1}, {0, 0}, {0, 0}, {0, 0}},
128 {{0, 0}, {0, 0}, {1, 1}, {2, 2}}, // 14: 2 coincident
129 {{0, 0}, {1, 1}, {0, 0}, {2, 2}},
130 {{0, 0}, {1, 1}, {2, 2}, {0, 0}},
131 {{1, 1}, {0, 0}, {0, 0}, {2, 2}}, // 17:
132 {{1, 1}, {0, 0}, {2, 2}, {0, 0}},
133 {{1, 1}, {2, 2}, {0, 0}, {0, 0}},
134 {{1, 1}, {2, 2}, {3, 3}, {2, 2}}, // middle-last coincident
135 {{1, 1}, {2, 2}, {3, 3}, {3, 3}}, // middle-last coincident
136 {{1, 1}, {1, 1}, {2, 2}, {2, 2}}, // 2 pairs coincident
137 {{1, 1}, {2, 2}, {1, 1}, {2, 2}},
138 {{1, 1}, {2, 2}, {2, 2}, {1, 1}},
139 {{1, 1}, {1, 1}, {3, 3}, {3, 3}}, // first-middle middle-last coincident
140 {{1, 1}, {2, 2}, {3, 3}, {4, 4}}, // no coincident
141 {{1, 1}, {3, 3}, {2, 2}, {4, 4}},
142 {{1, 1}, {2, 2}, {4, 4}, {3, 3}},
143 {{1, 1}, {3, 3}, {4, 4}, {2, 2}},
144 {{1, 1}, {4, 4}, {2, 2}, {3, 3}},
145 {{1, 1}, {4, 4}, {3, 3}, {2, 2}},
146 {{2, 2}, {1, 1}, {3, 3}, {4, 4}},
147 {{2, 2}, {1, 1}, {4, 4}, {3, 3}},
148 {{2, 2}, {3, 3}, {1, 1}, {4, 4}},
149 {{2, 2}, {3, 3}, {4, 4}, {1, 1}},
150 {{2, 2}, {4, 4}, {1, 1}, {3, 3}},
151 {{2, 2}, {4, 4}, {3, 3}, {1, 1}},
152};
153
154const size_t lines_count = sizeof(lines) / sizeof(lines[0]);
155
156// 'not a line' tries to fool the line detection code
157const Cubic notLines[] = {
158 {{0, 0}, {0, 0}, {0, 1}, {1, 0}},
159 {{0, 0}, {0, 1}, {0, 0}, {1, 0}},
160 {{0, 0}, {0, 1}, {1, 0}, {0, 0}},
161 {{0, 1}, {0, 0}, {0, 0}, {1, 0}},
162 {{0, 1}, {0, 0}, {1, 0}, {0, 0}},
163 {{0, 1}, {1, 0}, {0, 0}, {0, 0}},
164};
165
166const size_t notLines_count = sizeof(notLines) / sizeof(notLines[0]);
167
168static const double E = PointEpsilon * 2;
169static const double F = PointEpsilon * 3;
170static const double H = PointEpsilon * 4;
171static const double J = PointEpsilon * 5;
172static const double K = PointEpsilon * 8; // INVESTIGATE: why are larger multiples necessary?
173
174const Cubic modEpsilonLines[] = {
175 {{0, E}, {0, 0}, {0, 0}, {1, 0}}, // horizontal
176 {{0, 0}, {0, E}, {1, 0}, {0, 0}},
177 {{0, 0}, {1, 0}, {0, E}, {0, 0}},
178 {{1, 0}, {0, 0}, {0, 0}, {0, E}},
179 {{1, E}, {2, 0}, {3, 0}, {4, 0}},
180 {{E, 0}, {0, 0}, {0, 0}, {0, 1}}, // vertical
181 {{0, 0}, {E, 0}, {0, 1}, {0, 0}},
182 {{0, 0}, {0, 1}, {E, 0}, {0, 0}},
183 {{0, 1}, {0, 0}, {0, 0}, {E, 0}},
184 {{E, 1}, {0, 2}, {0, 3}, {0, 4}},
185 {{E, 0}, {0, 0}, {0, 0}, {1, 1}}, // 3 coincident
186 {{0, 0}, {E, 0}, {1, 1}, {0, 0}},
187 {{0, 0}, {1, 1}, {E, 0}, {0, 0}},
188 {{1, 1}, {0, 0}, {0, 0}, {E, 0}},
189 {{0, E}, {0, 0}, {1, 1}, {2, 2}}, // 2 coincident
190 {{0, 0}, {1, 1}, {0, E}, {2, 2}},
191 {{0, 0}, {1, 1}, {2, 2}, {0, E}},
192 {{1, 1}, {0, E}, {0, 0}, {2, 2}},
193 {{1, 1}, {0, E}, {2, 2}, {0, 0}},
194 {{1, 1}, {2, 2}, {E, 0}, {0, 0}},
195 {{1, 1}, {2, 2+E}, {3, 3}, {2, 2}}, // middle-last coincident
196 {{1, 1}, {2+E, 2}, {3, 3}, {3, 3}}, // middle-last coincident
197 {{1, 1}, {1, 1}, {2, 2}, {2+E, 2}}, // 2 pairs coincident
198 {{1, 1}, {2, 2}, {1, 1}, {2+E, 2}},
199 {{1, 1}, {2, 2}, {2, 2+E}, {1, 1}},
200 {{1, 1}, {1, 1+E}, {3, 3}, {3, 3}}, // first-middle middle-last coincident
201 {{1, 1}, {2+E, 2}, {3, 3}, {4, 4}}, // no coincident
202 {{1, 1}, {3, 3}, {2, 2}, {4, 4+F}}, // INVESTIGATE: why the epsilon is bigger
203 {{1, 1+F}, {2, 2}, {4, 4}, {3, 3}}, // INVESTIGATE: why the epsilon is bigger
204 {{1, 1}, {3, 3}, {4, 4+E}, {2, 2}},
205 {{1, 1}, {4, 4}, {2, 2}, {3, 3+E}},
206 {{1, 1}, {4, 4}, {3, 3}, {2+E, 2}},
207 {{2, 2}, {1, 1}, {3+E, 3}, {4, 4}},
208 {{2, 2}, {1+E, 1}, {4, 4}, {3, 3}},
209 {{2, 2+E}, {3, 3}, {1, 1}, {4, 4}},
210 {{2+E, 2}, {3, 3}, {4, 4}, {1, 1}},
211 {{2, 2}, {4+E, 4}, {1, 1}, {3, 3}},
212 {{2, 2}, {4, 4}, {3, 3}, {1, 1+E}},
213};
214
215const size_t modEpsilonLines_count = sizeof(modEpsilonLines) / sizeof(modEpsilonLines[0]);
216
217static const double D = PointEpsilon / 2;
218static const double G = PointEpsilon / 3;
219
220const Cubic lessEpsilonLines[] = {
221 {{0, D}, {0, 0}, {0, 0}, {1, 0}}, // horizontal
222 {{0, 0}, {0, D}, {1, 0}, {0, 0}},
223 {{0, 0}, {1, 0}, {0, D}, {0, 0}},
224 {{1, 0}, {0, 0}, {0, 0}, {0, D}},
225 {{1, D}, {2, 0}, {3, 0}, {4, 0}},
226 {{D, 0}, {0, 0}, {0, 0}, {0, 1}}, // vertical
227 {{0, 0}, {D, 0}, {0, 1}, {0, 0}},
228 {{0, 0}, {0, 1}, {D, 0}, {0, 0}},
229 {{0, 1}, {0, 0}, {0, 0}, {D, 0}},
230 {{D, 1}, {0, 2}, {0, 3}, {0, 4}},
231 {{D, 0}, {0, 0}, {0, 0}, {1, 1}}, // 3 coincident
232 {{0, 0}, {D, 0}, {1, 1}, {0, 0}},
233 {{0, 0}, {1, 1}, {D, 0}, {0, 0}},
234 {{1, 1}, {0, 0}, {0, 0}, {D, 0}},
235 {{0, D}, {0, 0}, {1, 1}, {2, 2}}, // 2 coincident
236 {{0, 0}, {1, 1}, {0, D}, {2, 2}},
237 {{0, 0}, {1, 1}, {2, 2}, {0, D}},
238 {{1, 1}, {0, D}, {0, 0}, {2, 2}},
239 {{1, 1}, {0, D}, {2, 2}, {0, 0}},
240 {{1, 1}, {2, 2}, {D, 0}, {0, 0}},
241 {{1, 1}, {2, 2+D}, {3, 3}, {2, 2}}, // middle-last coincident
242 {{1, 1}, {2+D, 2}, {3, 3}, {3, 3}}, // middle-last coincident
243 {{1, 1}, {1, 1}, {2, 2}, {2+D, 2}}, // 2 pairs coincident
244 {{1, 1}, {2, 2}, {1, 1}, {2+D, 2}},
245 {{1, 1}, {2, 2}, {2, 2+D}, {1, 1}},
246 {{1, 1}, {1, 1+D}, {3, 3}, {3, 3}}, // first-middle middle-last coincident
247 {{1, 1}, {2+D/2, 2}, {3, 3}, {4, 4}}, // no coincident (FIXME: N as opposed to N/2 failed)
248 {{1, 1}, {3, 3}, {2, 2}, {4, 4+D}},
249 {{1, 1+D}, {2, 2}, {4, 4}, {3, 3}},
250 {{1, 1}, {3, 3}, {4, 4+D}, {2, 2}},
251 {{1, 1}, {4, 4}, {2, 2}, {3, 3+D}},
252 {{1, 1}, {4, 4}, {3, 3}, {2+G, 2}}, // INVESTIGATE: why the epsilon is smaller
253 {{2, 2}, {1, 1}, {3+D, 3}, {4, 4}},
254 {{2, 2}, {1+D, 1}, {4, 4}, {3, 3}},
255 {{2, 2+D}, {3, 3}, {1, 1}, {4, 4}},
256 {{2+G, 2}, {3, 3}, {4, 4}, {1, 1}}, // INVESTIGATE: why the epsilon is smaller
257 {{2, 2}, {4+D, 4}, {1, 1}, {3, 3}},
258 {{2, 2}, {4, 4}, {3, 3}, {1, 1+D}},
259};
260
261const size_t lessEpsilonLines_count = sizeof(lessEpsilonLines) / sizeof(lessEpsilonLines[0]);
262
263static const double N = -PointEpsilon / 2;
264static const double M = -PointEpsilon / 3;
265
266const Cubic negEpsilonLines[] = {
267 {{0, N}, {0, 0}, {0, 0}, {1, 0}}, // horizontal
268 {{0, 0}, {0, N}, {1, 0}, {0, 0}},
269 {{0, 0}, {1, 0}, {0, N}, {0, 0}},
270 {{1, 0}, {0, 0}, {0, 0}, {0, N}},
271 {{1, N}, {2, 0}, {3, 0}, {4, 0}},
272 {{N, 0}, {0, 0}, {0, 0}, {0, 1}}, // vertical
273 {{0, 0}, {N, 0}, {0, 1}, {0, 0}},
274 {{0, 0}, {0, 1}, {N, 0}, {0, 0}},
275 {{0, 1}, {0, 0}, {0, 0}, {N, 0}},
276 {{N, 1}, {0, 2}, {0, 3}, {0, 4}},
277 {{N, 0}, {0, 0}, {0, 0}, {1, 1}}, // 3 coincident
278 {{0, 0}, {N, 0}, {1, 1}, {0, 0}},
279 {{0, 0}, {1, 1}, {N, 0}, {0, 0}},
280 {{1, 1}, {0, 0}, {0, 0}, {N, 0}},
281 {{0, N}, {0, 0}, {1, 1}, {2, 2}}, // 2 coincident
282 {{0, 0}, {1, 1}, {0, N}, {2, 2}},
283 {{0, 0}, {1, 1}, {2, 2}, {0, N}},
284 {{1, 1}, {0, N}, {0, 0}, {2, 2}},
285 {{1, 1}, {0, N}, {2, 2}, {0, 0}},
286 {{1, 1}, {2, 2}, {N, 0}, {0, 0}},
287 {{1, 1}, {2, 2+N}, {3, 3}, {2, 2}}, // middle-last coincident
288 {{1, 1}, {2+N, 2}, {3, 3}, {3, 3}}, // middle-last coincident
289 {{1, 1}, {1, 1}, {2, 2}, {2+N, 2}}, // 2 pairs coincident
290 {{1, 1}, {2, 2}, {1, 1}, {2+N, 2}},
291 {{1, 1}, {2, 2}, {2, 2+N}, {1, 1}},
292 {{1, 1}, {1, 1+N}, {3, 3}, {3, 3}}, // first-middle middle-last coincident
293 {{1, 1}, {2+N/2, 2}, {3, 3}, {4, 4}}, // no coincident (FIXME: N as opposed to N/2 failed)
294 {{1, 1}, {3, 3}, {2, 2}, {4, 4+N}},
295 {{1, 1+N}, {2, 2}, {4, 4}, {3, 3}},
296 {{1, 1}, {3, 3}, {4, 4+N}, {2, 2}},
297 {{1, 1}, {4, 4}, {2, 2}, {3, 3+N}},
298 {{1, 1}, {4, 4}, {3, 3}, {2+M, 2}}, // INVESTIGATE: why the epsilon is smaller
299 {{2, 2}, {1, 1}, {3+N, 3}, {4, 4}},
300 {{2, 2}, {1+N, 1}, {4, 4}, {3, 3}},
301 {{2, 2+N}, {3, 3}, {1, 1}, {4, 4}},
302 {{2+M, 2}, {3, 3}, {4, 4}, {1, 1}}, // INVESTIGATE: why the epsilon is smaller
303 {{2, 2}, {4+N, 4}, {1, 1}, {3, 3}},
304 {{2, 2}, {4, 4}, {3, 3}, {1, 1+N}},
305};
306
307const size_t negEpsilonLines_count = sizeof(negEpsilonLines) / sizeof(negEpsilonLines[0]);
308
309const Quadratic quadraticLines[] = {
310 {{0, 0}, {0, 0}, {1, 0}},
311 {{0, 0}, {1, 0}, {0, 0}},
312 {{1, 0}, {0, 0}, {0, 0}},
313 {{1, 0}, {2, 0}, {3, 0}},
314 {{0, 0}, {0, 0}, {0, 1}},
315 {{0, 0}, {0, 1}, {0, 0}},
316 {{0, 1}, {0, 0}, {0, 0}},
317 {{0, 1}, {0, 2}, {0, 3}},
318 {{0, 0}, {0, 0}, {1, 1}},
319 {{0, 0}, {1, 1}, {0, 0}},
320 {{1, 1}, {0, 0}, {0, 0}},
321 {{1, 1}, {2, 2}, {3, 3}},
322 {{1, 1}, {3, 3}, {3, 3}},
323 {{1, 1}, {1, 1}, {2, 2}},
324 {{1, 1}, {2, 2}, {1, 1}},
325 {{1, 1}, {1, 1}, {3, 3}},
326 {{1, 1}, {2, 2}, {4, 4}}, // no coincident
327 {{1, 1}, {3, 3}, {4, 4}},
328 {{1, 1}, {3, 3}, {2, 2}},
329 {{1, 1}, {4, 4}, {2, 2}},
330 {{1, 1}, {4, 4}, {3, 3}},
331 {{2, 2}, {1, 1}, {3, 3}},
332 {{2, 2}, {1, 1}, {4, 4}},
333 {{2, 2}, {3, 3}, {1, 1}},
334 {{2, 2}, {3, 3}, {4, 4}},
335 {{2, 2}, {4, 4}, {1, 1}},
336 {{2, 2}, {4, 4}, {3, 3}},
337};
338
339const size_t quadraticLines_count = sizeof(quadraticLines) / sizeof(quadraticLines[0]);
340
341const Quadratic quadraticModEpsilonLines[] = {
342 {{0, F}, {0, 0}, {1, 0}},
343 {{0, 0}, {1, 0}, {0, F}},
344 {{1, 0}, {0, F}, {0, 0}},
345 {{1, H}, {2, 0}, {3, 0}},
346 {{F, 0}, {0, 0}, {0, 1}},
347 {{0, 0}, {0, 1}, {F, 0}},
348 {{0, 1}, {F, 0}, {0, 0}},
349 {{H, 1}, {0, 2}, {0, 3}},
350 {{0, F}, {0, 0}, {1, 1}},
351 {{0, 0}, {1, 1}, {F, 0}},
352 {{1, 1}, {F, 0}, {0, 0}},
353 {{1, 1+J}, {2, 2}, {3, 3}},
354 {{1, 1}, {3, 3}, {3+F, 3}},
355 {{1, 1}, {1+F, 1}, {2, 2}},
356 {{1, 1}, {2, 2}, {1, 1+F}},
357 {{1, 1}, {1, 1+F}, {3, 3}},
358 {{1+H, 1}, {2, 2}, {4, 4}}, // no coincident
359 {{1, 1+K}, {3, 3}, {4, 4}},
360 {{1, 1}, {3+F, 3}, {2, 2}},
361 {{1, 1}, {4, 4+F}, {2, 2}},
362 {{1, 1}, {4, 4}, {3+F, 3}},
363 {{2, 2}, {1, 1}, {3, 3+F}},
364 {{2+F, 2}, {1, 1}, {4, 4}},
365 {{2, 2+F}, {3, 3}, {1, 1}},
366 {{2, 2}, {3+F, 3}, {4, 4}},
367 {{2, 2}, {4, 4+F}, {1, 1}},
368 {{2, 2}, {4, 4}, {3+F, 3}},
369};
370
371const size_t quadraticModEpsilonLines_count = sizeof(quadraticModEpsilonLines) / sizeof(quadraticModEpsilonLines[0]);
372
373