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

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2008 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  714025       // values from Numerical Recipes in C, 1988
  26 #define A    4096
  27 #define C  150889
  28
  29 void
  30 gc_set_seed(int seed)
  31 {
  32   last_val = ((unsigned int) seed) % M;
  33 }
  34
  35 float
  36 gc_uniform_deviate(void)
  37 {
  38   last_val = (last_val * A + C) % M;
  39   return (float) last_val / (float) M;
  40 }