diff options
Diffstat (limited to 'coreutils/factor.c')
| -rw-r--r-- | coreutils/factor.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/coreutils/factor.c b/coreutils/factor.c index 8ca3de71a..4e521dbc0 100644 --- a/coreutils/factor.c +++ b/coreutils/factor.c | |||
| @@ -54,7 +54,7 @@ static inline half_t isqrt(wide_t N) | |||
| 54 | // if (N == 0) | 54 | // if (N == 0) |
| 55 | // return 0; | 55 | // return 0; |
| 56 | 56 | ||
| 57 | /* First approximation x > sqrt(N) - half as many bits: | 57 | /* First approximation of x+1 > sqrt(N) - all-ones, half as many bits: |
| 58 | * 1xxxxx -> 111 (six bits to three) | 58 | * 1xxxxx -> 111 (six bits to three) |
| 59 | * 01xxxx -> 111 | 59 | * 01xxxx -> 111 |
| 60 | * 001xxx -> 011 | 60 | * 001xxx -> 011 |
| @@ -125,7 +125,7 @@ static NOINLINE void factorize(wide_t N) | |||
| 125 | if (factor >= max_factor) | 125 | if (factor >= max_factor) |
| 126 | break; | 126 | break; |
| 127 | factor += 2; | 127 | factor += 2; |
| 128 | /* Rudimentary wheel sieving: skip multiples of 3: | 128 | /* Rudimentary wheel sieving: skip multiples of 3, 5 and 7: |
| 129 | * Every third odd number is divisible by three and thus isn't a prime: | 129 | * Every third odd number is divisible by three and thus isn't a prime: |
| 130 | * 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47... | 130 | * 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47... |
| 131 | * ^ ^ ^ ^ ^ ^ ^ _ ^ ^ _ ^ ^ ^ ^ | 131 | * ^ ^ ^ ^ ^ ^ ^ _ ^ ^ _ ^ ^ ^ ^ |