]> matita.cs.unibo.it Git - helm.git/blobdiff - helm/software/components/ng_paramodulation/paramod.ml
proof reconstruction almost OK
[helm.git] / helm / software / components / ng_paramodulation / paramod.ml
index f30a9fcb6ba34949c3ed99f1753a6e6607c820ed..f0cd03758594d92411284fd0e1945d86d21a3967 100644 (file)
@@ -1,4 +1,9 @@
-let nparamod metasenv subst context t table =
+let debug s =
+  () (* prerr_endline s *)
+;;
+
+let nparamod rdb metasenv subst context t table =
+  let nb_iter = ref 200 in
   prerr_endline "========================================";
   let module C = struct
     let metasenv = metasenv
@@ -12,48 +17,190 @@ let nparamod metasenv subst context t table =
   let module IDX = Index.Index(B) in
   let module Sup = Superposition.Superposition(B) in
   let module Utils = FoUtils.Utils(B) in
-(*
-  let test_unification _ = function
-    | Terms.Node [_; _; lhs; rhs] ->
-       prerr_endline "Unification test :";
-       prerr_endline (Pp.pp_foterm lhs);
-       prerr_endline (Pp.pp_foterm rhs);
-        FU.unification [] [] lhs rhs
-    | _ -> assert false
-  in
-  let subst,vars = test_unification [] res in
-    prerr_endline "Result :";
-    prerr_endline (Pp.pp_foterm res);
-    prerr_endline "Substitution :";
-    prerr_endline (Pp.pp_substitution subst)
-*)
-  let mk_clause maxvar t =
-    let ty = B.embed t in
-    let proof = B.embed (NCic.Rel ~-1) in
-    Utils.mk_unit_clause maxvar ty proof 
-  in
-  let clause, maxvar = mk_clause 0 t in
-  prerr_endline "Input clause :";
-  prerr_endline (Pp.pp_unit_clause clause);
-  let bag = Utils.empty_bag in
-  let active_clauses, maxvar = 
+    
+  let module OrderedPassives =
+      struct
+       type t = B.t Terms.passive_clause
+           
+       let compare = Utils.compare_passive_clauses
+      end
+  in
+  let module PassiveSet = Set.Make(OrderedPassives)
+  in
+  let add_passive_clause passives cl =
+    PassiveSet.add (Utils.mk_passive_clause cl) passives
+  in
+    (* TODO : fairness condition *)
+  let select passives =
+    if PassiveSet.is_empty passives then None
+    else let cl = PassiveSet.min_elt passives in
+      Some (snd cl,PassiveSet.remove cl passives)
+  in
+  let rec given_clause bag maxvar actives      
+      passives g_actives g_passives =
+    
+    decr nb_iter; if !nb_iter = 0 then 
+      (*(prerr_endline "Bag :"; prerr_endline (Pp.pp_bag bag);
+      prerr_endline "Active table :"; 
+       (List.iter (fun x -> prerr_endline (Pp.pp_unit_clause x))
+         (fst actives));*)
+    raise (Failure "Timeout !");
+
+
+      
+    (* superposition left, simplifications on goals *)
+      debug "infer_left step...";
+      let bag, maxvar, g_actives, g_passives =
+       match select g_passives with
+      | None -> bag, maxvar, g_actives, g_passives
+      | Some (g_current, g_passives) -> 
+         debug ("Selected goal : " ^ Pp.pp_unit_clause g_current);
+          let bag, g_current = 
+            Sup.simplify_goal maxvar (snd actives) bag g_current
+          in
+          let bag, maxvar, new_goals = 
+            Sup.infer_left bag maxvar g_current actives 
+          in
+         let new_goals = List.fold_left add_passive_clause
+           PassiveSet.empty new_goals
+         in
+          bag, maxvar, g_current::g_actives,
+         (PassiveSet.union new_goals g_passives)
+    in
+      prerr_endline
+       (Printf.sprintf "Number of active goals : %d"
+          (List.length g_actives));
+      prerr_endline
+       (Printf.sprintf "Number of passive goals : %d"
+          (PassiveSet.cardinal g_passives));
+  
+      (* forward step *)
+      (* e = select P           *
+       * e' = demod A e         *
+       * A' = demod [e'] A      *
+       * A'' = A' + e'          *
+       * e'' = fresh e'         *
+       * new = supright e'' A'' *
+       * new'= demod A'' new    *
+       * P' = P + new'          *)
+      debug "Forward infer step...";
+    let bag, maxvar, actives, passives, g_passives =
+      let rec aux_simplify passives =
+       match select passives with
+         | None -> assert false
+         | Some (current, passives) -> 
+             debug ("Selected fact : " ^ Pp.pp_unit_clause current);
+              match Sup.keep_simplified current actives bag with
+             (* match Sup.one_pass_simplification current actives bag with *)
+               | None -> aux_simplify passives
+               | Some x -> x,passives
+      in
+      let (current, bag, actives),passives = aux_simplify passives
+      in                   
+       debug ("Fact after simplification :"
+              ^ Pp.pp_unit_clause current);
+       let bag, maxvar, actives, new_clauses = 
+          Sup.infer_right bag maxvar current actives 
+       in
+  debug "Demodulating goals with actives...";
+         (* keep goals demodulated w.r.t. actives and check if solved *)
+       let bag, g_actives = 
+         List.fold_left 
+            (fun (bag,acc) c -> 
+               let bag, c = Sup.simplify_goal maxvar (snd actives) bag c in
+                bag, c::acc) 
+            (bag,[]) g_actives 
+       in
+       let ctable = IDX.index_unit_clause IDX.DT.empty current in
+       let bag, maxvar, new_goals = 
+         List.fold_left 
+           (fun (bag,m,acc) g -> 
+              let bag, m, ng = Sup.infer_left bag maxvar g
+                ([current],ctable) in
+                bag,m,ng@acc) 
+           (bag,maxvar,[]) g_actives 
+       in
+       let new_clauses = List.fold_left add_passive_clause
+         PassiveSet.empty new_clauses in
+       let new_goals = List.fold_left add_passive_clause
+         PassiveSet.empty new_goals in
+         bag, maxvar, actives,
+      PassiveSet.union new_clauses passives,
+      PassiveSet.union new_goals g_passives
+    in
+      prerr_endline
+       (Printf.sprintf "Number of actives : %d" (List.length (fst actives)));
+      prerr_endline
+       (Printf.sprintf "Number of passives : %d"
+          (PassiveSet.cardinal passives));
+      given_clause bag maxvar actives passives g_actives g_passives
+  in
+
+  let mk_clause bag maxvar (t,ty) =
+    let (proof,ty) = B.saturate t ty in
+    let c, maxvar = Utils.mk_unit_clause maxvar ty proof in
+    let bag, c = Utils.add_to_bag bag c in
+    bag, maxvar, c
+  in
+  let bag, maxvar, goal = mk_clause Terms.M.empty 0 t in
+  let g_actives = [] in
+  let g_passives = PassiveSet.singleton (Utils.mk_passive_clause goal) in
+  let passives, bag, maxvar = 
     List.fold_left
