helm/software/components/extlib/hTopoSort.ml

   1 (* Copyright (C) 2006, HELM Team.
   2  *
   3  * This file is part of HELM, an Hypertextual, Electronic
   4  * Library of Mathematics, developed at the Computer Science
   5  * Department, University of Bologna, Italy.
   6  *
   7  * HELM is free software; you can redistribute it and/or
   8  * modify it under the terms of the GNU General Public License
   9  * as published by the Free Software Foundation; either version 2
  10  * of the License, or (at your option) any later version.
  11  *
  12  * HELM is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with HELM; if not, write to the Free Software
  19  * Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  20  * MA  02111-1307, USA.
  21  *
  22  * For details, see the HELM World-Wide-Web page,
  23  * http://helm.cs.unibo.it/
  24  *)
  25
  26 module Make (OT:Map.OrderedType) =
  27     struct
  28
  29       module M = Map.Make(OT);;
  30
  31       let topological_sort l find_deps =
  32         let find_deps m x =
  33           let deps = find_deps x in
  34           M.add x deps m
  35         in
  36         (* build the partial order relation *)
  37         let m = List.fold_left (fun m i -> find_deps m i) M.empty l in
  38         let m = (* keep only deps inside l *)
  39           List.fold_left
  40             (fun m' i ->
  41               M.add i (List.filter (fun x -> List.mem x l) (M.find i m)) m')
  42             M.empty l
  43         in
  44         let m = M.map (fun x -> Some x) m in
  45         (* utils *)
  46         let keys m = M.fold (fun i _ acc -> i::acc) m [] in
  47         let split l m = List.filter (fun i -> M.find i m = Some []) l in
  48         let purge l m =
  49           M.mapi
  50             (fun k v -> if List.mem k l then None else
  51                match v with
  52                | None -> None
  53                | Some ll -> Some (List.filter (fun i -> not (List.mem i l)) ll))
  54             m
  55         in
  56         let rec aux m res =
  57             let keys = keys m in
  58             let ok = split keys m in
  59             let m = purge ok m in
  60             let res = ok @ res in
  61             if ok = [] then res else aux m res
  62         in
  63         let rc = List.rev (aux m []) in
  64         rc
  65       ;;
  66
  67     end
  68 ;;
  69