diff --git a/external/random/jl_random.c b/external/random/jl_random.c index 5ac06cd59259b..3dc841c1f514c 100644 --- a/external/random/jl_random.c +++ b/external/random/jl_random.c @@ -2,40 +2,5 @@ random numbers */ -#include "jl_random.h" #include "dsfmt-2.1/dSFMT.c" #include "randmtzig.c" - -#define RANDN_RESET -99999999 - -static double randn_bm_next = RANDN_RESET; - -void randn_bm_reset(void); -double randn_bm(void); - -void randn_bm_reset(void) -{ - randn_bm_next = RANDN_RESET; -} - -// Box Muller -double randn_bm(void) -{ - double s, vre, vim, ure, uim; - - if (randn_bm_next != RANDN_RESET) { - s = randn_bm_next; - randn_bm_next = RANDN_RESET; - return s; - } - do { - ure = dsfmt_gv_genrand_close1_open2(); - uim = dsfmt_gv_genrand_close1_open2(); - vre = 2*ure - 3; - vim = 2*uim - 3; - s = vre*vre + vim*vim; - } while (s >= 1); - s = sqrt(-2*log(s)/s); - randn_bm_next = s * vre; - return s * vim; -} diff --git a/external/random/jl_random.h b/external/random/jl_random.h deleted file mode 100644 index 59ac597797c68..0000000000000 --- a/external/random/jl_random.h +++ /dev/null @@ -1,18 +0,0 @@ -#ifndef __JL_RANDOM_H_ -#define __JL_RANDOM_H_ - -#include -#include -#include -#include -#include - -double rand_double(); -float rand_float(); -double randn(); -void randomize(); -uint32_t genrand_int32(); -void randomseed32(uint32_t s); -void randomseed64(uint64_t s); - -#endif diff --git a/external/random/randmtzig.c b/external/random/randmtzig.c index e698492daa8cd..b11004d9e07a0 100644 --- a/external/random/randmtzig.c +++ b/external/random/randmtzig.c @@ -49,8 +49,6 @@ #include #include -#include "jl_random.h" - typedef int randmtzig_idx_type; typedef signed char randmtzig_int8_t; typedef unsigned char randmtzig_uint8_t;