src/Pure/proofterm.ML
author wenzelm
Tue Jul 04 19:49:51 2006 +0200 (2006-07-04)
changeset 20000 2fa767ab91aa
parent 19908 f035261fb5b9
child 20071 8f3e1ddb50e6
permissions -rw-r--r--
added map_proof_terms_option;
tuned generalize, instantiate;
berghofe@11519
     1
(*  Title:      Pure/proofterm.ML
berghofe@11519
     2
    ID:         $Id$
wenzelm@11540
     3
    Author:     Stefan Berghofer, TU Muenchen
berghofe@11519
     4
wenzelm@11540
     5
LF style proof terms.
berghofe@11519
     6
*)
berghofe@11519
     7
berghofe@11615
     8
infix 8 % %% %>;
berghofe@11519
     9
berghofe@11519
    10
signature BASIC_PROOFTERM =
berghofe@11519
    11
sig
wenzelm@11543
    12
  val proofs: int ref
berghofe@11519
    13
berghofe@11519
    14
  datatype proof =
berghofe@11519
    15
     PBound of int
berghofe@11519
    16
   | Abst of string * typ option * proof
berghofe@11519
    17
   | AbsP of string * term option * proof
wenzelm@17756
    18
   | % of proof * term option
wenzelm@17756
    19
   | %% of proof * proof
berghofe@11519
    20
   | Hyp of term
berghofe@11519
    21
   | PThm of (string * (string * string list) list) * proof * term * typ list option
berghofe@11519
    22
   | PAxm of string * term * typ list option
berghofe@11519
    23
   | Oracle of string * term * typ list option
berghofe@17017
    24
   | MinProof of ((string * term) * proof) list * (string * term) list * (string * term) list;
berghofe@11519
    25
berghofe@11615
    26
  val %> : proof * term -> proof
berghofe@11519
    27
end;
berghofe@11519
    28
berghofe@11519
    29
signature PROOFTERM =
berghofe@11519
    30
sig
berghofe@11519
    31
  include BASIC_PROOFTERM
berghofe@11519
    32
berghofe@11519
    33
  val infer_derivs : (proof -> proof -> proof) -> bool * proof -> bool * proof -> bool * proof
berghofe@11519
    34
  val infer_derivs' : (proof -> proof) -> (bool * proof -> bool * proof)
berghofe@11519
    35
berghofe@11519
    36
  (** primitive operations **)
berghofe@11519
    37
  val proof_combt : proof * term list -> proof
berghofe@11519
    38
  val proof_combt' : proof * term option list -> proof
berghofe@11519
    39
  val proof_combP : proof * proof list -> proof
berghofe@11519
    40
  val strip_combt : proof -> proof * term option list
berghofe@11519
    41
  val strip_combP : proof -> proof * proof list
berghofe@11519
    42
  val strip_thm : proof -> proof
wenzelm@20000
    43
  val map_proof_terms_option : (term -> term option) -> (typ -> typ option) -> proof -> proof
berghofe@11519
    44
  val map_proof_terms : (term -> term) -> (typ -> typ) -> proof -> proof
berghofe@11519
    45
  val fold_proof_terms : (term * 'a -> 'a) -> (typ * 'a -> 'a) -> 'a * proof -> 'a
berghofe@11519
    46
  val add_prf_names : string list * proof -> string list
berghofe@11519
    47
  val add_prf_tfree_names : string list * proof -> string list
berghofe@11519
    48
  val add_prf_tvar_ixns : indexname list * proof -> indexname list
berghofe@12868
    49
  val maxidx_of_proof : proof -> int
berghofe@13744
    50
  val size_of_proof : proof -> int
berghofe@12907
    51
  val change_type : typ list option -> proof -> proof
berghofe@11519
    52
  val prf_abstract_over : term -> proof -> proof
berghofe@11519
    53
  val prf_incr_bv : int -> int -> int -> int -> proof -> proof
berghofe@11519
    54
  val incr_pboundvars : int -> int -> proof -> proof
berghofe@11519
    55
  val prf_loose_bvar1 : proof -> int -> bool
berghofe@11519
    56
  val prf_loose_Pbvar1 : proof -> int -> bool
berghofe@12279
    57
  val prf_add_loose_bnos : int -> int -> proof ->
berghofe@12279
    58
    int list * int list -> int list * int list
berghofe@11519
    59
  val norm_proof : Envir.env -> proof -> proof
berghofe@11519
    60
  val norm_proof' : Envir.env -> proof -> proof
berghofe@11519
    61
  val prf_subst_bounds : term list -> proof -> proof
berghofe@11519
    62
  val prf_subst_pbounds : proof list -> proof -> proof
berghofe@11519
    63
  val freeze_thaw_prf : proof -> proof * (proof -> proof)
berghofe@17017
    64
  val proof_of_min_axm : string * term -> proof
berghofe@17017
    65
  val proof_of_min_thm : (string * term) * proof -> proof
berghofe@11519
    66
berghofe@17017
    67
  val thms_of_proof : proof -> (term * proof) list Symtab.table ->
berghofe@11519
    68
    (term * proof) list Symtab.table
berghofe@19357
    69
  val thms_of_proof' : proof -> (term * proof) list Symtab.table ->
berghofe@19357
    70
    (term * proof) list Symtab.table
berghofe@17017
    71
  val axms_of_proof : proof -> proof Symtab.table -> proof Symtab.table
berghofe@17017
    72
  val oracles_of_proof : (string * term) list -> proof -> (string * term) list
berghofe@11519
    73
berghofe@11519
    74
  (** proof terms for specific inference rules **)
berghofe@11519
    75
  val implies_intr_proof : term -> proof -> proof
berghofe@11519
    76
  val forall_intr_proof : term -> string -> proof -> proof
berghofe@15797
    77
  val varify_proof : term -> (string * sort) list -> proof -> proof
berghofe@11519
    78
  val freezeT : term -> proof -> proof
berghofe@11519
    79
  val rotate_proof : term list -> term -> int -> proof -> proof
berghofe@11519
    80
  val permute_prems_prf : term list -> int -> int -> proof -> proof
wenzelm@19908
    81
  val generalize: string list * string list -> int -> proof -> proof
wenzelm@16880
    82
  val instantiate : ((indexname * sort) * typ) list * ((indexname * typ) * term) list
wenzelm@16880
    83
    -> proof -> proof
berghofe@11519
    84
  val lift_proof : term -> int -> term -> proof -> proof
berghofe@11519
    85
  val assumption_proof : term list -> term -> int -> proof -> proof
wenzelm@18485
    86
  val bicompose_proof : bool -> term list -> term list -> term list -> term option ->
berghofe@11519
    87
    int -> proof -> proof -> proof
berghofe@11519
    88
  val equality_axms : (string * term) list
berghofe@11519
    89
  val reflexive_axm : proof
berghofe@11519
    90
  val symmetric_axm : proof
berghofe@11519
    91
  val transitive_axm : proof
berghofe@11519
    92
  val equal_intr_axm : proof
berghofe@11519
    93
  val equal_elim_axm : proof
berghofe@11519
    94
  val abstract_rule_axm : proof
berghofe@11519
    95
  val combination_axm : proof
berghofe@11519
    96
  val reflexive : proof
berghofe@11519
    97
  val symmetric : proof -> proof
berghofe@11519
    98
  val transitive : term -> typ -> proof -> proof -> proof
berghofe@11519
    99
  val abstract_rule : term -> string -> proof -> proof
berghofe@11519
   100
  val combination : term -> term -> term -> term -> typ -> proof -> proof -> proof
berghofe@11519
   101
  val equal_intr : term -> term -> proof -> proof -> proof
berghofe@11519
   102
  val equal_elim : term -> term -> proof -> proof -> proof
berghofe@11519
   103
  val axm_proof : string -> term -> proof
berghofe@11519
   104
  val oracle_proof : string -> term -> proof
wenzelm@16458
   105
  val thm_proof : theory -> string * (string * string list) list ->
berghofe@11519
   106
    term list -> term -> proof -> proof
wenzelm@12923
   107
  val get_name_tags : term list -> term -> proof -> string * (string * string list) list
berghofe@11519
   108
berghofe@11519
   109
  (** rewriting on proof terms **)
berghofe@12233
   110
  val add_prf_rrules : (proof * proof) list -> theory -> theory
berghofe@12233
   111
  val add_prf_rprocs : (string * (Term.typ list -> proof -> proof option)) list ->
berghofe@12233
   112
    theory -> theory
wenzelm@17203
   113
  val rewrite_proof : theory -> (proof * proof) list *
berghofe@11519
   114
    (string * (typ list -> proof -> proof option)) list -> proof -> proof
berghofe@11615
   115
  val rewrite_proof_notypes : (proof * proof) list *
berghofe@11615
   116
    (string * (typ list -> proof -> proof option)) list -> proof -> proof
wenzelm@16536
   117
  val init_data: theory -> theory
berghofe@11519
   118
  
berghofe@11519
   119
end
berghofe@11519
   120
berghofe@11519
   121
structure Proofterm : PROOFTERM =
berghofe@11519
   122
struct
berghofe@11519
   123
berghofe@11715
   124
open Envir;
berghofe@11715
   125
berghofe@11519
   126
datatype proof =
berghofe@11519
   127
   PBound of int
berghofe@11519
   128
 | Abst of string * typ option * proof
berghofe@11519
   129
 | AbsP of string * term option * proof
wenzelm@12497
   130
 | op % of proof * term option
wenzelm@12497
   131
 | op %% of proof * proof
berghofe@11519
   132
 | Hyp of term
berghofe@11519
   133
 | PThm of (string * (string * string list) list) * proof * term * typ list option
berghofe@11519
   134
 | PAxm of string * term * typ list option
berghofe@11519
   135
 | Oracle of string * term * typ list option
berghofe@17017
   136
 | MinProof of ((string * term) * proof) list * (string * term) list * (string * term) list;
berghofe@11519
   137
berghofe@17017
   138
fun proof_of_min_axm (s, prop) = PAxm (s, prop, NONE);
berghofe@17017
   139
fun proof_of_min_thm ((s, prop), prf) = PThm ((s, []), prf, prop, NONE);
berghofe@17017
   140
berghofe@17017
   141
val string_term_ord = prod_ord fast_string_ord Term.fast_term_ord;
berghofe@17017
   142
berghofe@17017
   143
fun oracles_of_proof oras prf =
berghofe@11519
   144
  let
berghofe@17017
   145
    fun oras_of (Abst (_, _, prf)) = oras_of prf
berghofe@17017
   146
      | oras_of (AbsP (_, _, prf)) = oras_of prf
berghofe@17017
   147
      | oras_of (prf % _) = oras_of prf
berghofe@17017
   148
      | oras_of (prf1 %% prf2) = oras_of prf1 #> oras_of prf2
berghofe@17017
   149
      | oras_of (PThm ((name, _), prf, prop, _)) = (fn tabs as (thms, oras) =>
wenzelm@17412
   150
          case Symtab.lookup thms name of
wenzelm@17412
   151
            NONE => oras_of prf (Symtab.update (name, [prop]) thms, oras)
wenzelm@17492
   152
          | SOME ps => if member (op =) ps prop then tabs else
wenzelm@17412
   153
              oras_of prf (Symtab.update (name, prop::ps) thms, oras))
berghofe@17017
   154
      | oras_of (Oracle (s, prop, _)) =
berghofe@17017
   155
          apsnd (OrdList.insert string_term_ord (s, prop))
berghofe@17017
   156
      | oras_of (MinProof (thms, _, oras)) =
berghofe@17017
   157
          apsnd (OrdList.union string_term_ord oras) #>
berghofe@17017
   158
          fold (oras_of o proof_of_min_thm) thms
berghofe@17017
   159
      | oras_of _ = I
berghofe@11519
   160
  in
berghofe@17017
   161
    snd (oras_of prf (Symtab.empty, oras))
berghofe@11519
   162
  end;
berghofe@11519
   163
berghofe@17017
   164
fun thms_of_proof (Abst (_, _, prf)) = thms_of_proof prf
berghofe@17017
   165
  | thms_of_proof (AbsP (_, _, prf)) = thms_of_proof prf
berghofe@17017
   166
  | thms_of_proof (prf1 %% prf2) = thms_of_proof prf1 #> thms_of_proof prf2
