3 * Copyright 2008,2009 Free Software Foundation, Inc.
5 * This file is part of GNU Radio
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)
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.
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.
21 #include <gcell/spu/gc_random.h>
23 static int last_val = 0;
25 # define M 259200 // values from Numerical Recipes in C, 1988
32 last_val = ((unsigned int) seed) % M;
36 * Return a uniformly distributed value in the range [0, 1.0)
37 * (Linear congruential generator. YMMV. Caveat emptor.)
40 gc_uniform_deviate(void)
42 last_val = (last_val * A + C) % M;
43 return (float) last_val * (1.0f / (float) M);