Source file `owl_stats_prng.ml`

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
# 1 "src/owl/stats/owl_stats_prng.ml"
(*
 * OWL - an OCaml numerical library for scientific computing
 * Copyright (c) 2016-2018 Liang Wang <liang.wang@cl.cam.ac.uk>
 *)


type state


external sfmt_seed : int -> unit = "owl_sfmt_seed"


external rand_int : unit -> int = "owl_sfmt_rand_int"


external ziggurat_init : unit -> unit = "owl_ziggurat_init"


external rand_exp : unit -> float = "owl_std_exponential_rvs"


external rand_gaussian : unit -> float = "owl_std_gaussian_rvs"


let self_init () =
  Owl_base_stats_prng.self_init ();
  let seed = Random.int 65535 in
  sfmt_seed seed;
  ziggurat_init ()


let init seed =
  Owl_base_stats_prng.init seed;
  let seed = Random.int 65535 in
  sfmt_seed seed;
  ziggurat_init ()