2 ||M|| This file is part of HELM, an Hypertextual, Electronic
3 ||A|| Library of Mathematics, developed at the Computer Science
4 ||T|| Department, University of Bologna, Italy.
6 ||T|| HELM is free software; you can redistribute it and/or
7 ||A|| modify it under the terms of the GNU General Public License
8 \ / version 2 or (at your option) any later version.
9 \ / This software is distributed as is, NO WARRANTY.
10 V_______________________________________________________________ *)
12 (* $Id: orderings.ml 9869 2009-06-11 22:52:38Z denes $ *)
14 let debug s = prerr_endline s ;;
17 module Paramod (B : Terms.Blob) = struct
18 exception Failure of string * B.t Terms.bag * int * int
19 type bag = B.t Terms.bag * int
21 module FU = FoUnif.Founif(B)
22 module IDX = Index.Index(B)
23 module Sup = Superposition.Superposition(B)
24 module Utils = FoUtils.Utils(B)
25 module WeightOrderedPassives =
27 type t = B.t Terms.passive_clause
28 let compare = Utils.compare_passive_clauses_weight
31 module AgeOrderedPassives =
33 type t = B.t Terms.passive_clause
34 let compare = Utils.compare_passive_clauses_age
37 module WeightPassiveSet = Set.Make(WeightOrderedPassives)
38 module AgePassiveSet = Set.Make(AgeOrderedPassives)
40 let add_passive_clause ?(no_weight=false) (passives_w,passives_a) cl =
41 let cl = if no_weight then (0,cl)
42 else Utils.mk_passive_clause cl in
43 WeightPassiveSet.add cl passives_w, AgePassiveSet.add cl passives_a
46 let remove_passive_clause (passives_w,passives_a) cl =
47 let passives_w = WeightPassiveSet.remove cl passives_w in
48 let passives_a = AgePassiveSet.remove cl passives_a in
52 let add_passive_clauses ?(no_weight=false)
53 (passives_w,passives_a) new_clauses =
54 let new_clauses_w,new_clauses_a =
55 List.fold_left (add_passive_clause ~no_weight)
56 (WeightPassiveSet.empty,AgePassiveSet.empty) new_clauses
58 (WeightPassiveSet.union new_clauses_w passives_w,
59 AgePassiveSet.union new_clauses_a passives_a)
62 let is_passive_set_empty (passives_w,passives_a) =
63 if (WeightPassiveSet.is_empty passives_w) then begin
64 assert (AgePassiveSet.is_empty passives_a); true
66 assert (not (AgePassiveSet.is_empty passives_a)); false
70 let passive_set_cardinal (passives_w,_) = WeightPassiveSet.cardinal passives_w
72 let passive_empty_set =
73 (WeightPassiveSet.empty,AgePassiveSet.empty)
76 let pick_min_passive ~use_age (passives_w,passives_a) =
77 if use_age then AgePassiveSet.min_elt passives_a
78 else WeightPassiveSet.min_elt passives_w
81 let mk_clause bag maxvar (t,ty) =
82 let (proof,ty) = B.saturate t ty in
83 let c, maxvar = Utils.mk_unit_clause maxvar ty proof in
84 let bag, c = Utils.add_to_bag bag c in
88 let mk_passive (bag,maxvar) = mk_clause bag maxvar;;
89 let mk_goal (bag,maxvar) = mk_clause bag maxvar;;
91 (* TODO : global age over facts and goals (without comparing weights) *)
92 let select ~use_age passives g_passives =
93 if is_passive_set_empty passives then begin
94 assert (not (is_passive_set_empty g_passives));
95 let g_cl = pick_min_passive ~use_age:use_age g_passives in
96 (true,snd g_cl,passives,remove_passive_clause g_passives g_cl)
98 else let cl = pick_min_passive ~use_age:use_age passives in
99 if is_passive_set_empty g_passives then
100 (false,snd cl,remove_passive_clause passives cl,g_passives)
102 let g_cl = pick_min_passive ~use_age:use_age g_passives in
103 if (fst cl <= fst g_cl) then
104 (false,snd cl,remove_passive_clause passives cl,g_passives)
106 (true,snd g_cl,passives,remove_passive_clause g_passives g_cl)
109 let backward_infer_step bag maxvar actives passives
110 g_actives g_passives g_current =
111 (* superposition left, simplifications on goals *)
112 debug "infer_left step...";
113 let bag, maxvar, new_goals =
114 Sup.infer_left bag maxvar g_current actives
116 debug "Performed infer_left step";
117 bag, maxvar, actives, passives, g_current::g_actives,
118 (add_passive_clauses g_passives new_goals)
121 let forward_infer_step bag maxvar actives passives g_actives
127 * A' = demod [e'] A *
130 * new = supright e'' A'' *
131 * new'= demod A'' new *
133 debug "Forward infer step...";
134 let bag, maxvar, actives, new_clauses =
135 Sup.infer_right bag maxvar current actives
137 debug "Demodulating goals with actives...";
138 (* keep goals demodulated w.r.t. actives and check if solved *)
143 Sup.simplify_goal ~no_demod:false maxvar (snd actives) bag acc c
146 | Some (bag,c) -> bag,c::acc)
149 let ctable = IDX.index_unit_clause IDX.DT.empty current in
150 let bag, maxvar, new_goals =
152 (fun (bag,m,acc) g ->
153 let bag, m, ng = Sup.infer_left bag m g
154 ([current],ctable) in
156 (bag,maxvar,[]) g_actives
158 bag, maxvar, actives,
159 add_passive_clauses passives new_clauses, g_actives,
160 add_passive_clauses g_passives new_goals
163 let rec given_clause ~noinfer
164 bag maxvar iterno max_steps timeout
165 actives passives g_actives g_passives
167 let iterno = iterno + 1 in
168 if iterno = max_steps then
169 raise (Failure ("No iterations left !",bag,maxvar,iterno));
170 (* timeout check: gettimeofday called only if timeout set *)
171 if timeout <> None &&
173 | None -> assert false
174 | Some timeout -> Unix.gettimeofday () > timeout) then
176 prerr_endline("Last chance: "^string_of_float (Unix.gettimeofday()));
177 given_clause ~noinfer:true bag maxvar iterno max_steps None
178 actives passives g_actives g_passives
179 (*raise (Failure ("Timeout !",bag,maxvar,iterno))*)
183 let use_age = iterno mod 10 = 0 in
185 let rec aux_select bag passives g_passives =
187 is_passive_set_empty passives &&
188 is_passive_set_empty g_passives then
191 ("Last chance: all is indexed " ^ string_of_float
192 (Unix.gettimeofday()));
193 let maxgoals = 100 in
197 ignore(Sup.simplify_goal ~no_demod:true
198 maxvar (snd actives) bag acc x)
203 raise (Failure (("Last chance: failed over " ^
204 string_of_int maxgoals^ " goal " ^
205 string_of_float (Unix.gettimeofday())),bag,maxvar,0));
207 let backward,current,passives,g_passives =
208 select ~use_age:false passives g_passives
212 if noinfer then Some (bag,current)
215 ~no_demod:false maxvar (snd actives) bag g_actives current
217 | None -> aux_select bag passives g_passives
218 | Some (bag,g_current) ->
220 let g_actives = g_current :: g_actives in
221 bag,maxvar,actives,passives,g_actives,g_passives
223 backward_infer_step bag maxvar actives passives
224 g_actives g_passives g_current
226 (* debug ("Selected fact : " ^ Pp.pp_unit_clause current); *)
228 if noinfer then bag, Some (current,actives)
229 else Sup.keep_simplified current actives bag maxvar with
230 (* match Sup.one_pass_simplification current actives bag maxvar with*)
231 | bag,None -> aux_select bag passives g_passives
232 | bag,Some (current,actives) ->
235 current::fst actives,
236 IDX.index_unit_clause (snd actives) current
238 bag,maxvar,actives,passives,g_actives,g_passives
240 forward_infer_step bag maxvar actives passives
241 g_actives g_passives current
244 (* prerr_endline "Bag :"; prerr_endline (Pp.pp_bag bag);
245 prerr_endline "Active table :";
246 (List.iter (fun x -> prerr_endline (Pp.pp_unit_clause x))
249 let bag,maxvar,actives,passives,g_actives,g_passives =
250 aux_select bag passives g_passives
253 (Printf.sprintf "Number of active goals : %d"
254 (List.length g_actives));
256 (Printf.sprintf "Number of passive goals : %d"
257 (passive_set_cardinal g_passives));
259 (Printf.sprintf "Number of actives : %d" (List.length (fst actives)));
261 (Printf.sprintf "Number of passives : %d"
262 (passive_set_cardinal passives));
263 given_clause ~noinfer
264 bag maxvar iterno max_steps timeout
265 actives passives g_actives g_passives
268 let paramod ~max_steps ?timeout (bag,maxvar) ~g_passives ~passives =
269 let initial_timestamp = Unix.gettimeofday () in
271 add_passive_clauses ~no_weight:true passive_empty_set passives
274 add_passive_clauses ~no_weight:true passive_empty_set g_passives
276 let g_actives = [] in
277 let actives = [], IDX.DT.empty in
279 given_clause ~noinfer:false
280 bag maxvar 0 max_steps timeout actives passives g_actives g_passives
282 | Sup.Success (bag, _, (i,_,_,_)) ->
284 let rec traverse ongoal (accg,acce) i =
285 match Terms.M.find i bag with
286 | (id,_,_,Terms.Exact _),_ ->
287 if ongoal then [i],acce else
288 if (List.mem i acce) then accg,acce else accg,acce@[i]
289 | (_,_,_,Terms.Step (_,i1,i2,_,_,_)),_ ->
290 if (not ongoal) && (List.mem i acce) then accg,acce
293 traverse false (traverse ongoal (accg,acce) i1) i2
295 if ongoal then i::accg,acce else accg,i::acce
297 let gsteps,esteps = traverse true ([],[]) i in
298 (List.rev esteps)@gsteps
301 (Printf.sprintf "Found proof, %fs"
302 (Unix.gettimeofday() -. initial_timestamp));
304 prerr_endline "Proof:";
306 prerr_endline (Pp.pp_unit_clause (fst(Terms.M.find x bag)))) l;
309 | Failure (msg,_bag,_maxvar,iterno) ->
311 prerr_endline (Printf.sprintf "FAILURE in %d iterations" iterno);