]> matita.cs.unibo.it Git - helm.git/blobdiff - helm/software/components/ng_paramodulation/paramod.ml
Removed debug printing raising Failure
[helm.git] / helm / software / components / ng_paramodulation / paramod.ml
index a33f1121dd3ceed53364778945541580144cfd8e..ec098a61597142d0b603674c27d8cac95487dca4 100644 (file)
@@ -1,77 +1,19 @@
-(*let nparamod metasenv subst context t table =
-  prerr_endline "========================================";
-  let module C = struct
-    let metasenv = metasenv
-    let subst = subst
-    let context = context
-  end
-  in
-  let module B = NCicBlob.NCicBlob(C) in
-  let module Pp = Pp.Pp (B) in
-  let module FU = FoUnif.Founif(B) in
-  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)
-*)
-(*
+(* LOOPING : COL057-1.ma *)
 
-  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 = 
-    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 "========================================";
-;;
-*)
 let debug s =
-  prerr_endline s
+   () (* prerr_endline s *)
 ;;
 
-let nparamod metasenv subst context t table =
-  prerr_endline "========================================";
+let nparamod rdb metasenv subst context t table =
+  let max_nb_iter = 999999999 in
+  let amount_of_time = 300.0 in
   let module C = struct
     let metasenv = metasenv
     let subst = subst
     let context = context
   end
   in
+  let nb_iter = ref 0 in
   let module B = NCicBlob.NCicBlob(C) in
   let module Pp = Pp.Pp (B) in
   let module FU = FoUnif.Founif(B) in
@@ -91,95 +33,130 @@ let nparamod metasenv subst context t table =
   let add_passive_clause passives cl =
     PassiveSet.add (Utils.mk_passive_clause cl) passives
   in
+  let timeout = Unix.gettimeofday () +. amount_of_time in 
+
     (* TODO : fairness condition *)
-  let select passives =
-    if PassiveSet.is_empty passives then None
+  let select passives g_passives =
+    if PassiveSet.is_empty passives then begin
+      assert (not (PassiveSet.is_empty g_passives));
+      let g_cl = PassiveSet.min_elt g_passives in
+       (true,snd g_cl,passives,PassiveSet.remove g_cl g_passives)
+    end
     else let cl = PassiveSet.min_elt passives in
-      Some (snd cl,PassiveSet.remove cl passives)
+      if PassiveSet.is_empty g_passives then
+       (false,snd cl,PassiveSet.remove cl passives,g_passives)
+      else
+       let g_cl = PassiveSet.min_elt g_passives in
+         if (fst cl <= fst g_cl) then
+           (false,snd cl,PassiveSet.remove cl passives,g_passives)
+         else
+           (true,snd g_cl,passives,PassiveSet.remove g_cl g_passives)
   in
-  let rec given_clause bag maxvar actives passives g_actives g_passives =
+
+  let backward_infer_step bag maxvar actives passives g_actives g_passives g_current =
+    (* superposition left, simplifications on goals *)
+      debug "infer_left step...";
+    debug ("Selected goal : " ^ Pp.pp_unit_clause g_current);
+    let bag, g_current = 
+      Sup.simplify_goal maxvar (snd actives) bag g_current
+    in
+      debug "Simplified goal";
+      let bag, maxvar, new_goals = 
+        Sup.infer_left bag maxvar g_current actives 
+      in
+       debug "Performed infer_left step";
+       let new_goals = List.fold_left add_passive_clause
+         PassiveSet.empty new_goals
+       in
+          bag, maxvar, actives, passives, g_current::g_actives,
+    (PassiveSet.union new_goals g_passives)
+  in
+
+  let forward_infer_step bag maxvar actives passives g_actives
+                         g_passives current =
+    (* forward step *)
     
-    (* keep goals demodulated w.r.t. actives and check if solved *)
-    (* we may move this at the end of infer_right and simplify with
-     * just new_clauses *) 
-    let bag, g_actives = 
-      List.fold_left 
-        (fun (bag,acc) c -> 
-           let bag, c = Sup.backward_simplify maxvar (snd actives) bag c in
-             bag, c::acc) 
-        (bag,[]) g_actives 
+    (* 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...";
+    debug "Selected and simplified";
+    (* 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 m 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, g_actives,
+    PassiveSet.union new_goals g_passives
+  in
+
+  let rec given_clause bag maxvar actives passives g_actives g_passives =
+    (* 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)); *)
+    incr nb_iter; if !nb_iter = max_nb_iter then       
+      raise (Failure "No iterations left !");
+    if Unix.gettimeofday () > timeout then
+      raise (Failure "Timeout !");
+
+
+    let rec aux_select passives g_passives =
+      let backward,current,passives,g_passives = select passives g_passives in
+       if backward then
+         backward_infer_step bag maxvar actives passives
+                             g_actives g_passives current
+       else
+         (* debug ("Selected fact : " ^ Pp.pp_unit_clause current); *)
+          match Sup.keep_simplified current actives bag maxvar with
+       (*  match Sup.one_pass_simplification current actives bag maxvar with*)
+             | None -> aux_select passives g_passives
+             | Some x -> let (current, bag, actives) = x in
+                 forward_infer_step bag maxvar actives passives
+                                    g_actives g_passives current
     in