berghofe@17017
   167
  | thms_of_proof (prf % _) = thms_of_proof prf
berghofe@17017
   168
  | thms_of_proof (prf' as PThm ((s, _), prf, prop, _)) = (fn tab =>
wenzelm@17412
   169
      case Symtab.lookup tab s of
wenzelm@17412
   170
        NONE => thms_of_proof prf (Symtab.update (s, [(prop, prf')]) tab)
wenzelm@19304
   171
      | SOME ps => if exists (fn (p, _) => p = prop) ps then tab else
wenzelm@17412
   172
          thms_of_proof prf (Symtab.update (s, (prop, prf')::ps) tab))
berghofe@17017
   173
  | thms_of_proof (MinProof (prfs, _, _)) = fold (thms_of_proof o proof_of_min_thm) prfs
berghofe@17017
   174
  | thms_of_proof _ = I;
berghofe@11519
   175
berghofe@19357
   176
(* this version does not recursively descend into proofs of (named) theorems *)
berghofe@19357
   177
fun thms_of_proof' (Abst (_, _, prf)) = thms_of_proof' prf
berghofe@19357
   178
  | thms_of_proof' (AbsP (_, _, prf)) = thms_of_proof' prf
berghofe@19357
   179
  | thms_of_proof' (prf1 %% prf2) = thms_of_proof' prf1 #> thms_of_proof' prf2
berghofe@19357
   180
  | thms_of_proof' (prf % _) = thms_of_proof' prf
berghofe@19357
   181
  | thms_of_proof' (PThm (("", _), prf, prop, _)) = thms_of_proof' prf
berghofe@19357
   182
  | thms_of_proof' (prf' as PThm ((s, _), _, prop, _)) = (fn tab =>
berghofe@19357
   183
      case Symtab.lookup tab s of
berghofe@19357
   184
        NONE => Symtab.update (s, [(prop, prf')]) tab
berghofe@19357
   185
      | SOME ps => if exists (fn (p, _) => p = prop) ps then tab else
berghofe@19357
   186
          Symtab.update (s, (prop, prf')::ps) tab)
berghofe@19357
   187
  | thms_of_proof' (MinProof (prfs, _, _)) = fold (thms_of_proof' o proof_of_min_thm) prfs
berghofe@19357
   188
  | thms_of_proof' _ = I;
berghofe@19357
   189
berghofe@17017
   190
fun axms_of_proof (Abst (_, _, prf)) = axms_of_proof prf
berghofe@17017
   191
  | axms_of_proof (AbsP (_, _, prf)) = axms_of_proof prf
berghofe@17017
   192
  | axms_of_proof (prf1 %% prf2) = axms_of_proof prf1 #> axms_of_proof prf2
berghofe@17017
   193
  | axms_of_proof (prf % _) = axms_of_proof prf
wenzelm@17412
   194
  | axms_of_proof (prf as PAxm (s, _, _)) = Symtab.update (s, prf)
berghofe@17017
   195
  | axms_of_proof (MinProof (_, prfs, _)) = fold (axms_of_proof o proof_of_min_axm) prfs
berghofe@17017
   196
  | axms_of_proof _ = I;
berghofe@11519
   197
berghofe@11519
   198
(** collect all theorems, axioms and oracles **)
berghofe@11519
   199
berghofe@17017
   200
fun map3 f g h (thms, axms, oras) = (f thms, g axms, h oras);
berghofe@17017
   201
berghofe@17017
   202
fun mk_min_proof (Abst (_, _, prf)) = mk_min_proof prf
berghofe@17017
   203
  | mk_min_proof (AbsP (_, _, prf)) = mk_min_proof prf
berghofe@17017
   204
  | mk_min_proof (prf % _) = mk_min_proof prf
berghofe@17017
   205
  | mk_min_proof (prf1 %% prf2) = mk_min_proof prf1 #> mk_min_proof prf2
berghofe@17017
   206
  | mk_min_proof (PThm ((s, _), prf, prop, _)) =
berghofe@17017
   207
      map3 (OrdList.insert (string_term_ord o pairself fst) ((s, prop), prf)) I I
berghofe@17017
   208
  | mk_min_proof (PAxm (s, prop, _)) =
berghofe@17017
   209
      map3 I (OrdList.insert string_term_ord (s, prop)) I
berghofe@17017
   210
  | mk_min_proof (Oracle (s, prop, _)) =
berghofe@17017
   211
      map3 I I (OrdList.insert string_term_ord (s, prop))
berghofe@17017
   212
  | mk_min_proof (MinProof (thms, axms, oras)) =
berghofe@17017
   213
      map3 (OrdList.union (string_term_ord o pairself fst) thms)
berghofe@17017
   214
        (OrdList.union string_term_ord axms) (OrdList.union string_term_ord oras)
berghofe@17017
   215
  | mk_min_proof _ = I;
berghofe@11519
   216
berghofe@11519
   217
(** proof objects with different levels of detail **)
berghofe@11519
   218
wenzelm@11543
   219
val proofs = ref 2;
berghofe@11519
   220
wenzelm@11543
   221
fun err_illegal_level i =
wenzelm@11543
   222
  error ("Illegal level of detail for proof objects: " ^ string_of_int i);
berghofe@11519
   223
berghofe@11519
   224
fun if_ora b = if b then oracles_of_proof else K;
wenzelm@19222
   225
val min_proof = MinProof ([], [], []);
berghofe@11519
   226
berghofe@11519
   227
fun infer_derivs f (ora1, prf1) (ora2, prf2) =
wenzelm@19222
   228
  let val ora = ora1 orelse ora2 in
wenzelm@19222
   229
    (ora,
wenzelm@19222
   230
      case !proofs of
wenzelm@19222
   231
        2 => f prf1 prf2
wenzelm@19222
   232
      | 1 => MinProof (([], [], []) |> mk_min_proof prf1 |> mk_min_proof prf2)
wenzelm@19222
   233
      | 0 => if ora then MinProof ([], [], if_ora ora2 (if_ora ora1 [] prf1) prf2) else min_proof
wenzelm@19222
   234
      | i => err_illegal_level i)
wenzelm@19222
   235
  end;
berghofe@11519
   236
wenzelm@19222
   237
fun infer_derivs' f = infer_derivs (K f) (false, min_proof);
berghofe@11519
   238
skalberg@15531
   239
fun (prf %> t) = prf % SOME t;
berghofe@11519
   240
skalberg@15570
   241
val proof_combt = Library.foldl (op %>);
skalberg@15570
   242
val proof_combt' = Library.foldl (op %);
skalberg@15570
   243
val proof_combP = Library.foldl (op %%);
berghofe@11519
   244
berghofe@11519
   245
fun strip_combt prf = 
berghofe@11615
   246
    let fun stripc (prf % t, ts) = stripc (prf, t::ts)
berghofe@11519
   247
          | stripc  x =  x 
berghofe@11519
   248
    in  stripc (prf, [])  end;
berghofe@11519
   249
berghofe@11519
   250
fun strip_combP prf = 
berghofe@11615
   251
    let fun stripc (prf %% prf', prfs) = stripc (prf, prf'::prfs)
berghofe@11519
   252
          | stripc  x =  x
berghofe@11519
   253
    in  stripc (prf, [])  end;
berghofe@11519
   254
berghofe@11519
   255
fun strip_thm prf = (case strip_combt (fst (strip_combP prf)) of
berghofe@11519
   256
      (PThm (_, prf', _, _), _) => prf'
berghofe@11519
   257
    | _ => prf);
berghofe@11519
   258
skalberg@15574
   259
val mk_Abst = foldr (fn ((s, T:typ), prf) => Abst (s, NONE, prf));
skalberg@15531
   260
fun mk_AbsP (i, prf) = funpow i (fn prf => AbsP ("H", NONE, prf)) prf;
berghofe@11519
   261
wenzelm@20000
   262
fun apsome f NONE = raise SAME
wenzelm@20000
   263
  | apsome f (SOME x) = (case f x of NONE => raise SAME | some => some);
wenzelm@20000
   264
skalberg@15531
   265
fun apsome' f NONE = raise SAME
skalberg@15531
   266
  | apsome' f (SOME x) = SOME (f x);
berghofe@11715
   267
wenzelm@20000
   268
fun map_proof_terms_option f g =
wenzelm@20000
   269
  let
wenzelm@20000
   270
    fun map_typs (T :: Ts) =
wenzelm@20000
   271
          (case g T of
wenzelm@20000
   272
            NONE => T :: map_typs Ts
wenzelm@20000
   273
          | SOME T' => T' :: (map_typs Ts handle SAME => Ts))
wenzelm@20000
   274
      | map_typs [] = raise SAME;
wenzelm@20000
   275
wenzelm@20000
   276
    fun mapp (Abst (s, T, prf)) = (Abst (s, apsome g T, mapph prf)
wenzelm@20000
   277
          handle SAME => Abst (s, T, mapp prf))
wenzelm@20000
   278
      | mapp (AbsP (s, t, prf)) = (AbsP (s, apsome f t, mapph prf)
wenzelm@20000
   279
          handle SAME => AbsP (s, t, mapp prf))
wenzelm@20000
   280
      | mapp (prf % t) = (mapp prf % (apsome f t handle SAME => t)
wenzelm@20000
   281
          handle SAME => prf % apsome f t)
wenzelm@20000
   282
      | mapp (prf1 %% prf2) = (mapp prf1 %% mapph prf2
wenzelm@20000
   283
          handle SAME => prf1 %% mapp prf2)
wenzelm@20000
   284
      | mapp (PThm (a, prf, prop, SOME Ts)) =
wenzelm@20000
   285
          PThm (a, prf, prop, SOME (map_typs Ts))
wenzelm@20000
   286
      | mapp (PAxm (a, prop, SOME Ts)) =
wenzelm@20000
   287
          PAxm (a, prop, SOME (map_typs Ts))
wenzelm@20000
   288
      | mapp _ = raise SAME
wenzelm@20000
   289
    and mapph prf = (mapp prf handle SAME => prf)
wenzelm@20000
   290
wenzelm@20000
   291
  in mapph end;
wenzelm@20000
   292
berghofe@11715
   293
fun same f x =
berghofe@11715
   294
  let val x' = f x
berghofe@11715
   295
  in if x = x' then raise SAME else x' end;
berghofe@11715
   296
berghofe@11715
   297
fun map_proof_terms f g =
wenzelm@20000
   298
  map_proof_terms_option
wenzelm@20000
   299
   (fn t => SOME (same f t) handle SAME => NONE)
wenzelm@20000
   300
   (fn T => SOME (same g T) handle SAME => NONE);
berghofe@11519
   301
skalberg@15531
   302
fun fold_proof_terms f g (a, Abst (_, SOME T, prf)) = fold_proof_terms f g (g (T, a), prf)
skalberg@15531
   303
  | fold_proof_terms f g (a, Abst (_, NONE, prf)) = fold_proof_terms f g (a, prf)
skalberg@15531
   304
  | fold_proof_terms f g (a, AbsP (_, SOME t, prf)) = fold_proof_terms f g (f (t, a), prf)
skalberg@15531
   305
  | fold_proof_terms f g (a, AbsP (_, NONE, prf)) = fold_proof_terms f g (a, prf)
skalberg@15531
   306
  | fold_proof_terms f g (a, prf % SOME t) = f (t, fold_proof_terms f g (a, prf))
skalberg@15531
   307
  | fold_proof_terms f g (a, prf % NONE) = fold_proof_terms f g (a, prf)
berghofe@11615
   308
  | fold_proof_terms f g (a, prf1 %% prf2) = fold_proof_terms f g
berghofe@11519
   309
      (fold_proof_terms f g (a, prf1), prf2)
skalberg@15574
   310
  | fold_proof_terms _ g (a, PThm (_, _, _, SOME Ts)) = foldr g a Ts
skalberg@15574
   311
  | fold_proof_terms _ g (a, PAxm (_, prop, SOME Ts)) = foldr g a Ts
berghofe@11519
   312
  | fold_proof_terms _ _ (a, _) = a;
berghofe@11519
   313
berghofe@11519
   314
val add_prf_names = fold_proof_terms add_term_names ((uncurry K) o swap);
berghofe@11519
   315
val add_prf_tfree_names = fold_proof_terms add_term_tfree_names add_typ_tfree_names;
berghofe@11519
   316
val add_prf_tvar_ixns = fold_proof_terms add_term_tvar_ixns (add_typ_ixns o swap);
berghofe@11519
   317
berghofe@12868
   318
fun maxidx_of_proof prf = fold_proof_terms
berghofe@12868
   319
  (Int.max o apfst maxidx_of_term) (Int.max o apfst maxidx_of_typ) (~1, prf); 
berghofe@12868
   320
berghofe@13744
   321
fun size_of_proof (Abst (_, _, prf)) = 1 + size_of_proof prf
berghofe@13749
   322
  | size_of_proof (AbsP (_, t, prf)) = 1 + size_of_proof prf
berghofe@13744
   323
  | size_of_proof (prf1 %% prf2) = size_of_proof prf1 + size_of_proof prf2
berghofe@13749
   324
  | size_of_proof (prf % _) = 1 + size_of_proof prf
berghofe@13744
   325
  | size_of_proof _ = 1;
berghofe@13744
   326
berghofe@12907
   327
fun change_type opTs (PThm (name, prf, prop, _)) = PThm (name, prf, prop, opTs)
berghofe@12907
   328
  | change_type opTs (PAxm (name, prop, _)) = PAxm (name, prop, opTs)
berghofe@12907
   329
  | change_type opTs (Oracle (name, prop, _)) = Oracle (name, prop, opTs)
berghofe@12907
   330
  | change_type _ prf = prf;
berghofe@12907
   331
berghofe@11519
   332
berghofe@11519
   333
(***** utilities *****)
berghofe@11519
   334
berghofe@11519
   335
fun strip_abs (_::Ts) (Abs (_, _, t)) = strip_abs Ts t
berghofe@11519
   336
  | strip_abs _ t = t;
berghofe@11519
   337
skalberg@15570
   338
fun mk_abs Ts t = Library.foldl (fn (t', T) => Abs ("", T, t')) (t, Ts);
berghofe@11519
   339
berghofe@11519
   340
berghofe@11519
   341
(*Abstraction of a proof term over its occurrences of v, 
berghofe@11519
   342
    which must contain no loose bound variables.
berghofe@11519
   343
  The resulting proof term is ready to become the body of an Abst.*)
berghofe@11519
   344
berghofe@11519
   345
fun prf_abstract_over v =
berghofe@11519
   346
  let
berghofe@11715
   347
    fun abst' lev u = if v aconv u then Bound lev else
berghofe@11715
   348
      (case u of
berghofe@11715
   349
         Abs (a, T, t) => Abs (a, T, abst' (lev + 1) t)
berghofe@11715
   350
       | f $ t => (abst' lev f $ absth' lev t handle SAME => f $ abst' lev t)
berghofe@11715
   351
       | _ => raise SAME)
berghofe@11715
   352
    and absth' lev t = (abst' lev t handle SAME => t);
berghofe@11519
   353
berghofe@11715
   354
    fun abst lev (AbsP (a, t, prf)) =
berghofe@11715
   355
          (AbsP (a, apsome' (abst' lev) t, absth lev prf)
berghofe@11715
   356
           handle SAME => AbsP (a, t, abst lev prf))
berghofe@11715
   357
      | abst lev (Abst (a, T, prf)) = Abst (a, T, abst (lev + 1) prf)
berghofe@11715
   358
      | abst lev (prf1 %% prf2) = (abst lev prf1 %% absth lev prf2
berghofe@11715
   359
          handle SAME => prf1 %% abst lev prf2)
skalberg@15570
   360
      | abst lev (prf % t) = (abst lev prf % Option.map (absth' lev) t
berghofe@11715
   361
          handle SAME => prf % apsome' (abst' lev) t)
berghofe@11715
   362
      | abst _ _ = raise SAME
berghofe@11715
   363
    and absth lev prf = (abst lev prf handle SAME => prf)
berghofe@11519
   364
berghofe@11715
   365
  in absth 0 end;
berghofe@11519
   366
berghofe@11519
   367
berghofe@11519
   368
(*increments a proof term's non-local bound variables
berghofe@11519
   369
  required when moving a proof term within abstractions
berghofe@11519
   370
     inc is  increment for bound variables
berghofe@11519
   371
     lev is  level at which a bound variable is considered 'loose'*)
berghofe@11519
   372
berghofe@11519
   373
fun incr_bv' inct tlev t = incr_bv (inct, tlev, t);
berghofe@11519
   374
berghofe@11715
   375
fun prf_incr_bv' incP inct Plev tlev (PBound i) =
berghofe@11715
   376
      if i >= Plev then PBound (i+incP) else raise SAME 
berghofe@11715
   377
  | prf_incr_bv' incP inct Plev tlev (AbsP (a, t, body)) =
berghofe@11715
   378
      (AbsP (a, apsome' (same (incr_bv' inct tlev)) t,
berghofe@11715
   379
         prf_incr_bv incP inct (Plev+1) tlev body) handle SAME =>
berghofe@11715
   380
           AbsP (a, t, prf_incr_bv' incP inct (Plev+1) tlev body))
berghofe@11715
   381
  | prf_incr_bv' incP inct Plev tlev (Abst (a, T, body)) =
berghofe@11715
   382
      Abst (a, T, prf_incr_bv' incP inct Plev (tlev+1) body)
berghofe@11715
   383
  | prf_incr_bv' incP inct Plev tlev (prf %% prf') = 
berghofe@11715
   384
      (prf_incr_bv' incP inct Plev tlev prf %% prf_incr_bv incP inct Plev tlev prf'
berghofe@11715
   385
       handle SAME => prf %% prf_incr_bv' incP inct Plev tlev prf')
berghofe@11715
   386
  | prf_incr_bv' incP inct Plev tlev (prf % t) = 
skalberg@15570
   387
      (prf_incr_bv' incP inct Plev tlev prf % Option.map (incr_bv' inct tlev) t
berghofe@11715
   388
       handle SAME => prf % apsome' (same (incr_bv' inct tlev)) t)
berghofe@11715
   389
  | prf_incr_bv' _ _ _ _ _ = raise SAME
berghofe@11715
   390
and prf_incr_bv incP inct Plev tlev prf =
berghofe@11715
   391
      (prf_incr_bv' incP inct Plev tlev prf handle SAME => prf);
berghofe@11519
   392
berghofe@11519
   393
fun incr_pboundvars  0 0 prf = prf
berghofe@11519
   394
  | incr_pboundvars incP inct prf = prf_incr_bv incP inct 0 0 prf;
berghofe@11519
   395
berghofe@11519
   396
berghofe@11615
   397
fun prf_loose_bvar1 (prf1 %% prf2) k = prf_loose_bvar1 prf1 k orelse prf_loose_bvar1 prf2 k
skalberg@15531
   398
  | prf_loose_bvar1 (prf % SOME t) k = prf_loose_bvar1 prf k orelse loose_bvar1 (t, k)
skalberg@15531
   399
  | prf_loose_bvar1 (_ % NONE) _ = true
skalberg@15531
   400
  | prf_loose_bvar1 (AbsP (_, SOME t, prf)) k = loose_bvar1 (t, k) orelse prf_loose_bvar1 prf k
skalberg@15531
   401
  | prf_loose_bvar1 (AbsP (_, NONE, _)) k = true
berghofe@11519
   402
  | prf_loose_bvar1 (Abst (_, _, prf)) k = prf_loose_bvar1 prf (k+1)
berghofe@11519
   403
  | prf_loose_bvar1 _ _ = false;
berghofe@11519
   404
berghofe@11519
   405
fun prf_loose_Pbvar1 (PBound i) k = i = k
berghofe@11615
   406
  | prf_loose_Pbvar1 (prf1 %% prf2) k = prf_loose_Pbvar1 prf1 k orelse prf_loose_Pbvar1 prf2 k
berghofe@11615
   407
  | prf_loose_Pbvar1 (prf % _) k = prf_loose_Pbvar1 prf k
berghofe@11519
   408
  | prf_loose_Pbvar1 (AbsP (_, _, prf)) k = prf_loose_Pbvar1 prf (k+1)
berghofe@11519
   409
  | prf_loose_Pbvar1 (Abst (_, _, prf)) k = prf_loose_Pbvar1 prf k
berghofe@11519
   410
  | prf_loose_Pbvar1 _ _ = false;
berghofe@11519
   411
berghofe@12279
   412
fun prf_add_loose_bnos plev tlev (PBound i) (is, js) =
wenzelm@17492
   413
      if i < plev then (is, js) else (insert (op =) (i-plev) is, js)
berghofe@12279
   414
  | prf_add_loose_bnos plev tlev (prf1 %% prf2) p =
berghofe@12279
   415
      prf_add_loose_bnos plev tlev prf2
berghofe@12279
   416
        (prf_add_loose_bnos plev tlev prf1 p)
berghofe@12279
   417
  | prf_add_loose_bnos plev tlev (prf % opt) (is, js) =
berghofe@12279
   418
      prf_add_loose_bnos plev tlev prf (case opt of
wenzelm@17492
   419
          NONE => (is, insert (op =) ~1 js)
skalberg@15531
   420
        | SOME t => (is, add_loose_bnos (t, tlev, js)))
berghofe@12279
   421
  | prf_add_loose_bnos plev tlev (AbsP (_, opt, prf)) (is, js) =
berghofe@12279
   422
      prf_add_loose_bnos (plev+1) tlev prf (case opt of
wenzelm@17492
   423
          NONE => (is, insert (op =) ~1 js)
skalberg@15531
   424
        | SOME t => (is, add_loose_bnos (t, tlev, js)))
berghofe@12279
   425
  | prf_add_loose_bnos plev tlev (Abst (_, _, prf)) p =
berghofe@12279
   426
      prf_add_loose_bnos plev (tlev+1) prf p
berghofe@12279
   427
  | prf_add_loose_bnos _ _ _ _ = ([], []);
berghofe@12279
   428
berghofe@11519
   429
berghofe@11519
   430
(**** substitutions ****)
berghofe@11519
   431
berghofe@18316
   432
fun del_conflicting_tvars envT T = Term.instantiateT
wenzelm@19482
   433
  (map_filter (fn ixnS as (_, S) =>
berghofe@18316
   434
     (Type.lookup (envT, ixnS); NONE) handle TYPE _ =>
berghofe@18316
   435
        SOME (ixnS, TFree ("'dummy", S))) (typ_tvars T)) T;
berghofe@18316
   436
berghofe@18316
   437
fun del_conflicting_vars env t = Term.instantiate
wenzelm@19482
   438
  (map_filter (fn ixnS as (_, S) =>
berghofe@18316
   439
     (Type.lookup (type_env env, ixnS); NONE) handle TYPE _ =>
berghofe@18316
   440
        SOME (ixnS, TFree ("'dummy", S))) (term_tvars t),
wenzelm@19482
   441
   map_filter (fn Var (ixnT as (_, T)) =>
berghofe@18316
   442
     (Envir.lookup (env, ixnT); NONE) handle TYPE _ =>
berghofe@18316
   443
        SOME (ixnT, Free ("dummy", T))) (term_vars t)) t;
berghofe@18316
   444
berghofe@11519
   445
fun norm_proof env =
berghofe@11519
   446
  let
wenzelm@12497
   447
    val envT = type_env env;
berghofe@18316
   448
    fun msg s = warning ("type conflict in norm_proof:\n" ^ s);
berghofe@18316
   449
    fun htype f t = f env t handle TYPE (s, _, _) =>
berghofe@18316
   450
      (msg s; f env (del_conflicting_vars env t));
berghofe@18316
   451
    fun htypeT f T = f envT T handle TYPE (s, _, _) =>
berghofe@18316
   452
      (msg s; f envT (del_conflicting_tvars envT T));
berghofe@18316
   453
    fun htypeTs f Ts = f envT Ts handle TYPE (s, _, _) =>
berghofe@18316
   454
      (msg s; f envT (map (del_conflicting_tvars envT) Ts));
berghofe@18316
   455
    fun norm (Abst (s, T, prf)) = (Abst (s, apsome' (htypeT norm_type_same) T, normh prf)
berghofe@11519
   456
          handle SAME => Abst (s, T, norm prf))
berghofe@18316
   457
      | norm (AbsP (s, t, prf)) = (AbsP (s, apsome' (htype norm_term_same) t, normh prf)
berghofe@11519
   458
          handle SAME => AbsP (s, t, norm prf))
berghofe@18316
   459
      | norm (prf % t) = (norm prf % Option.map (htype norm_term) t
berghofe@18316
   460
          handle SAME => prf % apsome' (htype norm_term_same) t)
berghofe@11615
   461
      | norm (prf1 %% prf2) = (norm prf1 %% normh prf2
berghofe@11615
   462
          handle SAME => prf1 %% norm prf2)
berghofe@18316
   463
      | norm (PThm (s, prf, t, Ts)) = PThm (s, prf, t, apsome' (htypeTs norm_types_same) Ts)
berghofe@18316
   464
      | norm (PAxm (s, prop, Ts)) = PAxm (s, prop, apsome' (htypeTs norm_types_same) Ts)
berghofe@11519
   465
      | norm _ = raise SAME
berghofe@11519
   466
    and normh prf = (norm prf handle SAME => prf);
berghofe@11519
   467
  in normh end;
berghofe@11519
   468
berghofe@11519
   469
(***** Remove some types in proof term (to save space) *****)
berghofe@11519
   470
berghofe@11519
   471
fun remove_types (Abs (s, _, t)) = Abs (s, dummyT, remove_types t)
berghofe@11519
   472
  | remove_types (t $ u) = remove_types t $ remove_types u
berghofe@11519
   473
  | remove_types (Const (s, _)) = Const (s, dummyT)
berghofe@11519
   474
  | remove_types t = t;
berghofe@11519
   475
berghofe@11519
   476
fun remove_types_env (Envir.Envir {iTs, asol, maxidx}) =
berghofe@15797
   477
  Envir.Envir {iTs = iTs, asol = Vartab.map (apsnd remove_types) asol,
berghofe@15797
   478
    maxidx = maxidx};
berghofe@11519
   479
berghofe@11519
   480
fun norm_proof' env prf = norm_proof (remove_types_env env) prf;
berghofe@11519
   481
berghofe@11519
   482
(**** substitution of bound variables ****)
berghofe@11519
   483
berghofe@11519
   484
fun prf_subst_bounds args prf =
berghofe@11519
   485
  let
berghofe@11519
   486
    val n = length args;
berghofe@11519
   487
    fun subst' lev (Bound i) =
berghofe@11519
   488
         (if i<lev then raise SAME    (*var is locally bound*)
berghofe@11519
   489
          else  incr_boundvars lev (List.nth (args, i-lev))
berghofe@11519
   490
                  handle Subscript => Bound (i-n)  (*loose: change it*))
berghofe@11519
   491
      | subst' lev (Abs (a, T, body)) = Abs (a, T,  subst' (lev+1) body)
berghofe@11519
   492
      | subst' lev (f $ t) = (subst' lev f $ substh' lev t
berghofe@11519
   493
          handle SAME => f $ subst' lev t)
berghofe@11519
   494
      | subst' _ _ = raise SAME
berghofe@11519
   495
    and substh' lev t = (subst' lev t handle SAME => t);
berghofe@11519
   496
berghofe@11519
   497
    fun subst lev (AbsP (a, t, body)) = (AbsP (a, apsome' (subst' lev) t, substh lev body)
berghofe@11519
   498
          handle SAME => AbsP (a, t, subst lev body))
berghofe@11519
   499
      | subst lev (Abst (a, T, body)) = Abst (a, T, subst (lev+1) body)
berghofe@11615
   500
      | subst lev (prf %% prf') = (subst lev prf %% substh lev prf'
berghofe@11615
   501
          handle SAME => prf %% subst lev prf')
skalberg@15570
   502
      | subst lev (prf % t) = (subst lev prf % Option.map (substh' lev) t
berghofe@11615
   503
          handle SAME => prf % apsome' (subst' lev) t)
berghofe@11519
   504
      | subst _ _ = raise SAME
berghofe@11519
   505
    and substh lev prf = (subst lev prf handle SAME => prf)
berghofe@11519
   506
  in case args of [] => prf | _ => substh 0 prf end;
berghofe@11519
   507
berghofe@11519
   508
fun prf_subst_pbounds args prf =
berghofe@11519
   509
  let
berghofe@11519
   510
    val n = length args;
berghofe@11519
   511
    fun subst (PBound i) Plev tlev =
berghofe@11519
   512
 	 (if i < Plev then raise SAME    (*var is locally bound*)
berghofe@11519
   513
          else incr_pboundvars Plev tlev (List.nth (args, i-Plev))
berghofe@11519
   514
                 handle Subscript => PBound (i-n)  (*loose: change it*))
berghofe@11519
   515
      | subst (AbsP (a, t, body)) Plev tlev = AbsP (a, t, subst body (Plev+1) tlev)
berghofe@11519
   516
      | subst (Abst (a, T, body)) Plev tlev = Abst (a, T, subst body Plev (tlev+1))
berghofe@11615
   517
      | subst (prf %% prf') Plev tlev = (subst prf Plev tlev %% substh prf' Plev tlev
berghofe@11615
   518
          handle SAME => prf %% subst prf' Plev tlev)
berghofe@11615
   519
      | subst (prf % t) Plev tlev = subst prf Plev tlev % t
berghofe@11519
   520
      | subst  prf _ _ = raise SAME
berghofe@11519
   521
    and substh prf Plev tlev = (subst prf Plev tlev handle SAME => prf)
berghofe@11519
   522
  in case args of [] => prf | _ => substh prf 0 0 end;
berghofe@11519
   523
berghofe@11519
   524
berghofe@11519
   525
(**** Freezing and thawing of variables in proof terms ****)
berghofe@11519
   526
berghofe@11519
   527
fun frzT names =
haftmann@17325
   528
  map_type_tvar (fn (ixn, xs) => TFree ((the o AList.lookup (op =) names) ixn, xs));
berghofe@11519
   529
berghofe@11519
   530
fun thawT names =
haftmann@17325
   531
  map_type_tfree (fn (s, xs) => case AList.lookup (op =) names s of
skalberg@15531
   532
      NONE => TFree (s, xs)
skalberg@15531
   533
    | SOME ixn => TVar (ixn, xs));
berghofe@11519
   534
berghofe@11519
   535
fun freeze names names' (t $ u) =
berghofe@11519
   536
      freeze names names' t $ freeze names names' u
berghofe@11519
   537
  | freeze names names' (Abs (s, T, t)) =
berghofe@11519
   538
      Abs (s, frzT names' T, freeze names names' t)
berghofe@11519
   539
  | freeze names names' (Const (s, T)) = Const (s, frzT names' T)
berghofe@11519
   540
  | freeze names names' (Free (s, T)) = Free (s, frzT names' T)
berghofe@11519
   541
  | freeze names names' (Var (ixn, T)) =
haftmann@17325
   542
      Free ((the o AList.lookup (op =) names) ixn, frzT names' T)
berghofe@11519
   543
  | freeze names names' t = t;
berghofe@11519
   544
berghofe@11519
   545
fun thaw names names' (t $ u) =
berghofe@11519
   546
      thaw names names' t $ thaw names names' u
berghofe@11519
   547
  | thaw names names' (Abs (s, T, t)) =
berghofe@11519
   548
      Abs (s, thawT names' T, thaw names names' t)
berghofe@11519
   549
  | thaw names names' (Const (s, T)) = Const (s, thawT names' T)
berghofe@11519
   550
  | thaw names names' (Free (s, T)) = 
berghofe@11519
   551
      let val T' = thawT names' T
haftmann@17325
   552
      in case AList.lookup (op =) names s of
skalberg@15531
   553
          NONE => Free (s, T')
skalberg@15531
   554
        | SOME ixn => Var (ixn, T')
berghofe@11519
   555
      end
berghofe@11519
   556
  | thaw names names' (Var (ixn, T)) = Var (ixn, thawT names' T)
berghofe@11519
   557
  | thaw names names' t = t;
berghofe@11519
   558
berghofe@11519
   559
fun freeze_thaw_prf prf =
berghofe@11519
   560
  let
berghofe@11519
   561
    val (fs, Tfs, vs, Tvs) = fold_proof_terms
berghofe@11519
   562
      (fn (t, (fs, Tfs, vs, Tvs)) =>
berghofe@11519
   563
         (add_term_frees (t, fs), add_term_tfree_names (t, Tfs),
berghofe@11519
   564
          add_term_vars (t, vs), add_term_tvar_ixns (t, Tvs)))
berghofe@11519
   565
      (fn (T, (fs, Tfs, vs, Tvs)) =>
berghofe@11519
   566
         (fs, add_typ_tfree_names (T, Tfs),
berghofe@11519
   567
          vs, add_typ_ixns (Tvs, T)))
berghofe@11519
   568
            (([], [], [], []), prf);
berghofe@11519
   569
    val fs' = map (fst o dest_Free) fs;
berghofe@11519
   570
    val vs' = map (fst o dest_Var) vs;
berghofe@11519
   571
    val names = vs' ~~ variantlist (map fst vs', fs');
berghofe@11519
   572
    val names' = Tvs ~~ variantlist (map fst Tvs, Tfs);
berghofe@11519
   573
    val rnames = map swap names;
berghofe@11519
   574
    val rnames' = map swap names';
berghofe@11519
   575
  in
berghofe@11519
   576
    (map_proof_terms (freeze names names') (frzT names') prf,
berghofe@11519
   577
     map_proof_terms (thaw rnames rnames') (thawT rnames'))
berghofe@11519
   578
  end;
berghofe@11519
   579
berghofe@11519
   580
berghofe@11519
   581
(***** implication introduction *****)
berghofe@11519
   582
berghofe@11519
   583
fun implies_intr_proof h prf =
berghofe@11519
   584
  let
berghofe@11715
   585
    fun abshyp i (Hyp t) = if h aconv t then PBound i else raise SAME
berghofe@11519
   586
      | abshyp i (Abst (s, T, prf)) = Abst (s, T, abshyp i prf)
berghofe@11519
   587
      | abshyp i (AbsP (s, t, prf)) = AbsP (s, t, abshyp (i+1) prf)
berghofe@11615
   588
      | abshyp i (prf % t) = abshyp i prf % t
berghofe@11715
   589
      | abshyp i (prf1 %% prf2) = (abshyp i prf1 %% abshyph i prf2
berghofe@11715
   590
          handle SAME => prf1 %% abshyp i prf2)
berghofe@11715
   591
      | abshyp _ _ = raise SAME
berghofe@11715
   592
    and abshyph i prf = (abshyp i prf handle SAME => prf)
berghofe@11519
   593
  in
skalberg@15531
   594
    AbsP ("H", NONE (*h*), abshyph 0 prf)
berghofe@11519
   595
  end;
berghofe@11519
   596
berghofe@11519
   597
berghofe@11519
   598
(***** forall introduction *****)
berghofe@11519
   599
skalberg@15531
   600
fun forall_intr_proof x a prf = Abst (a, NONE, prf_abstract_over x prf);
berghofe@11519
   601
berghofe@11519
   602
berghofe@11519
   603
(***** varify *****)
berghofe@11519
   604
berghofe@11519
   605
fun varify_proof t fixed prf =
berghofe@11519
   606
  let
wenzelm@19304
   607
    val fs = Term.fold_types (Term.fold_atyps
wenzelm@19304
   608
      (fn TFree v => if member (op =) fixed v then I else insert (op =) v | _ => I)) t [];
berghofe@11519
   609
    val ixns = add_term_tvar_ixns (t, []);
berghofe@15797
   610
    val fmap = fs ~~ variantlist (map fst fs, map #1 ixns)
berghofe@11519
   611
    fun thaw (f as (a, S)) =
haftmann@17314
   612
      (case AList.lookup (op =) fmap f of
skalberg@15531
   613
        NONE => TFree f
skalberg@15531
   614
      | SOME b => TVar ((b, 0), S));
berghofe@11519
   615
  in map_proof_terms (map_term_types (map_type_tfree thaw)) (map_type_tfree thaw) prf
berghofe@11519
   616
  end;
berghofe@11519
   617
berghofe@11519
   618
berghofe@11519
   619
local
berghofe@11519
   620
berghofe@11519
   621
fun new_name (ix, (pairs,used)) =
berghofe@11519
   622
  let val v = variant used (string_of_indexname ix)
berghofe@11519
   623
  in  ((ix, v) :: pairs, v :: used)  end;
berghofe@11519
   624
haftmann@17325
   625
fun freeze_one alist (ix, sort) = (case AList.lookup (op =) alist ix of
skalberg@15531
   626
    NONE => TVar (ix, sort)
skalberg@15531
   627
  | SOME name => TFree (name, sort));
berghofe@11519
   628
berghofe@11519
   629
in
berghofe@11519
   630
berghofe@11519
   631
fun freezeT t prf =
berghofe@11519
   632
  let
berghofe@11519
   633
    val used = it_term_types add_typ_tfree_names (t, [])
berghofe@11519
   634
    and tvars = map #1 (it_term_types add_typ_tvars (t, []));
skalberg@15574
   635
    val (alist, _) = foldr new_name ([], used) tvars;
berghofe@11519
   636
  in
berghofe@11519
   637
    (case alist of
berghofe@11519
   638
      [] => prf (*nothing to do!*)
berghofe@11519
   639
    | _ =>
berghofe@11519
   640
      let val frzT = map_type_tvar (freeze_one alist)
berghofe@11519
   641
      in map_proof_terms (map_term_types frzT) frzT prf end)
berghofe@11519
   642
  end;
berghofe@11519
   643
berghofe@11519
   644
end;
berghofe@11519
   645
berghofe@11519
   646
berghofe@11519
   647
(***** rotate assumptions *****)
berghofe@11519
   648
berghofe@11519
   649
fun rotate_proof Bs Bi m prf =
berghofe@11519
   650
  let
berghofe@11519
   651
    val params = Term.strip_all_vars Bi;
berghofe@11519
   652
    val asms = Logic.strip_imp_prems (Term.strip_all_body Bi);
berghofe@11519
   653
    val i = length asms;
berghofe@11519
   654
    val j = length Bs;
berghofe@11519
   655
  in
berghofe@11519
   656
    mk_AbsP (j+1, proof_combP (prf, map PBound
skalberg@15574
   657
      (j downto 1) @ [mk_Abst (mk_AbsP (i,
berghofe@11519
   658
        proof_combP (proof_combt (PBound i, map Bound ((length params - 1) downto 0)),
skalberg@15574
   659
          map PBound (((i-m-1) downto 0) @ ((i-1) downto (i-m)))))) params]))
berghofe@11519
   660
  end;
berghofe@11519
   661
berghofe@11519
   662
berghofe@11519
   663
(***** permute premises *****)
berghofe@11519
   664
berghofe@11519
   665
fun permute_prems_prf prems j k prf =
berghofe@11519
   666
  let val n = length prems
berghofe@11519
   667
  in mk_AbsP (n, proof_combP (prf,
berghofe@11519
   668
    map PBound ((n-1 downto n-j) @ (k-1 downto 0) @ (n-j-1 downto k))))
berghofe@11519
   669
  end;
berghofe@11519
   670
berghofe@11519
   671
wenzelm@19908
   672
(***** generalization *****)
wenzelm@19908
   673
wenzelm@20000
   674
fun generalize (tfrees, frees) idx =
wenzelm@20000
   675
  map_proof_terms_option
wenzelm@20000
   676
    (Term.generalize_option (tfrees, frees) idx)
wenzelm@20000
   677
    (Term.generalizeT_option tfrees idx);
wenzelm@19908
   678
wenzelm@19908
   679
berghofe@11519
   680
(***** instantiation *****)
berghofe@11519
   681
wenzelm@20000
   682
fun instantiate (instT, inst) =
wenzelm@20000
   683
  map_proof_terms_option
wenzelm@20000
   684
    (Term.instantiate_option (instT, map (apsnd remove_types) inst))
wenzelm@20000
   685
    (Term.instantiateT_option instT);
berghofe@11519
   686
berghofe@11519
   687
berghofe@11519
   688
(***** lifting *****)
berghofe@11519
   689
berghofe@11519
   690
fun lift_proof Bi inc prop prf =
berghofe@11519
   691
  let
berghofe@11519
   692
    fun lift'' Us Ts t = strip_abs Ts (Logic.incr_indexes (Us, inc) (mk_abs Ts t));
berghofe@11519
   693
berghofe@11715
   694
    fun lift' Us Ts (Abst (s, T, prf)) =
wenzelm@16880
   695
          (Abst (s, apsome' (same (Logic.incr_tvar inc)) T, lifth' Us (dummyT::Ts) prf)
berghofe@11715
   696
           handle SAME => Abst (s, T, lift' Us (dummyT::Ts) prf))
berghofe@11715
   697
      | lift' Us Ts (AbsP (s, t, prf)) =
berghofe@11715
   698
          (AbsP (s, apsome' (same (lift'' Us Ts)) t, lifth' Us Ts prf)
berghofe@11715
   699
           handle SAME => AbsP (s, t, lift' Us Ts prf))
skalberg@15570
   700
      | lift' Us Ts (prf % t) = (lift' Us Ts prf % Option.map (lift'' Us Ts) t
berghofe@11715
   701
          handle SAME => prf % apsome' (same (lift'' Us Ts)) t)
berghofe@11715
   702
      | lift' Us Ts (prf1 %% prf2) = (lift' Us Ts prf1 %% lifth' Us Ts prf2
berghofe@11715
   703
          handle SAME => prf1 %% lift' Us Ts prf2)
berghofe@11715
   704
      | lift' _ _ (PThm (s, prf, prop, Ts)) =
wenzelm@16880
   705
          PThm (s, prf, prop, apsome' (same (map (Logic.incr_tvar inc))) Ts)
berghofe@11715
   706
      | lift' _ _ (PAxm (s, prop, Ts)) =
wenzelm@16880
   707
          PAxm (s, prop, apsome' (same (map (Logic.incr_tvar inc))) Ts)
berghofe@11715
   708
      | lift' _ _ _ = raise SAME
berghofe@11715
   709
    and lifth' Us Ts prf = (lift' Us Ts prf handle SAME => prf);
berghofe@11519
   710
wenzelm@18030
   711
    val ps = map (Logic.lift_all inc Bi) (Logic.strip_imp_prems prop);
berghofe@11519
   712
    val k = length ps;
berghofe@11519
   713
berghofe@11519
   714
    fun mk_app (b, (i, j, prf)) = 
berghofe@11615
   715
          if b then (i-1, j, prf %% PBound i) else (i, j-1, prf %> Bound j);
berghofe@11519
   716
berghofe@11519
   717
    fun lift Us bs i j (Const ("==>", _) $ A $ B) =
skalberg@15531
   718
	    AbsP ("H", NONE (*A*), lift Us (true::bs) (i+1) j B)
berghofe@11519
   719
      | lift Us bs i j (Const ("all", _) $ Abs (a, T, t)) = 
skalberg@15531
   720
	    Abst (a, NONE (*T*), lift (T::Us) (false::bs) i (j+1) t)
berghofe@11715
   721
      | lift Us bs i j _ = proof_combP (lifth' (rev Us) [] prf,
skalberg@15574
   722
            map (fn k => (#3 (foldr mk_app (i-1, j-1, PBound k) bs)))
berghofe@11519
   723
              (i + k - 1 downto i));
berghofe@11519
   724
  in
berghofe@11519
   725
    mk_AbsP (k, lift [] [] 0 0 Bi)
berghofe@11519
   726
  end;
berghofe@11519
   727
berghofe@11519
   728
berghofe@11519
   729
(***** proof by assumption *****)
berghofe@11519
   730
skalberg@15531
   731
fun mk_asm_prf (Const ("==>", _) $ A $ B) i = AbsP ("H", NONE (*A*), mk_asm_prf B (i+1))
skalberg@15531
   732
  | mk_asm_prf (Const ("all", _) $ Abs (a, T, t)) i = Abst (a, NONE (*T*), mk_asm_prf t i)
berghofe@11519
   733
  | mk_asm_prf _ i = PBound i;
berghofe@11519
   734
berghofe@11519
   735
fun assumption_proof Bs Bi n prf =
berghofe@11519
   736
  mk_AbsP (length Bs, proof_combP (prf,
berghofe@11519
   737
    map PBound (length Bs - 1 downto 0) @ [mk_asm_prf Bi (~n)]));
berghofe@11519
   738
berghofe@11519
   739
berghofe@11519
   740
(***** Composition of object rule with proof state *****)
berghofe@11519
   741
berghofe@11519
   742
fun flatten_params_proof i j n (Const ("==>", _) $ A $ B, k) =
skalberg@15531
   743
      AbsP ("H", NONE (*A*), flatten_params_proof (i+1) j n (B, k))
berghofe@11519
   744
  | flatten_params_proof i j n (Const ("all", _) $ Abs (a, T, t), k) =
skalberg@15531
   745
      Abst (a, NONE (*T*), flatten_params_proof i (j+1) n (t, k))
berghofe@11519
   746
  | flatten_params_proof i j n (_, k) = proof_combP (proof_combt (PBound (k+i),
wenzelm@19304
   747
      map Bound (j-1 downto 0)), map PBound (remove (op =) (i-n) (i-1 downto 0)));
berghofe@11519
   748
wenzelm@18485
   749
fun bicompose_proof flatten Bs oldAs newAs A n rprf sprf =
berghofe@11519
   750
  let
berghofe@11519
   751
    val la = length newAs;
berghofe@11519
   752
    val lb = length Bs;
berghofe@11519
   753
  in
berghofe@11519
   754
    mk_AbsP (lb+la, proof_combP (sprf,
berghofe@11615
   755
      map PBound (lb + la - 1 downto la)) %%
wenzelm@18485
   756
        proof_combP (rprf, (if n>0 then [mk_asm_prf (the A) (~n)] else []) @
wenzelm@18485
   757
          map (if flatten then flatten_params_proof 0 0 n else PBound o snd)
wenzelm@18485
   758
            (oldAs ~~ (la - 1 downto 0))))
berghofe@11519
   759
  end;
berghofe@11519
   760
berghofe@11519
   761
berghofe@11519
   762
(***** axioms for equality *****)
berghofe@11519
   763
wenzelm@14854
   764
val aT = TFree ("'a", []);
wenzelm@14854
   765
val bT = TFree ("'b", []);
berghofe@11519
   766
val x = Free ("x", aT);
berghofe@11519
   767
val y = Free ("y", aT);
berghofe@11519
   768
val z = Free ("z", aT);
berghofe@11519
   769
val A = Free ("A", propT);
berghofe@11519
   770
val B = Free ("B", propT);
berghofe@11519
   771
val f = Free ("f", aT --> bT);
berghofe@11519
   772
val g = Free ("g", aT --> bT);
berghofe@11519
   773
berghofe@11519
   774
local open Logic in
berghofe@11519
   775
berghofe@11519
   776
val equality_axms =
berghofe@11519
   777
  [("reflexive", mk_equals (x, x)),
berghofe@11519
   778
   ("symmetric", mk_implies (mk_equals (x, y), mk_equals (y, x))),
berghofe@11519
   779
   ("transitive", list_implies ([mk_equals (x, y), mk_equals (y, z)], mk_equals (x, z))),
berghofe@11519
   780
   ("equal_intr", list_implies ([mk_implies (A, B), mk_implies (B, A)], mk_equals (A, B))),
berghofe@11519
   781
   ("equal_elim", list_implies ([mk_equals (A, B), A], B)),
berghofe@11519
   782
   ("abstract_rule", Logic.mk_implies
berghofe@11519
   783
      (all aT $ Abs ("x", aT, equals bT $ (f $ Bound 0) $ (g $ Bound 0)),
berghofe@11519
   784
       equals (aT --> bT) $
berghofe@11519
   785
         Abs ("x", aT, f $ Bound 0) $ Abs ("x", aT, g $ Bound 0))),
berghofe@11519
   786
   ("combination", Logic.list_implies
berghofe@11519
   787
      ([Logic.mk_equals (f, g), Logic.mk_equals (x, y)],
berghofe@11519
   788
       Logic.mk_equals (f $ x, g $ y)))];
berghofe@11519
   789
berghofe@11519
   790
val [reflexive_axm, symmetric_axm, transitive_axm, equal_intr_axm,
berghofe@11519
   791
  equal_elim_axm, abstract_rule_axm, combination_axm] =
skalberg@15531
   792
    map (fn (s, t) => PAxm ("ProtoPure." ^ s, varify t, NONE)) equality_axms;
berghofe@11519
   793
berghofe@11519
   794
end;
berghofe@11519
   795
skalberg@15531
   796
val reflexive = reflexive_axm % NONE;
berghofe@11519
   797
berghofe@11615
   798
fun symmetric (prf as PAxm ("ProtoPure.reflexive", _, _) % _) = prf
skalberg@15531
   799
  | symmetric prf = symmetric_axm % NONE % NONE %% prf;
berghofe@11519
   800
berghofe@11615
   801
fun transitive _ _ (PAxm ("ProtoPure.reflexive", _, _) % _) prf2 = prf2
berghofe@11615
   802
  | transitive _ _ prf1 (PAxm ("ProtoPure.reflexive", _, _) % _) = prf1
berghofe@11519
   803
  | transitive u (Type ("prop", [])) prf1 prf2 =
skalberg@15531
   804
      transitive_axm % NONE % SOME (remove_types u) % NONE %% prf1 %% prf2
berghofe@11519
   805
  | transitive u T prf1 prf2 =
skalberg@15531
   806
      transitive_axm % NONE % NONE % NONE %% prf1 %% prf2;
berghofe@11519
   807
berghofe@11519
   808
fun abstract_rule x a prf =
skalberg@15531
   809
  abstract_rule_axm % NONE % NONE %% forall_intr_proof x a prf;
berghofe@11519
   810
berghofe@11615
   811
fun check_comb (PAxm ("ProtoPure.combination", _, _) % f % g % _ % _ %% prf %% _) =
wenzelm@19502
   812
      is_some f orelse check_comb prf
berghofe@11615
   813
  | check_comb (PAxm ("ProtoPure.transitive", _, _) % _ % _ % _ %% prf1 %% prf2) =
berghofe@11519
   814
      check_comb prf1 andalso check_comb prf2
berghofe@11615
   815
  | check_comb (PAxm ("ProtoPure.symmetric", _, _) % _ % _ %% prf) = check_comb prf
berghofe@11519
   816
  | check_comb _ = false;
berghofe@11519
   817
berghofe@11519
   818
fun combination f g t u (Type (_, [T, U])) prf1 prf2 =
berghofe@11519
   819
  let
berghofe@11519
   820
    val f = Envir.beta_norm f;
berghofe@11519
   821
    val g = Envir.beta_norm g;
berghofe@11519
   822
    val prf =  if check_comb prf1 then
skalberg@15531
   823
        combination_axm % NONE % NONE
berghofe@11519
   824
      else (case prf1 of
berghofe@11615
   825
          PAxm ("ProtoPure.reflexive", _, _) % _ =>
skalberg@15531
   826
            combination_axm %> remove_types f % NONE
berghofe@11615
   827
        | _ => combination_axm %> remove_types f %> remove_types g)
berghofe@11519
   828
  in
berghofe@11519
   829
    (case T of
berghofe@11615
   830
       Type ("fun", _) => prf %
berghofe@11519
   831
         (case head_of f of
skalberg@15531
   832
            Abs _ => SOME (remove_types t)
skalberg@15531
   833
          | Var _ => SOME (remove_types t)
skalberg@15531
   834
          | _ => NONE) %
berghofe@11519
   835
         (case head_of g of
skalberg@15531
   836
            Abs _ => SOME (remove_types u)
skalberg@15531
   837
          | Var _ => SOME (remove_types u)
skalberg@15531
   838
          | _ => NONE) %% prf1 %% prf2
skalberg@15531
   839
     | _ => prf % NONE % NONE %% prf1 %% prf2)
berghofe@11519
   840
  end;
berghofe@11519
   841
berghofe@11519
   842
fun equal_intr A B prf1 prf2 =
berghofe@11615
   843
  equal_intr_axm %> remove_types A %> remove_types B %% prf1 %% prf2;
berghofe@11519
   844
berghofe@11519
   845
fun equal_elim A B prf1 prf2 =
berghofe@11615
   846
  equal_elim_axm %> remove_types A %> remove_types B %% prf1 %% prf2;
berghofe@11519
   847
berghofe@11519
   848
berghofe@11519
   849
(***** axioms and theorems *****)
berghofe@11519
   850
wenzelm@17492
   851
fun vars_of t = rev (fold_aterms (fn v as Var _ => insert (op =) v | _ => I) t []);
berghofe@11519
   852
berghofe@11519
   853
fun test_args _ [] = true
berghofe@11519
   854
  | test_args is (Bound i :: ts) =
wenzelm@17492
   855
      not (member (op =) is i) andalso test_args (i :: is) ts
berghofe@11519
   856
  | test_args _ _ = false;
berghofe@11519
   857
berghofe@11519
   858
fun is_fun (Type ("fun", _)) = true
berghofe@11519
   859
  | is_fun (TVar _) = true
berghofe@11519
   860
  | is_fun _ = false;
berghofe@11519
   861
berghofe@11519
   862
fun add_funvars Ts (vs, t) =
berghofe@11519
   863
  if is_fun (fastype_of1 (Ts, t)) then
wenzelm@19482
   864
    vs union map_filter (fn Var (ixn, T) =>
skalberg@15531
   865
      if is_fun T then SOME ixn else NONE | _ => NONE) (vars_of t)
berghofe@11519
   866
  else vs;
berghofe@11519
   867
berghofe@11519
   868
fun add_npvars q p Ts (vs, Const ("==>", _) $ t $ u) =
berghofe@11519
   869
      add_npvars q p Ts (add_npvars q (not p) Ts (vs, t), u)
berghofe@11519
   870
  | add_npvars q p Ts (vs, Const ("all", Type (_, [Type (_, [T, _]), _])) $ t) =
berghofe@11519
   871
      add_npvars q p Ts (vs, if p andalso q then betapply (t, Var (("",0), T)) else t)
berghofe@12041
   872
  | add_npvars q p Ts (vs, Abs (_, T, t)) = add_npvars q p (T::Ts) (vs, t)
berghofe@12041
   873
  | add_npvars _ _ Ts (vs, t) = add_npvars' Ts (vs, t)
berghofe@12041
   874
and add_npvars' Ts (vs, t) = (case strip_comb t of
berghofe@11519
   875
    (Var (ixn, _), ts) => if test_args [] ts then vs
haftmann@17314
   876
      else Library.foldl (add_npvars' Ts)
haftmann@17314
   877
        (AList.update (op =) (ixn,
haftmann@17314
   878
          Library.foldl (add_funvars Ts) ((these ooo AList.lookup) (op =) vs ixn, ts)) vs, ts)
skalberg@15570
   879
  | (Abs (_, T, u), ts) => Library.foldl (add_npvars' (T::Ts)) (vs, u :: ts)
skalberg@15570
   880
  | (_, ts) => Library.foldl (add_npvars' Ts) (vs, ts));
berghofe@11519
   881
berghofe@11519
   882
fun prop_vars (Const ("==>", _) $ P $ Q) = prop_vars P union prop_vars Q
berghofe@11519
   883
  | prop_vars (Const ("all", _) $ Abs (_, _, t)) = prop_vars t
berghofe@11519
   884
  | prop_vars t = (case strip_comb t of
berghofe@11519
   885
      (Var (ixn, _), _) => [ixn] | _ => []);
berghofe@11519
   886
berghofe@11519
   887
fun is_proj t =
berghofe@11519
   888
  let
berghofe@11519
   889
    fun is_p i t = (case strip_comb t of
berghofe@11519
   890
        (Bound j, []) => false
berghofe@11519
   891
      | (Bound j, ts) => j >= i orelse exists (is_p i) ts
berghofe@11519
   892
      | (Abs (_, _, u), _) => is_p (i+1) u
berghofe@11519
   893
      | (_, ts) => exists (is_p i) ts)
berghofe@11519
   894
  in (case strip_abs_body t of
berghofe@11519
   895
        Bound _ => true
berghofe@11519
   896
      | t' => is_p 0 t')
berghofe@11519
   897
  end;
berghofe@11519
   898
berghofe@11519
   899
fun needed_vars prop = 
skalberg@15570
   900
  Library.foldl op union ([], map op ins (add_npvars true true [] ([], prop))) union
berghofe@11519
   901
  prop_vars prop;
berghofe@11519
   902
berghofe@11519
   903
fun gen_axm_proof c name prop =
berghofe@11519
   904
  let
berghofe@11519
   905
    val nvs = needed_vars prop;
berghofe@11519
   906
    val args = map (fn (v as Var (ixn, _)) =>
wenzelm@17492
   907
        if member (op =) nvs ixn then SOME v else NONE) (vars_of prop) @
wenzelm@16983
   908
      map SOME (sort Term.term_ord (term_frees prop));
berghofe@11519
   909
  in
skalberg@15531
   910
    proof_combt' (c (name, prop, NONE), args)
berghofe@11519
   911
  end;
berghofe@11519
   912
berghofe@11519
   913
val axm_proof = gen_axm_proof PAxm;
berghofe@17017
   914
berghofe@17017
   915
val dummy = Const (Term.dummy_patternN, dummyT);
berghofe@17017
   916
berghofe@17017
   917
fun oracle_proof name prop =
berghofe@17017
   918
  if !proofs = 0 then Oracle (name, dummy, NONE)
berghofe@17017
   919
  else gen_axm_proof Oracle name prop;
berghofe@11519
   920
wenzelm@17492
   921
fun shrink_proof thy =
wenzelm@17492
   922
  let
wenzelm@17492
   923
    val compress_typ = Compress.typ thy;
wenzelm@17492
   924
    val compress_term = Compress.term thy;
wenzelm@17492
   925
  
wenzelm@17492
   926
    fun shrink ls lev (prf as Abst (a, T, body)) =
wenzelm@17492
   927
          let val (b, is, ch, body') = shrink ls (lev+1) body
wenzelm@17492
   928
          in (b, is, ch, if ch then Abst (a, Option.map compress_typ T, body') else prf) end
wenzelm@17492
   929
      | shrink ls lev (prf as AbsP (a, t, body)) =
wenzelm@17492
   930
          let val (b, is, ch, body') = shrink (lev::ls) lev body
wenzelm@19482
   931
          in (b orelse member (op =) is 0, map_filter (fn 0 => NONE | i => SOME (i-1)) is,
wenzelm@17492
   932
            ch, if ch then AbsP (a, Option.map compress_term t, body') else prf)
wenzelm@17492
   933
          end
wenzelm@17492
   934
      | shrink ls lev prf =
wenzelm@17492
   935
          let val (is, ch, _, prf') = shrink' ls lev [] [] prf
wenzelm@17492
   936
          in (false, is, ch, prf') end
wenzelm@17492
   937
    and shrink' ls lev ts prfs (prf as prf1 %% prf2) =
wenzelm@17492
   938
          let
wenzelm@17492
   939
            val p as (_, is', ch', prf') = shrink ls lev prf2;
wenzelm@17492
   940
            val (is, ch, ts', prf'') = shrink' ls lev ts (p::prfs) prf1
wenzelm@17492
   941
          in (is union is', ch orelse ch', ts',
wenzelm@17492
   942
              if ch orelse ch' then prf'' %% prf' else prf)
wenzelm@17492
   943
          end
wenzelm@17492
   944
      | shrink' ls lev ts prfs (prf as prf1 % t) =
wenzelm@17492
   945
          let val (is, ch, (ch', t')::ts', prf') = shrink' ls lev (t::ts) prfs prf1
wenzelm@17492
   946
          in (is, ch orelse ch', ts',
wenzelm@17492
   947
              if ch orelse ch' then prf' % Option.map compress_term t' else prf) end
wenzelm@17492
   948
      | shrink' ls lev ts prfs (prf as PBound i) =
wenzelm@17492
   949
          (if exists (fn SOME (Bound j) => lev-j <= List.nth (ls, i) | _ => true) ts
haftmann@18928
   950
             orelse has_duplicates (op =)
haftmann@18928
   951
               (Library.foldl (fn (js, SOME (Bound j)) => j :: js | (js, _) => js) ([], ts))
wenzelm@17492
   952
             orelse exists #1 prfs then [i] else [], false, map (pair false) ts, prf)
wenzelm@17492
   953
      | shrink' ls lev ts prfs (Hyp t) = ([], false, map (pair false) ts, Hyp (compress_term t))
wenzelm@17492
   954
      | shrink' ls lev ts prfs (prf as MinProof _) =
wenzelm@17492
   955
          ([], false, map (pair false) ts, prf)
wenzelm@17492
   956
      | shrink' ls lev ts prfs prf =
wenzelm@17492
   957
          let
wenzelm@17492
   958
            val prop = (case prf of PThm (_, _, prop, _) => prop | PAxm (_, prop, _) => prop
wenzelm@17492
   959
              | Oracle (_, prop, _) => prop | _ => error "shrink: proof not in normal form");
wenzelm@17492
   960
            val vs = vars_of prop;
wenzelm@19012
   961
            val (ts', ts'') = chop (length vs) ts;
wenzelm@17492
   962
            val insts = Library.take (length ts', map (fst o dest_Var) vs) ~~ ts';
wenzelm@17492
   963
            val nvs = Library.foldl (fn (ixns', (ixn, ixns)) =>
wenzelm@17492
   964
              insert (op =) ixn (case AList.lookup (op =) insts ixn of
wenzelm@17492
   965
                  SOME (SOME t) => if is_proj t then ixns union ixns' else ixns'
wenzelm@17492
   966
                | _ => ixns union ixns'))
wenzelm@17492
   967
                  (needed prop ts'' prfs, add_npvars false true [] ([], prop));
wenzelm@17492
   968
            val insts' = map
wenzelm@17492
   969
              (fn (ixn, x as SOME _) => if member (op =) nvs ixn then (false, x) else (true, NONE)
wenzelm@17492
   970
                | (_, x) => (false, x)) insts
wenzelm@17492
   971
          in ([], false, insts' @ map (pair false) ts'', prf) end
wenzelm@17492
   972
    and needed (Const ("==>", _) $ t $ u) ts ((b, _, _, _)::prfs) =
wenzelm@17492
   973
          (if b then map (fst o dest_Var) (vars_of t) else []) union needed u ts prfs
wenzelm@17492
   974
      | needed (Var (ixn, _)) (_::_) _ = [ixn]
wenzelm@17492
   975
      | needed _ _ _ = [];
wenzelm@17492
   976
  in shrink end;
berghofe@11519
   977
berghofe@11519
   978
berghofe@11519
   979
(**** Simple first order matching functions for terms and proofs ****)
berghofe@11519
   980
berghofe@11519
   981
exception PMatch;
berghofe@11519
   982
berghofe@11519
   983
(** see pattern.ML **)
berghofe@11519
   984
skalberg@15570
   985
fun flt (i: int) = List.filter (fn n => n < i);
berghofe@12279
   986
berghofe@12279
   987
fun fomatch Ts tymatch j =
berghofe@11519
   988
  let
berghofe@11519
   989
    fun mtch (instsp as (tyinsts, insts)) = fn
berghofe@11519
   990
        (Var (ixn, T), t)  =>
berghofe@12279
   991
          if j>0 andalso not (null (flt j (loose_bnos t)))
berghofe@12279
   992
          then raise PMatch
berghofe@12279
   993
          else (tymatch (tyinsts, fn () => (T, fastype_of1 (Ts, t))),
berghofe@12279
   994
            (ixn, t) :: insts)
berghofe@11519
   995
      | (Free (a, T), Free (b, U)) =>
berghofe@12279
   996
	  if a=b then (tymatch (tyinsts, K (T, U)), insts) else raise PMatch
berghofe@11519
   997
      | (Const (a, T), Const (b, U))  =>
berghofe@12279
   998
	  if a=b then (tymatch (tyinsts, K (T, U)), insts) else raise PMatch
berghofe@11519
   999
      | (f $ t, g $ u) => mtch (mtch instsp (f, g)) (t, u)
berghofe@12279
  1000
      | (Bound i, Bound j) => if i=j then instsp else raise PMatch
berghofe@11519
  1001
      | _ => raise PMatch
berghofe@11519
  1002
  in mtch end;
berghofe@11519
  1003
berghofe@12279
  1004
fun match_proof Ts tymatch =
berghofe@11519
  1005
  let
skalberg@15531
  1006
    fun optmatch _ inst (NONE, _) = inst
skalberg@15531
  1007
      | optmatch _ _ (SOME _, NONE) = raise PMatch
skalberg@15531
  1008
      | optmatch mtch inst (SOME x, SOME y) = mtch inst (x, y)
berghofe@12279
  1009
berghofe@12279
  1010
    fun matcht Ts j (pinst, tinst) (t, u) =
berghofe@12279
  1011
      (pinst, fomatch Ts tymatch j tinst (t, Envir.beta_norm u));
berghofe@12279
  1012
    fun matchT (pinst, (tyinsts, insts)) p =
berghofe@12279
  1013
      (pinst, (tymatch (tyinsts, K p), insts));
skalberg@15570
  1014
    fun matchTs inst (Ts, Us) = Library.foldl (uncurry matchT) (inst, Ts ~~ Us);
berghofe@12279
  1015
berghofe@12279
  1016
    fun mtch Ts i j (pinst, tinst) (Hyp (Var (ixn, _)), prf) =
berghofe@12279
  1017
          if i = 0 andalso j = 0 then ((ixn, prf) :: pinst, tinst)
berghofe@12279
  1018
          else (case apfst (flt i) (apsnd (flt j)
berghofe@12279
  1019
                  (prf_add_loose_bnos 0 0 prf ([], []))) of
berghofe@12279
  1020
              ([], []) => ((ixn, incr_pboundvars (~i) (~j) prf) :: pinst, tinst)
berghofe@12279
  1021
            | ([], _) => if j = 0 then
berghofe@12279
  1022
                   ((ixn, incr_pboundvars (~i) (~j) prf) :: pinst, tinst)
berghofe@12279
  1023
                 else raise PMatch
berghofe@12279
  1024
            | _ => raise PMatch)
berghofe@12279
  1025
      | mtch Ts i j inst (prf1 % opt1, prf2 % opt2) =
berghofe@12279
  1026
          optmatch (matcht Ts j) (mtch Ts i j inst (prf1, prf2)) (opt1, opt2)
berghofe@12279
  1027
      | mtch Ts i j inst (prf1 %% prf2, prf1' %% prf2') =
berghofe@12279
  1028
          mtch Ts i j (mtch Ts i j inst (prf1, prf1')) (prf2, prf2')
berghofe@12279
  1029
      | mtch Ts i j inst (Abst (_, opT, prf1), Abst (_, opU, prf2)) =
wenzelm@18485
  1030
          mtch (the_default dummyT opU :: Ts) i (j+1)
berghofe@12279
  1031
            (optmatch matchT inst (opT, opU)) (prf1, prf2)
berghofe@12279
  1032
      | mtch Ts i j inst (prf1, Abst (_, opU, prf2)) =
wenzelm@18485
  1033
          mtch (the_default dummyT opU :: Ts) i (j+1) inst
berghofe@12279
  1034
            (incr_pboundvars 0 1 prf1 %> Bound 0, prf2)
berghofe@12279
  1035
      | mtch Ts i j inst (AbsP (_, opt, prf1), AbsP (_, opu, prf2)) =
berghofe@12279
  1036
          mtch Ts (i+1) j (optmatch (matcht Ts j) inst (opt, opu)) (prf1, prf2)
berghofe@12279
  1037
      | mtch Ts i j inst (prf1, AbsP (_, _, prf2)) =
berghofe@12279
  1038
          mtch Ts (i+1) j inst (incr_pboundvars 1 0 prf1 %% PBound 0, prf2)
berghofe@12279
  1039
      | mtch Ts i j inst (PThm ((name1, _), _, prop1, opTs),
berghofe@12279
  1040
            PThm ((name2, _), _, prop2, opUs)) =
berghofe@11519
  1041
          if name1=name2 andalso prop1=prop2 then
berghofe@12279
  1042
            optmatch matchTs inst (opTs, opUs)
berghofe@11519
  1043
          else raise PMatch
berghofe@12279
  1044
      | mtch Ts i j inst (PAxm (s1, _, opTs), PAxm (s2, _, opUs)) =
berghofe@12279
  1045
          if s1=s2 then optmatch matchTs inst (opTs, opUs)
berghofe@11519
  1046
          else raise PMatch
berghofe@12279
  1047
      | mtch _ _ _ inst (PBound i, PBound j) = if i = j then inst else raise PMatch
berghofe@12279
  1048
      | mtch _ _ _ _ _ = raise PMatch
berghofe@12279
  1049
  in mtch Ts 0 0 end;
berghofe@11519
  1050
berghofe@11519
  1051
fun prf_subst (pinst, (tyinsts, insts)) =
berghofe@11519
  1052
  let
berghofe@15797
  1053
    val substT = Envir.typ_subst_TVars tyinsts;
berghofe@11519
  1054
haftmann@17325
  1055
    fun subst' lev (t as Var (ixn, _)) = (case AList.lookup (op =) insts ixn of
skalberg@15531
  1056
          NONE => t
skalberg@15531
  1057
        | SOME u => incr_boundvars lev u)
berghofe@11519
  1058
      | subst' lev (Const (s, T)) = Const (s, substT T)
berghofe@11519
  1059
      | subst' lev (Free (s, T)) = Free (s, substT T)
berghofe@11519
  1060
      | subst' lev (Abs (a, T, body)) = Abs (a, substT T, subst' (lev+1) body)
berghofe@11519
  1061
      | subst' lev (f $ t) = subst' lev f $ subst' lev t
berghofe@11519
  1062
      | subst' _ t = t;
berghofe@11519
  1063
berghofe@11519
  1064
    fun subst plev tlev (AbsP (a, t, body)) =
skalberg@15570
  1065
          AbsP (a, Option.map (subst' tlev) t, subst (plev+1) tlev body)
berghofe@11519
  1066
      | subst plev tlev (Abst (a, T, body)) =
skalberg@15570
  1067
          Abst (a, Option.map substT T, subst plev (tlev+1) body)
berghofe@11615
  1068
      | subst plev tlev (prf %% prf') = subst plev tlev prf %% subst plev tlev prf'
skalberg@15570
  1069
      | subst plev tlev (prf % t) = subst plev tlev prf % Option.map (subst' tlev) t
haftmann@17325
  1070
      | subst plev tlev (prf as Hyp (Var (ixn, _))) = (case AList.lookup (op =) pinst ixn of
skalberg@15531
  1071
          NONE => prf
skalberg@15531
  1072
        | SOME prf' => incr_pboundvars plev tlev prf')
berghofe@11519
  1073
      | subst _ _ (PThm (id, prf, prop, Ts)) =
skalberg@15570
  1074
          PThm (id, prf, prop, Option.map (map substT) Ts)
berghofe@11519
  1075
      | subst _ _ (PAxm (id, prop, Ts)) =
skalberg@15570
  1076
          PAxm (id, prop, Option.map (map substT) Ts)
berghofe@11519
  1077
      | subst _ _ t = t
berghofe@11519
  1078
  in subst 0 0 end;
berghofe@11519
  1079
berghofe@12871
  1080
(*A fast unification filter: true unless the two terms cannot be unified. 
berghofe@12871
  1081
  Terms must be NORMAL.  Treats all Vars as distinct. *)
berghofe@12871
  1082
fun could_unify prf1 prf2 =
berghofe@12871
  1083
  let
berghofe@12871
  1084
    fun matchrands (prf1 %% prf2) (prf1' %% prf2') =
berghofe@12871
  1085
          could_unify prf2 prf2' andalso matchrands prf1 prf1'
skalberg@15531
  1086
      | matchrands (prf % SOME t) (prf' % SOME t') =
berghofe@12871
  1087
          Term.could_unify (t, t') andalso matchrands prf prf'
berghofe@12871
  1088
      | matchrands (prf % _) (prf' % _) = matchrands prf prf'
berghofe@12871
  1089
      | matchrands _ _ = true
berghofe@12871
  1090
berghofe@12871
  1091
    fun head_of (prf %% _) = head_of prf
berghofe@12871
  1092
      | head_of (prf % _) = head_of prf
berghofe@12871
  1093
      | head_of prf = prf
berghofe@12871
  1094
berghofe@12871
  1095
  in case (head_of prf1, head_of prf2) of
berghofe@12871
  1096
        (_, Hyp (Var _)) => true
berghofe@12871
  1097
      | (Hyp (Var _), _) => true
berghofe@12871
  1098
      | (PThm ((a, _), _, propa, _), PThm ((b, _), _, propb, _)) =>
berghofe@12871
  1099
          a = b andalso propa = propb andalso matchrands prf1 prf2
berghofe@12871
  1100
      | (PAxm (a, _, _), PAxm (b, _, _)) => a = b andalso matchrands prf1 prf2
berghofe@12871
  1101
      | (PBound i, PBound j) =>  i = j andalso matchrands prf1 prf2
berghofe@12871
  1102
      | (AbsP _, _) =>  true   (*because of possible eta equality*)
berghofe@12871
  1103
      | (Abst _, _) =>  true
berghofe@12871
  1104
      | (_, AbsP _) =>  true
berghofe@12871
  1105
      | (_, Abst _) =>  true
berghofe@12871
  1106
      | _ => false
berghofe@12871
  1107
  end;
berghofe@12871
  1108
berghofe@11519
  1109
(**** rewriting on proof terms ****)
berghofe@11519
  1110
berghofe@13102
  1111
val skel0 = PBound 0;
berghofe@13102
  1112
berghofe@12279
  1113
fun rewrite_prf tymatch (rules, procs) prf =
berghofe@11519
  1114
  let
skalberg@15531
  1115
    fun rew _ (Abst (_, _, body) % SOME t) = SOME (prf_subst_bounds [t] body, skel0)
skalberg@15531
  1116
      | rew _ (AbsP (_, _, body) %% prf) = SOME (prf_subst_pbounds [prf] body, skel0)
berghofe@11519
  1117
      | rew Ts prf = (case get_first (fn (_, r) => r Ts prf) procs of
skalberg@15531
  1118
          SOME prf' => SOME (prf', skel0)
skalberg@15531
  1119
        | NONE => get_first (fn (prf1, prf2) => SOME (prf_subst
berghofe@13102
  1120
            (match_proof Ts tymatch ([], (Vartab.empty, [])) (prf1, prf)) prf2, prf2)
skalberg@15570
  1121
               handle PMatch => NONE) (List.filter (could_unify prf o fst) rules));
berghofe@11519
  1122
berghofe@11615
  1123
    fun rew0 Ts (prf as AbsP (_, _, prf' %% PBound 0)) =
berghofe@11519
  1124
          if prf_loose_Pbvar1 prf' 0 then rew Ts prf
berghofe@11519
  1125
          else
berghofe@11519
  1126
            let val prf'' = incr_pboundvars (~1) 0 prf'
wenzelm@19502
  1127
            in SOME (the_default (prf'', skel0) (rew Ts prf'')) end
skalberg@15531
  1128
      | rew0 Ts (prf as Abst (_, _, prf' % SOME (Bound 0))) =
berghofe@11519
  1129
          if prf_loose_bvar1 prf' 0 then rew Ts prf
berghofe@11519
  1130
          else
berghofe@11519
  1131
            let val prf'' = incr_pboundvars 0 (~1) prf'
wenzelm@19502
  1132
            in SOME (the_default (prf'', skel0) (rew Ts prf'')) end
berghofe@11519
  1133
      | rew0 Ts prf = rew Ts prf;
berghofe@11519
  1134
skalberg@15531
  1135
    fun rew1 _ (Hyp (Var _)) _ = NONE
berghofe@13102
  1136
      | rew1 Ts skel prf = (case rew2 Ts skel prf of
skalberg@15531
  1137
          SOME prf1 => (case rew0 Ts prf1 of
wenzelm@19502
  1138
              SOME (prf2, skel') => SOME (the_default prf2 (rew1 Ts skel' prf2))
skalberg@15531
  1139
            | NONE => SOME prf1)
skalberg@15531
  1140
        | NONE => (case rew0 Ts prf of
wenzelm@19502
  1141
              SOME (prf1, skel') => SOME (the_default prf1 (rew1 Ts skel' prf1))
skalberg@15531
  1142
            | NONE => NONE))
berghofe@11519
  1143
skalberg@15531
  1144
    and rew2 Ts skel (prf % SOME t) = (case prf of
berghofe@11519
  1145
            Abst (_, _, body) =>
berghofe@11519
  1146
              let val prf' = prf_subst_bounds [t] body
wenzelm@19502
  1147
              in SOME (the_default prf' (rew2 Ts skel0 prf')) end
berghofe@13102
  1148
          | _ => (case rew1 Ts (case skel of skel' % _ => skel' | _ => skel0) prf of
skalberg@15531
  1149
              SOME prf' => SOME (prf' % SOME t)
skalberg@15531
  1150
            | NONE => NONE))
skalberg@15570
  1151
      | rew2 Ts skel (prf % NONE) = Option.map (fn prf' => prf' % NONE)
berghofe@13102
  1152
          (rew1 Ts (case skel of skel' % _ => skel' | _ => skel0) prf)
berghofe@13102
  1153
      | rew2 Ts skel (prf1 %% prf2) = (case prf1 of
berghofe@11519
  1154
            AbsP (_, _, body) =>
berghofe@11519
  1155
              let val prf' = prf_subst_pbounds [prf2] body
wenzelm@19502
  1156
              in SOME (the_default prf' (rew2 Ts skel0 prf')) end
berghofe@13102
  1157
          | _ =>
berghofe@13102
  1158
            let val (skel1, skel2) = (case skel of
berghofe@13102
  1159
                skel1 %% skel2 => (skel1, skel2)
berghofe@13102
  1160
              | _ => (skel0, skel0))
berghofe@13102
  1161
            in case rew1 Ts skel1 prf1 of
skalberg@15531
  1162
                SOME prf1' => (case rew1 Ts skel2 prf2 of
skalberg@15531
  1163
                    SOME prf2' => SOME (prf1' %% prf2')
skalberg@15531
  1164
                  | NONE => SOME (prf1' %% prf2))
skalberg@15531
  1165
              | NONE => (case rew1 Ts skel2 prf2 of
skalberg@15531
  1166
                    SOME prf2' => SOME (prf1 %% prf2')
skalberg@15531
  1167
                  | NONE => NONE)
berghofe@13102
  1168
            end)
wenzelm@19502
  1169
      | rew2 Ts skel (Abst (s, T, prf)) = (case rew1 (the_default dummyT T :: Ts)
berghofe@13102
  1170
              (case skel of Abst (_, _, skel') => skel' | _ => skel0) prf of
skalberg@15531
  1171
            SOME prf' => SOME (Abst (s, T, prf'))
skalberg@15531
  1172
          | NONE => NONE)
berghofe@13102
  1173
      | rew2 Ts skel (AbsP (s, t, prf)) = (case rew1 Ts
berghofe@13102
  1174
              (case skel of AbsP (_, _, skel') => skel' | _ => skel0) prf of
skalberg@15531
  1175
            SOME prf' => SOME (AbsP (s, t, prf'))
skalberg@15531
  1176
          | NONE => NONE)
skalberg@15531
  1177
      | rew2 _ _ _ = NONE
berghofe@11519
  1178
wenzelm@19502
  1179
  in the_default prf (rew1 [] skel0 prf) end;
berghofe@11519
  1180
wenzelm@17203
  1181
fun rewrite_proof thy = rewrite_prf (fn (tyenv, f) =>
wenzelm@17203
  1182
  Sign.typ_match thy (f ()) tyenv handle Type.TYPE_MATCH => raise PMatch);
berghofe@11519
  1183
berghofe@11715
  1184
fun rewrite_proof_notypes rews = rewrite_prf fst rews;
berghofe@11615
  1185
wenzelm@16940
  1186
berghofe@11519
  1187
(**** theory data ****)
berghofe@11519
  1188
wenzelm@16458
  1189
structure ProofData = TheoryDataFun
wenzelm@16458
  1190
(struct
berghofe@11519
  1191
  val name = "Pure/proof";
berghofe@11519
  1192
  type T = ((proof * proof) list *
berghofe@12233
  1193
    (string * (typ list -> proof -> proof option)) list);
berghofe@11519
  1194
berghofe@12233
  1195
  val empty = ([], []);
berghofe@12233
  1196
  val copy = I;
wenzelm@16458
  1197
  val extend = I;
wenzelm@16458
  1198
  fun merge _ ((rules1, procs1), (rules2, procs2)) =
wenzelm@12293
  1199
    (merge_lists rules1 rules2, merge_alists procs1 procs2);
berghofe@11519
  1200
  fun print _ _ = ();
wenzelm@16458
  1201
end);
berghofe@11519
  1202
wenzelm@16536
  1203
val init_data = ProofData.init;
berghofe@11519
  1204
berghofe@12233
  1205
fun add_prf_rrules rs thy =
berghofe@11519
  1206
  let val r = ProofData.get thy
berghofe@12233
  1207
  in ProofData.put (rs @ fst r, snd r) thy end;
berghofe@11519
  1208
berghofe@12233
  1209
fun add_prf_rprocs ps thy =
berghofe@11519
  1210
  let val r = ProofData.get thy
berghofe@12233
  1211
  in ProofData.put (fst r, ps @ snd r) thy end;
berghofe@11519
  1212
wenzelm@16458
  1213
fun thm_proof thy (name, tags) hyps prop prf =
berghofe@11519
  1214
  let
wenzelm@12923
  1215
    val prop = Logic.list_implies (hyps, prop);
berghofe@11519
  1216
    val nvs = needed_vars prop;
berghofe@11519
  1217
    val args = map (fn (v as Var (ixn, _)) =>
wenzelm@17492
  1218
        if member (op =) nvs ixn then SOME v else NONE) (vars_of prop) @
wenzelm@16983
  1219
      map SOME (sort Term.term_ord (term_frees prop));
wenzelm@11543
  1220
    val opt_prf = if ! proofs = 2 then
wenzelm@17492
  1221
        #4 (shrink_proof thy [] 0 (rewrite_prf fst (ProofData.get thy)
skalberg@15574
  1222
          (foldr (uncurry implies_intr_proof) prf hyps)))
berghofe@17017
  1223
      else MinProof (mk_min_proof prf ([], [], []));
berghofe@12233
  1224
    val head = (case strip_combt (fst (strip_combP prf)) of
skalberg@15531
  1225
        (PThm ((old_name, _), prf', prop', NONE), args') =>
berghofe@11519
  1226
          if (old_name="" orelse old_name=name) andalso
berghofe@11519
  1227
             prop = prop' andalso args = args' then
skalberg@15531
  1228
            PThm ((name, tags), prf', prop, NONE)
berghofe@11519
  1229
          else
skalberg@15531
  1230
            PThm ((name, tags), opt_prf, prop, NONE)
skalberg@15531
  1231
      | _ => PThm ((name, tags), opt_prf, prop, NONE))
berghofe@11519
  1232
  in
wenzelm@12923
  1233
    proof_combP (proof_combt' (head, args), map Hyp hyps)
berghofe@11519
  1234
  end;
berghofe@11519
  1235
wenzelm@12923
  1236
fun get_name_tags hyps prop prf =
wenzelm@12923
  1237
  let val prop = Logic.list_implies (hyps, prop) in
wenzelm@12923
  1238
    (case strip_combt (fst (strip_combP prf)) of
berghofe@11519
  1239
      (PThm ((name, tags), _, prop', _), _) =>
berghofe@11519
  1240
        if prop=prop' then (name, tags) else ("", [])
berghofe@11519
  1241
    | (PAxm (name, prop', _), _) =>
berghofe@11519
  1242
        if prop=prop' then (name, []) else ("", [])
wenzelm@12923
  1243
    | _ => ("", []))
wenzelm@12923
  1244
  end;
berghofe@11519
  1245
berghofe@11519
  1246
end;
berghofe@11519
  1247
berghofe@11519
  1248
structure BasicProofterm : BASIC_PROOFTERM = Proofterm;
berghofe@11519
  1249
open BasicProofterm;