git.gag.com Git - debian/gnuradio/blob - gcell/lib/runtime/spu/gc_random.c

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2008,2009 Free Software Foundation, Inc.
   4  *
   5  * This file is part of GNU Radio
   6  *
   7  * GNU Radio is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU General Public License as published by
   9  * the Free Software Foundation; either version 3, or (at your option)
  10  * any later version.
  11  *
  12  * GNU Radio is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License along
  18  * with this program; if not, write to the Free Software Foundation, Inc.,
  19  * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
  20  */
  21 #include <gcell/spu/gc_random.h>
  22
  23 static int last_val = 0;
  24
  25 # define M  259200      // values from Numerical Recipes in C, 1988
  26 # define A    7141
  27 # define C   54773
  28
  29 void
  30 gc_set_seed(int seed)
  31 {
  32   last_val = ((unsigned int) seed) % M;
  33 }
  34
  35 /*
  36  * Return a uniformly distributed value in the range [0, 1.0)
  37  * (Linear congruential generator. YMMV. Caveat emptor.)
  38  */
  39 float
  40 gc_uniform_deviate(void)
  41 {
  42   last_val = (last_val * A + C) % M;
  43   return (float) last_val * (1.0f / (float) M);
  44 }