Source file `owl_ppl.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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
# 1 "src/owl/ppl/owl_ppl.ml"
(*
 * OWL - OCaml Scientific and Engineering Computing
 * Copyright (c) 2016-2018 Liang Wang <liang.wang@cl.cam.ac.uk>
 *)

(** Experimental: Probabilistic Programming *)

open Owl_types

module Make (A : Stats_Dist) = struct

  include Owl_distribution.Make (A)

  module Lazy = Owl_lazy.Make (A)


  (* graph manipulation *)

  (* Draw random variables from different distributions. *)
(*
  let uniform ~a ~b =
    let draw_samples args =
      let a = Lazy.value_to_arr args.(0) in
      let b = Lazy.value_to_arr args.(1) in
      (*let s = A.shape (to_arr args.(2)) in*)
      let n = 0 in
      let t = Uniform.make ~a ~b in
      Lazy.arr_to_value (Uniform.sample t n)
    in
    let shape_holder = Lazy.var_arr ~name:"" in (* FIXME *)
    map ~name:"uniform" draw_samples [|a;b;shape_holder|]

  let gaussian ~mu ~sigma =
    let draw_samples args =
      let mu = Lazy.value_to_arr args.(0) in
      let sigma = Lazy.value_to_arr args.(1) in
      (*let s = A.shape (to_arr args.(2)) in*)
      let n = 0 in
      let t = Gaussian.make ~mu ~sigma in
      Lazy.arr_to_value (Gaussian.sample t n)
    in
    let shape_holder = variable () in (* FIXME *)
    map ~name:"gaussian" draw_samples [|mu;sigma;shape_holder|]


  let sample x s =
    invalidate x;
    Array.iter (fun n ->
      assign_arr n (A.empty s)
    ) (get_by_name x "variable");
    eval x;
    x

  let infer x = ()
*)

  (* Mathematical operators *)


end