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
# 1 "src/owl/stats/owl_stats_prng.ml"
(*
 * OWL - OCaml Scientific and Engineering Computing
 * Copyright (c) 2016-2020 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 ()