src/integer/random/cl_I_random.cc

   1 // random_I().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_integer.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_random_impl.h"
  13 #include "cl_I.h"
  14 #include "cl_DS.h"
  15
  16 const cl_I random_I (cl_random_state& randomstate, const cl_I& n)
  17 {
  18         CL_ALLOCA_STACK;
  19         var const uintD* n_MSDptr;
  20         var uintC n_len;
  21         var const uintD* n_LSDptr;
  22         I_to_NDS_nocopy(n, n_MSDptr=,n_len=,n_LSDptr=,cl_false,); // Digit sequence >0 zu n
  23         var uintD* MSDptr;
  24         var uintC len = n_len + ceiling(16,intDsize); // 16 Bits mehr
  25         // neue UDS mit len Zufallsdigits bilden:
  26         num_stack_alloc(len,MSDptr=,);
  27         random_UDS(randomstate,MSDptr,len);
  28         // und durch n dividieren:
  29         var DS q;
  30         var DS r;
  31         UDS_divide(MSDptr,len,MSDptr mspop len, n_MSDptr,n_len,n_LSDptr, &q,&r);
  32         // Rest in Integer umwandeln:
  33         return NUDS_to_I(r.MSDptr,r.len);
  34 }