let int_of_string s =
  try
    int_of_string s
  with e ->
      failwith (Printf.sprintf  "int_of_string: error with [%s]" (String.escaped s))

let on_file filename = 
  let ic = open_in filename in
  let count = ref 0 in
  let mins = ref [||] in
  let maxs = ref [||] in
  (try
      let line = input_line ic in
      incr count;
      let line = String2.split_simplify line ' ' in
      mins := Array.create (List.length line) 0.;
      maxs := Array.create (List.length line) 0.;
      let mins = !mins in
      let maxs = !maxs in
      List2.iteri (fun i x ->
          let x = float_of_string x in
          mins.(i) <- x;
          maxs.(i) <- x) line;
      while true do
        let line = input_line ic in
        incr count;
        let line = String2.split_simplify line ' ' in
        List2.iteri (fun i x ->
            let x = float_of_string x in
            mins.(i) <- min x mins.(i);
            maxs.(i) <- max x maxs.(i)) line;
      done
      
    with End_of_file -> ());
  !count, !mins, !maxs

    
let value_per_time x_of_y filename int_table element =
  Printf.printf "%s...\n%!" filename;
  let oc = open_out (Printf.sprintf  "%s.txt" filename) in
  let maximum = ref int_table.(0) in
  let nrounds = Array.length int_table in  
  for round = 0 to nrounds - 1 do
    maximum := max !maximum int_table.(round);
    output_string oc (Printf.sprintf "%d %s\n" round
      (x_of_y int_table.(round)))
  done;
  close_out oc;

  let oc = open_out (Printf.sprintf "%s.plot" filename) in
  output_string oc (Printf.sprintf "
set out '%s.eps'
set terminal postscript eps enhanced \"Helvetica\" 14 
set out '%s.eps'
set yrange [0:%s+1]
set xrange [0:%d]
set ylabel 'Number of %s'
set xlabel 'Time (minutes)'

set data style lines
      
plot \"%s.txt\" using ($1):($2) title '%s' 
  " filename filename (x_of_y !maximum) nrounds element filename element);
  close_out oc



  let npoints = 1000
  
let cdf_of_list x_of_y filename int_list options xlabel ylabel end_plot =
  Printf.printf "%s...\n%!" filename;
  match int_list with
    [] -> Printf.printf "  not enough data\n"
  | init :: _ ->
      let int_list = List.sort compare int_list in
      let nints = List.length int_list in
      
      let oc = open_out (Printf.sprintf  "%s.txt" filename) in
      let maximum = ref init in
      List2.iteri (fun i x ->
          maximum := max !maximum  x;
          output_string oc (Printf.sprintf "%d %s\n" i 
            (x_of_y x))
      ) int_list;
      close_out oc;

      let every = 
        if nints > npoints then 
          Printf.sprintf " every %d " (nints / npoints) else "" in
      
    let oc = open_out (Printf.sprintf "%s.plot" filename) in
  output_string oc (Printf.sprintf "
set out '%s.eps'
set terminal postscript eps enhanced \"Helvetica\" 14 
set out '%s.eps'
set yrange [1:%s+1]
set xrange [0:%d]
set ylabel '%s'
set xlabel '%s'
set data style lines
      
%s

plot \"%s.txt\" %s using ($1):($2) %s
  " filename filename (x_of_y !maximum) (nints+2) ylabel xlabel options filename every end_plot);
  close_out oc

      
let for_every every filename =
  let ic = open_in filename in
  let nlines = ref 0 in
  let line = ref "" in
  try
    while true do
      line := input_line ic;
      incr nlines
    done
  with End_of_file ->
      close_in ic;
      if !nlines > 0 then        
        let x = (!nlines - 1) mod every in
        if x > 0 then
        let oc = open_out_gen [Open_append] 0o644 filename in       
        for i = x to every-1 do 
          Printf.printf "adding line\n";
          Printf.fprintf oc "%s\n" !line
        done;
        close_out oc