File: eratosthenes.ml

package info (click to toggle)
js-of-ocaml 5.9.1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 32,020 kB
  • sloc: ml: 91,250; javascript: 57,289; ansic: 315; makefile: 271; lisp: 23; sh: 6; perl: 4
file content (200 lines) | stat: -rw-r--r-- 5,544 bytes parent folder | download
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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
(*
Copyright (c) 2015, KC Sivaramakrishnan <sk826@cl.cam.ac.uk>

Permission to use, copy, modify, and/or distribute this software for any
purpose with or without fee is hereby granted, provided that the above
copyright notice and this permission notice appear in all copies.

THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
*)

open Effect
(** Message-passing parallel prime number generation using Sieve of Eratosthenes **)

open Effect.Deep

(* A message is either a [Stop] signal or a [Candidate] prime number *)
type message =
  | Stop
  | Candidate of int

let string_of_msg = function
  | Stop -> "Stop"
  | Candidate i -> Printf.sprintf "%d" i

type pid = int
(** Process primitives **)

type _ Effect.t += Spawn : (pid -> unit) -> pid Effect.t

let spawn p = perform (Spawn p)

type _ Effect.t += Yield : unit Effect.t

let yield () = perform Yield

(** Communication primitives **)
type _ Effect.t += Send : pid * message -> unit Effect.t

let send pid data =
  perform (Send (pid, data));
  yield ()

type _ Effect.t += Recv : pid -> message option Effect.t

let rec recv pid =
  match perform (Recv pid) with
  | Some m -> m
  | None ->
      yield ();
      recv pid

(** A mailbox is indexed by process ids (PIDs), each process has its own message queue **)
module Mailbox = struct
  module Make (Ord : Map.OrderedType) = struct
    include Map.Make (Ord)

    let empty = empty

    let lookup key mb = try Some (find key mb) with Not_found -> None

    let pop key mb =
      ( (match lookup key mb with
        | Some msg_q -> if Queue.is_empty msg_q then None else Some (Queue.pop msg_q)
        | None -> None)
      , mb )

    let push key msg mb =
      match lookup key mb with
      | Some msg_q ->
          Queue.push msg msg_q;
          mb
      | None ->
          let msg_q = Queue.create () in
          Queue.push msg msg_q;
          add key msg_q mb
  end
end

(** Communication handler **)
let mailbox f =
  let module Mailbox = Mailbox.Make (struct
    type t = pid

    let compare = compare
  end) in
  let mailbox = ref Mailbox.empty in
  let lookup pid =
    let msg, mb = Mailbox.pop pid !mailbox in
    mailbox := mb;
    msg
  in
  try_with
    f
    ()
    { effc =
        (fun (type a) (e : a Effect.t) ->
          match e with
          | Send (pid, msg) ->
              Some
                (fun (k : (a, _) continuation) ->
                  mailbox := Mailbox.push pid msg !mailbox;
                  continue k ())
          | Recv who ->
              Some
                (fun k ->
                  let msg = lookup who in
                  continue k msg)
          | _ -> None)
    }

(** Process handler
    Slightly modified version of sched.ml **)
let run main () =
  let run_q = Queue.create () in
  let enqueue k = Queue.push k run_q in
  let dequeue () = if Queue.is_empty run_q then () else (Queue.pop run_q) () in
  let pid = ref (-1) in
  let rec spawn f =
    pid := 1 + !pid;
    match_with
      f
      !pid
      { retc = (fun () -> dequeue ())
      ; exnc = (fun e -> raise e)
      ; effc =
          (fun (type a) (e : a Effect.t) ->
            match e with
            | Yield ->
                Some
                  (fun (k : (a, _) continuation) ->
                    enqueue (fun () -> continue k ());
                    dequeue ())
            | Spawn p ->
                Some
                  (fun k ->
                    enqueue (fun () -> continue k !pid);
                    spawn p)
            | _ -> None)
      }
  in
  spawn main

let fromSome = function
  | Some x -> x
  | _ -> failwith "Attempt to unwrap None"

(** The prime number generator **)
let rec generator n : pid -> unit =
 fun _ ->
  let first = spawn sieve in
  (* Spawn first sieve *)
  Printf.printf "Primes in [2..%d]: " n;
  for i = 2 to n do
    send first (Candidate i) (* Send candidate prime to first sieve *)
  done;
  send first Stop;
  (* Stop the pipeline *)
  Printf.printf "\n"

and sieve : pid -> unit =
 fun mypid ->
  match recv mypid with
  | Candidate myprime ->
      let _ = Printf.printf "%d " myprime in
      let succ = ref None in
      let rec loop () =
        let msg = recv mypid in
        match msg with
        | Candidate prime when prime mod myprime <> 0 ->
            let succ_pid =
              if !succ = None
              then (
                let pid = spawn sieve in
                (* Create a successor process *)
                succ := Some pid;
                pid)
              else fromSome !succ
            in
            send succ_pid (Candidate prime);
            (* Send candidate prime to successor process *)
            loop ()
        | Stop when !succ <> None -> send (fromSome !succ) Stop (* Forward stop command *)
        | Stop -> ()
        | _ -> loop ()
      in
      loop ()
  | _ -> ()

(* Run application *)

let%expect_test _ =
  let _ = mailbox (run (generator 202)) in
  [%expect
    {| Primes in [2..202]: 2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199 |}]