-      (fun (cl,maxvar) t -> 
-         let c, m = mk_clause maxvar t in
-         c::cl, m)
-      ([],maxvar) table 
-  in
-  let table =  
-    List.fold_left IDX.index_unit_clause IDX.DT.empty active_clauses
-  in
-  prerr_endline "Active table:";
-  List.iter (fun uc -> prerr_endline (Pp.pp_unit_clause uc)) active_clauses;
-  let bag, maxvar, _, newclauses = 
-    Sup.infer_right bag maxvar clause (active_clauses, table)
-  in
-  prerr_endline "Output clauses :";
-  List.iter (fun c -> prerr_endline (Pp.pp_unit_clause c)) newclauses;
-  prerr_endline "Proofs: ";
-  prerr_endline (Pp.pp_bag bag);
-  prerr_endline "========================================";
+      (fun (cl, bag, maxvar) t -> 
+         let bag, maxvar, c = mk_clause bag maxvar t in
+         (add_passive_clause cl c), bag, maxvar)
+      (PassiveSet.empty, bag, maxvar) table 
+  in
+  let actives = [], IDX.DT.empty in
+  try given_clause bag maxvar actives passives g_actives g_passives
+  with 
+  | Sup.Success (bag, _, (i,_,_,_)) ->
+      let l =
+        let module S  = 
+          HTopoSort.Make(struct type t=int let compare=Pervasives.compare end)
+        in
+        let module C : Set.S with type elt = int = 
+          Set.Make(struct type t=int let compare=Pervasives.compare end)
+        in
+        let all id = 
+          let rec traverse ongoal (accg,acce) i =
+            match Terms.M.find i bag with
+            | (_,_,_,Terms.Exact _) -> accg, acce
+            | (_,_,_,Terms.Step (_,i1,i2,_,_,_)) -> 
+               let accg, acce = 
+                 if ongoal then C.add i1 accg, acce
+                 else accg, C.add i1 acce
+               in
+               let acce = C.add i2 acce in
+                 traverse false (traverse ongoal (accg,acce) i1) i2    
+          in
+            traverse true (C.empty,C.empty) id
+        in
+        let esteps = 
+          S.topological_sort (C.elements (snd (all i))) 
+           (fun i -> C.elements (snd (all i)))
+        in
+        let gsteps = 
+          S.topological_sort (C.elements (fst (all i))) 
+           (fun i -> C.elements (fst (all i)))
+        in
+        let gsteps = List.rev gsteps in
+        esteps@gsteps
+      in
+      prerr_endline "YES!";
+      prerr_endline "Proof:"; 
+      List.iter (fun x -> 
+              prerr_endline (Pp.pp_unit_clause (Terms.M.find x bag))) l;
+      let proofterm = B.mk_proof bag i l in
+      prerr_endline
+       (NCicPp.ppterm ~metasenv:C.metasenv ~subst:C.subst ~context:C.context
+         proofterm); 
+      let _metasenv, _subst, _proofterm, _prooftype = 
+        NCicRefiner.typeof rdb C.metasenv C.subst C.context proofterm None
+      in
+      prerr_endline "REFINED!";
+      ()
+  | Failure _ -> prerr_endline "FAILURE";
 ;;