src/base/digitseq/cl_2DS.h

   1 // Digit sequence 2-adic arithmetic
   2
   3 #ifndef _CL_2DS_H
   4 #define _CL_2DS_H
   5
   6 // div2adic(a_len,a_LSDptr,b_len,b_LSDptr,dest_LSDptr);
   7 // dividiert die UDS a_LSDptr[-a_len..-1] mod 2^(intDsize*b_len)
   8 // durch die ungerade UDS b_LSDptr[-b_len..-1] mod 2^(intDsize*b_len)
   9 // (a_len >= b_len > 0) und liefert
  10 // den Quotienten q als UDS dest_LSDptr[-b_len..-1] mod 2^(intDsize*b_len) und
  11 // den "Rest" (a-b*q)/2^(intDsize*b_len) als UDS dest_LSDptr[-a_len..-b_len-1].
  12 // Falls a_len > b_len, wird b implizit als durch Nullen fortgesetzt angenommen.
  13   extern void div2adic (uintC a_len, const uintD* a_LSDptr, uintC b_len, const uintD* b_LSDptr, uintD* dest_LSDptr);
  14
  15 // div2adic(len,a_LSDptr,b_LSDptr,dest_LSDptr);
  16 // dividiert die UDS a_LSDptr[-len..-1] mod 2^(intDsize*len)
  17 // durch die ungerade UDS b_LSDptr[-len..-1] mod 2^(intDsize*len) (len>0) und
  18 // liefert den Quotienten als UDS dest_LSDptr[-len..-1] mod 2^(intDsize*len).
  19   inline void div2adic (uintC len, const uintD* a_LSDptr, const uintD* b_LSDptr, uintD* dest_LSDptr)
  20   { div2adic(len,a_LSDptr,len,b_LSDptr,dest_LSDptr); }
  21
  22 // recip2adic(len,a_LSDptr,dest_LSDptr);
  23 // bildet den Kehrwert der ungeraden UDS a_LSDptr[-len..-1] mod 2^(intDsize*len)
  24 // (len>0) und liefert sie als UDS dest_LSDptr[-len..-1] mod 2^(intDsize*len).
  25   extern void recip2adic (uintC len, const uintD* a_LSDptr, uintD* dest_LSDptr);
  26
  27 #endif /* _CL_2DS_H */