-      
-      (* backward step : superposition left, simplifications on goals *)
-      debug "Backward infer 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.backward_simplify 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)
+
+    let bag,maxvar,actives,passives,g_actives,g_passives =      
+      aux_select passives g_passives
     in
-      prerr_endline
+      assert (PassiveSet.cardinal g_passives < 2);
+      debug
        (Printf.sprintf "Number of active goals : %d"
           (List.length g_actives));
-      prerr_endline
+      debug
        (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 =
-      match select passives with
-      | None -> bag, maxvar, actives, passives, g_passives
-      | Some (current, passives) -> 
-         debug ("Selected fact : " ^ Pp.pp_unit_clause current);
-          match Sup.forward_simplify (snd actives) bag current with
-          | None -> debug ("Discarded fact");
-             bag, maxvar, actives, passives, g_passives
-          | Some (bag, current) ->
-             debug ("Fact after simplification :" ^ Pp.pp_unit_clause current);
-              let bag, maxvar, actives, new_clauses = 
-                Sup.infer_right bag maxvar current 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
+      debug
        (Printf.sprintf "Number of actives : %d" (List.length (fst actives)));
-      prerr_endline
+      debug
        (Printf.sprintf "Number of passives : %d"
           (PassiveSet.cardinal passives));
       given_clause bag maxvar actives passives g_actives g_passives
@@ -203,8 +180,53 @@ let nparamod metasenv subst context t table =
   in
   let actives = [], IDX.DT.empty in
   try given_clause bag maxvar actives passives g_actives g_passives
-  with Sup.Success (bag, _, mp) ->
-    prerr_endline "YES!";
-    prerr_endline "Meeting point :"; prerr_endline (Pp.pp_unit_clause mp);
-    (* prerr_endline "Bag :"; prerr_endline (Pp.pp_bag bag) *)
+  with 
+  | Sup.Success (bag, _, (i,_,_,_)) ->
+      let l =
+       let rec traverse ongoal (accg,acce) i =
+         match Terms.M.find i bag with
+           | (id,_,_,Terms.Exact _) ->
+               if ongoal then [i],acce else
+                 if (List.mem i acce) then accg,acce else accg,acce@[i]
+           | (_,_,_,Terms.Step (_,i1,i2,_,_,_)) ->
+               if (not ongoal) && (List.mem i acce) then accg,acce
+               else
+                  let accg,acce = 
+                   traverse false (traverse ongoal (accg,acce) i1) i2
+                 in
+                   if ongoal then i::accg,acce else accg,i::acce
+        in
+       let gsteps,esteps = traverse true ([],[]) i in
+         (List.rev esteps)@gsteps
+      in
+       prerr_endline (Printf.sprintf "Found proof in %d iterations, %fs"
+                        !nb_iter
+                        (Unix.gettimeofday() -. timeout +. amount_of_time));
+      (* prerr_endline "Proof:"; 
+      List.iter (fun x -> prerr_endline (string_of_int x);
+              prerr_endline (Pp.pp_unit_clause (Terms.M.find x bag))) l;*)
+      let proofterm = B.mk_proof bag i l in
+       prerr_endline (Printf.sprintf "Got proof term, %fs"
+                        (Unix.gettimeofday() -. timeout +. amount_of_time));
+      let metasenv, proofterm = 
+        let rec aux k metasenv = function
+          | NCic.Meta _ as t -> metasenv, t
+          | NCic.Implicit _ -> 
+              let metasenv,i,_,_=NCicMetaSubst.mk_meta metasenv context `Term in
+              metasenv, NCic.Meta (i,(k,NCic.Irl (List.length context)))
+          | t -> NCicUntrusted.map_term_fold_a 
+                  (fun _ k -> k+1) k aux metasenv t
+        in
+         aux 0 metasenv proofterm
+      in
+      let metasenv, subst, proofterm, _prooftype = 
+        NCicRefiner.typeof 
+          (rdb#set_coerc_db NCicCoercion.empty_db) 
+          metasenv subst context proofterm None
+      in
+      proofterm, metasenv, subst
+  | Failure s -> 
+      prerr_endline s;
+      prerr_endline
+      (Printf.sprintf "FAILURE in %d iterations" !nb_iter); assert false
 ;;