4 #include "base/cl_sysdep.h"
7 #include "numtheory/cl_IF.h"
12 #include "base/cl_low.h"
16 uint32 cl_trialdivision (uint32 n, uint32 d1, uint32 d2)
18 var uintL i = cl_small_prime_table_search(d1);
19 var const uint16 * ptr = &cl_small_prime_table[i];
20 var const uint16 * ptr_limit = &cl_small_prime_table[cl_small_prime_table_search(d2+1)];
21 for ( ; ptr < ptr_limit; ptr++) {
22 var uint32 prime = *ptr;
24 r = n % prime; // or: divu_3232_3232(n,prime,,r=);