let make ~pos ~neg = 
  let pos = Array.of_stream pos
  and neg = Array.of_stream neg in
  Array.sort (Fn.flip compare) pos ;
  Array.sort (Fn.flip compare) neg ;
  let sorted_elements = 
    Stream.merge
      ~cmp:(Fn.flip compare)
      (Array.to_stream pos |! Stream.map ~f:(fun x -> x, `pos))
      (Array.to_stream neg |! Stream.map ~f:(fun x -> x, `neg))
  and initial = {
    tp = 0 ;
    tn = Array.length neg ;
    fp = 0 ;
    fn = Array.length pos
  } 
  in
    Stream.append
      (Stream.singleton (Float.infinity, initial))
      (Stream.unfold
         initial
         (fun accu -> 
           if Stream.is_empty sorted_elements then None
           else match Stream.next sorted_elements with
             | Some (x, `pos) -> 
               let next = { accu with tp = accu.tp + 1 ; fn = accu.fn - 1 } in
               Some ((x, next), next)
             | Some (x, `neg) ->
               let next = { accu with fp = accu.fp + 1 ; tn = accu.tn - 1 } in
               Some ((x, next), next)
             | None -> None))