84 lines
1.6 KiB
C
84 lines
1.6 KiB
C
/*
|
|
|
|
random.c - Glue code for BreezySLAM pseudorandom number generator
|
|
|
|
Copyright (C) 2014 Simon D. Levy
|
|
|
|
This code is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU Lesser General Public License as
|
|
published by the Free Software Foundation, either version 3 of the
|
|
License, or (at your option) any later version.
|
|
|
|
This code is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU Lesser General Public License
|
|
along with this code. If not, see <http:#www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#include "ziggurat.h"
|
|
#include "random.h"
|
|
|
|
#include <stdlib.h>
|
|
#include <string.h>
|
|
|
|
typedef struct random_t
|
|
{
|
|
float fn[128];
|
|
uint32_t kn[128];
|
|
float wn[128];
|
|
uint32_t seed;
|
|
|
|
} random_t;
|
|
|
|
size_t random_size(void)
|
|
{
|
|
return sizeof(random_t);
|
|
}
|
|
|
|
|
|
void * random_new(int seed)
|
|
{
|
|
random_t * r = (random_t *)malloc(sizeof(random_t));
|
|
|
|
random_init(r, seed);
|
|
|
|
return r;
|
|
|
|
}
|
|
|
|
|
|
void random_init(void * v, int seed)
|
|
{
|
|
random_t * r = (random_t *)v;
|
|
|
|
r->seed = seed;
|
|
|
|
r4_nor_setup (r->kn, r->fn, r->wn );
|
|
}
|
|
|
|
|
|
double random_normal(void * v, double mu, double sigma)
|
|
{
|
|
random_t * r = (random_t *)v;
|
|
|
|
return mu + sigma * r4_nor ( &r->seed, r->kn, r->fn, r->wn );
|
|
}
|
|
|
|
void random_free(void * v)
|
|
{
|
|
free(v);
|
|
}
|
|
|
|
void * random_copy(void * v)
|
|
{
|
|
random_t * r = (random_t *)malloc(sizeof(random_t));
|
|
|
|
memcpy(r, v, sizeof(random_t));
|
|
|
|
return r;
|
|
}
|
|
|