src/FOLP/simp.ML
author wenzelm
Sat Aug 29 12:01:25 2009 +0200 (2009-08-29)
changeset 32449 696d64ed85da
parent 32091 30e2ffbba718
child 32740 9dd0a2f83429
permissions -rw-r--r--
eliminated hard tabs;
wenzelm@29265
     1
(*  Title:      FOLP/simp.ML
clasohm@0
     2
    Author:     Tobias Nipkow
clasohm@0
     3
    Copyright   1993  University of Cambridge
clasohm@0
     4
clasohm@0
     5
FOLP version of...
clasohm@0
     6
clasohm@0
     7
Generic simplifier, suitable for most logics.  (from Provers)
clasohm@0
     8
clasohm@0
     9
This version allows instantiation of Vars in the subgoal, since the proof
clasohm@0
    10
term must change.
clasohm@0
    11
*)
clasohm@0
    12
clasohm@0
    13
signature SIMP_DATA =
clasohm@0
    14
sig
clasohm@0
    15
  val case_splits  : (thm * string) list
clasohm@0
    16
  val dest_red     : term -> term * term * term
clasohm@0
    17
  val mk_rew_rules : thm -> thm list
clasohm@0
    18
  val norm_thms    : (thm*thm) list (* [(?x>>norm(?x), norm(?x)>>?x), ...] *)
clasohm@0
    19
  val red1         : thm        (*  ?P>>?Q  ==>  ?P  ==>  ?Q  *)
clasohm@0
    20
  val red2         : thm        (*  ?P>>?Q  ==>  ?Q  ==>  ?P  *)
clasohm@0
    21
  val refl_thms    : thm list
clasohm@0
    22
  val subst_thms   : thm list   (* [ ?a>>?b ==> ?P(?a) ==> ?P(?b), ...] *)
clasohm@0
    23
  val trans_thms   : thm list
clasohm@0
    24
end;
clasohm@0
    25
clasohm@0
    26
clasohm@0
    27
infix 4 addrews addcongs delrews delcongs setauto;
clasohm@0
    28
clasohm@0
    29
signature SIMP =
clasohm@0
    30
sig
clasohm@0
    31
  type simpset
clasohm@0
    32
  val empty_ss  : simpset
clasohm@0
    33
  val addcongs  : simpset * thm list -> simpset
clasohm@0
    34
  val addrews   : simpset * thm list -> simpset
clasohm@0
    35
  val delcongs  : simpset * thm list -> simpset
clasohm@0
    36
  val delrews   : simpset * thm list -> simpset
clasohm@0
    37
  val dest_ss   : simpset -> thm list * thm list
clasohm@0
    38
  val print_ss  : simpset -> unit
clasohm@0
    39
  val setauto   : simpset * (int -> tactic) -> simpset
clasohm@0
    40
  val ASM_SIMP_CASE_TAC : simpset -> int -> tactic
clasohm@0
    41
  val ASM_SIMP_TAC      : simpset -> int -> tactic
clasohm@0
    42
  val CASE_TAC          : simpset -> int -> tactic
clasohm@0
    43
  val SIMP_CASE2_TAC    : simpset -> int -> tactic
clasohm@0
    44
  val SIMP_THM          : simpset -> thm -> thm
clasohm@0
    45
  val SIMP_TAC          : simpset -> int -> tactic
clasohm@0
    46
  val SIMP_CASE_TAC     : simpset -> int -> tactic
clasohm@0
    47
  val mk_congs          : theory -> string list -> thm list
clasohm@0
    48
  val mk_typed_congs    : theory -> (string * string) list -> thm list
clasohm@0
    49
(* temporarily disabled:
clasohm@0
    50
  val extract_free_congs        : unit -> thm list
clasohm@0
    51
*)
clasohm@0
    52
  val tracing   : bool ref
clasohm@0
    53
end;
clasohm@0
    54
wenzelm@32449
    55
functor SimpFun (Simp_data: SIMP_DATA) : SIMP =
clasohm@0
    56
struct
clasohm@0
    57
wenzelm@19805
    58
local open Simp_data in
clasohm@0
    59
clasohm@0
    60
(*For taking apart reductions into left, right hand sides*)
clasohm@0
    61
val lhs_of = #2 o dest_red;
clasohm@0
    62
val rhs_of = #3 o dest_red;
clasohm@0
    63
clasohm@0
    64
(*** Indexing and filtering of theorems ***)
clasohm@0
    65
wenzelm@22360
    66
fun eq_brl ((b1 : bool, th1), (b2, th2)) = b1 = b2 andalso Thm.eq_thm_prop (th1, th2);
clasohm@0
    67
clasohm@0
    68
(*insert a thm in a discrimination net by its lhs*)
clasohm@0
    69
fun lhs_insert_thm (th,net) =
wenzelm@16800
    70
    Net.insert_term eq_brl (lhs_of (concl_of th), (false,th)) net
clasohm@0
    71
    handle  Net.INSERT => net;
clasohm@0
    72
clasohm@0
    73
(*match subgoal i against possible theorems in the net.
clasohm@0
    74
  Similar to match_from_nat_tac, but the net does not contain numbers;
clasohm@0
    75
  rewrite rules are not ordered.*)
clasohm@0
    76
fun net_tac net =
wenzelm@32449
    77
  SUBGOAL(fn (prem,i) =>
wenzelm@19805
    78
          resolve_tac (Net.unify_term net (Logic.strip_assums_concl prem)) i);
clasohm@0
    79
