prime_sieve.c File Reference

#include <errno.h>
#include <limits.h>
#include <stdbool.h>
#include <stdio.h>
#include <stdlib.h>
#include <derive-c/macros/iterators.h>
#include <derive-c/structures/vector/template.h>

Include dependency graph for prime_sieve.c:

Go to the source code of this file.

Macros
#define	SELF   sieve_vec
#define	T   bool

Functions
size_t	sqrt_size_t (size_t n)
void	display (sieve_vec const *sieve)
void	compute (sieve_vec *sieve)
int	main ()

Macro Definition Documentation

SELF

#define SELF   sieve_vec

Definition at line 14 of file prime_sieve.c.

T

#define T   bool

Definition at line 15 of file prime_sieve.c.

Function Documentation

compute()

void compute

(

sieve_vec *

sieve

)

Examples

complex/prime_sieve.c.

Definition at line 52 of file prime_sieve.c.

                               {
    size_t size = sieve_vec_size(sieve);
    size_t sqrt = sqrt_size_t(size);
 
    for (size_t factor = 2; factor < sqrt; factor++) {
        for (size_t index = factor * 2; index < size; index += factor) {
            *sieve_vec_write(sieve, index) = true;
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

display()

void display

(

sieve_vec const *

sieve

)

Examples

complex/prime_sieve.c.

Definition at line 40 of file prime_sieve.c.

                                     {
    sieve_vec_iter_const iter = sieve_vec_get_iter_const(sieve);
    sieve_vec_iter_const_next(&iter); // skip 0
    sieve_vec_iter_const_next(&iter); // skip 1
 
    ITER_ENUMERATE_LOOP(sieve_vec_iter_const, iter, bool const*, is_not_prime, size_t, index) {
        if (!*is_not_prime) {
            printf("%zu is prime\n", sieve_vec_iter_const_position(&iter) - 1);
        }
    }
}

Here is the caller graph for this function:

main()

int main

(

)

Definition at line 63 of file prime_sieve.c.

           {
    size_t up_to = 23;
    printf("Listing primes up to: %zu\n", up_to);
 
    sieve_vec values = sieve_vec_new_with_defaults(up_to, false);
    compute(&values);
    display(&values);
    sieve_vec_delete(&values);
 
    return 0;
}

Here is the call graph for this function:

sqrt_size_t()

size_t sqrt_size_t

(

size_t

n

)

Examples

complex/prime_sieve.c.

Definition at line 18 of file prime_sieve.c.

                             {
    if (n == 0 || n == 1) {
        return n;
    }
    size_t left = 1, right = n, mid, result = 0;
    while (left <= right) {
        mid = left + (right - left) / 2;
        size_t square = mid * mid;
 
        if (square == n) {
            return mid;
        } else if (square < n) {
            result = mid; // Store the last valid result
            left = mid + 1;
        } else {
            right = mid - 1;
        }
    }
 
    return result;
}

Here is the caller graph for this function: