-let main () =
- let pool = Matrix.create ~rows:5 ~cols:5 ~data:() in
- Matrix.iter pool ~f:(
- fun ~row ~col ~data:() -> printf "R: %d, K: %d\n" row col
- )
+module Conway : CELL = struct
+ type t = D | A
+
+ let of_int = function
+ | 0 -> D
+ | 1 -> A
+ | _ -> assert false
+
+ let to_int = function
+ | D -> 0
+ | A -> 1
+
+ let to_string = function
+ | D -> " "
+ | A -> "o"
+
+ let create () =
+ Random.int 2 |> of_int
+
+ let state = to_int
+
+ let react t ~states =
+ let live_neighbors = List.fold_left states ~init:0 ~f:(+) in
+ match t with
+ | A when live_neighbors < 2 -> D
+ | A when live_neighbors < 4 -> A
+ | A when live_neighbors > 3 -> D
+ | D when live_neighbors = 3 -> A
+ | A -> A
+ | D -> D
+end
+
+
+let main rs ks () =
+ Random.self_init ();
+ let grid = Matrix.create ~rs ~ks ~data:() |> Matrix.map ~f:Conway.create in
+ Matrix.print grid ~to_string:Conway.to_string;
+ print_endline (String.make 80 '-');
+ let grid =
+ Matrix.mapi grid ~f:(fun ~r ~k ~data:cell ->
+ let neighbors = Matrix.get_neighbors grid ~r ~k in
+ Conway.react cell ~states:(List.map neighbors ~f:Conway.state)
+ )
+ in
+ Matrix.print grid ~to_string:Conway.to_string
+
+
+let spec =
+ let summary = "Polymorphic Cellular Automata" in
+ let spec =
+ let open Command.Spec in
+ empty
+ +> flag "-rows" (optional_with_default 5 int) ~doc:"Height"
+ +> flag "-cols" (optional_with_default 5 int) ~doc:"Width"
+ in
+ Command.basic ~summary spec main