moewsmoews_sieve.h

Go to the documentation of this file.

 
#ifndef POMYANG_KPARENT_INC_MOEWSMOEWS_SIEVE_H_
#define POMYANG_KPARENT_INC_MOEWSMOEWS_SIEVE_H_
 
#include <stdbool.h>
#include <stdint.h>
#include <stdlib.h>
#include <sys/types.h>
 
typedef struct {
    size_t bound;
    size_t seg_len;
    size_t buf_len;
    size_t primes_len;
    uint32_t *primes;
} sieve_config_t;
 
typedef struct {
    const sieve_config_t *cfg;
    size_t seg_start;
    uint64_t *sigma_buf;
    uint64_t *q_buf;
    ssize_t last_sieve_standard;
    bool squared;
    bool *is_prime_buf;
} sieve_worker_t;
 
void moews_sieve_odd_standard(sieve_worker_t *worker, uint64_t seg_start);
 
void moews_sieve_odd_squared(sieve_worker_t *worker, uint64_t seg_start);
 
uint64_t moews_lookup_sigma_m(const sieve_worker_t *worker, uint64_t m);
 
uint64_t moews_lookup_sigma_m_squared(const sieve_worker_t *worker, uint64_t m);
 
sieve_config_t *moews_init_sieve(size_t bound, size_t seg_len);
 
void moews_destroy_sieve(sieve_config_t *cfg);
 
sieve_worker_t *moews_init_worker(const sieve_config_t *cfg);
 
void destroy_worker(sieve_worker_t *worker);
 
bool moews_check_prime(const sieve_worker_t *worker, uint64_t m);
 
size_t moews_estimate_heap_usage(size_t bound, size_t seg_len, size_t num_workers);
 
#endif  // POMYANG_KPARENT_INC_MOEWSMOEWS_SIEVE_H_