clasohm@0
    80
(*match subgoal i against possible theorems indexed by lhs in the net*)
clasohm@0
    81
fun lhs_net_tac net =
wenzelm@32449
    82
  SUBGOAL(fn (prem,i) =>
clasohm@1459
    83
          biresolve_tac (Net.unify_term net
wenzelm@19805
    84
                       (lhs_of (Logic.strip_assums_concl prem))) i);
clasohm@0
    85
skalberg@15570
    86
fun nth_subgoal i thm = List.nth(prems_of thm,i-1);
clasohm@0
    87
wenzelm@19805
    88
fun goal_concl i thm = Logic.strip_assums_concl (nth_subgoal i thm);
clasohm@0
    89
clasohm@0
    90
fun lhs_of_eq i thm = lhs_of(goal_concl i thm)
clasohm@0
    91
and rhs_of_eq i thm = rhs_of(goal_concl i thm);
clasohm@0
    92
clasohm@0
    93
fun var_lhs(thm,i) =
clasohm@0
    94
let fun var(Var _) = true
clasohm@0
    95
      | var(Abs(_,_,t)) = var t
clasohm@0
    96
      | var(f$_) = var f
clasohm@0
    97
      | var _ = false;
clasohm@0
    98
in var(lhs_of_eq i thm) end;
clasohm@0
    99
clasohm@0
   100
fun contains_op opns =
clasohm@0
   101
    let fun contains(Const(s,_)) = s mem opns |
clasohm@0
   102
            contains(s$t) = contains s orelse contains t |
clasohm@0
   103
            contains(Abs(_,_,t)) = contains t |
clasohm@0
   104
            contains _ = false;
clasohm@0
   105
    in contains end;
clasohm@0
   106
clasohm@0
   107
fun may_match(match_ops,i) = contains_op match_ops o lhs_of_eq i;
clasohm@0
   108
clasohm@0
   109
val (normI_thms,normE_thms) = split_list norm_thms;
clasohm@0
   110
clasohm@0
   111
(*Get the norm constants from norm_thms*)
clasohm@0
   112
val norms =
wenzelm@32449
   113
  let fun norm thm =
clasohm@0
   114
      case lhs_of(concl_of thm) of
clasohm@1459
   115
          Const(n,_)$_ => n
wenzelm@32091
   116
        | _ => error "No constant in lhs of a norm_thm"
clasohm@0
   117
  in map norm normE_thms end;
clasohm@0
   118
clasohm@0
   119
fun lhs_is_NORM(thm,i) = case lhs_of_eq i thm of
clasohm@1459
   120
        Const(s,_)$_ => s mem norms | _ => false;
clasohm@0
   121
clasohm@0
   122
val refl_tac = resolve_tac refl_thms;
clasohm@0
   123
clasohm@0
   124
fun find_res thms thm =
wenzelm@32091
   125
    let fun find [] = error "Check Simp_Data"
wenzelm@6969
   126
          | find(th::thms) = thm RS th handle THM _ => find thms
clasohm@0
   127
    in find thms end;
clasohm@0
   128
clasohm@0
   129
val mk_trans = find_res trans_thms;
clasohm@0
   130
clasohm@0
   131
fun mk_trans2 thm =
clasohm@0
   132
let fun mk[] = error"Check transitivity"
wenzelm@6969
   133
      | mk(t::ts) = (thm RSN (2,t))  handle THM _  => mk ts
clasohm@0
   134
in mk trans_thms end;
clasohm@0
   135
clasohm@0
   136
(*Applies tactic and returns the first resulting state, FAILS if none!*)
wenzelm@4271
   137
