src/HOL/Tools/Sledgehammer/sledgehammer_mash.ML
author blanchet
Mon, 23 Jul 2012 15:32:30 +0200
changeset 49448 9e9b6e363859
parent 49423 5493e67982ee
child 49449 aaaec69db3db
permissions -rw-r--r--
don't relearn old facts in Isar mode
     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     -> ('a * thm) list -> (('a * thm) * real) list * ('a * thm) 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_get ctxt =
   555   Synchronized.change_result global_state (load ctxt #> `snd)
   556 
   557 fun mash_unlearn ctxt =
   558   Synchronized.change global_state (fn _ =>
   559       (mash_CLEAR ctxt; wipe_out (mash_state_path ()); (true, empty_state)))
   560 
   561 end
   562 
   563 fun mash_could_suggest_facts () = mash_home () <> ""
   564 fun mash_can_suggest_facts ctxt = not (Graph.is_empty (#fact_G (mash_get ctxt)))
   565 
   566 fun num_keys keys = Graph.Keys.fold (K (Integer.add 1)) keys 0
   567 
   568 fun maximal_in_graph fact_G facts =
   569   let
   570     val facts = [] |> fold (cons o nickname_of o snd) facts
   571     val tab = Symtab.empty |> fold (fn name => Symtab.default (name, ())) facts
   572     fun insert_new seen name =
   573       not (Symtab.defined seen name) ? insert (op =) name
   574     fun find_maxes _ (maxs, []) = map snd maxs
   575       | find_maxes seen (maxs, new :: news) =
   576         find_maxes
   577             (seen |> num_keys (Graph.imm_succs fact_G new) > 1
   578                      ? Symtab.default (new, ()))
   579             (if Symtab.defined tab new then
   580                let
   581                  val newp = Graph.all_preds fact_G [new]
   582                  fun is_ancestor x yp = member (op =) yp x
   583                  val maxs =
   584                    maxs |> filter (fn (_, max) => not (is_ancestor max newp))
   585                in
   586                  if exists (is_ancestor new o fst) maxs then
   587                    (maxs, news)
   588                  else
   589                    ((newp, new)
   590                     :: filter_out (fn (_, max) => is_ancestor max newp) maxs,
   591                     news)
   592                end
   593              else
   594                (maxs, Graph.Keys.fold (insert_new seen)
   595                                       (Graph.imm_preds fact_G new) news))
   596   in find_maxes Symtab.empty ([], Graph.maximals fact_G) end
   597 
   598 (* Generate more suggestions than requested, because some might be thrown out
   599    later for various reasons and "meshing" gives better results with some
   600    slack. *)
   601 fun max_suggs_of max_facts = max_facts + Int.min (200, max_facts)
   602 
   603 fun is_fact_in_graph fact_G (_, th) =
   604   can (Graph.get_node fact_G) (nickname_of th)
   605 
   606 fun mash_suggested_facts ctxt ({overlord, ...} : params) prover max_facts hyp_ts
   607                          concl_t facts =
   608   let
   609     val thy = Proof_Context.theory_of ctxt
   610     val fact_G = #fact_G (mash_get ctxt)
   611     val parents = maximal_in_graph fact_G facts
   612     val feats = features_of ctxt prover thy (Local, General) (concl_t :: hyp_ts)
   613     val suggs =
   614       if Graph.is_empty fact_G then []
   615       else mash_QUERY ctxt overlord (max_suggs_of max_facts) (parents, feats)
   616     val selected =
   617       facts |> suggested_facts suggs
   618             (* The weights currently returned by "mash.py" are too extreme to
   619                make any sense. *)
   620             |> map fst |> weight_mepo_facts
   621     val unknown = facts |> filter_out (is_fact_in_graph fact_G)
   622   in (selected, unknown) end
   623 
   624 fun add_to_fact_graph ctxt (name, parents, feats, deps) (adds, graph) =
   625   let
   626     fun maybe_add_from from (accum as (parents, graph)) =
   627       try_graph ctxt "updating graph" accum (fn () =>
   628           (from :: parents, Graph.add_edge_acyclic (from, name) graph))
   629     val graph = graph |> Graph.default_node (name, ())
   630     val (parents, graph) = ([], graph) |> fold maybe_add_from parents
   631     val (deps, _) = ([], graph) |> fold maybe_add_from deps
   632   in ((name, parents, feats, deps) :: adds, graph) end
   633 
   634 val learn_timeout_slack = 2.0
   635 
   636 fun launch_thread timeout task =
   637   let
   638     val hard_timeout = time_mult learn_timeout_slack timeout
   639     val birth_time = Time.now ()
   640     val death_time = Time.+ (birth_time, hard_timeout)
   641     val desc = ("machine learner for Sledgehammer", "")
   642   in Async_Manager.launch MaShN birth_time death_time desc task end
   643 
   644 fun freshish_name () =
   645   Date.fmt ".%Y_%m_%d_%H_%M_%S__" (Date.fromTimeLocal (Time.now ())) ^
   646   serial_string ()
   647 
   648 fun mash_learn_proof ctxt ({overlord, timeout, ...} : params) prover t facts
   649                      used_ths =
   650   if is_smt_prover ctxt prover then
   651     ()
   652   else
   653     launch_thread timeout (fn () =>
   654         let
   655           val thy = Proof_Context.theory_of ctxt
   656           val name = freshish_name ()
   657           val feats = features_of ctxt prover thy (Local, General) [t]
   658           val deps = used_ths |> map nickname_of
   659           val {fact_G} = mash_get ctxt
   660           val parents = maximal_in_graph fact_G facts
   661         in
   662           mash_ADD ctxt overlord [(name, parents, feats, deps)]; (true, "")
   663         end)
   664 
   665 fun sendback sub =
   666   Markup.markup Isabelle_Markup.sendback (sledgehammerN ^ " " ^ sub)
   667 
   668 val commit_timeout = seconds 30.0
   669 
   670 (* The timeout is understood in a very slack fashion. *)
   671 fun mash_learn_facts ctxt (params as {debug, verbose, overlord, ...}) prover
   672                      auto_level atp learn_timeout facts =
   673   let
   674     val timer = Timer.startRealTimer ()
   675     fun next_commit_time () =
   676       Time.+ (Timer.checkRealTimer timer, commit_timeout)
   677     val {fact_G} = mash_get ctxt
   678     val (old_facts, new_facts) =
   679       facts |> List.partition (is_fact_in_graph fact_G)
   680             ||> sort (thm_ord o pairself snd)
   681   in
   682     if null new_facts andalso (not atp orelse null old_facts) then
   683       if auto_level < 2 then
   684         "No new " ^ (if atp then "ATP" else "Isar") ^ " proofs to learn." ^
   685         (if auto_level = 0 andalso not atp then
   686            "\n\nHint: Try " ^ sendback learn_atpN ^ " to learn from ATP proofs."
   687          else
   688            "")
   689       else
   690         ""
   691     else
   692       let
   693         val all_names =
   694           facts |> map snd
   695                 |> filter_out is_likely_tautology_or_too_meta
   696                 |> map (rpair () o nickname_of)
   697                 |> Symtab.make
   698         val deps_of =
   699           if atp then
   700             atp_dependencies_of ctxt params prover auto_level facts all_names
   701           else
   702             isar_dependencies_of all_names
   703         fun do_commit [] [] state = state
   704           | do_commit adds reps {fact_G} =
   705             let
   706               val (adds, fact_G) =
   707                 ([], fact_G) |> fold (add_to_fact_graph ctxt) adds
   708             in
   709               mash_ADD ctxt overlord (rev adds);
   710               mash_REPROVE ctxt overlord reps;
   711               {fact_G = fact_G}
   712             end
   713         fun commit last adds reps =
   714           (if debug andalso auto_level = 0 then
   715              Output.urgent_message "Committing..."
   716            else
   717              ();
   718            mash_map ctxt (do_commit (rev adds) reps);
   719            if not last andalso auto_level = 0 then
   720              let val num_proofs = length adds + length reps in
   721                "Learned " ^ string_of_int num_proofs ^ " " ^
   722                (if atp then "ATP" else "Isar") ^ " proof" ^
   723                plural_s num_proofs ^ " in the last " ^
   724                string_from_time commit_timeout ^ "."
   725                |> Output.urgent_message
   726              end
   727            else
   728              ())
   729         fun learn_new_fact _ (accum as (_, (_, _, _, true))) = accum
   730           | learn_new_fact ((_, stature), th)
   731                            (adds, (parents, n, next_commit, _)) =
   732             let
   733               val name = nickname_of th
   734               val feats =
   735                 features_of ctxt prover (theory_of_thm th) stature [prop_of th]
   736               val deps = deps_of th |> these
   737               val n = n |> not (null deps) ? Integer.add 1
   738               val adds = (name, parents, feats, deps) :: adds
   739               val (adds, next_commit) =
   740                 if Time.> (Timer.checkRealTimer timer, next_commit) then
   741                   (commit false adds []; ([], next_commit_time ()))
   742                 else
   743                   (adds, next_commit)
   744               val timed_out = Time.> (Timer.checkRealTimer timer, learn_timeout)
   745             in (adds, ([name], n, next_commit, timed_out)) end
   746         val n =
   747           if null new_facts then
   748             0
   749           else
   750             let
   751               val last_th = new_facts |> List.last |> snd
   752               (* crude approximation *)
   753               val ancestors =
   754                 old_facts
   755                 |> filter (fn (_, th) => thm_ord (th, last_th) <> GREATER)
   756               val parents = maximal_in_graph fact_G ancestors
   757               val (adds, (_, n, _, _)) =
   758                 ([], (parents, 0, next_commit_time (), false))
   759                 |> fold learn_new_fact new_facts
   760             in commit true adds []; n end
   761         fun relearn_old_fact _ (accum as (_, (_, _, true))) = accum
   762           | relearn_old_fact (_, th) (reps, (n, next_commit, _)) =
   763             let
   764               val name = nickname_of th
   765               val (n, reps) =
   766                 case deps_of th of
   767                   SOME deps => (n + 1, (name, deps) :: reps)
   768                 | NONE => (n, reps)
   769               val (reps, next_commit) =
   770                 if Time.> (Timer.checkRealTimer timer, next_commit) then
   771                   (commit false [] reps; ([], next_commit_time ()))
   772                 else
   773                   (reps, next_commit)
   774               val timed_out = Time.> (Timer.checkRealTimer timer, learn_timeout)
   775             in (reps, (n, next_commit, timed_out)) end
   776         val n =
   777           if not atp orelse null old_facts then
   778             n
   779           else
   780             let
   781               fun priority_of (_, th) =
   782                 random_range 0 (1000 * max_dependencies)
   783                 - 500 * (th |> isar_dependencies_of all_names
   784                             |> Option.map length
   785                             |> the_default max_dependencies)
   786               val old_facts =
   787                 old_facts |> map (`priority_of)
   788                           |> sort (int_ord o pairself fst)
   789                           |> map snd
   790               val (reps, (n, _, _)) =
   791                 ([], (n, next_commit_time (), false))
   792                 |> fold relearn_old_fact old_facts
   793             in commit true [] reps; n end
   794       in
   795         if verbose orelse auto_level < 2 then
   796           "Learned " ^ string_of_int n ^ " nontrivial " ^
   797           (if atp then "ATP" else "Isar") ^ " proof" ^ plural_s n ^
   798           (if verbose then
   799              " in " ^ string_from_time (Timer.checkRealTimer timer)
   800            else
   801              "") ^ "."
   802         else
   803           ""
   804       end
   805   end
   806 
   807 fun mash_learn ctxt (params as {provers, timeout, ...}) fact_override chained
   808                atp =
   809   let
   810     val css = Sledgehammer_Fact.clasimpset_rule_table_of ctxt
   811     val ctxt = ctxt |> Config.put instantiate_inducts false
   812     val facts =
   813       nearly_all_facts ctxt false fact_override Symtab.empty css chained []
   814                        @{prop True}
   815     val num_facts = length facts
   816     val prover = hd provers
   817     fun learn auto_level atp =
   818       mash_learn_facts ctxt params prover auto_level atp infinite_timeout facts
   819       |> Output.urgent_message
   820   in
   821     (if atp then
   822        ("MaShing through " ^ string_of_int num_facts ^ " fact" ^
   823         plural_s num_facts ^ " for ATP proofs (" ^ quote prover ^ " timeout: " ^
   824         string_from_time timeout ^ ").\n\nCollecting Isar proofs first..."
   825         |> Output.urgent_message;
   826         learn 1 false;
   827         "Now collecting ATP proofs. This may take several hours. You can \
   828         \safely stop the learning process at any point."
   829         |> Output.urgent_message;
   830         learn 0 true)
   831      else
   832        ("MaShing through " ^ string_of_int num_facts ^ " fact" ^
   833         plural_s num_facts ^ " for Isar proofs..."
   834         |> Output.urgent_message;
   835         learn 0 false))
   836   end
   837 
   838 (* The threshold should be large enough so that MaSh doesn't kick in for Auto
   839    Sledgehammer and Try. *)
   840 val min_secs_for_learning = 15
   841 
   842 fun relevant_facts ctxt (params as {learn, fact_filter, timeout, ...}) prover
   843         max_facts ({add, only, ...} : fact_override) hyp_ts concl_t facts =
   844   if not (subset (op =) (the_list fact_filter, fact_filters)) then
   845     error ("Unknown fact filter: " ^ quote (the fact_filter) ^ ".")
   846   else if only then
   847     facts
   848   else if max_facts <= 0 orelse null facts then
   849     []
   850   else
   851     let
   852       fun maybe_learn () =
   853         if learn andalso not (Async_Manager.has_running_threads MaShN) andalso
   854            Time.toSeconds timeout >= min_secs_for_learning then
   855           let val timeout = time_mult learn_timeout_slack timeout in
   856             launch_thread timeout
   857                 (fn () => (true, mash_learn_facts ctxt params prover 2 false
   858                                                   timeout facts))
   859           end
   860         else
   861           ()
   862       val fact_filter =
   863         case fact_filter of
   864           SOME ff => (() |> ff <> mepoN ? maybe_learn; ff)
   865         | NONE =>
   866           if is_smt_prover ctxt prover then
   867             mepoN
   868           else if mash_could_suggest_facts () then
   869             (maybe_learn ();
   870              if mash_can_suggest_facts ctxt then meshN else mepoN)
   871           else
   872             mepoN
   873       val add_ths = Attrib.eval_thms ctxt add
   874       fun prepend_facts ths accepts =
   875         ((facts |> filter (member Thm.eq_thm_prop ths o snd)) @
   876          (accepts |> filter_out (member Thm.eq_thm_prop ths o snd)))
   877         |> take max_facts
   878       fun mepo () =
   879         facts |> mepo_suggested_facts ctxt params prover max_facts NONE hyp_ts
   880                                       concl_t
   881               |> weight_mepo_facts
   882       fun mash () =
   883         mash_suggested_facts ctxt params prover max_facts hyp_ts concl_t facts
   884       val mess =
   885         [] |> (if fact_filter <> mashN then cons (mepo (), []) else I)
   886            |> (if fact_filter <> mepoN then cons (mash ()) else I)
   887     in
   888       mesh_facts max_facts mess
   889       |> not (null add_ths) ? prepend_facts add_ths
   890     end
   891 
   892 fun kill_learners () = Async_Manager.kill_threads MaShN "learner"
   893 fun running_learners () = Async_Manager.running_threads MaShN "learner"
   894 
   895 end;