X-Git-Url: https://git.xandkar.net/?a=blobdiff_plain;f=polymorphic-life%2F001%2Fsrc%2Fpolymorphic_life.ml;h=e08dc2acd8978d16bd06b0ddd15b2596c1379239;hb=a1665c9210867ffc55a4a609885d1d4cba840dbd;hp=9a9163e1c1be2f6cb0a0132730fc075f71b19c89;hpb=0d6f78335b05b14736c63234d7f15a9cf2420b1a;p=cellular-automata.git

diff --git a/polymorphic-life/001/src/polymorphic_life.ml b/polymorphic-life/001/src/polymorphic_life.ml
index 9a9163e..e08dc2a 100644
--- a/polymorphic-life/001/src/polymorphic_life.ml
+++ b/polymorphic-life/001/src/polymorphic_life.ml
@@ -2,54 +2,118 @@ open Core.Std
 
 
 module type MATRIX = sig
+  module Point : sig
+    type t = {r : int; k : int}
+  end
+
   type 'a t
 
-  val create : rows:int -> cols:int -> data:'a -> 'a t
+  val create : rs:int -> ks:int -> 'a -> 'a t
+
+  val get_neighbors : 'a t -> Point.t -> 'a list
 
-  val get : 'a t -> row:int -> col:int -> 'a
+  val map : 'a t -> f:('a -> 'b) -> 'b t
 
-  val set : 'a t -> row:int -> col:int -> data:'a -> unit
+  val mapi : 'a t -> f:(Point.t -> 'a -> 'b) -> 'b t
 
-  val map : 'a t -> f:(row:int -> col:int -> data:'a -> 'b) -> 'b t
+  val iter : 'a t -> f:(Point.t -> 'a -> unit) -> unit
 
-  val iter : 'a t -> f:(row:int -> col:int -> data:'a -> unit) -> unit
+  val print : 'a t -> to_string:('a -> string) -> unit
 end
 
 module Matrix : MATRIX = struct
+  module Point = struct
+    type t = {r : int; k : int}
+
+    let (+) p p' =
+      { r = p.r + p'.r
+      ; k = p.k + p'.k
+      }
+  end
+
+  module Direction = struct
+    type t = NW | N | NE
+           | W  |     E
+           | SW | S | SE
+
+    let all = [ NW ; N ; NE
+              ; W  ;     E
+              ; SW ; S ; SE
+              ]
+
+    let to_offset =
+      let open Point in
+      function
+      | NW -> {r = -1; k = -1}
+      | N  -> {r = -1; k =  0}
+      | NE -> {r = -1; k =  1}
+      | W  -> {r =  0; k = -1}
+      | E  -> {r =  0; k =  1}
+      | SW -> {r =  1; k = -1}
+      | S  -> {r =  1; k =  0}
+      | SE -> {r =  1; k =  1}
+  end
+
   type 'a t = 'a array array
 
-  let create ~rows ~cols ~data =
-    Array.make_matrix ~dimx:rows ~dimy:cols data
+  let create ~rs ~ks x =
+    Array.make_matrix ~dimx:rs ~dimy:ks x
 
   let iter t ~f =
     Array.iteri t ~f:(
-      fun row cols ->
-        Array.iteri cols ~f:(
-          fun col data ->
-            f ~row ~col ~data
+      fun r ks ->
+        Array.iteri ks ~f:(
+          fun k x ->
+            f {Point.r; Point.k} x
         )
     )
 
+  let print t ~to_string =
+    Array.iter t ~f:(
+      fun r ->
+        Array.iter r ~f:(fun x -> printf "%s" (to_string x));
+        print_newline ()
+    )
+
   let map t ~f =
+    Array.map t ~f:(Array.map ~f:(fun x -> f x))
+
+  let mapi t ~f =
     Array.mapi t ~f:(
-      fun row cols ->
-        Array.mapi cols ~f:(
-          fun col data ->
-            f ~row ~col ~data
+      fun r ks ->
+        Array.mapi ks ~f:(
+          fun k x ->
+            f {Point.r; Point.k} x
         )
     )
 
-  let get t ~row ~col =
-    t.(row).(col)
+  let get t {Point.r; Point.k} =
+    t.(r).(k)
 
-  let set t ~row ~col ~data =
-    t.(row).(col) <- data
+  let is_within_bounds t {Point.r; Point.k} =
+    match t with
+    | [||] -> assert false
+    | t ->
+      r >= 0 && r < Array.length t &&
+      k >= 0 && k < Array.length t.(0)
+
+  let neighborhood t point =
+    List.map Direction.all ~f:Direction.to_offset
+    |> List.map ~f:(fun offset_point -> Point.(point + offset_point))
+    |> List.filter ~f:(is_within_bounds t)
+
+  let get_neighbors t point =
+    List.map (neighborhood t point) ~f:(get t)
 end
 
 
 module type CELL = sig
   type t
 
+  val create : unit -> t
+
+  val to_string : t -> string
+
   val state : t -> int
 
   val react : t -> states:int list -> t
@@ -59,10 +123,24 @@ end
 module Conway : CELL = struct
   type t = D | A
 
-  let state = function
+  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
@@ -70,15 +148,45 @@ module Conway : CELL = struct
     | A when live_neighbors < 4 -> A
     | A when live_neighbors > 3 -> D
     | D when live_neighbors = 3 -> A
-    | t -> t
+    | A -> A
+    | D -> D
 end
 
 
+type opt = { interval : Time.Span.t
+           ; bar      : string
+           }
+
+
+let rec loop opt grid =
+  print_endline opt.bar;
+  Matrix.print grid ~to_string:Conway.to_string;
+  print_endline opt.bar;
+
+  let grid =
+    Matrix.mapi grid ~f:(fun point cell ->
+      let neighbors = Matrix.get_neighbors grid point in
+      Conway.react cell ~states:(List.map neighbors ~f:Conway.state)
+    )
+  in
+  Time.pause opt.interval;
+  loop opt grid
+
+
 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
-  )
+  Random.self_init ();
+  let rs, ks = Or_error.ok_exn Linux_ext.get_terminal_size () in
+  Matrix.create ~rs:(rs - 3) ~ks ()
+  |> Matrix.map ~f:Conway.create
+  |> loop { interval = Time.Span.of_float 0.1
+          ; bar      = String.make ks '-'
+          }
+
+
+let spec =
+  let summary = "Polymorphic Cellular Automata" in
+  let spec = Command.Spec.empty in
+  Command.basic ~summary spec main
 
 
-let () = main ()
+let () = Command.run spec