fun one_result(tac,thm) = case Seq.pull(tac thm) of
skalberg@15531
   138
        SOME(thm',_) => thm'
skalberg@15531
   139
      | NONE => raise THM("Simplifier: could not continue", 0, [thm]);
clasohm@0
   140
clasohm@0
   141
fun res1(thm,thms,i) = one_result(resolve_tac thms i,thm);
clasohm@0
   142
clasohm@0
   143
clasohm@0
   144
(**** Adding "NORM" tags ****)
clasohm@0
   145
clasohm@0
   146
(*get name of the constant from conclusion of a congruence rule*)
wenzelm@32449
   147
fun cong_const cong =
clasohm@0
   148
    case head_of (lhs_of (concl_of cong)) of
clasohm@1459
   149
        Const(c,_) => c
clasohm@1459
   150
      | _ => ""                 (*a placeholder distinct from const names*);
clasohm@0
   151
clasohm@0
   152
(*true if the term is an atomic proposition (no ==> signs) *)
wenzelm@19805
   153
val atomic = null o Logic.strip_assums_hyp;
clasohm@0
   154
clasohm@0
   155
(*ccs contains the names of the constants possessing congruence rules*)
clasohm@0
   156
fun add_hidden_vars ccs =
haftmann@21078
   157
  let fun add_hvars tm hvars = case tm of
wenzelm@29265
   158
              Abs(_,_,body) => OldTerm.add_term_vars(body,hvars)
wenzelm@32449
   159
            | _$_ => let val (f,args) = strip_comb tm
clasohm@1459
   160
                     in case f of
wenzelm@32449
   161
                            Const(c,T) =>
haftmann@21078
   162
                                if member (op =) ccs c
haftmann@21078
   163
                                then fold_rev add_hvars args hvars
wenzelm@29265
   164
                                else OldTerm.add_term_vars (tm, hvars)
wenzelm@29265
   165
                          | _ => OldTerm.add_term_vars (tm, hvars)
clasohm@1459
   166
                     end
clasohm@1459
   167
            | _ => hvars;
clasohm@0
   168
  in add_hvars end;
clasohm@0
   169
clasohm@0
   170
fun add_new_asm_vars new_asms =
haftmann@21078
   171
    let fun itf (tm, at) vars =
wenzelm@29265
   172
                if at then vars else OldTerm.add_term_vars(tm,vars)
clasohm@1459
   173
        fun add_list(tm,al,vars) = let val (_,tml) = strip_comb tm
clasohm@1459
   174
                in if length(tml)=length(al)
haftmann@21078
   175
                   then fold_rev itf (tml ~~ al) vars
clasohm@1459
   176
                   else vars
clasohm@1459
   177
                end
clasohm@1459
   178
        fun add_vars (tm,vars) = case tm of
clasohm@1459
   179
                  Abs (_,_,body) => add_vars(body,vars)
clasohm@1459
   180
                | r$s => (case head_of tm of
haftmann@17325
   181
                          Const(c,T) => (case AList.lookup (op =) new_asms c of
skalberg@15531
   182
                                  NONE => add_vars(r,add_vars(s,vars))
skalberg@15531
   183
                                | SOME(al) => add_list(tm,al,vars))
clasohm@1459
   184
                        | _ => add_vars(r,add_vars(s,vars)))
clasohm@1459
   185
                | _ => vars
clasohm@0
   186
    in add_vars end;
clasohm@0
   187
clasohm@0
   188
clasohm@0
   189
fun add_norms(congs,ccs,new_asms) thm =
clasohm@0
   190
let val thm' = mk_trans2 thm;
clasohm@0
   191
(* thm': [?z -> l; Prems; r -> ?t] ==> ?z -> ?t *)
clasohm@0
   192
    val nops = nprems_of thm'
clasohm@0
   193
    val lhs = rhs_of_eq 1 thm'
clasohm@0
   194
    val rhs = lhs_of_eq nops thm'
clasohm@0
   195
    val asms = tl(rev(tl(prems_of thm')))
haftmann@21078
   196
    val hvars = fold_rev (add_hidden_vars ccs) (lhs::rhs::asms) []
clasohm@0
   197
    val hvars = add_new_asm_vars new_asms (rhs,hvars)
haftmann@21078
   198
    fun it_asms asm hvars =
clasohm@1459
   199
        if atomic asm then add_new_asm_vars new_asms (asm,hvars)
wenzelm@29265
   200
        else OldTerm.add_term_frees(asm,hvars)
haftmann@21078
   201
    val hvars = fold_rev it_asms asms hvars
clasohm@0
   202
    val hvs = map (#1 o dest_Var) hvars
clasohm@0
   203
    val refl1_tac = refl_tac 1
paulson@3537
   204
    fun norm_step_tac st = st |>
wenzelm@32449
   205
         (case head_of(rhs_of_eq 1 st) of
wenzelm@32449
   206
            Var(ixn,_) => if ixn mem hvs then refl1_tac
wenzelm@32449
   207
                          else resolve_tac normI_thms 1 ORELSE refl1_tac
wenzelm@32449
   208
          | Const _ => resolve_tac normI_thms 1 ORELSE
wenzelm@32449
   209
                       resolve_tac congs 1 ORELSE refl1_tac
wenzelm@32449
   210
          | Free _ => resolve_tac congs 1 ORELSE refl1_tac
wenzelm@32449
   211
          | _ => refl1_tac)
paulson@3537
   212
    val add_norm_tac = DEPTH_FIRST (has_fewer_prems nops) norm_step_tac
skalberg@15531
   213
    val SOME(thm'',_) = Seq.pull(add_norm_tac thm')
clasohm@0
   214
in thm'' end;
clasohm@0
   215
clasohm@0
   216
fun add_norm_tags congs =
clasohm@0
   217
    let val ccs = map cong_const congs
skalberg@15570
   218
        val new_asms = List.filter (exists not o #2)
clasohm@1459
   219
                (ccs ~~ (map (map atomic o prems_of) congs));
clasohm@0
   220
    in add_norms(congs,ccs,new_asms) end;
clasohm@0
   221
clasohm@0
   222
fun normed_rews congs =
wenzelm@19925
   223
  let val add_norms = add_norm_tags congs in
wenzelm@21287
   224
    fn thm => Variable.tradeT
wenzelm@21287
   225
      (K (map (add_norms o mk_trans) o maps mk_rew_rules)) (Variable.thm_context thm) [thm]
clasohm@0
   226
  end;
clasohm@0
   227
clasohm@1459
   228
fun NORM norm_lhs_tac = EVERY'[rtac red2 , norm_lhs_tac, refl_tac];
clasohm@0
   229
clasohm@0
   230
val trans_norms = map mk_trans normE_thms;
clasohm@0
   231
clasohm@0
   232
clasohm@0
   233
(* SIMPSET *)
clasohm@0
   234
clasohm@0
   235
datatype simpset =
clasohm@1459
   236
        SS of {auto_tac: int -> tactic,
clasohm@1459
   237
               congs: thm list,
clasohm@1459
   238
               cong_net: thm Net.net,
clasohm@1459
   239
               mk_simps: thm -> thm list,
clasohm@1459
   240
               simps: (thm * thm list) list,
clasohm@1459
   241
               simp_net: thm Net.net}
clasohm@0
   242
clasohm@0
   243
val empty_ss = SS{auto_tac= K no_tac, congs=[], cong_net=Net.empty,
clasohm@1459
   244
                  mk_simps=normed_rews[], simps=[], simp_net=Net.empty};
clasohm@0
   245
clasohm@0
   246
(** Insertion of congruences and rewrites **)
clasohm@0
   247
clasohm@0
   248
(*insert a thm in a thm net*)
wenzelm@32449
   249
fun insert_thm_warn th net =
wenzelm@22360
   250
  Net.insert_term Thm.eq_thm_prop (concl_of th, th) net
wenzelm@32449
   251
  handle Net.INSERT =>
wenzelm@32091
   252
    (writeln ("Duplicate rewrite or congruence rule:\n" ^
wenzelm@32091
   253
        Display.string_of_thm_without_context th); net);
clasohm@0
   254
haftmann@21078
   255
val insert_thms = fold_rev insert_thm_warn;
clasohm@0
   256
clasohm@0
   257
fun addrew(SS{auto_tac,congs,cong_net,mk_simps,simps,simp_net}, thm) =
clasohm@0
   258
let val thms = mk_simps thm
clasohm@0
   259
in SS{auto_tac=auto_tac,congs=congs, cong_net=cong_net, mk_simps=mk_simps,
haftmann@21078
   260
      simps = (thm,thms)::simps, simp_net = insert_thms thms simp_net}
clasohm@0
   261
end;
clasohm@0
   262
skalberg@15570
   263
val op addrews = Library.foldl addrew;
clasohm@0
   264
clasohm@0
   265
fun op addcongs(SS{auto_tac,congs,cong_net,mk_simps,simps,simp_net}, thms) =
clasohm@0
   266
let val congs' = thms @ congs;
clasohm@0
   267
in SS{auto_tac=auto_tac, congs= congs',
haftmann@21078
   268
      cong_net= insert_thms (map mk_trans thms) cong_net,
clasohm@0
   269
      mk_simps= normed_rews congs', simps=simps, simp_net=simp_net}
clasohm@0
   270
end;
clasohm@0
   271
clasohm@0
   272
(** Deletion of congruences and rewrites **)
clasohm@0
   273
clasohm@0
   274
(*delete a thm from a thm net*)
wenzelm@32449
   275
fun delete_thm_warn th net =
wenzelm@22360
   276
  Net.delete_term Thm.eq_thm_prop (concl_of th, th) net
wenzelm@32449
   277
  handle Net.DELETE =>
wenzelm@32091
   278
    (writeln ("No such rewrite or congruence rule:\n" ^
wenzelm@32091
   279
        Display.string_of_thm_without_context th); net);
clasohm@0
   280
haftmann@21078
   281
val delete_thms = fold_rev delete_thm_warn;
clasohm@0
   282
clasohm@0
   283
fun op delcongs(SS{auto_tac,congs,cong_net,mk_simps,simps,simp_net}, thms) =
wenzelm@22360
   284
let val congs' = fold (remove Thm.eq_thm_prop) thms congs
clasohm@0
   285
in SS{auto_tac=auto_tac, congs= congs',
haftmann@21078
   286
      cong_net= delete_thms (map mk_trans thms) cong_net,
clasohm@0
   287
      mk_simps= normed_rews congs', simps=simps, simp_net=simp_net}
clasohm@0
   288
end;
clasohm@0
   289
clasohm@0
   290
fun delrew(SS{auto_tac,congs,cong_net,mk_simps,simps,simp_net}, thm) =
clasohm@0
   291
let fun find((p as (th,ths))::ps',ps) =
wenzelm@22360
   292
          if Thm.eq_thm_prop(thm,th) then (ths,ps@ps') else find(ps',p::ps)
wenzelm@32091
   293
      | find([],simps') =
wenzelm@32091
   294
          (writeln ("No such rewrite or congruence rule:\n" ^
wenzelm@32091
   295
              Display.string_of_thm_without_context thm); ([], simps'))
clasohm@0
   296
    val (thms,simps') = find(simps,[])
clasohm@0
   297
in SS{auto_tac=auto_tac, congs=congs, cong_net=cong_net, mk_simps=mk_simps,
haftmann@21078
   298
      simps = simps', simp_net = delete_thms thms simp_net }
clasohm@0
   299
end;
clasohm@0
   300
skalberg@15570
   301
val op delrews = Library.foldl delrew;
clasohm@0
   302
clasohm@0
   303
clasohm@0
   304
fun op setauto(SS{congs,cong_net,mk_simps,simps,simp_net,...}, auto_tac) =
clasohm@0
   305
    SS{auto_tac=auto_tac, congs=congs, cong_net=cong_net, mk_simps=mk_simps,
clasohm@0
   306
       simps=simps, simp_net=simp_net};
clasohm@0
   307
clasohm@0
   308
clasohm@0
   309
(** Inspection of a simpset **)
clasohm@0
   310
clasohm@0
   311
fun dest_ss(SS{congs,simps,...}) = (congs, map #1 simps);
clasohm@0
   312
clasohm@0
   313
fun print_ss(SS{congs,simps,...}) =
wenzelm@32091
   314
  writeln (cat_lines
wenzelm@32091
   315
   (["Congruences:"] @ map Display.string_of_thm_without_context congs @
wenzelm@32091
   316
    ["Rewrite Rules:"] @ map (Display.string_of_thm_without_context o #1) simps));
clasohm@0
   317
clasohm@0
   318
clasohm@0
   319
(* Rewriting with conditionals *)
clasohm@0
   320
clasohm@0
   321
val (case_thms,case_consts) = split_list case_splits;
clasohm@0
   322
val case_rews = map mk_trans case_thms;
clasohm@0
   323
clasohm@0
   324
fun if_rewritable ifc i thm =
clasohm@0
   325
    let val tm = goal_concl i thm
clasohm@1459
   326
        fun nobound(Abs(_,_,tm),j,k) = nobound(tm,j,k+1)
clasohm@1459
   327
          | nobound(s$t,j,k) = nobound(s,j,k) andalso nobound(t,j,k)
clasohm@1459
   328
          | nobound(Bound n,j,k) = n < k orelse k+j <= n
clasohm@1459
   329
          | nobound(_) = true;
clasohm@1459
   330
        fun check_args(al,j) = forall (fn t => nobound(t,j,0)) al
clasohm@1459
   331
        fun find_if(Abs(_,_,tm),j) = find_if(tm,j+1)
clasohm@1459
   332
          | find_if(tm as s$t,j) = let val (f,al) = strip_comb tm in
clasohm@1459
   333
                case f of Const(c,_) => if c=ifc then check_args(al,j)
clasohm@1459
   334
                        else find_if(s,j) orelse find_if(t,j)
clasohm@1459
   335
                | _ => find_if(s,j) orelse find_if(t,j) end
clasohm@1459
   336
          | find_if(_) = false;
clasohm@0
   337
    in find_if(tm,0) end;
clasohm@0
   338
clasohm@0
   339
fun IF1_TAC cong_tac i =
wenzelm@32449
   340
    let fun seq_try (ifth::ifths,ifc::ifcs) thm =
paulson@1512
   341
                (COND (if_rewritable ifc i) (DETERM(rtac ifth i))
paulson@1512
   342
                        (seq_try(ifths,ifcs))) thm
paulson@1512
   343
              | seq_try([],_) thm = no_tac thm
paulson@1512
   344
        and try_rew thm = (seq_try(case_rews,case_consts) ORELSE one_subt) thm
clasohm@1459
   345
        and one_subt thm =
clasohm@1459
   346
                let val test = has_fewer_prems (nprems_of thm + 1)
wenzelm@32449
   347
                    fun loop thm =
wenzelm@32449
   348
                        COND test no_tac
paulson@1512
   349
                          ((try_rew THEN DEPTH_FIRST test (refl_tac i))
wenzelm@32449
   350
                           ORELSE (refl_tac i THEN loop)) thm
paulson@1512
   351
                in (cong_tac THEN loop) thm end
paulson@1512
   352
    in COND (may_match(case_consts,i)) try_rew no_tac end;
clasohm@0
   353
clasohm@0
   354
fun CASE_TAC (SS{cong_net,...}) i =
clasohm@0
   355
let val cong_tac = net_tac cong_net i
clasohm@0
   356
in NORM (IF1_TAC cong_tac) i end;
clasohm@0
   357
clasohm@0
   358
(* Rewriting Automaton *)
clasohm@0
   359
clasohm@0
   360
datatype cntrl = STOP | MK_EQ | ASMS of int | SIMP_LHS | REW | REFL | TRUE
clasohm@1459
   361
               | PROVE | POP_CS | POP_ARTR | IF;
wenzelm@22578
   362
clasohm@0
   363
fun simp_refl([],_,ss) = ss
clasohm@0
   364
  | simp_refl(a'::ns,a,ss) = if a'=a then simp_refl(ns,a,SIMP_LHS::REFL::ss)
clasohm@1459
   365
        else simp_refl(ns,a,ASMS(a)::SIMP_LHS::REFL::POP_ARTR::ss);
clasohm@0
   366
clasohm@0
   367
(** Tracing **)
clasohm@0
   368
clasohm@0
   369
val tracing = ref false;
clasohm@0
   370
clasohm@0
   371
(*Replace parameters by Free variables in P*)
clasohm@0
   372
fun variants_abs ([],P) = P
clasohm@0
   373
  | variants_abs ((a,T)::aTs, P) =
wenzelm@20194
   374
      variants_abs (aTs, #2 (Syntax.variant_abs(a,T,P)));
clasohm@0
   375
clasohm@0
   376
(*Select subgoal i from proof state; substitute parameters, for printing*)
clasohm@0
   377
fun prepare_goal i st =
clasohm@0
   378
    let val subgi = nth_subgoal i st
wenzelm@19805
   379
        val params = rev (Logic.strip_params subgi)
wenzelm@19805
   380
    in variants_abs (params, Logic.strip_assums_concl subgi) end;
clasohm@0
   381
clasohm@0
   382
(*print lhs of conclusion of subgoal i*)
clasohm@0
   383
fun pr_goal_lhs i st =
wenzelm@32449
   384
    writeln (Syntax.string_of_term_global (Thm.theory_of_thm st)
clasohm@1459
   385
             (lhs_of (prepare_goal i st)));
clasohm@0
   386
clasohm@0
   387
(*print conclusion of subgoal i*)
clasohm@0
   388
fun pr_goal_concl i st =
wenzelm@32449
   389
    writeln (Syntax.string_of_term_global (Thm.theory_of_thm st) (prepare_goal i st))
clasohm@0
   390
clasohm@0
   391
(*print subgoals i to j (inclusive)*)
clasohm@0
   392
fun pr_goals (i,j) st =
clasohm@0
   393
    if i>j then ()
clasohm@0
   394
    else (pr_goal_concl i st;  pr_goals (i+1,j) st);
clasohm@0
   395
clasohm@0
   396
(*Print rewrite for tracing; i=subgoal#, n=number of new subgoals,
clasohm@0
   397
  thm=old state, thm'=new state *)
clasohm@0
   398
fun pr_rew (i,n,thm,thm',not_asms) =
clasohm@0
   399
    if !tracing
clasohm@0
   400
    then (if not_asms then () else writeln"Assumption used in";
clasohm@0
   401
          pr_goal_lhs i thm; writeln"->"; pr_goal_lhs (i+n) thm';
clasohm@1459
   402
          if n>0 then (writeln"Conditions:"; pr_goals (i, i+n-1) thm')
clasohm@0
   403
          else ();
clasohm@0
   404
          writeln"" )
clasohm@0
   405
    else ();
clasohm@0
   406
clasohm@0
   407
(* Skip the first n hyps of a goal, and return the rest in generalized form *)
clasohm@0
   408
fun strip_varify(Const("==>", _) $ H $ B, n, vs) =
clasohm@1459
   409
        if n=0 then subst_bounds(vs,H)::strip_varify(B,0,vs)
clasohm@1459
   410
        else strip_varify(B,n-1,vs)
clasohm@0
   411
  | strip_varify(Const("all",_)$Abs(_,T,t), n, vs) =
clasohm@1459
   412
        strip_varify(t,n,Var(("?",length vs),T)::vs)
clasohm@0
   413
  | strip_varify  _  = [];
clasohm@0
   414
clasohm@0
   415
fun execute(ss,if_fl,auto_tac,cong_tac,net,i,thm) = let
clasohm@0
   416
clasohm@0
   417
fun simp_lhs(thm,ss,anet,ats,cs) =
clasohm@0
   418
    if var_lhs(thm,i) then (ss,thm,anet,ats,cs) else
clasohm@0
   419
    if lhs_is_NORM(thm,i) then (ss, res1(thm,trans_norms,i), anet,ats,cs)
wenzelm@4271
   420
    else case Seq.pull(cong_tac i thm) of
skalberg@15531
   421
            SOME(thm',_) =>
clasohm@1459
   422
                    let val ps = prems_of thm and ps' = prems_of thm';
clasohm@1459
   423
                        val n = length(ps')-length(ps);
wenzelm@19805
   424
                        val a = length(Logic.strip_assums_hyp(List.nth(ps,i-1)))
wenzelm@19805
   425
                        val l = map (fn p => length(Logic.strip_assums_hyp(p)))
skalberg@15570
   426
                                    (Library.take(n,Library.drop(i-1,ps')));
clasohm@1459
   427
                    in (simp_refl(rev(l),a,REW::ss),thm',anet,ats,cs) end
skalberg@15531
   428
          | NONE => (REW::ss,thm,anet,ats,cs);
clasohm@0
   429
clasohm@0
   430
(*NB: the "Adding rewrites:" trace will look strange because assumptions
clasohm@0
   431
      are represented by rules, generalized over their parameters*)
clasohm@0
   432
fun add_asms(ss,thm,a,anet,ats,cs) =
clasohm@0
   433
    let val As = strip_varify(nth_subgoal i thm, a, []);
wenzelm@22596
   434
        val thms = map (trivial o cterm_of(Thm.theory_of_thm thm)) As;
skalberg@15570
   435
        val new_rws = List.concat(map mk_rew_rules thms);
skalberg@15570
   436
        val rwrls = map mk_trans (List.concat(map mk_rew_rules thms));
wenzelm@30190
   437
        val anet' = List.foldr lhs_insert_thm anet rwrls
clasohm@0
   438
    in  if !tracing andalso not(null new_rws)
wenzelm@32091
   439
        then writeln (cat_lines
wenzelm@32091
   440
          ("Adding rewrites:" :: map Display.string_of_thm_without_context new_rws))
clasohm@1459
   441
        else ();
wenzelm@32449
   442
        (ss,thm,anet',anet::ats,cs)
clasohm@0
   443
    end;
clasohm@0
   444
wenzelm@4271
   445
fun rew(seq,thm,ss,anet,ats,cs, more) = case Seq.pull seq of
skalberg@15531
   446
      SOME(thm',seq') =>
clasohm@1459
   447
            let val n = (nprems_of thm') - (nprems_of thm)
clasohm@1459
   448
            in pr_rew(i,n,thm,thm',more);
clasohm@1459
   449
               if n=0 then (SIMP_LHS::ss, thm', anet, ats, cs)
clasohm@1459
   450
               else ((replicate n PROVE) @ (POP_CS::SIMP_LHS::ss),
clasohm@1459
   451
                     thm', anet, ats, (ss,thm,anet,ats,seq',more)::cs)
clasohm@1459
   452
            end
skalberg@15531
   453
    | NONE => if more
paulson@1512
   454
            then rew((lhs_net_tac anet i THEN assume_tac i) thm,
clasohm@1459
   455
                     thm,ss,anet,ats,cs,false)
clasohm@1459
   456
            else (ss,thm,anet,ats,cs);
clasohm@0
   457
clasohm@0
   458
fun try_true(thm,ss,anet,ats,cs) =
wenzelm@4271
   459
    case Seq.pull(auto_tac i thm) of
skalberg@15531
   460
      SOME(thm',_) => (ss,thm',anet,ats,cs)
skalberg@15531
   461
    | NONE => let val (ss0,thm0,anet0,ats0,seq,more)::cs0 = cs
clasohm@1459
   462
              in if !tracing
clasohm@1459
   463
                 then (writeln"*** Failed to prove precondition. Normal form:";
clasohm@1459
   464
                       pr_goal_concl i thm;  writeln"")
clasohm@1459
   465
                 else ();
clasohm@1459
   466
                 rew(seq,thm0,ss0,anet0,ats0,cs0,more)
clasohm@1459
   467
              end;
clasohm@0
   468
clasohm@0
   469
fun if_exp(thm,ss,anet,ats,cs) =
wenzelm@4271
   470
        case Seq.pull (IF1_TAC (cong_tac i) i thm) of
skalberg@15531
   471
                SOME(thm',_) => (SIMP_LHS::IF::ss,thm',anet,ats,cs)
skalberg@15531
   472
              | NONE => (ss,thm,anet,ats,cs);
clasohm@0
   473
clasohm@0
   474
fun step(s::ss, thm, anet, ats, cs) = case s of
clasohm@1459
   475
          MK_EQ => (ss, res1(thm,[red2],i), anet, ats, cs)
clasohm@1459
   476
        | ASMS(a) => add_asms(ss,thm,a,anet,ats,cs)
clasohm@1459
   477
        | SIMP_LHS => simp_lhs(thm,ss,anet,ats,cs)
paulson@1512
   478
        | REW => rew(net_tac net i thm,thm,ss,anet,ats,cs,true)
clasohm@1459
   479
        | REFL => (ss, res1(thm,refl_thms,i), anet, ats, cs)
clasohm@1459
   480
        | TRUE => try_true(res1(thm,refl_thms,i),ss,anet,ats,cs)
clasohm@1459
   481
        | PROVE => (if if_fl then MK_EQ::SIMP_LHS::IF::TRUE::ss
clasohm@1459
   482
                    else MK_EQ::SIMP_LHS::TRUE::ss, thm, anet, ats, cs)
clasohm@1459
   483
        | POP_ARTR => (ss,thm,hd ats,tl ats,cs)
clasohm@1459
   484
        | POP_CS => (ss,thm,anet,ats,tl cs)
clasohm@1459
   485
        | IF => if_exp(thm,ss,anet,ats,cs);
clasohm@0
   486
clasohm@0
   487
fun exec(state as (s::ss, thm, _, _, _)) =
clasohm@1459
   488
        if s=STOP then thm else exec(step(state));
clasohm@0
   489
clasohm@0
   490
in exec(ss, thm, Net.empty, [], []) end;
clasohm@0
   491
clasohm@0
   492
clasohm@0
   493
fun EXEC_TAC(ss,fl) (SS{auto_tac,cong_net,simp_net,...}) =
clasohm@0
   494
let val cong_tac = net_tac cong_net
wenzelm@32449
   495
in fn i =>
paulson@1512
   496
    (fn thm =>
wenzelm@4271
   497
     if i <= 0 orelse nprems_of thm < i then Seq.empty
wenzelm@4271
   498
     else Seq.single(execute(ss,fl,auto_tac,cong_tac,simp_net,i,thm)))
paulson@1512
   499
    THEN TRY(auto_tac i)
clasohm@0
   500
end;
clasohm@0
   501
clasohm@0
   502
val SIMP_TAC = EXEC_TAC([MK_EQ,SIMP_LHS,REFL,STOP],false);
clasohm@0
   503
val SIMP_CASE_TAC = EXEC_TAC([MK_EQ,SIMP_LHS,IF,REFL,STOP],false);
clasohm@0
   504
clasohm@0
   505
val ASM_SIMP_TAC = EXEC_TAC([ASMS(0),MK_EQ,SIMP_LHS,REFL,STOP],false);
clasohm@0
   506
val ASM_SIMP_CASE_TAC = EXEC_TAC([ASMS(0),MK_EQ,SIMP_LHS,IF,REFL,STOP],false);
clasohm@0
   507
clasohm@0
   508
val SIMP_CASE2_TAC = EXEC_TAC([MK_EQ,SIMP_LHS,IF,REFL,STOP],true);
clasohm@0
   509
clasohm@0
   510
fun REWRITE (ss,fl) (SS{auto_tac,cong_net,simp_net,...}) =
clasohm@0
   511
let val cong_tac = net_tac cong_net
clasohm@0
   512
in fn thm => let val state = thm RSN (2,red1)
clasohm@1459
   513
             in execute(ss,fl,auto_tac,cong_tac,simp_net,1,state) end
clasohm@0
   514
end;
clasohm@0
   515
clasohm@0
   516
val SIMP_THM = REWRITE ([ASMS(0),SIMP_LHS,IF,REFL,STOP],false);
clasohm@0
   517
clasohm@0
   518
clasohm@0
   519
(* Compute Congruence rules for individual constants using the substition
clasohm@0
   520
   rules *)
clasohm@0
   521
clasohm@0
   522
val subst_thms = map standard subst_thms;
clasohm@0
   523
clasohm@0
   524
clasohm@0
   525
fun exp_app(0,t) = t
clasohm@0
   526
  | exp_app(i,t) = exp_app(i-1,t $ Bound (i-1));
clasohm@0
   527
clasohm@0
   528
fun exp_abs(Type("fun",[T1,T2]),t,i) =
clasohm@1459
   529
        Abs("x"^string_of_int i,T1,exp_abs(T2,t,i+1))
clasohm@0
   530
  | exp_abs(T,t,i) = exp_app(i,t);
clasohm@0
   531
clasohm@0
   532
fun eta_Var(ixn,T) = exp_abs(T,Var(ixn,T),0);
clasohm@0
   533
clasohm@0
   534
clasohm@0
   535
fun Pinst(f,fT,(eq,eqT),k,i,T,yik,Ts) =
clasohm@0
   536
let fun xn_list(x,n) =
clasohm@1459
   537
        let val ixs = map (fn i => (x^(radixstring(26,"a",i)),0)) (0 upto n);
skalberg@15570
   538
        in ListPair.map eta_Var (ixs, Library.take(n+1,Ts)) end
clasohm@0
   539
    val lhs = list_comb(f,xn_list("X",k-1))
clasohm@0
   540
    val rhs = list_comb(f,xn_list("X",i-1) @ [Bound 0] @ yik)
clasohm@0
   541
in Abs("", T, Const(eq,[fT,fT]--->eqT) $ lhs $ rhs) end;
clasohm@0
   542
wenzelm@16931
   543
fun find_subst sg T =
clasohm@0
   544
let fun find (thm::thms) =
clasohm@1459
   545
        let val (Const(_,cT), va, vb) = dest_red(hd(prems_of thm));
wenzelm@29265
   546
            val [P] = OldTerm.add_term_vars(concl_of thm,[]) \\ [va,vb]
clasohm@1459
   547
            val eqT::_ = binder_types cT
wenzelm@16931
   548
        in if Sign.typ_instance sg (T,eqT) then SOME(thm,va,vb,P)
clasohm@1459
   549
           else find thms
clasohm@1459
   550
        end
skalberg@15531
   551
      | find [] = NONE
clasohm@0
   552
in find subst_thms end;
clasohm@0
   553
clasohm@0
   554
fun mk_cong sg (f,aTs,rT) (refl,eq) =
wenzelm@16931
   555
let val k = length aTs;
clasohm@0
   556
    fun ri((subst,va as Var(_,Ta),vb as Var(_,Tb),P),i,si,T,yik) =
clasohm@1459
   557
        let val ca = cterm_of sg va
clasohm@1459
   558
            and cx = cterm_of sg (eta_Var(("X"^si,0),T))
clasohm@1459
   559
            val cb = cterm_of sg vb
clasohm@1459
   560
            and cy = cterm_of sg (eta_Var(("Y"^si,0),T))
clasohm@1459
   561
            val cP = cterm_of sg P
clasohm@1459
   562
            and cp = cterm_of sg (Pinst(f,rT,eq,k,i,T,yik,aTs))
clasohm@1459
   563
        in cterm_instantiate [(ca,cx),(cb,cy),(cP,cp)] subst end;
clasohm@0
   564
    fun mk(c,T::Ts,i,yik) =
clasohm@1459
   565
        let val si = radixstring(26,"a",i)
wenzelm@16931
   566
        in case find_subst sg T of
skalberg@15531
   567
             NONE => mk(c,Ts,i-1,eta_Var(("X"^si,0),T)::yik)
skalberg@15531
   568
           | SOME s => let val c' = c RSN (2,ri(s,i,si,T,yik))
clasohm@1459
   569
                       in mk(c',Ts,i-1,eta_Var(("Y"^si,0),T)::yik) end
clasohm@1459
   570
        end
clasohm@0
   571
      | mk(c,[],_,_) = c;
clasohm@0
   572
in mk(refl,rev aTs,k-1,[]) end;
clasohm@0
   573
clasohm@0
   574
fun mk_cong_type sg (f,T) =
clasohm@0
   575
let val (aTs,rT) = strip_type T;
clasohm@0
   576
    fun find_refl(r::rs) =
clasohm@1459
   577
        let val (Const(eq,eqT),_,_) = dest_red(concl_of r)
wenzelm@16931
   578
        in if Sign.typ_instance sg (rT, hd(binder_types eqT))
skalberg@15531
   579
           then SOME(r,(eq,body_type eqT)) else find_refl rs
clasohm@1459
   580
        end
skalberg@15531
   581
      | find_refl([]) = NONE;
clasohm@0
   582
in case find_refl refl_thms of
skalberg@15531
   583
     NONE => []  |  SOME(refl) => [mk_cong sg (f,aTs,rT) refl]
clasohm@0
   584
end;
clasohm@0
   585
clasohm@0
   586
fun mk_cong_thy thy f =
wenzelm@22578
   587
let val T = case Sign.const_type thy f of
skalberg@15531
   588
                NONE => error(f^" not declared") | SOME(T) => T;
wenzelm@16876
   589
    val T' = Logic.incr_tvar 9 T;
wenzelm@22578
   590
in mk_cong_type thy (Const(f,T'),T') end;
clasohm@0
   591
skalberg@15570
   592
fun mk_congs thy = List.concat o map (mk_cong_thy thy);
clasohm@0
   593
clasohm@0
   594
fun mk_typed_congs thy =
wenzelm@22675
   595
let
wenzelm@22675
   596
  fun readfT(f,s) =
wenzelm@22675
   597
    let
wenzelm@24707
   598
      val T = Logic.incr_tvar 9 (Syntax.read_typ_global thy s);
wenzelm@22675
   599
      val t = case Sign.const_type thy f of
wenzelm@22675
   600
                  SOME(_) => Const(f,T) | NONE => Free(f,T)
wenzelm@22675
   601
    in (t,T) end
wenzelm@22578
   602
in List.concat o map (mk_cong_type thy o readfT) end;
clasohm@0
   603
wenzelm@22675
   604
end;
wenzelm@22675
   605
end;