src/HOL/Tools/res_atp.ML
author quigley
Thu Mar 31 20:12:54 2005 +0200 (2005-03-31 ago)
changeset 15644 f2ef8c258fa4
parent 15608 f161fa6f8fd5
child 15653 3549ff7158f3
permissions -rw-r--r--
*** empty log message ***
paulson@15608
     1
(*  Author: Jia Meng, Cambridge University Computer Laboratory
paulson@15608
     2
    ID: $Id$
paulson@15608
     3
    Copyright 2004 University of Cambridge
paulson@15347
     4
paulson@15347
     5
ATPs with TPTP format input.
paulson@15347
     6
*)
paulson@15452
     7
paulson@15452
     8
(*Jia: changed: isar_atp now processes entire proof context.  fetch thms from delta simpset/claset*)
quigley@15644
     9
(*Claire: changed: added actual watcher calls *)
paulson@15452
    10
paulson@15347
    11
signature RES_ATP = 
paulson@15347
    12
sig
paulson@15608
    13
val trace_res : bool ref
paulson@15608
    14
val axiom_file : Path.T
paulson@15608
    15
val hyps_file : Path.T
paulson@15608
    16
val isar_atp : ProofContext.context * Thm.thm -> unit
quigley@15644
    17
(*val prob_file : Path.T*)
quigley@15644
    18
(*val atp_ax_tac : Thm.thm list -> int -> Tactical.tactic*)
quigley@15644
    19
(*val atp_tac : int -> Tactical.tactic*)
paulson@15608
    20
val debug: bool ref
paulson@15347
    21
paulson@15347
    22
end;
paulson@15347
    23
paulson@15347
    24
structure ResAtp : RES_ATP =
paulson@15608
    25
paulson@15347
    26
struct
paulson@15347
    27
quigley@15644
    28
quigley@15644
    29
paulson@15608
    30
(* used for debug *)
paulson@15608
    31
val debug = ref false;
paulson@15452
    32
paulson@15608
    33
fun debug_tac tac = (warning "testing";tac);
paulson@15608
    34
(* default value is false *)
paulson@15347
    35
paulson@15608
    36
val trace_res = ref false;
paulson@15347
    37
paulson@15608
    38
val skolem_tac = skolemize_tac;
paulson@15608
    39
quigley@15644
    40
val num_of_clauses = ref 0;
quigley@15644
    41
val clause_arr = Array.array(3500, ("empty", 0));
quigley@15644
    42
paulson@15347
    43
paulson@15608
    44
val atomize_tac =
paulson@15608
    45
    SUBGOAL
paulson@15608
    46
     (fn (prop,_) =>
paulson@15608
    47
	 let val ts = Logic.strip_assums_hyp prop
paulson@15608
    48
	 in EVERY1 
paulson@15608
    49
		[METAHYPS
paulson@15608
    50
		     (fn hyps => (cut_facts_tac (map (ObjectLogic.atomize_thm o forall_intr_vars) hyps) 1)),
paulson@15608
    51
	  REPEAT_DETERM_N (length ts) o (etac thin_rl)]
paulson@15608
    52
     end);
paulson@15347
    53
paulson@15608
    54
(* temporarily use these files, which will be loaded by Vampire *)
quigley@15644
    55
val file_id_num =ref 0;
paulson@15608
    56
quigley@15644
    57
fun new_prob_file () =  (file_id_num := (!file_id_num) + 1;"prob"^(string_of_int (!file_id_num)));
paulson@15608
    58
paulson@15347
    59
quigley@15644
    60
val axiom_file = File.tmp_path (Path.basic "axioms");
quigley@15644
    61
val clasimp_file = File.tmp_path (Path.basic "clasimp");
quigley@15644
    62
val hyps_file = File.tmp_path (Path.basic "hyps");
quigley@15644
    63
val prob_file = File.tmp_path (Path.basic "prob");
quigley@15644
    64
val dummy_tac = PRIMITIVE(fn thm => thm );
quigley@15644
    65
quigley@15644
    66
 
quigley@15644
    67
fun concat_with_and [] str = str
quigley@15644
    68
|   concat_with_and (x::[]) str = str^" ("^x^")"
quigley@15644
    69
