src/HOL/Tools/Sledgehammer/sledgehammer_mash.ML
author blanchet
Mon, 23 Jul 2012 15:32:30 +0200
changeset 49450 f30eb5eb7927
parent 49449 aaaec69db3db
child 49451 72a31418ff8d
permissions -rw-r--r--
include unknown local facts in MaSh
     1 (*  Title:      HOL/Tools/Sledgehammer/sledgehammer_mash.ML
     2     Author:     Jasmin Blanchette, TU Muenchen
     3 
     4 Sledgehammer's machine-learning-based relevance filter (MaSh).
     5 *)
     6 
     7 signature SLEDGEHAMMER_MASH =
     8 sig
     9   type stature = ATP_Problem_Generate.stature
    10   type fact = Sledgehammer_Fact.fact
    11   type fact_override = Sledgehammer_Fact.fact_override
    12   type params = Sledgehammer_Provers.params
    13   type relevance_fudge = Sledgehammer_Provers.relevance_fudge
    14   type prover_result = Sledgehammer_Provers.prover_result
    15 
    16   val trace : bool Config.T
    17   val MaShN : string
    18   val mepoN : string
    19   val mashN : string
    20   val meshN : string
    21   val unlearnN : string
    22   val learn_isarN : string
    23   val learn_atpN : string
    24   val relearn_isarN : string
    25   val relearn_atpN : string
    26   val fact_filters : string list
    27   val escape_meta : string -> string
    28   val escape_metas : string list -> string
    29   val unescape_meta : string -> string
    30   val unescape_metas : string -> string list
    31   val extract_query : string -> string * (string * real) list
    32   val nickname_of : thm -> string
    33   val suggested_facts :
    34     (string * 'a) list -> ('b * thm) list -> (('b * thm) * 'a) list
    35   val mesh_facts :
    36     int -> ((('a * thm) * real) list * ('a * thm) list) list -> ('a * thm) list
    37   val theory_ord : theory * theory -> order
    38   val thm_ord : thm * thm -> order
    39   val goal_of_thm : theory -> thm -> thm
    40   val run_prover_for_mash :
    41     Proof.context -> params -> string -> fact list -> thm -> prover_result
    42   val features_of :
    43     Proof.context -> string -> theory -> stature -> term list -> string list
    44   val isar_dependencies_of : unit Symtab.table -> thm -> string list option
    45   val atp_dependencies_of :
    46     Proof.context -> params -> string -> int -> fact list -> unit Symtab.table
    47     -> thm -> string list option
    48   val mash_CLEAR : Proof.context -> unit
    49   val mash_ADD :
    50     Proof.context -> bool
    51     -> (string * string list * string list * string list) list -> unit
    52   val mash_REPROVE :
    53     Proof.context -> bool -> (string * string list) list -> unit
    54   val mash_QUERY :
    55     Proof.context -> bool -> int -> string list * string list
    56     -> (string * real) list
    57   val mash_unlearn : Proof.context -> unit
    58   val mash_could_suggest_facts : unit -> bool
    59   val mash_can_suggest_facts : Proof.context -> bool
    60   val mash_suggested_facts :
    61     Proof.context -> params -> string -> int -> term list -> term
    62     -> fact list -> (fact * real) list * fact list
    63   val mash_learn_proof :
    64     Proof.context -> params -> string -> term -> ('a * thm) list -> thm list
    65     -> unit
    66   val mash_learn :
    67     Proof.context -> params -> fact_override -> thm list -> bool -> unit
    68   val relevant_facts :
    69     Proof.context -> params -> string -> int -> fact_override -> term list
    70     -> term -> fact list -> fact list
    71   val kill_learners : unit -> unit
    72   val running_learners : unit -> unit
    73 end;
    74 
    75 structure Sledgehammer_MaSh : SLEDGEHAMMER_MASH =
    76 struct
    77 
    78 open ATP_Util
    79 open ATP_Problem_Generate
    80 open Sledgehammer_Util
    81 open Sledgehammer_Fact
    82 open Sledgehammer_Provers
    83 open Sledgehammer_Minimize
    84 open Sledgehammer_MePo
    85 
    86 val trace =
    87   Attrib.setup_config_bool @{binding sledgehammer_mash_trace} (K false)
    88 fun trace_msg ctxt msg = if Config.get ctxt trace then tracing (msg ()) else ()
    89 
    90 val MaShN = "MaSh"
    91 
    92 val mepoN = "mepo"
    93 val mashN = "mash"
    94 val meshN = "mesh"
    95 
    96 val fact_filters = [meshN, mepoN, mashN]
    97 
    98 val unlearnN = "unlearn"
    99 val learn_isarN = "learn_isar"
   100 val learn_atpN = "learn_atp"
   101 val relearn_isarN = "relearn_isar"
   102 val relearn_atpN = "relearn_atp"
   103 
   104 fun mash_home () = getenv "MASH_HOME"
   105 fun mash_model_dir () =
   106   getenv "ISABELLE_HOME_USER" ^ "/mash"
   107   |> tap (Isabelle_System.mkdir o Path.explode)
   108 val mash_state_dir = mash_model_dir
   109 fun mash_state_path () = mash_state_dir () ^ "/state" |> Path.explode
   110 
   111 
   112 (*** Isabelle helpers ***)
   113 
   114 fun meta_char c =
   115   if Char.isAlphaNum c orelse c = #"_" orelse c = #"." orelse c = #"(" orelse
   116      c = #")" orelse c = #"," then
   117     String.str c
   118   else
   119     (* fixed width, in case more digits follow *)
   120     "%" ^ stringN_of_int 3 (Char.ord c)
   121 
   122 fun unmeta_chars accum [] = String.implode (rev accum)
   123   | unmeta_chars accum (#"%" :: d1 :: d2 :: d3 :: cs) =
   124     (case Int.fromString (String.implode [d1, d2, d3]) of
   125        SOME n => unmeta_chars (Char.chr n :: accum) cs
   126      | NONE => "" (* error *))
   127   | unmeta_chars _ (#"%" :: _) = "" (* error *)
   128   | unmeta_chars accum (c :: cs) = unmeta_chars (c :: accum) cs
   129 
   130 val escape_meta = String.translate meta_char
   131 val escape_metas = map escape_meta #> space_implode " "
   132 val unescape_meta = String.explode #> unmeta_chars []
   133 val unescape_metas =
   134   space_explode " " #> filter_out (curry (op =) "") #> map unescape_meta
   135 
   136 fun extract_node line =
   137   case space_explode ":" line of
   138     [name, parents] => (unescape_meta name, unescape_metas parents)
   139   | _ => ("", [])
   140 
   141 fun extract_suggestion sugg =
   142   case space_explode "=" sugg of
   143     [name, weight] =>
   144     SOME (unescape_meta name, Real.fromString weight |> the_default 0.0)
   145   | _ => NONE
   146 
   147 fun extract_query line =
   148   case space_explode ":" line of
   149     [goal, suggs] =>
   150     (unescape_meta goal,
   151      map_filter extract_suggestion (space_explode " " suggs))
   152   | _ => ("", [])
   153 
   154 fun parent_of_local_thm th =
   155   let
   156     val thy = th |> Thm.theory_of_thm
   157     val facts = thy |> Global_Theory.facts_of
   158     val space = facts |> Facts.space_of
   159     fun id_of s = #id (Name_Space.the_entry space s)
   160     fun max_id (s', _) (s, id) =
   161       let val id' = id_of s' in if id > id' then (s, id) else (s', id') end
   162   in ("", ~1) |> Facts.fold_static max_id facts |> fst end
   163 
   164 val local_prefix = "local" ^ Long_Name.separator
   165 
   166 fun nickname_of th =
   167   if Thm.has_name_hint th then
   168     let val hint = Thm.get_name_hint th in
   169       (* FIXME: There must be a better way to detect local facts. *)
   170       case try (unprefix local_prefix) hint of
   171         SOME suf =>
   172         parent_of_local_thm th ^ Long_Name.separator ^ Long_Name.separator ^ suf
   173       | NONE => hint
   174     end
   175   else
   176     backquote_thm th
   177 
   178 fun suggested_facts suggs facts =
   179   let
   180     fun add_fact (fact as (_, th)) = Symtab.default (nickname_of th, fact)
   181     val tab = Symtab.empty |> fold add_fact facts
   182     fun find_sugg (name, weight) =
   183       Symtab.lookup tab name |> Option.map (rpair weight)
   184   in map_filter find_sugg suggs end
   185 
   186 fun sum_avg [] = 0
   187   | sum_avg xs =
   188     Real.ceil (100000000.0 * fold (curry (op +)) xs 0.0) div length xs
   189 
   190 fun normalize_scores [] = []
   191   | normalize_scores ((fact, score) :: tail) =
   192     (fact, 1.0) :: map (apsnd (curry Real.* (1.0 / score))) tail
   193 
   194 fun mesh_facts max_facts [(sels, unks)] =
   195     map fst (take max_facts sels) @ take (max_facts - length sels) unks
   196   | mesh_facts max_facts mess =
   197     let
   198       val mess = mess |> map (apfst (normalize_scores #> `length))
   199       val fact_eq = Thm.eq_thm o pairself snd
   200       fun score_at sels = try (nth sels) #> Option.map snd
   201       fun score_in fact ((sel_len, sels), unks) =
   202         case find_index (curry fact_eq fact o fst) sels of
   203           ~1 => (case find_index (curry fact_eq fact) unks of
   204                    ~1 => score_at sels sel_len
   205                  | _ => NONE)
   206         | rank => score_at sels rank
   207       fun weight_of fact = mess |> map_filter (score_in fact) |> sum_avg
   208       val facts =
   209         fold (union fact_eq o map fst o take max_facts o snd o fst) mess []
   210     in
   211       facts |> map (`weight_of) |> sort (int_ord o swap o pairself fst)
   212             |> map snd |> take max_facts
   213     end
   214 
   215 val thy_feature_name_of = prefix "y"
   216 val const_name_of = prefix "c"
   217 val type_name_of = prefix "t"
   218 val class_name_of = prefix "s"
   219 
   220 fun theory_ord p =
   221   if Theory.eq_thy p then
   222     EQUAL
   223   else if Theory.subthy p then
   224     LESS
   225   else if Theory.subthy (swap p) then
   226     GREATER
   227   else case int_ord (pairself (length o Theory.ancestors_of) p) of
   228     EQUAL => string_ord (pairself Context.theory_name p)
   229   | order => order
   230 
   231 val thm_ord = theory_ord o pairself theory_of_thm
   232 
   233 val freezeT = Type.legacy_freeze_type
   234 
   235 fun freeze (t $ u) = freeze t $ freeze u
   236   | freeze (Abs (s, T, t)) = Abs (s, freezeT T, freeze t)
   237   | freeze (Var ((s, _), T)) = Free (s, freezeT T)
   238   | freeze (Const (s, T)) = Const (s, freezeT T)
   239   | freeze (Free (s, T)) = Free (s, freezeT T)
   240   | freeze t = t
   241 
   242 fun goal_of_thm thy = prop_of #> freeze #> cterm_of thy #> Goal.init
   243 
   244 fun run_prover_for_mash ctxt params prover facts goal =
   245   let
   246     val problem =
   247       {state = Proof.init ctxt, goal = goal, subgoal = 1, subgoal_count = 1,
   248        facts = facts |> map (apfst (apfst (fn name => name ())))
   249                      |> map Untranslated_Fact}
   250   in
   251     get_minimizing_prover ctxt MaSh (K (K ())) prover params (K (K (K "")))
   252                           problem
   253   end
   254 
   255 val bad_types = [@{type_name prop}, @{type_name bool}, @{type_name fun}]
   256 
   257 val logical_consts =
   258   [@{const_name prop}, @{const_name Pure.conjunction}] @ atp_logical_consts
   259 
   260 fun interesting_terms_types_and_classes ctxt prover term_max_depth
   261                                         type_max_depth ts =
   262   let
   263     fun is_bad_const (x as (s, _)) args =
   264       member (op =) logical_consts s orelse
   265       fst (is_built_in_const_for_prover ctxt prover x args)
   266     fun add_classes @{sort type} = I
   267       | add_classes S = union (op =) (map class_name_of S)
   268     fun do_add_type (Type (s, Ts)) =
   269         (not (member (op =) bad_types s) ? insert (op =) (type_name_of s))
   270         #> fold do_add_type Ts
   271       | do_add_type (TFree (_, S)) = add_classes S
   272       | do_add_type (TVar (_, S)) = add_classes S
   273     fun add_type T = type_max_depth >= 0 ? do_add_type T
   274     fun mk_app s args =
   275       if member (op <>) args "" then s ^ "(" ^ space_implode "," args ^ ")"
   276       else s
   277     fun patternify ~1 _ = ""
   278       | patternify depth t =
   279         case strip_comb t of
   280           (Const (x as (s, _)), args) =>
   281           if is_bad_const x args then ""
   282           else mk_app (const_name_of s) (map (patternify (depth - 1)) args)
   283         | _ => ""
   284     fun add_pattern depth t =
   285       case patternify depth t of "" => I | s => insert (op =) s
   286     fun add_term_patterns ~1 _ = I
   287       | add_term_patterns depth t =
   288         add_pattern depth t #> add_term_patterns (depth - 1) t
   289     val add_term = add_term_patterns term_max_depth
   290     fun add_patterns t =
   291       let val (head, args) = strip_comb t in
   292         (case head of
   293            Const (_, T) => add_term t #> add_type T
   294          | Free (_, T) => add_type T
   295          | Var (_, T) => add_type T
   296          | Abs (_, T, body) => add_type T #> add_patterns body
   297          | _ => I)
   298         #> fold add_patterns args
   299       end
   300   in [] |> fold add_patterns ts end
   301 
   302 fun is_exists (s, _) = (s = @{const_name Ex} orelse s = @{const_name Ex1})
   303 
   304 val term_max_depth = 1
   305 val type_max_depth = 1
   306 
   307 (* TODO: Generate type classes for types? *)
   308 fun features_of ctxt prover thy (scope, status) ts =
   309   thy_feature_name_of (Context.theory_name thy) ::
   310   interesting_terms_types_and_classes ctxt prover term_max_depth type_max_depth
   311                                       ts
   312   |> forall is_lambda_free ts ? cons "no_lams"
   313   |> forall (not o exists_Const is_exists) ts ? cons "no_skos"
   314   |> scope <> Global ? cons "local"
   315   |> (case status of
   316         General => I
   317       | Induction => cons "induction"
   318       | Intro => cons "intro"
   319       | Inductive => cons "inductive"
   320       | Elim => cons "elim"
   321       | Simp => cons "simp"
   322       | Def => cons "def")
   323 
   324 (* Too many dependencies is a sign that a decision procedure is at work. There
   325    isn't much too learn from such proofs. *)
   326 val max_dependencies = 10
   327 val atp_dependency_default_max_fact = 50
   328 
   329 fun trim_dependencies deps =
   330   if length deps <= max_dependencies then SOME deps else NONE
   331 
   332 fun isar_dependencies_of all_names =
   333   thms_in_proof (SOME all_names) #> trim_dependencies
   334 
   335 fun atp_dependencies_of ctxt (params as {verbose, max_facts, ...}) prover
   336                         auto_level facts all_names th =
   337   case isar_dependencies_of all_names th of
   338     SOME [] => NONE
   339   | isar_deps =>
   340     let
   341       val thy = Proof_Context.theory_of ctxt
   342       val goal = goal_of_thm thy th
   343       val (_, hyp_ts, concl_t) = ATP_Util.strip_subgoal ctxt goal 1
   344       val facts = facts |> filter (fn (_, th') => thm_ord (th', th) = LESS)
   345       fun fix_name ((_, stature), th) = ((fn () => nickname_of th, stature), th)
   346       fun is_dep dep (_, th) = nickname_of th = dep
   347       fun add_isar_dep facts dep accum =
   348         if exists (is_dep dep) accum then
   349           accum
   350         else case find_first (is_dep dep) facts of
   351           SOME ((name, status), th) => accum @ [((name, status), th)]
   352         | NONE => accum (* shouldn't happen *)
   353       val facts =
   354         facts |> mepo_suggested_facts ctxt params prover
   355                      (max_facts |> the_default atp_dependency_default_max_fact)
   356                      NONE hyp_ts concl_t
   357               |> fold (add_isar_dep facts) (these isar_deps)
   358               |> map fix_name
   359     in
   360       if verbose andalso auto_level = 0 then
   361         let val num_facts = length facts in
   362           "MaSh: " ^ quote prover ^ " on " ^ quote (nickname_of th) ^
   363           " with " ^ string_of_int num_facts ^ " fact" ^ plural_s num_facts ^
   364           "."
   365           |> Output.urgent_message
   366         end
   367       else
   368         ();
   369       case run_prover_for_mash ctxt params prover facts goal of
   370         {outcome = NONE, used_facts, ...} =>
   371         (if verbose andalso auto_level = 0 then
   372            let val num_facts = length used_facts in
   373              "Found proof with " ^ string_of_int num_facts ^ " fact" ^
   374              plural_s num_facts ^ "."
   375              |> Output.urgent_message
   376            end
   377          else
   378            ();
   379          used_facts |> map fst |> trim_dependencies)
   380       | _ => NONE
   381     end
   382 
   383 
   384 (*** Low-level communication with MaSh ***)
   385 
   386 (* more friendly than "try o File.rm" for those who keep the files open in their
   387    text editor *)
   388 fun wipe_out file = File.write file ""
   389 
   390 fun write_file (xs, f) file =
   391   let val path = Path.explode file in
   392     wipe_out path;
   393     xs |> chunk_list 500
   394        |> List.app (File.append path o space_implode "" o map f)
   395   end
   396 
   397 fun run_mash_tool ctxt overlord save max_suggs write_cmds read_suggs =
   398   let
   399     val (temp_dir, serial) =
   400       if overlord then (getenv "ISABELLE_HOME_USER", "")
   401       else (getenv "ISABELLE_TMP", serial_string ())
   402     val log_file = if overlord then temp_dir ^ "/mash_log" else "/dev/null"
   403     val err_file = temp_dir ^ "/mash_err" ^ serial
   404     val sugg_file = temp_dir ^ "/mash_suggs" ^ serial
   405     val cmd_file = temp_dir ^ "/mash_commands" ^ serial
   406     val core =
   407       "--inputFile " ^ cmd_file ^ " --predictions " ^ sugg_file ^
   408       " --numberOfPredictions " ^ string_of_int max_suggs ^
   409       (if save then " --saveModel" else "")
   410     val command =
   411       mash_home () ^ "/mash --quiet --outputDir " ^ mash_model_dir () ^
   412       " --log " ^ log_file ^ " " ^ core ^ " >& " ^ err_file
   413   in
   414     write_file ([], K "") sugg_file;
   415     write_file write_cmds cmd_file;
   416     trace_msg ctxt (fn () => "Running " ^ command);
   417     Isabelle_System.bash command;
   418     read_suggs (fn () => try File.read_lines (Path.explode sugg_file) |> these)
   419     |> tap (fn _ => trace_msg ctxt (fn () =>
   420            case try File.read (Path.explode err_file) of
   421              NONE => "Done"
   422            | SOME "" => "Done"
   423            | SOME s => "Error: " ^ elide_string 1000 s))
   424     |> not overlord
   425        ? tap (fn _ => List.app (wipe_out o Path.explode)
   426                                [err_file, sugg_file, cmd_file])
   427   end
   428 
   429 fun str_of_add (name, parents, feats, deps) =
   430   "! " ^ escape_meta name ^ ": " ^ escape_metas parents ^ "; " ^
   431   escape_metas feats ^ "; " ^ escape_metas deps ^ "\n"
   432 
   433 fun str_of_reprove (name, deps) =
   434   "p " ^ escape_meta name ^ ": " ^ escape_metas deps ^ "\n"
   435 
   436 fun str_of_query (parents, feats) =
   437   "? " ^ escape_metas parents ^ "; " ^ escape_metas feats ^ "\n"
   438 
   439 fun mash_CLEAR ctxt =
   440   let val path = mash_model_dir () |> Path.explode in
   441     trace_msg ctxt (K "MaSh CLEAR");
   442     File.fold_dir (fn file => fn _ =>
   443                       try File.rm (Path.append path (Path.basic file)))
   444                   path NONE;
   445     ()
   446   end
   447 
   448 fun mash_ADD _ _ [] = ()
   449   | mash_ADD ctxt overlord adds =
   450     (trace_msg ctxt (fn () => "MaSh ADD " ^
   451          elide_string 1000 (space_implode " " (map #1 adds)));
   452      run_mash_tool ctxt overlord true 0 (adds, str_of_add) (K ()))
   453 
   454 fun mash_REPROVE _ _ [] = ()
   455   | mash_REPROVE ctxt overlord reps =
   456     (trace_msg ctxt (fn () => "MaSh REPROVE " ^
   457          elide_string 1000 (space_implode " " (map #1 reps)));
   458      run_mash_tool ctxt overlord true 0 (reps, str_of_reprove) (K ()))
   459 
   460 fun mash_QUERY ctxt overlord max_suggs (query as (_, feats)) =
   461   (trace_msg ctxt (fn () => "MaSh QUERY " ^ space_implode " " feats);
   462    run_mash_tool ctxt overlord false max_suggs
   463        ([query], str_of_query)
   464        (fn suggs =>
   465            case suggs () of
   466              [] => []
   467            | suggs => snd (extract_query (List.last suggs)))
   468    handle List.Empty => [])
   469 
   470 
   471 (*** High-level communication with MaSh ***)
   472 
   473 fun try_graph ctxt when def f =
   474   f ()
   475   handle Graph.CYCLES (cycle :: _) =>
   476          (trace_msg ctxt (fn () =>
   477               "Cycle involving " ^ commas cycle ^ " when " ^ when); def)
   478        | Graph.DUP name =>
   479          (trace_msg ctxt (fn () =>
   480               "Duplicate fact " ^ quote name ^ " when " ^ when); def)
   481        | Graph.UNDEF name =>
   482          (trace_msg ctxt (fn () =>
   483               "Unknown fact " ^ quote name ^ " when " ^ when); def)
   484        | exn =>
   485          if Exn.is_interrupt exn then
   486            reraise exn
   487          else
   488            (trace_msg ctxt (fn () =>
   489                 "Internal error when " ^ when ^ ":\n" ^
   490                 ML_Compiler.exn_message exn); def)
   491 
   492 fun graph_info G =
   493   string_of_int (length (Graph.keys G)) ^ " node(s), " ^
   494   string_of_int (fold (Integer.add o length o snd) (Graph.dest G) 0) ^
   495   " edge(s), " ^
   496   string_of_int (length (Graph.minimals G)) ^ " minimal, " ^
   497   string_of_int (length (Graph.maximals G)) ^ " maximal"
   498 
   499 type mash_state = {fact_G : unit Graph.T}
   500 
   501 val empty_state = {fact_G = Graph.empty}
   502 
   503 local
   504 
   505 val version = "*** MaSh 0.0 ***"
   506 
   507 fun load _ (state as (true, _)) = state
   508   | load ctxt _ =
   509     let val path = mash_state_path () in
   510       (true,
   511        case try File.read_lines path of
   512          SOME (version' :: node_lines) =>
   513          let
   514            fun add_edge_to name parent =
   515              Graph.default_node (parent, ()) #> Graph.add_edge (parent, name)
   516            fun add_node line =
   517              case extract_node line of
   518                ("", _) => I (* shouldn't happen *)
   519              | (name, parents) =>
   520                Graph.default_node (name, ()) #> fold (add_edge_to name) parents
   521            val fact_G =
   522              try_graph ctxt "loading state" Graph.empty (fn () =>
   523                  Graph.empty |> version' = version ? fold add_node node_lines)
   524          in
   525            trace_msg ctxt (fn () =>
   526                "Loaded fact graph (" ^ graph_info fact_G ^ ")");
   527            {fact_G = fact_G}
   528          end
   529        | _ => empty_state)
   530     end
   531 
   532 fun save ctxt {fact_G} =
   533   let
   534     val path = mash_state_path ()
   535     fun fact_line_for name parents =
   536       escape_meta name ^ ": " ^ escape_metas parents
   537     val append_fact = File.append path o suffix "\n" oo fact_line_for
   538     fun append_entry (name, ((), (parents, _))) () =
   539       append_fact name (Graph.Keys.dest parents)
   540   in
   541     File.write path (version ^ "\n");
   542     Graph.fold append_entry fact_G ();
   543     trace_msg ctxt (fn () => "Saved fact graph (" ^ graph_info fact_G ^ ")")
   544   end
   545 
   546 val global_state =
   547   Synchronized.var "Sledgehammer_MaSh.global_state" (false, empty_state)
   548 
   549 in
   550 
   551 fun mash_map ctxt f =
   552   Synchronized.change global_state (load ctxt ##> (f #> tap (save ctxt)))
   553 
   554 fun mash_peek ctxt f =
   555   Synchronized.change_result global_state (load ctxt #> `snd #>> f)
   556 
   557 fun mash_get ctxt =
   558   Synchronized.change_result global_state (load ctxt #> `snd)
   559 
   560 fun mash_unlearn ctxt =
   561   Synchronized.change global_state (fn _ =>
   562       (mash_CLEAR ctxt; wipe_out (mash_state_path ()); (true, empty_state)))
   563 
   564 end
   565 
   566 fun mash_could_suggest_facts () = mash_home () <> ""
   567 fun mash_can_suggest_facts ctxt = not (Graph.is_empty (#fact_G (mash_get ctxt)))
   568 
   569 fun num_keys keys = Graph.Keys.fold (K (Integer.add 1)) keys 0
   570 
   571 fun maximal_in_graph fact_G facts =
   572   let
   573     val facts = [] |> fold (cons o nickname_of o snd) facts
   574     val tab = Symtab.empty |> fold (fn name => Symtab.default (name, ())) facts
   575     fun insert_new seen name =
   576       not (Symtab.defined seen name) ? insert (op =) name
   577     fun find_maxes _ (maxs, []) = map snd maxs
   578       | find_maxes seen (maxs, new :: news) =
   579         find_maxes
   580             (seen |> num_keys (Graph.imm_succs fact_G new) > 1
   581                      ? Symtab.default (new, ()))
   582             (if Symtab.defined tab new then
   583                let
   584                  val newp = Graph.all_preds fact_G [new]
   585                  fun is_ancestor x yp = member (op =) yp x
   586                  val maxs =
   587                    maxs |> filter (fn (_, max) => not (is_ancestor max newp))
   588                in
   589                  if exists (is_ancestor new o fst) maxs then
   590                    (maxs, news)
   591                  else
   592                    ((newp, new)
   593                     :: filter_out (fn (_, max) => is_ancestor max newp) maxs,
   594                     news)
   595                end
   596              else
   597                (maxs, Graph.Keys.fold (insert_new seen)
   598                                       (Graph.imm_preds fact_G new) news))
   599   in find_maxes Symtab.empty ([], Graph.maximals fact_G) end
   600 
   601 (* Generate more suggestions than requested, because some might be thrown out
   602    later for various reasons and "meshing" gives better results with some
   603    slack. *)
   604 fun max_suggs_of max_facts = max_facts + Int.min (50, max_facts)
   605 
   606 fun is_fact_in_graph fact_G (_, th) =
   607   can (Graph.get_node fact_G) (nickname_of th)
   608 
   609 fun interleave [] ys = ys
   610   | interleave xs [] = xs
   611   | interleave (x :: xs) (y :: ys) = x :: y :: interleave xs ys
   612 
   613 fun mash_suggested_facts ctxt ({overlord, ...} : params) prover max_facts hyp_ts
   614                          concl_t facts =
   615   let
   616     val thy = Proof_Context.theory_of ctxt
   617     val (fact_G, suggs) =
   618       mash_peek ctxt (fn {fact_G} =>
   619           if Graph.is_empty fact_G then
   620             (fact_G, [])
   621           else
   622             let
   623               val parents = maximal_in_graph fact_G facts
   624               val feats =
   625                 features_of ctxt prover thy (Local, General) (concl_t :: hyp_ts)
   626             in
   627               (fact_G, mash_QUERY ctxt overlord (max_suggs_of max_facts)
   628                                   (parents, feats))
   629             end)
   630     val sels =
   631       facts |> suggested_facts suggs
   632             (* The weights currently returned by "mash.py" are too extreme to
   633                make any sense. *)
   634             |> map fst
   635     val (unk_global, unk_local) =
   636       facts |> filter_out (is_fact_in_graph fact_G)
   637             |> List.partition (fn ((_, (loc, _)), _) => loc = Global)
   638   in (interleave unk_local sels |> weight_mepo_facts, unk_global) end
   639 
   640 fun add_to_fact_graph ctxt (name, parents, feats, deps) (adds, graph) =
   641   let
   642     fun maybe_add_from from (accum as (parents, graph)) =
   643       try_graph ctxt "updating graph" accum (fn () =>
   644           (from :: parents, Graph.add_edge_acyclic (from, name) graph))
   645     val graph = graph |> Graph.default_node (name, ())
   646     val (parents, graph) = ([], graph) |> fold maybe_add_from parents
   647     val (deps, _) = ([], graph) |> fold maybe_add_from deps
   648   in ((name, parents, feats, deps) :: adds, graph) end
   649 
   650 val learn_timeout_slack = 2.0
   651 
   652 fun launch_thread timeout task =
   653   let
   654     val hard_timeout = time_mult learn_timeout_slack timeout
   655     val birth_time = Time.now ()
   656     val death_time = Time.+ (birth_time, hard_timeout)
   657     val desc = ("machine learner for Sledgehammer", "")
   658   in Async_Manager.launch MaShN birth_time death_time desc task end
   659 
   660 fun freshish_name () =
   661   Date.fmt ".%Y_%m_%d_%H_%M_%S__" (Date.fromTimeLocal (Time.now ())) ^
   662   serial_string ()
   663 
   664 fun mash_learn_proof ctxt ({overlord, timeout, ...} : params) prover t facts
   665                      used_ths =
   666   if is_smt_prover ctxt prover then
   667     ()
   668   else
   669     launch_thread timeout (fn () =>
   670         let
   671           val thy = Proof_Context.theory_of ctxt
   672           val name = freshish_name ()
   673           val feats = features_of ctxt prover thy (Local, General) [t]
   674           val deps = used_ths |> map nickname_of
   675         in
   676           mash_peek ctxt (fn {fact_G} =>
   677               let val parents = maximal_in_graph fact_G facts in
   678                 mash_ADD ctxt overlord [(name, parents, feats, deps)]
   679               end);
   680           (true, "")
   681         end)
   682 
   683 fun sendback sub =
   684   Markup.markup Isabelle_Markup.sendback (sledgehammerN ^ " " ^ sub)
   685 
   686 val commit_timeout = seconds 30.0
   687 
   688 (* The timeout is understood in a very slack fashion. *)
   689 fun mash_learn_facts ctxt (params as {debug, verbose, overlord, ...}) prover
   690                      auto_level atp learn_timeout facts =
   691   let
   692     val timer = Timer.startRealTimer ()
   693     fun next_commit_time () =
   694       Time.+ (Timer.checkRealTimer timer, commit_timeout)
   695     val {fact_G} = mash_get ctxt
   696     val (old_facts, new_facts) =
   697       facts |> List.partition (is_fact_in_graph fact_G)
   698             ||> sort (thm_ord o pairself snd)
   699   in
   700     if null new_facts andalso (not atp orelse null old_facts) then
   701       if auto_level < 2 then
   702         "No new " ^ (if atp then "ATP" else "Isar") ^ " proofs to learn." ^
   703         (if auto_level = 0 andalso not atp then
   704            "\n\nHint: Try " ^ sendback learn_atpN ^ " to learn from ATP proofs."
   705          else
   706            "")
   707       else
   708         ""
   709     else
   710       let
   711         val all_names =
   712           facts |> map snd
   713                 |> filter_out is_likely_tautology_or_too_meta
   714                 |> map (rpair () o nickname_of)
   715                 |> Symtab.make
   716         val deps_of =
   717           if atp then
   718             atp_dependencies_of ctxt params prover auto_level facts all_names
   719           else
   720             isar_dependencies_of all_names
   721         fun do_commit [] [] state = state
   722           | do_commit adds reps {fact_G} =
   723             let
   724               val (adds, fact_G) =
   725                 ([], fact_G) |> fold (add_to_fact_graph ctxt) adds
   726             in
   727               mash_ADD ctxt overlord (rev adds);
   728               mash_REPROVE ctxt overlord reps;
   729               {fact_G = fact_G}
   730             end
   731         fun commit last adds reps =
   732           (if debug andalso auto_level = 0 then
   733              Output.urgent_message "Committing..."
   734            else
   735              ();
   736            mash_map ctxt (do_commit (rev adds) reps);
   737            if not last andalso auto_level = 0 then
   738              let val num_proofs = length adds + length reps in
   739                "Learned " ^ string_of_int num_proofs ^ " " ^
   740                (if atp then "ATP" else "Isar") ^ " proof" ^
   741                plural_s num_proofs ^ " in the last " ^
   742                string_from_time commit_timeout ^ "."
   743                |> Output.urgent_message
   744              end
   745            else
   746              ())
   747         fun learn_new_fact _ (accum as (_, (_, _, _, true))) = accum
   748           | learn_new_fact ((_, stature), th)
   749                            (adds, (parents, n, next_commit, _)) =
   750             let
   751               val name = nickname_of th
   752               val feats =
   753                 features_of ctxt prover (theory_of_thm th) stature [prop_of th]
   754               val deps = deps_of th |> these
   755               val n = n |> not (null deps) ? Integer.add 1
   756               val adds = (name, parents, feats, deps) :: adds
   757               val (adds, next_commit) =
   758                 if Time.> (Timer.checkRealTimer timer, next_commit) then
   759                   (commit false adds []; ([], next_commit_time ()))
   760                 else
   761                   (adds, next_commit)
   762               val timed_out = Time.> (Timer.checkRealTimer timer, learn_timeout)
   763             in (adds, ([name], n, next_commit, timed_out)) end
   764         val n =
   765           if null new_facts then
   766             0
   767           else
   768             let
   769               val last_th = new_facts |> List.last |> snd
   770               (* crude approximation *)
   771               val ancestors =
   772                 old_facts
   773                 |> filter (fn (_, th) => thm_ord (th, last_th) <> GREATER)
   774               val parents = maximal_in_graph fact_G ancestors
   775               val (adds, (_, n, _, _)) =
   776                 ([], (parents, 0, next_commit_time (), false))
   777                 |> fold learn_new_fact new_facts
   778             in commit true adds []; n end
   779         fun relearn_old_fact _ (accum as (_, (_, _, true))) = accum
   780           | relearn_old_fact (_, th) (reps, (n, next_commit, _)) =
   781             let
   782               val name = nickname_of th
   783               val (n, reps) =
   784                 case deps_of th of
   785                   SOME deps => (n + 1, (name, deps) :: reps)
   786                 | NONE => (n, reps)
   787               val (reps, next_commit) =
   788                 if Time.> (Timer.checkRealTimer timer, next_commit) then
   789                   (commit false [] reps; ([], next_commit_time ()))
   790                 else
   791                   (reps, next_commit)
   792               val timed_out = Time.> (Timer.checkRealTimer timer, learn_timeout)
   793             in (reps, (n, next_commit, timed_out)) end
   794         val n =
   795           if not atp orelse null old_facts then
   796             n
   797           else
   798             let
   799               fun priority_of (_, th) =
   800                 random_range 0 (1000 * max_dependencies)
   801                 - 500 * (th |> isar_dependencies_of all_names
   802                             |> Option.map length
   803                             |> the_default max_dependencies)
   804               val old_facts =
   805                 old_facts |> map (`priority_of)
   806                           |> sort (int_ord o pairself fst)
   807                           |> map snd
   808               val (reps, (n, _, _)) =
   809                 ([], (n, next_commit_time (), false))
   810                 |> fold relearn_old_fact old_facts
   811             in commit true [] reps; n end
   812       in
   813         if verbose orelse auto_level < 2 then
   814           "Learned " ^ string_of_int n ^ " nontrivial " ^
   815           (if atp then "ATP" else "Isar") ^ " proof" ^ plural_s n ^
   816           (if verbose then
   817              " in " ^ string_from_time (Timer.checkRealTimer timer)
   818            else
   819              "") ^ "."
   820         else
   821           ""
   822       end
   823   end
   824 
   825 fun mash_learn ctxt (params as {provers, timeout, ...}) fact_override chained
   826                atp =
   827   let
   828     val css = Sledgehammer_Fact.clasimpset_rule_table_of ctxt
   829     val ctxt = ctxt |> Config.put instantiate_inducts false
   830     val facts =
   831       nearly_all_facts ctxt false fact_override Symtab.empty css chained []
   832                        @{prop True}
   833     val num_facts = length facts
   834     val prover = hd provers
   835     fun learn auto_level atp =
   836       mash_learn_facts ctxt params prover auto_level atp infinite_timeout facts
   837       |> Output.urgent_message
   838   in
   839     (if atp then
   840        ("MaShing through " ^ string_of_int num_facts ^ " fact" ^
   841         plural_s num_facts ^ " for ATP proofs (" ^ quote prover ^ " timeout: " ^
   842         string_from_time timeout ^ ").\n\nCollecting Isar proofs first..."
   843         |> Output.urgent_message;
   844         learn 1 false;
   845         "Now collecting ATP proofs. This may take several hours. You can \
   846         \safely stop the learning process at any point."
   847         |> Output.urgent_message;
   848         learn 0 true)
   849      else
   850        ("MaShing through " ^ string_of_int num_facts ^ " fact" ^
   851         plural_s num_facts ^ " for Isar proofs..."
   852         |> Output.urgent_message;
   853         learn 0 false))
   854   end
   855 
   856 (* The threshold should be large enough so that MaSh doesn't kick in for Auto
   857    Sledgehammer and Try. *)
   858 val min_secs_for_learning = 15
   859 
   860 fun relevant_facts ctxt (params as {learn, fact_filter, timeout, ...}) prover
   861         max_facts ({add, only, ...} : fact_override) hyp_ts concl_t facts =
   862   if not (subset (op =) (the_list fact_filter, fact_filters)) then
   863     error ("Unknown fact filter: " ^ quote (the fact_filter) ^ ".")
   864   else if only then
   865     facts
   866   else if max_facts <= 0 orelse null facts then
   867     []
   868   else
   869     let
   870       fun maybe_learn () =
   871         if learn andalso not (Async_Manager.has_running_threads MaShN) andalso
   872            Time.toSeconds timeout >= min_secs_for_learning then
   873           let val timeout = time_mult learn_timeout_slack timeout in
   874             launch_thread timeout
   875                 (fn () => (true, mash_learn_facts ctxt params prover 2 false
   876                                                   timeout facts))
   877           end
   878         else
   879           ()
   880       val fact_filter =
   881         case fact_filter of
   882           SOME ff => (() |> ff <> mepoN ? maybe_learn; ff)
   883         | NONE =>
   884           if is_smt_prover ctxt prover then
   885             mepoN
   886           else if mash_could_suggest_facts () then
   887             (maybe_learn ();
   888              if mash_can_suggest_facts ctxt then meshN else mepoN)
   889           else
   890             mepoN
   891       val add_ths = Attrib.eval_thms ctxt add
   892       fun prepend_facts ths accepts =
   893         ((facts |> filter (member Thm.eq_thm_prop ths o snd)) @
   894          (accepts |> filter_out (member Thm.eq_thm_prop ths o snd)))
   895         |> take max_facts
   896       fun mepo () =
   897         facts |> mepo_suggested_facts ctxt params prover max_facts NONE hyp_ts
   898                                       concl_t
   899               |> weight_mepo_facts
   900       fun mash () =
   901         mash_suggested_facts ctxt params prover max_facts hyp_ts concl_t facts
   902       val mess =
   903         [] |> (if fact_filter <> mashN then cons (mepo (), []) else I)
   904            |> (if fact_filter <> mepoN then cons (mash ()) else I)
   905     in
   906       mesh_facts max_facts mess
   907       |> not (null add_ths) ? prepend_facts add_ths
   908     end
   909 
   910 fun kill_learners () = Async_Manager.kill_threads MaShN "learner"
   911 fun running_learners () = Async_Manager.running_threads MaShN "learner"
   912 
   913 end;