Source file xoshiro256plusplus_bindings.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
module LowLevel = struct
  type t = (int64, Bigarray.int64_elt, Bigarray.c_layout) Bigarray.Array1.t

  let of_int64_array a =
    if Array.length a <> 4 then
      invalid_arg "Xoshiro256plusplus.LowLevel.of_int64_array";
    Bigarray.(Array1.of_array int64 c_layout) a

  let to_int64_array ba =
    Array.init 4 (Bigarray.Array1.get ba)

  external next : t -> int64 = "caml_x256pp_next"
  external jump : t -> unit = "caml_x256pp_jump"
  external long_jump : t -> unit = "caml_x256pp_long_jump"
end

include MakeRandom.Full30Init64(struct
    type state

    external bits : state -> int = "caml_x256pp_bits"
    external new_state : unit -> state = "caml_x256pp_new_state"
    external assign : state -> state -> unit = "caml_x256pp_assign"

    let init_size = 4
    external init : state -> int64 array -> unit = "caml_x256pp_init"

    let default_seed = 135801055
  end)

(* Note to self [from Niols]: One can see that the basic functions (eg. [val
   bits: unit -> int]) in fact hide an OCaml default state which will be passed
   to C every time, which means that, in order to run a [unit -> int], we have
   to copy and unbox a state. It is then tempting to duplicate the code, to
   write C function using a C state, and to bind these functions to OCaml. I
   have tried, and it really increases the complexity of the code for...
   nothing. like at all. like it is even possibly slower than before. *)