|   concat_with_and (x::xs) str = (concat_with_and xs (str^"("^x^")"^" & "))
quigley@15644
    70
paulson@15347
    71
paulson@15347
    72
paulson@15347
    73
(**** for Isabelle/ML interface  ****)
paulson@15347
    74
quigley@15644
    75
fun is_proof_char ch = ((33 <= (ord ch)) andalso ((ord ch ) <= 126) andalso (not ((ord ch ) = 63))) orelse (ch = " ")
paulson@15608
    76
quigley@15644
    77
fun proofstring x = let val exp = explode x 
quigley@15644
    78
                    in
quigley@15644
    79
                        List.filter (is_proof_char ) exp
quigley@15644
    80
                    end
paulson@15608
    81
paulson@15608
    82
paulson@15452
    83
quigley@15644
    84
(*
quigley@15644
    85
fun call_atp_tac thms n = (tptp_inputs thms ; dummy_tac);
paulson@15347
    86
quigley@15644
    87
*)
paulson@15347
    88
paulson@15347
    89
(**** For running in Isar ****)
paulson@15347
    90
paulson@15608
    91
(* same function as that in res_axioms.ML *)
paulson@15608
    92
fun repeat_RS thm1 thm2 =
paulson@15608
    93
    let val thm1' =  thm1 RS thm2 handle THM _ => thm1
paulson@15608
    94
    in
