src/integer/bitwise/cl_I_mkfx.cc

   1 // mkf_extract().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_I_byte.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_integer.h"
  13 #include "cl_I.h"
  14 #include "cl_DS.h"
  15
  16 const cl_I mkf_extract (const cl_I& x, uintL p, uintL q)
  17 {
  18       CL_ALLOCA_STACK;
  19       var const uintD* MSDptr;
  20       var uintC len;
  21       var const uintD* LSDptr;
  22       I_to_NDS_nocopy(x, MSDptr=,len=,LSDptr=,cl_true, { return 0; } ); // NDS zu x bilden
  23       // MSDptr erhöhen und len erniedrigen, so daß len = ceiling(q/intDsize) wird:
  24       { var uintL 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 ((uintL)len - qD); // MSDptr um len-qD Digits erhöhen
  28         len = qD; // len um len-qD erniedrigen
  29       }
  30       // Platz (len Digits) für die neue UDS bereitstellen:
  31       var uintD* newMSDptr;
  32       num_stack_alloc_1((uintL)len, newMSDptr = ,); // Platz belegen
  33       {var uintL pD = p/intDsize; // floor(p/intDsize), paßt in ein uintC
  34        // Kopiere len-pD Digits aus der DS zu x heraus:
  35        var uintD* midptr = copy_loop_msp(MSDptr,newMSDptr,len-(uintC)pD);
  36        // Lösche p-intDsize*floor(p/intDsize) Bits im Digit unterhalb von midptr:
  37        {var uintL p_D = p%intDsize;
  38         if (!(p_D==0)) { lspref(midptr,0) &= minus_bit(p_D); }
  39        }
  40        // Lösche pD Digits darüber:
  41        clear_loop_msp(midptr,pD);
  42       }
  43       // Lösche intDsize*ceiling(q/intDsize)-q Bits im ersten Digit:
  44       {var uintL q_D = q%intDsize;
  45        if (!(q_D==0))
  46          { mspref(newMSDptr,0) &= (uintD)((1L<<q_D)-1); } // intDsize-q_D Bits löschen
  47       }
  48       // Jetzt enthält die UDS newMSDptr/len/.. die extrahierten Bits.
  49       return UDS_to_I(newMSDptr,len);
  50 }