sasa.ml

(* Time-stamp: <modified the 13/03/2019 (at 18:04) by Erwan Jahier> *)

(* XXX Je pourrais utiliser Lwt pour rendre step non-bloquant, ce qui
   permettrait d'accelerer la simu sur les machines qui ont plusieurs
   coeurs

  step : action -> (string -> value) Lwt.t ;
*)

open Algo
       
let (update_env_with_init : Env.t -> Process.t list -> Algo.neighbor list list -> Env.t) =
  fun e pl neighbors ->
    let (aux: Env.t -> Process.t -> Algo.neighbor list -> Env.t) =
      fun e p nl ->
        List.fold_left
          (fun e (n,_t) -> Env.set e p.pid n (p.init nl n))
          e
          p.variables
    in
    List.fold_left2 aux e pl neighbors
    
let (get_neighors: Process.t -> Topology.neighbor list) =
  fun p ->
    let id = p.Process.pid in
    let idl = try Hashtbl.find Topology.node_succ id with Not_found -> [] in
    List.map
      (fun id ->
         let node =
           try Hashtbl.find Topology.node_info id with Not_found -> assert false
         in
         let algo_id = Filename.chop_suffix node.file ".cmxs" in
          {
            Topology.n_id = node.id;
            Topology.n_vars = Algo.get_vars algo_id;
          }
      )
      idl

let (dump_process: Process.t * Topology.neighbor list -> unit) =
  fun (p,nl) ->
    let pvars = StringOf.algo_vars p.variables in
    let neighbors = List.map StringOf.topology_neighbor nl in
    Printf.printf "process %s\n\tvars:%s\n\tneighors: \n\t\t%s\n" p.pid  pvars
      (String.concat "\n\t\t" neighbors)

let (update_env: Env.t -> Process.t * Algo.local_env -> Env.t) =
  fun e (p, lenv) -> 
    List.fold_left
      (fun e (n,_) -> Env.set e p.pid n (lenv n))
      e
      p.variables

open Process
open SasArg
    
let (to_algo_neighbor: Env.t -> Topology.neighbor -> Algo.neighbor) =
  fun e n ->
    { 
      lenv= Env.get e n.Topology.n_id;
      n_vars = n.Topology.n_vars
    }

exception Silent of int


let (print_step : int -> int -> Env.t -> Process.t list -> string -> string -> unit) =
  fun n i e pl activate_val enable_val ->
    if SasArg.args.rif then (
      if SasArg.args.demon = Demon.Custom then (
        (* in custom mode, to be able to talk with lurette, this should not be 
           printed on stdout
        *)
        Printf.eprintf "\n#step %s\n" (string_of_int (n-i+1)) ;
        Printf.eprintf "%s #outs " activate_val; flush stderr
      ) else (
        Printf.printf "\n#step %s\n" (string_of_int (n-i+1)) ;
        Printf.printf "%s #outs " activate_val; flush stdout
      );
      Printf.printf "%s %s\n" (StringOf.env_rif e pl)  enable_val;
      flush stdout
    )
    else (
      Printf.eprintf "step %s: %s %s\n"
        (string_of_int (n-i+1)) (StringOf.env e pl) activate_val;
      flush stderr
    )
          

let rec (simu: int -> int -> Process.t list -> string -> 
         (Process.t * Topology.neighbor list) list -> Env.t -> unit) =
  fun n i pl activate_val pl_n e ->
    let custom = args.demon = Demon.Custom in
    let all = List.fold_left
        (fun acc (p,nl) ->
           let nl4algo = List.map (to_algo_neighbor e) nl in
           let lenv = Env.get e p.pid in
           let al = p.enable nl4algo lenv in
           let al = List.map (fun a -> p,nl,a) al in
           al::acc)
        [] pl_n
    in
    assert (List.length pl = List.length all);
    let all = List.rev all in 
    let enab_ll =
      if not custom then [] else
        List.map2
          (fun p al ->
             let al = List.map (fun (_,_,a) -> a) al in
             List.map (fun a_static -> List.mem a_static al) p.actions)
          pl
          all
    in
    let enable_val =
      String.concat " " (List.map (fun b -> if b then "t" else "f")
                           (List.flatten enab_ll))
    in
    if (List.flatten all = []) then (
      print_step n i e pl activate_val enable_val; 
      raise (Silent (n-i+1))
    );
    print_step n i e pl activate_val enable_val;
    let next_activate_val, pnal =
      Demon.f (SasArg.args.verbose > 1) args.demon pl all enab_ll
    in

    (* Do the steps *)
    let lenv_list =
      List.map (fun (p,nl,a) ->
          let nl4algo = List.map (to_algo_neighbor e) nl in
          let lenv = Env.get e p.pid in
          p, p.step nl4algo lenv a)
        pnal
    in
    (* update the env *)
    let ne = List.fold_left update_env e lenv_list in

    if i > 0 then simu n (i-1) pl next_activate_val pl_n ne else (
      if SasArg.args.rif then (
        print_string "q\n"; flush stdout
      ))

let () =
  ( try SasArg.parse Sys.argv;
	 with
	     Failure(e) ->
	       output_string stdout e;
	       flush stdout ;
	       exit 2
	   | e ->
	     output_string stdout (Printexc.to_string e);
	     flush stdout;
	     exit 2
  );
  let dot_file = SasArg.args.topo in
  let nl = Topology.read dot_file in
  let nstrl = List.map (fun n -> n.Topology.id) nl in
  let nstr = String.concat "," nstrl in
  try
    Algo.verbose_level := SasArg.args.verbose;
    Random.self_init();
    if !Algo.verbose_level > 0 then Printf.printf "nodes: %s\nedges:\n" nstr;
    let e = Env.init () in
    let pl = List.map (Process.make (args.demon=Custom)) nl in
    let neighors = List.map get_neighors pl in
    let algo_neighors = List.map (List.map (to_algo_neighbor e)) neighors in
    let e = update_env_with_init e pl algo_neighors in
    let pl_n = List.combine pl neighors in
    if !Algo.verbose_level > 0 then List.iter dump_process pl_n;
    let n = SasArg.args.length in
    if SasArg.args.rif then (
      Printf.printf "#inputs %s\n"
        (if SasArg.args.demon = Demon.Custom then (
            let f p = List.map
                (fun a -> "\""^p.pid ^(if a="" then "" else "_")^a^ "\":bool")
                p.actions
            in
            String.concat " " (List.flatten (List.map f pl))
      ) else "");
      Printf.printf "#outputs %s\n" (StringOf.env_rif_decl pl);
      flush stdout
    );
    if SasArg.args.ifi then (
      List.iter
        (fun p -> List.iter
            (fun a -> ignore (RifRead.bool (args.verbose > 1) p a)) p.actions)
        pl;
      Printf.eprintf "Ignoring the first vectors of sasa inputs\n"; flush stderr;
    );
    simu n n pl "" pl_n e
  with
  | Dynlink.Error e -> Printf.printf "Error: %s\n" (Dynlink.error_message e)
  | Failure msg -> Printf.printf "Error: %s\n" msg
  | Silent i ->
    let str = if SasArg.args.rif then "#" else "" in
     Printf.printf "%sThis algo is silent after %i steps\n" str i ;
    flush stdout;
    if SasArg.args.rif && args.demon = Custom then (
      print_string "q\n"; flush stdout
    )