primal-sieve 0.2.0

A high performance prime sieve.
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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
// automatically generated
use wheel::{WheelInit, Wheel, WheelElem};

#[derive(Debug)]
pub struct Wheel30;
impl Wheel for Wheel30 {
    #[inline(always)]
    fn modulo(&self) -> usize { MODULO }

    #[inline(always)]
    fn size(&self) -> usize { SIZE }

    #[inline(always)]
    fn wheel(&self) -> &'static [WheelElem] { WHEEL }

    #[inline(always)]
    fn init(&self) -> &'static [WheelInit] { INIT }

    #[inline(always)]
    unsafe fn hardcoded_sieve(&self,
                              bytes: &mut [u8], si_: &mut usize, wi_: &mut usize, prime: usize) {
        hardcoded_sieve(bytes, si_, wi_, prime)
    }
}

pub const SIZE: usize = 8;

pub const MODULO: usize = 30;

const INIT: &'static [WheelInit; 30] = &[
    WheelInit { next_mult_factor: 1, wheel_index: 0 }, // 0
    WheelInit { next_mult_factor: 0, wheel_index: 0 }, // 1
    WheelInit { next_mult_factor: 5, wheel_index: 1 }, // 2
    WheelInit { next_mult_factor: 4, wheel_index: 1 }, // 3
    WheelInit { next_mult_factor: 3, wheel_index: 1 }, // 4
    WheelInit { next_mult_factor: 2, wheel_index: 1 }, // 5
    WheelInit { next_mult_factor: 1, wheel_index: 1 }, // 6
    WheelInit { next_mult_factor: 0, wheel_index: 1 }, // 7
    WheelInit { next_mult_factor: 3, wheel_index: 2 }, // 8
    WheelInit { next_mult_factor: 2, wheel_index: 2 }, // 9
    WheelInit { next_mult_factor: 1, wheel_index: 2 }, // 10
    WheelInit { next_mult_factor: 0, wheel_index: 2 }, // 11
    WheelInit { next_mult_factor: 1, wheel_index: 3 }, // 12
    WheelInit { next_mult_factor: 0, wheel_index: 3 }, // 13
    WheelInit { next_mult_factor: 3, wheel_index: 4 }, // 14
    WheelInit { next_mult_factor: 2, wheel_index: 4 }, // 15
    WheelInit { next_mult_factor: 1, wheel_index: 4 }, // 16
    WheelInit { next_mult_factor: 0, wheel_index: 4 }, // 17
    WheelInit { next_mult_factor: 1, wheel_index: 5 }, // 18
    WheelInit { next_mult_factor: 0, wheel_index: 5 }, // 19
    WheelInit { next_mult_factor: 3, wheel_index: 6 }, // 20
    WheelInit { next_mult_factor: 2, wheel_index: 6 }, // 21
    WheelInit { next_mult_factor: 1, wheel_index: 6 }, // 22
    WheelInit { next_mult_factor: 0, wheel_index: 6 }, // 23
    WheelInit { next_mult_factor: 5, wheel_index: 7 }, // 24
    WheelInit { next_mult_factor: 4, wheel_index: 7 }, // 25
    WheelInit { next_mult_factor: 3, wheel_index: 7 }, // 26
    WheelInit { next_mult_factor: 2, wheel_index: 7 }, // 27
    WheelInit { next_mult_factor: 1, wheel_index: 7 }, // 28
    WheelInit { next_mult_factor: 0, wheel_index: 7 }, // 29
];
const WHEEL: &'static [WheelElem; 64] = &[
    // remainder 1
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 6, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 4, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 2, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 4, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 2, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 4, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 6, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 2, correction: 1, next: -7 },
    // remainder 7
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 4, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 4, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 4, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 6, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 6, correction: 1, next: -7 },
    // remainder 11
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 2, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 2, correction: 0, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 6, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 6, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 4, correction: 2, next: -7 },
    // remainder 13
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 4, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 6, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 6, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 4, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 2, correction: 1, next: -7 },
    // remainder 17
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 4, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 6, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 6, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 4, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 4, correction: 2, next: -7 },
    // remainder 19
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 6, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 6, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 2, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 4, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 2, correction: 2, next: -7 },
    // remainder 23
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 6, correction: 5, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 6, correction: 5, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 4, correction: 3, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 4, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 4, correction: 3, next: -7 },
    // remainder 29
    WheelElem { unset_bit: 1u8 << 0u8, next_mult_factor: 2, correction: 1, next: 1 },
    WheelElem { unset_bit: 1u8 << 7u8, next_mult_factor: 6, correction: 6, next: 1 },
    WheelElem { unset_bit: 1u8 << 6u8, next_mult_factor: 4, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 5u8, next_mult_factor: 2, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 4u8, next_mult_factor: 4, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 3u8, next_mult_factor: 2, correction: 2, next: 1 },
    WheelElem { unset_bit: 1u8 << 2u8, next_mult_factor: 4, correction: 4, next: 1 },
    WheelElem { unset_bit: 1u8 << 1u8, next_mult_factor: 6, correction: 6, next: -7 },
];
pub unsafe fn hardcoded_sieve(bytes: &mut [u8], si_: &mut usize, wi_: &mut usize, prime: usize) {
    let bytes = bytes;
    let start = bytes.as_mut_ptr();
    let len = bytes.len() as isize;
    let largest_step = (28 * prime + 23) as isize;
    let loop_len = len - largest_step;
    let loop_end = start.offset(loop_len);
    let end = start.offset(len);
    let si = *si_ as isize;
    let mut p = start.offset(si);
    let mut wi = *wi_;
    let prime_ = prime as isize;

    'outer: loop {
    match wi {
        0...7 => { // 30 * x + 1
            loop {
             'label7: loop {
              'label6: loop {
               'label5: loop {
                'label4: loop {
                 'label3: loop {
                  'label2: loop {
                   'label1: loop {
                    'label0: loop {
                     match wi {
                      0 => break 'label0,
                      1 => break 'label1,
                      2 => break 'label2,
                      3 => break 'label3,
                      4 => break 'label4,
                      5 => break 'label5,
                      6 => break 'label6,
                      _ => break 'label7,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 1;
                        *p.offset(prime_ * 6 + 0) |= 2;
                        *p.offset(prime_ * 10 + 0) |= 4;
                        *p.offset(prime_ * 12 + 0) |= 8;
                        *p.offset(prime_ * 16 + 0) |= 16;
                        *p.offset(prime_ * 18 + 0) |= 32;
                        *p.offset(prime_ * 22 + 0) |= 64;
                        *p.offset(prime_ * 28 + 0) |= 128;

                        p = p.offset(prime_ * 30 + 1)
                    }
                    if p >= end { wi = 0; break 'outer; }
                    *p |= 1; p = p.offset(prime_ * 6 + 0);
                    break 'label1
                   }
                   if p >= end { wi = 1; break 'outer; }
                   *p |= 2; p = p.offset(prime_ * 4 + 0);
                   break 'label2
                  }
                  if p >= end { wi = 2; break 'outer; }
                  *p |= 4; p = p.offset(prime_ * 2 + 0);
                  break 'label3
                 }
                 if p >= end { wi = 3; break 'outer; }
                 *p |= 8; p = p.offset(prime_ * 4 + 0);
                 break 'label4
                }
                if p >= end { wi = 4; break 'outer; }
                *p |= 16; p = p.offset(prime_ * 2 + 0);
                break 'label5
               }
               if p >= end { wi = 5; break 'outer; }
               *p |= 32; p = p.offset(prime_ * 4 + 0);
               break 'label6
              }
              if p >= end { wi = 6; break 'outer; }
              *p |= 64; p = p.offset(prime_ * 6 + 0);
              break 'label7
             }
             if p >= end { wi = 7; break 'outer; }
             *p |= 128; p = p.offset(prime_ * 2 + 1);
             wi = 0
            }
        }
        8...15 => { // 30 * x + 7
            loop {
             'label15: loop {
              'label14: loop {
               'label13: loop {
                'label12: loop {
                 'label11: loop {
                  'label10: loop {
                   'label9: loop {
                    'label8: loop {
                     match wi {
                      8 => break 'label8,
                      9 => break 'label9,
                      10 => break 'label10,
                      11 => break 'label11,
                      12 => break 'label12,
                      13 => break 'label13,
                      14 => break 'label14,
                      _ => break 'label15,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 32;
                        *p.offset(prime_ * 4 + 1) |= 16;
                        *p.offset(prime_ * 6 + 2) |= 1;
                        *p.offset(prime_ * 10 + 2) |= 128;
                        *p.offset(prime_ * 12 + 3) |= 8;
                        *p.offset(prime_ * 16 + 4) |= 4;
                        *p.offset(prime_ * 22 + 5) |= 64;
                        *p.offset(prime_ * 24 + 6) |= 2;

                        p = p.offset(prime_ * 30 + 7)
                    }
                    if p >= end { wi = 8; break 'outer; }
                    *p |= 32; p = p.offset(prime_ * 4 + 1);
                    break 'label9
                   }
                   if p >= end { wi = 9; break 'outer; }
                   *p |= 16; p = p.offset(prime_ * 2 + 1);
                   break 'label10
                  }
                  if p >= end { wi = 10; break 'outer; }
                  *p |= 1; p = p.offset(prime_ * 4 + 0);
                  break 'label11
                 }
                 if p >= end { wi = 11; break 'outer; }
                 *p |= 128; p = p.offset(prime_ * 2 + 1);
                 break 'label12
                }
                if p >= end { wi = 12; break 'outer; }
                *p |= 8; p = p.offset(prime_ * 4 + 1);
                break 'label13
               }
               if p >= end { wi = 13; break 'outer; }
               *p |= 4; p = p.offset(prime_ * 6 + 1);
               break 'label14
              }
              if p >= end { wi = 14; break 'outer; }
              *p |= 64; p = p.offset(prime_ * 2 + 1);
              break 'label15
             }
             if p >= end { wi = 15; break 'outer; }
             *p |= 2; p = p.offset(prime_ * 6 + 1);
             wi = 8
            }
        }
        16...23 => { // 30 * x + 11
            loop {
             'label23: loop {
              'label22: loop {
               'label21: loop {
                'label20: loop {
                 'label19: loop {
                  'label18: loop {
                   'label17: loop {
                    'label16: loop {
                     match wi {
                      16 => break 'label16,
                      17 => break 'label17,
                      18 => break 'label18,
                      19 => break 'label19,
                      20 => break 'label20,
                      21 => break 'label21,
                      22 => break 'label22,
                      _ => break 'label23,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 1;
                        *p.offset(prime_ * 2 + 0) |= 64;
                        *p.offset(prime_ * 6 + 2) |= 2;
                        *p.offset(prime_ * 8 + 2) |= 128;
                        *p.offset(prime_ * 12 + 4) |= 8;
                        *p.offset(prime_ * 18 + 6) |= 32;
                        *p.offset(prime_ * 20 + 7) |= 4;
                        *p.offset(prime_ * 26 + 9) |= 16;

                        p = p.offset(prime_ * 30 + 11)
                    }
                    if p >= end { wi = 16; break 'outer; }
                    *p |= 1; p = p.offset(prime_ * 2 + 0);
                    break 'label17
                   }
                   if p >= end { wi = 17; break 'outer; }
                   *p |= 64; p = p.offset(prime_ * 4 + 2);
                   break 'label18
                  }
                  if p >= end { wi = 18; break 'outer; }
                  *p |= 2; p = p.offset(prime_ * 2 + 0);
                  break 'label19
                 }
                 if p >= end { wi = 19; break 'outer; }
                 *p |= 128; p = p.offset(prime_ * 4 + 2);
                 break 'label20
                }
                if p >= end { wi = 20; break 'outer; }
                *p |= 8; p = p.offset(prime_ * 6 + 2);
                break 'label21
               }
               if p >= end { wi = 21; break 'outer; }
               *p |= 32; p = p.offset(prime_ * 2 + 1);
               break 'label22
              }
              if p >= end { wi = 22; break 'outer; }
              *p |= 4; p = p.offset(prime_ * 6 + 2);
              break 'label23
             }
             if p >= end { wi = 23; break 'outer; }
             *p |= 16; p = p.offset(prime_ * 4 + 2);
             wi = 16
            }
        }
        24...31 => { // 30 * x + 13
            loop {
             'label31: loop {
              'label30: loop {
               'label29: loop {
                'label28: loop {
                 'label27: loop {
                  'label26: loop {
                   'label25: loop {
                    'label24: loop {
                     match wi {
                      24 => break 'label24,
                      25 => break 'label25,
                      26 => break 'label26,
                      27 => break 'label27,
                      28 => break 'label28,
                      29 => break 'label29,
                      30 => break 'label30,
                      _ => break 'label31,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 32;
                        *p.offset(prime_ * 4 + 2) |= 4;
                        *p.offset(prime_ * 6 + 3) |= 2;
                        *p.offset(prime_ * 10 + 4) |= 128;
                        *p.offset(prime_ * 16 + 7) |= 16;
                        *p.offset(prime_ * 18 + 8) |= 8;
                        *p.offset(prime_ * 24 + 11) |= 1;
                        *p.offset(prime_ * 28 + 12) |= 64;

                        p = p.offset(prime_ * 30 + 13)
                    }
                    if p >= end { wi = 24; break 'outer; }
                    *p |= 32; p = p.offset(prime_ * 4 + 2);
                    break 'label25
                   }
                   if p >= end { wi = 25; break 'outer; }
                   *p |= 4; p = p.offset(prime_ * 2 + 1);
                   break 'label26
                  }
                  if p >= end { wi = 26; break 'outer; }
                  *p |= 2; p = p.offset(prime_ * 4 + 1);
                  break 'label27
                 }
                 if p >= end { wi = 27; break 'outer; }
                 *p |= 128; p = p.offset(prime_ * 6 + 3);
                 break 'label28
                }
                if p >= end { wi = 28; break 'outer; }
                *p |= 16; p = p.offset(prime_ * 2 + 1);
                break 'label29
               }
               if p >= end { wi = 29; break 'outer; }
               *p |= 8; p = p.offset(prime_ * 6 + 3);
               break 'label30
              }
              if p >= end { wi = 30; break 'outer; }
              *p |= 1; p = p.offset(prime_ * 4 + 1);
              break 'label31
             }
             if p >= end { wi = 31; break 'outer; }
             *p |= 64; p = p.offset(prime_ * 2 + 1);
             wi = 24
            }
        }
        32...39 => { // 30 * x + 17
            loop {
             'label39: loop {
              'label38: loop {
               'label37: loop {
                'label36: loop {
                 'label35: loop {
                  'label34: loop {
                   'label33: loop {
                    'label32: loop {
                     match wi {
                      32 => break 'label32,
                      33 => break 'label33,
                      34 => break 'label34,
                      35 => break 'label35,
                      36 => break 'label36,
                      37 => break 'label37,
                      38 => break 'label38,
                      _ => break 'label39,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 32;
                        *p.offset(prime_ * 2 + 1) |= 64;
                        *p.offset(prime_ * 6 + 4) |= 1;
                        *p.offset(prime_ * 12 + 7) |= 8;
                        *p.offset(prime_ * 14 + 8) |= 16;
                        *p.offset(prime_ * 20 + 11) |= 128;
                        *p.offset(prime_ * 24 + 14) |= 2;
                        *p.offset(prime_ * 26 + 15) |= 4;

                        p = p.offset(prime_ * 30 + 17)
                    }
                    if p >= end { wi = 32; break 'outer; }
                    *p |= 32; p = p.offset(prime_ * 2 + 1);
                    break 'label33
                   }
                   if p >= end { wi = 33; break 'outer; }
                   *p |= 64; p = p.offset(prime_ * 4 + 3);
                   break 'label34
                  }
                  if p >= end { wi = 34; break 'outer; }
                  *p |= 1; p = p.offset(prime_ * 6 + 3);
                  break 'label35
                 }
                 if p >= end { wi = 35; break 'outer; }
                 *p |= 8; p = p.offset(prime_ * 2 + 1);
                 break 'label36
                }
                if p >= end { wi = 36; break 'outer; }
                *p |= 16; p = p.offset(prime_ * 6 + 3);
                break 'label37
               }
               if p >= end { wi = 37; break 'outer; }
               *p |= 128; p = p.offset(prime_ * 4 + 3);
               break 'label38
              }
              if p >= end { wi = 38; break 'outer; }
              *p |= 2; p = p.offset(prime_ * 2 + 1);
              break 'label39
             }
             if p >= end { wi = 39; break 'outer; }
             *p |= 4; p = p.offset(prime_ * 4 + 2);
             wi = 32
            }
        }
        40...47 => { // 30 * x + 19
            loop {
             'label47: loop {
              'label46: loop {
               'label45: loop {
                'label44: loop {
                 'label43: loop {
                  'label42: loop {
                   'label41: loop {
                    'label40: loop {
                     match wi {
                      40 => break 'label40,
                      41 => break 'label41,
                      42 => break 'label42,
                      43 => break 'label43,
                      44 => break 'label44,
                      45 => break 'label45,
                      46 => break 'label46,
                      _ => break 'label47,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 1;
                        *p.offset(prime_ * 4 + 2) |= 16;
                        *p.offset(prime_ * 10 + 6) |= 4;
                        *p.offset(prime_ * 12 + 7) |= 32;
                        *p.offset(prime_ * 18 + 11) |= 8;
                        *p.offset(prime_ * 22 + 13) |= 128;
                        *p.offset(prime_ * 24 + 15) |= 2;
                        *p.offset(prime_ * 28 + 17) |= 64;

                        p = p.offset(prime_ * 30 + 19)
                    }
                    if p >= end { wi = 40; break 'outer; }
                    *p |= 1; p = p.offset(prime_ * 4 + 2);
                    break 'label41
                   }
                   if p >= end { wi = 41; break 'outer; }
                   *p |= 16; p = p.offset(prime_ * 6 + 4);
                   break 'label42
                  }
                  if p >= end { wi = 42; break 'outer; }
                  *p |= 4; p = p.offset(prime_ * 2 + 1);
                  break 'label43
                 }
                 if p >= end { wi = 43; break 'outer; }
                 *p |= 32; p = p.offset(prime_ * 6 + 4);
                 break 'label44
                }
                if p >= end { wi = 44; break 'outer; }
                *p |= 8; p = p.offset(prime_ * 4 + 2);
                break 'label45
               }
               if p >= end { wi = 45; break 'outer; }
               *p |= 128; p = p.offset(prime_ * 2 + 2);
               break 'label46
              }
              if p >= end { wi = 46; break 'outer; }
              *p |= 2; p = p.offset(prime_ * 4 + 2);
              break 'label47
             }
             if p >= end { wi = 47; break 'outer; }
             *p |= 64; p = p.offset(prime_ * 2 + 2);
             wi = 40
            }
        }
        48...55 => { // 30 * x + 23
            loop {
             'label55: loop {
              'label54: loop {
               'label53: loop {
                'label52: loop {
                 'label51: loop {
                  'label50: loop {
                   'label49: loop {
                    'label48: loop {
                     match wi {
                      48 => break 'label48,
                      49 => break 'label49,
                      50 => break 'label50,
                      51 => break 'label51,
                      52 => break 'label52,
                      53 => break 'label53,
                      54 => break 'label54,
                      _ => break 'label55,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 32;
                        *p.offset(prime_ * 6 + 5) |= 2;
                        *p.offset(prime_ * 8 + 6) |= 64;
                        *p.offset(prime_ * 14 + 11) |= 4;
                        *p.offset(prime_ * 18 + 14) |= 8;
                        *p.offset(prime_ * 20 + 15) |= 128;
                        *p.offset(prime_ * 24 + 19) |= 1;
                        *p.offset(prime_ * 26 + 20) |= 16;

                        p = p.offset(prime_ * 30 + 23)
                    }
                    if p >= end { wi = 48; break 'outer; }
                    *p |= 32; p = p.offset(prime_ * 6 + 5);
                    break 'label49
                   }
                   if p >= end { wi = 49; break 'outer; }
                   *p |= 2; p = p.offset(prime_ * 2 + 1);
                   break 'label50
                  }
                  if p >= end { wi = 50; break 'outer; }
                  *p |= 64; p = p.offset(prime_ * 6 + 5);
                  break 'label51
                 }
                 if p >= end { wi = 51; break 'outer; }
                 *p |= 4; p = p.offset(prime_ * 4 + 3);
                 break 'label52
                }
                if p >= end { wi = 52; break 'outer; }
                *p |= 8; p = p.offset(prime_ * 2 + 1);
                break 'label53
               }
               if p >= end { wi = 53; break 'outer; }
               *p |= 128; p = p.offset(prime_ * 4 + 4);
               break 'label54
              }
              if p >= end { wi = 54; break 'outer; }
              *p |= 1; p = p.offset(prime_ * 2 + 1);
              break 'label55
             }
             if p >= end { wi = 55; break 'outer; }
             *p |= 16; p = p.offset(prime_ * 4 + 3);
             wi = 48
            }
        }
        56...63 => { // 30 * x + 29
            loop {
             'label63: loop {
              'label62: loop {
               'label61: loop {
                'label60: loop {
                 'label59: loop {
                  'label58: loop {
                   'label57: loop {
                    'label56: loop {
                     match wi {
                      56 => break 'label56,
                      57 => break 'label57,
                      58 => break 'label58,
                      59 => break 'label59,
                      60 => break 'label60,
                      61 => break 'label61,
                      62 => break 'label62,
                      _ => break 'label63,
                     }
                    }
                    while p < loop_end {
                        p = ::b(p);
                        *p.offset(prime_ * 0 + 0) |= 1;
                        *p.offset(prime_ * 2 + 1) |= 128;
                        *p.offset(prime_ * 8 + 7) |= 64;
                        *p.offset(prime_ * 12 + 11) |= 32;
                        *p.offset(prime_ * 14 + 13) |= 16;
                        *p.offset(prime_ * 18 + 17) |= 8;
                        *p.offset(prime_ * 20 + 19) |= 4;
                        *p.offset(prime_ * 24 + 23) |= 2;

                        p = p.offset(prime_ * 30 + 29)
                    }
                    if p >= end { wi = 56; break 'outer; }
                    *p |= 1; p = p.offset(prime_ * 2 + 1);
                    break 'label57
                   }
                   if p >= end { wi = 57; break 'outer; }
                   *p |= 128; p = p.offset(prime_ * 6 + 6);
                   break 'label58
                  }
                  if p >= end { wi = 58; break 'outer; }
                  *p |= 64; p = p.offset(prime_ * 4 + 4);
                  break 'label59
                 }
                 if p >= end { wi = 59; break 'outer; }
                 *p |= 32; p = p.offset(prime_ * 2 + 2);
                 break 'label60
                }
                if p >= end { wi = 60; break 'outer; }
                *p |= 16; p = p.offset(prime_ * 4 + 4);
                break 'label61
               }
               if p >= end { wi = 61; break 'outer; }
               *p |= 8; p = p.offset(prime_ * 2 + 2);
               break 'label62
              }
              if p >= end { wi = 62; break 'outer; }
              *p |= 4; p = p.offset(prime_ * 4 + 4);
              break 'label63
             }
             if p >= end { wi = 63; break 'outer; }
             *p |= 2; p = p.offset(prime_ * 6 + 6);
             wi = 56
            }
        }
        _ => unreachable!("{}", wi),
    }
    }
    *si_ = (p as usize).wrapping_sub(end as usize);
    *wi_ = wi;
}