15 #error "intDsize too small for trialdivision!"
18 uint32 cl_trialdivision (const cl_I& n, uint32 d1, uint32 d2)
20 var uintL i = cl_small_prime_table_search(d1);
21 var const uint16 * ptr = &cl_small_prime_table[i];
22 var const uint16 * ptr_limit = &cl_small_prime_table[cl_small_prime_table_search(d2+1)];
25 var const uintD* n_MSDptr;
27 I_to_NDS_nocopy(n, n_MSDptr=,n_len=,,cl_false,);
28 if (mspref(n_MSDptr,0)==0) { msshrink(n_MSDptr); n_len--; }
29 // Make room for a quotient.
31 num_stack_alloc(n_len,q_MSDptr=,);
33 for ( ; ptr < ptr_limit; ptr++) {
34 var uint32 prime = *ptr;
35 var uintD r = divucopy_loop_msp(prime,n_MSDptr,q_MSDptr,n_len);