← Back to branch summary

~ubuntu-branches/ubuntu/trusty/gnuradio/trusty

« back to all changes in this revision

Viewing changes to gcell/lib/runtime/spu/gc_random.c

  • Committer: Bazaar Package Importer
  • Author(s): Kamal Mostafa
  • Date: 2010-03-13 07:46:01 UTC
  • mfrom: (2.1.2 sid)
  • Revision ID: james.westby@ubuntu.com-20100313074601-zjsa893a87bozyh7
Tags: 3.2.2.dfsg-1ubuntu1
https://launchpad.net/bugs/260406
* Fix build for Ubuntu lucid (LP: #260406)
  - add binary package dep for libusrp0, libusrp2-0: adduser
  - debian/rules clean: remove pre-built Qt moc files

Show diffs side-by-side

added added

removed removed

Lines of Context:
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
}