12 #include "cln/integer.h"
18 cl_boolean ldb_extract_test (const cl_I& x, uintC p, uintC q)
19 { var const uintD* MSDptr;
21 var const uintD* LSDptr;
22 I_to_NDS_nocopy(x, MSDptr=,len=,LSDptr=,cl_true, { return cl_false; } ); // NDS zu x bilden
23 // MSDptr erhöhen und len erniedrigen, so daß len = ceiling(q/intDsize) wird:
24 { var uintC qD = ceiling(q,intDsize); // ceiling(q/intDsize)
25 // wegen q<=l ist qD = ceiling(q/intDsize) <= ceiling((l+1)/intDsize) = len, also
26 // paßt qD ebenso wie len in ein uintC.
27 MSDptr = MSDptr mspop (len - qD); // MSDptr um len-qD Digits erhöhen
28 len = qD; // len um len-qD erniedrigen
30 // LSDptr und len um floor(p/intDsize) erniedrigen:
31 { var uintC pD = p/intDsize; // floor(p/intDsize)
32 LSDptr = LSDptr lspop pD;
35 // Jetzt enthält MSDptr/len/LSDptr genau die maßgeblichen Digits.
36 if (len==0) return cl_false; // len=0 -> keine Bits abzutesten
37 q = ((q-1)%intDsize); // q := intDsize - (intDsize*ceiling(q/intDsize) - q) - 1
38 p = p%intDsize; // p := p - intDsize*floor(p/intDsize)
39 // Jetzt ist 0 <= q < intDsize, 0 <= p < intDsize.
40 // Vom ersten Digit müssen die vorderen intDsize-1-q Bits unberücksichtigt bleiben.
41 // Ein AND 2^(q+1)-1 erreicht dies.
42 // Vom letzten Digit müssen die hinteren p Bits unberücksichtigt bleiben.
43 // Ein AND -2^p erreicht dies.
45 // 1 Digit maßgeblich, wird von beiden Seiten angeschnitten:
46 // Ein AND 2^(q+1)-2^p erreicht dies.
47 if (!(((uintD)(bitm(q+1)-bit(p)) & mspref(MSDptr,0)) == 0))
51 // mindestens 2 Digits. Teste erst die Randdigits, dann die inneren:
52 if (!(((msprefnext(MSDptr) & (uintD)(bitm(q+1)-1)) == 0) &&
53 ((lsprefnext(LSDptr) & (uintD)(minus_bit(p))) == 0)
56 len--; // die beiden Randdigits sind jetzt abgezogen.
57 if (DS_test_loop(MSDptr,len,LSDptr)) { return cl_true; } else { return cl_false; }