[BACK]Return to rand.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp

Diff for /OpenXM_contrib/gmp/Attic/rand.c between version 1.1 and 1.1.1.2

version 1.1, 2000/09/09 14:12:16 version 1.1.1.2, 2003/08/25 16:06:00
Line 1 
Line 1 
 /* gmp_randinit (state, algorithm, ...) -- Initialize a random state.  /* gmp_randinit (state, algorithm, ...) -- Initialize a random state.
   
 Copyright (C) 1999, 2000  Free Software Foundation, Inc.  Copyright 1999, 2000, 2001 Free Software Foundation, Inc.
   
 This file is part of the GNU MP Library.  This file is part of the GNU MP Library.
   
Line 19  along with the GNU MP Library; see the file COPYING.LI
Line 19  along with the GNU MP Library; see the file COPYING.LI
 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,  the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 MA 02111-1307, USA. */  MA 02111-1307, USA. */
   
   #include "config.h"
   
 #include <stdio.h> /* for NULL */  #include <stdio.h> /* for NULL */
 #if __STDC__  
 # include <stdarg.h>  #if HAVE_STDARG
   #include <stdarg.h>
 #else  #else
 # include <varargs.h>  #include <varargs.h>
 #endif  #endif
   
 #include "gmp.h"  #include "gmp.h"
 #include "gmp-impl.h"  #include "gmp-impl.h"
   
 /* Array of CL-schemes, ordered in increasing order of the first  
    member (the 'm2exp' value).  The end of the array is indicated with  
    an entry containing all zeros.  */  
   
 /* All multipliers are in the range 0.01*m and 0.99*m, and are  
 congruent to 5 (mod 8).  
 They all pass the spectral test with Vt >= 2^(30/t) and merit >= 1.  
 (Up to and including 196 bits, merit is >= 3.)  */  
   
 struct __gmp_rand_lc_scheme_struct  
 {  
   unsigned long int m2exp;      /* Modulus is 2 ^ m2exp. */  
   char *astr;                   /* Multiplier in string form. */  
   unsigned long int c;          /* Adder. */  
 };  
   
 struct __gmp_rand_lc_scheme_struct __gmp_rand_lc_scheme[] =  
 {  
   {32, "43840821",           1},  
   {33, "85943917",           1},  
   {34, "171799469",          1},  
   {35, "343825285",          1},  
   {36, "687285701",          1},  
   {37, "1374564613",         1},  
   {38, "2749193437",         1},  
   {39, "5497652029",         1},  
   {40, "10995212661",        1},  
   {56, "47988680294711517",  1},  
   {64, "13469374875402548381", 1},  
   {100, "203786806069096950756900463357", 1},  
   {128, "96573135900076068624591706046897650309", 1},  
   {156, "43051576988660538262511726153887323360449035333", 1},  
   {196, "1611627857640767981443524165616850972435303571524033586421", 1},  
   {200, "491824250216153841876046962368396460896019632211283945747141", 1},  
   {256, "79336254595106925775099152154558630917988041692672147726148065355845551082677", 1},  
   {0, NULL, 0}                  /* End of array. */  
 };  
   
 void  void
 #if __STDC__  #if HAVE_STDARG
 gmp_randinit (gmp_randstate_t rstate,  gmp_randinit (gmp_randstate_t rstate,
               gmp_randalg_t alg,                gmp_randalg_t alg,
               ...)                ...)
Line 78  gmp_randinit (va_alist)
Line 43  gmp_randinit (va_alist)
 #endif  #endif
 {  {
   va_list ap;    va_list ap;
 #if __STDC__  #if HAVE_STDARG
     va_start (ap, alg);
 #else  #else
   __gmp_randstate_struct *rstate;    __gmp_randstate_struct *rstate;
   gmp_randalg_t alg;    gmp_randalg_t alg;
 #endif  
   
 #if __STDC__  
   va_start (ap, alg);  
 #else  
   va_start (ap);    va_start (ap);
   
   rstate = va_arg (ap, __gmp_randstate_struct *);    rstate = va_arg (ap, __gmp_randstate_struct *);
   alg = va_arg (ap, gmp_randalg_t);    alg = va_arg (ap, gmp_randalg_t);
 #endif  #endif
   
   switch (alg)    switch (alg) {
     {    case GMP_RAND_ALG_LC:
     case GMP_RAND_ALG_LC:       /* Linear congruential.  */      if (! gmp_randinit_lc_2exp_size (rstate, va_arg (ap, unsigned long)))
       {        gmp_errno |= GMP_ERROR_INVALID_ARGUMENT;
         unsigned long int size;      break;
         struct __gmp_rand_lc_scheme_struct *sp;    default:
         mpz_t a;      gmp_errno |= GMP_ERROR_UNSUPPORTED_ARGUMENT;
       break;
     }
     va_end (ap);
   }
   
         size = va_arg (ap, unsigned long int);  
   
         /* Pick a scheme.  */  
         for (sp = __gmp_rand_lc_scheme; sp->m2exp != 0; sp++)  
           if (sp->m2exp / 2 >= size)  
             break;  
   
         if (sp->m2exp == 0)     /* Nothing big enough found.  */  
           {  
             gmp_errno |= GMP_ERROR_INVALID_ARGUMENT;  
             return;  
           }  
   
         /* Install scheme.  */  
         mpz_init_set_str (a, sp->astr, 0);  
         gmp_randinit_lc_2exp (rstate, a, sp->c, sp->m2exp);  
         mpz_clear (a);  
         break;  
       }  
   
 #if 0  #if 0
     case GMP_RAND_ALG_BBS:      /* Blum, Blum, and Shub. */      case GMP_RAND_ALG_BBS:      /* Blum, Blum, and Shub. */
       {        {
         mpz_t p, q;          mpz_t p, q;
         mpz_t ztmp;          mpz_t ztmp;
   
Line 134  gmp_randinit (va_alist)
Line 80  gmp_randinit (va_alist)
         /* These are around 128 bits. */          /* These are around 128 bits. */
         mpz_init_set_str (p, "148028650191182616877187862194899201391", 10);          mpz_init_set_str (p, "148028650191182616877187862194899201391", 10);
         mpz_init_set_str (q, "315270837425234199477225845240496832591", 10);          mpz_init_set_str (q, "315270837425234199477225845240496832591", 10);
   
         /* Allocate algorithm specific data. */          /* Allocate algorithm specific data. */
         rstate->data.bbs = (__gmp_rand_data_bbs *)          rstate->data.bbs = (__gmp_rand_data_bbs *)
           (*_mp_allocate_func) (sizeof (__gmp_rand_data_bbs));            (*__gmp_allocate_func) (sizeof (__gmp_rand_data_bbs));
   
         mpz_init (rstate->data.bbs->bi); /* The Blum integer. */          mpz_init (rstate->data.bbs->bi); /* The Blum integer. */
         mpz_mul (rstate->data.bbs->bi, p, q);          mpz_mul (rstate->data.bbs->bi, p, q);
Line 162  gmp_randinit (va_alist)
Line 108  gmp_randinit (va_alist)
         break;          break;
       }        }
 #endif /* 0 */  #endif /* 0 */
   
     default:                    /* Bad choice. */  
       gmp_errno |= GMP_ERROR_UNSUPPORTED_ARGUMENT;  
     }  
   
   va_end (ap);  
 }  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>