paulson@15608
    95
	if eq_thm(thm1,thm1') then thm1' else (repeat_RS thm1' thm2)
paulson@15608
    96
    end;
paulson@15608
    97
paulson@15608
    98
(* a special version of repeat_RS *)
paulson@15608
    99
fun repeat_someI_ex thm = repeat_RS thm someI_ex;
paulson@15608
   100
quigley@15644
   101
(*********************************************************************)
paulson@15608
   102
(* convert clauses from "assume" to conjecture. write to file "hyps" *)
quigley@15644
   103
(* hypotheses of the goal currently being proved                     *)
quigley@15644
   104
(*********************************************************************)
quigley@15644
   105
paulson@15608
   106
fun isar_atp_h thms =
quigley@15644
   107
        
paulson@15608
   108
    let val prems = map (skolemize o make_nnf o ObjectLogic.atomize_thm) thms
quigley@15644
   109
        val prems' = map repeat_someI_ex prems
quigley@15644
   110
        val prems'' = make_clauses prems'
quigley@15644
   111
        val prems''' = ResAxioms.rm_Eps [] prems''
quigley@15644
   112
        val clss = map ResClause.make_conjecture_clause prems'''
paulson@15608
   113
	val (tptp_clss,tfree_litss) = ResLib.unzip (map ResClause.clause2tptp clss) 
paulson@15608
   114
	val tfree_lits = ResLib.flat_noDup tfree_litss
paulson@15608
   115
	val tfree_clss = map ResClause.tfree_clause tfree_lits 
paulson@15608
   116
        val hypsfile = File.sysify_path hyps_file
paulson@15608
   117
	val out = TextIO.openOut(hypsfile)
paulson@15608
   118
    in
paulson@15608
   119
	((ResLib.writeln_strs out (tfree_clss @ tptp_clss); TextIO.closeOut out; if !trace_res then (warning hypsfile) else ());tfree_lits) 
paulson@15608
   120
    end;
paulson@15347
   121
quigley@15644
   122
quigley@15644
   123
(*********************************************************************)
quigley@15644
   124
(* write out a subgoal as tptp clauses to the file "probN"           *)
quigley@15644
   125
(* where N is the number of this subgoal                             *)
quigley@15644
   126
(*********************************************************************)
quigley@15644
   127
paulson@15608
   128
fun tptp_inputs_tfrees thms n tfrees = 
paulson@15608
   129
    let val clss = map (ResClause.make_conjecture_clause_thm) thms
paulson@15608
   130
	val (tptp_clss,tfree_litss) = ResLib.unzip (map ResClause.clause2tptp clss)
paulson@15608
   131
	val tfree_clss = map ResClause.tfree_clause ((ResLib.flat_noDup tfree_litss) \\ tfrees) 
paulson@15608
   132
        val probfile = (File.sysify_path prob_file) ^ "_" ^ (string_of_int n)
paulson@15608
   133
	val out = TextIO.openOut(probfile)
paulson@15608
   134
    in
paulson@15608
   135
	(ResLib.writeln_strs out (tfree_clss @ tptp_clss); TextIO.closeOut out; (if !trace_res then (warning probfile) else ()))
paulson@15608
   136
    end;
paulson@15452
   137
paulson@15608
   138
paulson@15608
   139
quigley@15644
   140
(*********************************************************************)
quigley@15644
   141
(* call SPASS with settings and problem file for the current subgoal *)
quigley@15644
   142
(* should be modified to allow other provers to be called            *)
quigley@15644
   143
(*********************************************************************)
paulson@15608
   144
quigley@15644
   145
fun call_resolve_tac thms sg_term (childin, childout,pid)  = let
quigley@15644
   146
                            val newprobfile = new_prob_file ()
quigley@15644
   147
                             val thmstring = concat_with_and (map string_of_thm thms) ""
quigley@15644
   148
                           val _ = warning ("thmstring in call_res is: "^thmstring)
quigley@15644
   149
                            val goalstr = Sign.string_of_term Mainsign sg_term 
quigley@15644
   150
                            val goalproofstring = proofstring goalstr
quigley@15644
   151
                               val no_returns =List.filter not_newline ( goalproofstring)
quigley@15644
   152
                            val goalstring = implode no_returns
quigley@15644
   153
                            val _ = warning ("goalstring in call_res is: "^goalstring)
quigley@15644
   154
        
quigley@15644
   155
                            val prob_file =File.tmp_path (Path.basic newprobfile); 
quigley@15644
   156
                             val clauses = make_clauses thms
quigley@15644
   157
                            (*val _ = tptp_inputs clauses prob_file*)
quigley@15644
   158
                            val thmstring = concat_with_and (map string_of_thm thms) ""
quigley@15644
   159
                           
quigley@15644
   160
                            val goalstr = Sign.string_of_term Mainsign sg_term 
quigley@15644
   161
                            val goalproofstring = proofstring goalstr
quigley@15644
   162
                             val no_returns =List.filter not_newline ( goalproofstring)
quigley@15644
   163
                            val goalstring = implode no_returns
paulson@15608
   164
quigley@15644
   165
                            val thmproofstring = proofstring ( thmstring)
quigley@15644
   166
                            val no_returns =List.filter   not_newline ( thmproofstring)
quigley@15644
   167
                            val thmstr = implode no_returns
quigley@15644
   168
                            
quigley@15644
   169
                            val prob_path = File.sysify_path prob_file
quigley@15644
   170
                            val outfile = TextIO.openOut("/home/clq20/Jia_Code/hellofile")
quigley@15644
   171
                            val _ = TextIO.output(outfile, "prob file path is "^prob_path^" thmstring is "^thmstr^" goalstring is "^goalstring);
quigley@15644
   172
                            val _ = TextIO.flushOut outfile;
quigley@15644
   173
                            val _ =  TextIO.closeOut outfile
quigley@15644
   174
                          in
quigley@15644
   175
                           (* without paramodulation *)
quigley@15644
   176
                           (*(warning ("goalstring in call_res_tac is: "^goalstring));
quigley@15644
   177
                           (warning ("prob path in cal_res_tac is: "^prob_path));
quigley@15644
   178
                            Watcher.callResProvers(childout,
quigley@15644
   179
                            [("spass",thmstr,goalstring,"/homes/clq20/bin/SPASS",  
quigley@15644
   180
                             "-FullRed=0%-Auto=0%-ISRe%-ISFc%-RTaut%-RFSub%-RBSub%-DocProof", 
quigley@15644
   181
                             prob_path)]);*)
quigley@15644
   182
                           (* with paramodulation *)
quigley@15644
   183
                           (*Watcher.callResProvers(childout,
quigley@15644
   184
                                  [("spass",thmstr,goalstring,"/homes/clq20/bin/SPASS",
quigley@15644
   185
                                  "-FullRed=0%-ISPm=1%-Split=0%-PObv=0%-DocProof", 
quigley@15644
   186
                                    prob_path)]); *)
quigley@15644
   187
                           Watcher.callResProvers(childout,
quigley@15644
   188
                           [("spass",thmstr,goalstring,"/homes/clq20/bin/SPASS", 
quigley@15644
   189
                           "-DocProof",  prob_path)]);
quigley@15644
   190
                           dummy_tac
quigley@15644
   191
                         end
paulson@15452
   192
quigley@15644
   193
(************************************************)
quigley@15644
   194
(* pass in subgoal as a term and watcher info   *)
quigley@15644
   195
(* process subgoal into skolemized, negated form*)
quigley@15644
   196
(* then call call_resolve_tac to send to ATP    *)
quigley@15644
   197
(************************************************)
paulson@15608
   198
quigley@15644
   199
fun resolve_tac sg_term  (childin, childout,pid) = 
quigley@15644
   200
   let val _ = warning ("in resolve_tac ")
quigley@15644
   201
   in
quigley@15644
   202
   SELECT_GOAL
quigley@15644
   203
  (EVERY1 [rtac ccontr,atomize_tac,skolemize_tac,  METAHYPS(fn negs => (warning ("calling call_resolve_tac next ");call_resolve_tac negs sg_term (childin, childout,pid)))])
quigley@15644
   204
   end;
paulson@15452
   205
webertj@15603
   206
paulson@15608
   207
paulson@15452
   208
quigley@15644
   209
(* Need to replace call_atp_tac_tfrees with call res_provers as it's the dummy one *)
quigley@15644
   210
quigley@15644
   211
(**********************************************************)
quigley@15644
   212
(* write out the current subgoal as a tptp file, probN,   *)
quigley@15644
   213
(* then call dummy_tac - should be call_res_tac           *)
quigley@15644
   214
(**********************************************************)
quigley@15644
   215
quigley@15644
   216
fun call_atp_tac_tfrees thms n tfrees sg_term (childin, childout,pid) = 
quigley@15644
   217
                                         (tptp_inputs_tfrees (make_clauses thms) n tfrees; 
quigley@15644
   218
                                          resolve_tac sg_term (childin, childout, pid);
quigley@15644
   219
  					  dummy_tac);
quigley@15644
   220
quigley@15644
   221
fun atp_tac_tfrees tfrees sg_term (childin, childout,pid)  n  = 
quigley@15644
   222
let val _ = (warning ("in atp_tac_tfrees "))
quigley@15644
   223
   in
quigley@15644
   224
SELECT_GOAL
quigley@15644
   225
  (EVERY1 [rtac ccontr,atomize_tac, skolemize_tac, 
quigley@15644
   226
  METAHYPS(fn negs => (call_atp_tac_tfrees (negs) n tfrees sg_term (childin, childout,pid)))]) n
quigley@15644
   227
end;
quigley@15644
   228
quigley@15644
   229
quigley@15644
   230
fun isar_atp_g tfrees sg_term (childin, childout, pid)  = 
quigley@15644
   231
                                        
quigley@15644
   232
(	(warning("in isar_atp_g"));atp_tac_tfrees tfrees sg_term (childin, childout, pid));
quigley@15644
   233
quigley@15644
   234
quigley@15644
   235
quigley@15644
   236
(**********************************************)
quigley@15644
   237
(* recursively call atp_tac_g on all subgoals *)
quigley@15644
   238
(* sg_term is the nth subgoal as a term - used*)
quigley@15644
   239
(* in proof reconstruction                    *)
quigley@15644
   240
(**********************************************)
quigley@15644
   241
quigley@15644
   242
fun isar_atp_goal' thm k n tfree_lits  (childin, childout, pid) = 
quigley@15644
   243
                  	if (k > n) 
quigley@15644
   244
                        then () 
quigley@15644
   245
	  		else 
quigley@15644
   246
                           (let val  prems = prems_of thm 
quigley@15644
   247
                                val sg_term = get_nth n prems
quigley@15644
   248
                            in   
quigley@15644
   249
                                 
quigley@15644
   250
                		(warning("in isar_atp_goal'"));
quigley@15644
   251
 				 isar_atp_g tfree_lits  sg_term (childin, childout, pid) k thm; 
quigley@15644
   252
                                 isar_atp_goal' thm (k+1) n tfree_lits  (childin, childout, pid) 
quigley@15644
   253
                            end);
quigley@15644
   254
quigley@15644
   255
quigley@15644
   256
fun isar_atp_goal thm n_subgoals tfree_lits   (childin, childout, pid) = (if (!debug) then warning (string_of_thm thm) else (isar_atp_goal' thm 1 n_subgoals tfree_lits  (childin, childout, pid) ));
quigley@15644
   257
quigley@15644
   258
(**************************************************)
quigley@15644
   259
(* convert clauses from "assume" to conjecture.   *)
quigley@15644
   260
(* i.e. apply make_clauses and then get tptp for  *)
quigley@15644
   261
(* any hypotheses in the goal                     *)
quigley@15644
   262
(* write to file "hyps"                           *)
quigley@15644
   263
(**************************************************)
quigley@15644
   264
quigley@15644
   265
quigley@15644
   266
fun isar_atp_aux thms thm n_subgoals  (childin, childout, pid) = 
paulson@15608
   267
    let val tfree_lits = isar_atp_h thms 
paulson@15608
   268
    in
quigley@15644
   269
	(warning("in isar_atp_aux"));isar_atp_goal thm n_subgoals tfree_lits   (childin, childout, pid)
paulson@15608
   270
    end;
quigley@15644
   271
quigley@15644
   272
(******************************************************************)
quigley@15644
   273
(* called in Isar automatically                                   *)
quigley@15644
   274
(* writes out the current clasimpset to a tptp file               *)
quigley@15644
   275
(* passes all subgoals on to isar_atp_aux for further processing  *)
quigley@15644
   276
(* turns off xsymbol at start of function, restoring it at end    *)
quigley@15644
   277
(******************************************************************)
paulson@15452
   278
paulson@15608
   279
fun isar_atp' (thms, thm) =
quigley@15644
   280
    let val _=  (print_mode := (Library.gen_rems (op =) (! print_mode, ["xsymbols", "symbols"])))
quigley@15644
   281
        val _= (warning ("in isar_atp'"))
quigley@15644
   282
        val prems = prems_of thm
quigley@15644
   283
        val thms_string =concat_with_and (map  string_of_thm thms) ""
quigley@15644
   284
        val thmstring = string_of_thm thm
quigley@15644
   285
        val prems_string = concat_with_and (map (Sign.string_of_term Mainsign) prems) ""
quigley@15644
   286
        (* set up variables for writing out the clasimps to a tptp file *)
quigley@15644
   287
        (* val _ = write_out_clasimp (File.sysify_path axiom_file)*)
quigley@15644
   288
        (* cq: add write out clasimps to file *)
quigley@15644
   289
        (* cq:create watcher and pass to isar_atp_aux *)                    
quigley@15644
   290
        val (childin,childout,pid) = Watcher.createWatcher(thm)
quigley@15644
   291
        val pidstring = string_of_int(Word.toInt (Word.fromLargeWord ( Posix.Process.pidToWord pid )))
paulson@15608
   292
    in
quigley@15644
   293
	case prems of [] => () 
quigley@15644
   294
		    | _ => ((warning ("initial thms: "^thms_string)); 
quigley@15644
   295
                           (warning ("initial thm: "^thmstring));
quigley@15644
   296
                           (warning ("subgoals: "^prems_string));
quigley@15644
   297
                           (warning ("pid: "^ pidstring))); 
quigley@15644
   298
                            isar_atp_aux thms thm (length prems) (childin, childout, pid) ;
quigley@15644
   299
                           
quigley@15644
   300
                           print_mode := (["xsymbols", "symbols"] @ ! print_mode)
paulson@15608
   301
    end;
paulson@15608
   302
paulson@15452
   303
paulson@15608
   304
paulson@15608
   305
paulson@15608
   306
local
paulson@15608
   307
paulson@15608
   308
fun get_thms_cs claset =
paulson@15608
   309
    let val clsset = rep_cs claset
paulson@15608
   310
	val safeEs = #safeEs clsset
paulson@15608
   311
	val safeIs = #safeIs clsset
paulson@15608
   312
	val hazEs = #hazEs clsset
paulson@15608
   313
	val hazIs = #hazIs clsset
paulson@15608
   314
    in
paulson@15608
   315
	safeEs @ safeIs @ hazEs @ hazIs
paulson@15608
   316
    end;
paulson@15608
   317
paulson@15608
   318
paulson@15452
   319
paulson@15608
   320
fun append_name name [] _ = []
paulson@15608
   321
  | append_name name (thm::thms) k = (Thm.name_thm ((name ^ "_" ^ (string_of_int k)),thm)) :: (append_name name thms (k+1));
paulson@15608
   322
paulson@15608
   323
fun append_names (name::names) (thms::thmss) =
paulson@15608
   324
    let val thms' = append_name name thms 0
paulson@15608
   325
    in
paulson@15608
   326
	thms'::(append_names names thmss)
paulson@15608
   327
    end;
paulson@15608
   328
paulson@15452
   329
paulson@15608
   330
fun get_thms_ss [] = []
paulson@15608
   331
  | get_thms_ss thms =
paulson@15608
   332
    let val names = map Thm.name_of_thm thms 
paulson@15608
   333
        val thms' = map (mksimps mksimps_pairs) thms
paulson@15608
   334
        val thms'' = append_names names thms'
paulson@15608
   335
    in
paulson@15608
   336
	ResLib.flat_noDup thms''
paulson@15608
   337
    end;
paulson@15608
   338
paulson@15452
   339
paulson@15608
   340
paulson@15608
   341
paulson@15608
   342
in
paulson@15608
   343
paulson@15452
   344
paulson@15608
   345
(* convert locally declared rules to axiom clauses *)
paulson@15608
   346
(* write axiom clauses to ax_file *)
quigley@15644
   347
(* what about clasimpset - it should already be in the ax file - perhaps append to ax file rather than just *)
quigley@15644
   348
(* write out ? Or keep as a separate file and then cat them all together in the watcher, like we do with the *)
quigley@15644
   349
(*claset file and prob file*)
paulson@15608
   350
fun isar_local_thms (delta_cs, delta_ss_thms) =
paulson@15608
   351
    let val thms_cs = get_thms_cs delta_cs
paulson@15608
   352
	val thms_ss = get_thms_ss delta_ss_thms
paulson@15608
   353
	val thms_clauses = ResLib.flat_noDup (map ResAxioms.clausify_axiom (thms_cs @ thms_ss))
paulson@15608
   354
	val clauses_strs = ResLib.flat_noDup (map ResClause.tptp_clause thms_clauses) (*string list*)
paulson@15608
   355
	val ax_file = File.sysify_path axiom_file
paulson@15608
   356
	val out = TextIO.openOut ax_file
paulson@15608
   357
    in
quigley@15644
   358
	(ResLib.writeln_strs out clauses_strs; (warning ("axiom file is: "^ax_file));TextIO.closeOut out)
paulson@15608
   359
    end;
paulson@15347
   360
paulson@15608
   361
paulson@15608
   362
paulson@15608
   363
paulson@15347
   364
paulson@15608
   365
(* called in Isar automatically *)
paulson@15608
   366
fun isar_atp (ctxt,thm) =
paulson@15608
   367
    let val prems = ProofContext.prems_of ctxt
quigley@15644
   368
        val d_cs = Classical.get_delta_claset ctxt 
quigley@15644
   369
        val d_ss_thms = Simplifier.get_delta_simpset ctxt
quigley@15644
   370
        val thmstring = string_of_thm thm
quigley@15644
   371
        val prem_no = length prems
quigley@15644
   372
        val prems_string = concat_with_and (map string_of_thm prems) ""
paulson@15608
   373
    in
quigley@15644
   374
         
quigley@15644
   375
          (warning ("initial thm in isar_atp: "^thmstring));
quigley@15644
   376
          (warning ("subgoals in isar_atp: "^prems_string));
quigley@15644
   377
    	   (warning ("number of subgoals in isar_atp: "^(string_of_int prem_no)));
quigley@15644
   378
          (isar_local_thms (d_cs,d_ss_thms); (warning("about to call isar_atp'"));
quigley@15644
   379
           isar_atp' (prems, thm))
paulson@15608
   380
    end;
paulson@15347
   381
paulson@15608
   382
end
paulson@15608
   383
paulson@15608
   384
paulson@15608
   385
paulson@15452
   386
paulson@15347
   387
end;
paulson@15347
   388
paulson@15347
   389
Proof.atp_hook := ResAtp.isar_atp;