src/HOL/Tools/SMT/smt_normalize.ML
author boehmes
Tue Oct 26 11:39:26 2010 +0200 (2010-10-26)
changeset 40161 539d07b00e5f
parent 39483 9f0e5684f04b
child 40162 7f58a9a843c2
permissions -rw-r--r--
keep track of theorems initially given to SMT (even if they are rewritten); provide interface to filter theorems necessary for SMT proofs
boehmes@36898
     1
(*  Title:      HOL/Tools/SMT/smt_normalize.ML
boehmes@36898
     2
    Author:     Sascha Boehme, TU Muenchen
boehmes@36898
     3
boehmes@36898
     4
Normalization steps on theorems required by SMT solvers:
boehmes@36898
     5
  * simplify trivial distincts (those with less than three elements),
boehmes@36898
     6
  * rewrite bool case expressions as if expressions,
boehmes@36898
     7
  * normalize numerals (e.g. replace negative numerals by negated positive
boehmes@36898
     8
    numerals),
boehmes@36898
     9
  * embed natural numbers into integers,
boehmes@36898
    10
  * add extra rules specifying types and constants which occur frequently,
boehmes@36898
    11
  * fully translate into object logic, add universal closure,
boehmes@39483
    12
  * monomorphize (create instances of schematic rules),
boehmes@36898
    13
  * lift lambda terms,
boehmes@36898
    14
  * make applications explicit for functions with varying number of arguments.
boehmes@39483
    15
  * add (hypothetical definitions for) missing datatype selectors,
boehmes@36898
    16
*)
boehmes@36898
    17
boehmes@36898
    18
signature SMT_NORMALIZE =
boehmes@36898
    19
sig
boehmes@40161
    20
  type extra_norm = (int * thm) list -> Proof.context ->
boehmes@40161
    21
    (int * thm) list * Proof.context
boehmes@40161
    22
  val normalize: extra_norm -> bool -> (int * thm) list -> Proof.context ->
boehmes@40161
    23
    (int * thm) list * Proof.context
boehmes@36899
    24
  val atomize_conv: Proof.context -> conv
boehmes@36898
    25
  val eta_expand_conv: (Proof.context -> conv) -> Proof.context -> conv
boehmes@36898
    26
end
boehmes@36898
    27
boehmes@36898
    28
structure SMT_Normalize: SMT_NORMALIZE =
boehmes@36898
    29
struct
boehmes@36898
    30
boehmes@36898
    31
infix 2 ??
boehmes@36898
    32
fun (test ?? f) x = if test x then f x else x
boehmes@36898
    33
boehmes@36898
    34
fun if_conv c cv1 cv2 ct = (if c (Thm.term_of ct) then cv1 else cv2) ct
boehmes@36898
    35
fun if_true_conv c cv = if_conv c cv Conv.all_conv
boehmes@36898
    36
boehmes@36898
    37
boehmes@36898
    38
boehmes@36898
    39
(* simplification of trivial distincts (distinct should have at least
boehmes@36898
    40
   three elements in the argument list) *)
boehmes@36898
    41
boehmes@36898
    42
local
boehmes@36898
    43
  fun is_trivial_distinct (Const (@{const_name distinct}, _) $ t) =
boehmes@36898
    44
        length (HOLogic.dest_list t) <= 2
boehmes@36898
    45
    | is_trivial_distinct _ = false
boehmes@36898
    46
boehmes@37786
    47
  val thms = map mk_meta_eq @{lemma
boehmes@37786
    48
    "distinct [] = True"
boehmes@37786
    49
    "distinct [x] = True"
boehmes@37786
    50
    "distinct [x, y] = (x ~= y)"
boehmes@36898
    51
    by simp_all}
boehmes@36898
    52
  fun distinct_conv _ =
wenzelm@36936
    53
    if_true_conv is_trivial_distinct (Conv.rewrs_conv thms)
boehmes@36898
    54
in
boehmes@36898
    55
fun trivial_distinct ctxt =
boehmes@40161
    56
  map (apsnd ((Term.exists_subterm is_trivial_distinct o Thm.prop_of) ??
boehmes@40161
    57
    Conv.fconv_rule (Conv.top_conv distinct_conv ctxt)))
boehmes@36898
    58
end
boehmes@36898
    59
boehmes@36898
    60
boehmes@36898
    61
boehmes@36898
    62
(* rewrite bool case expressions as if expressions *)
boehmes@36898
    63
boehmes@36898
    64
local
boehmes@36898
    65
  val is_bool_case = (fn
boehmes@36898
    66
      Const (@{const_name "bool.bool_case"}, _) $ _ $ _ $ _ => true
boehmes@36898
    67
    | _ => false)
boehmes@36898
    68
boehmes@37786
    69
  val thms = map mk_meta_eq @{lemma
boehmes@37786
    70
    "(case P of True => x | False => y) = (if P then x else y)"
boehmes@37786
    71
    "(case P of False => y | True => x) = (if P then x else y)"
boehmes@37786
    72
    by simp_all}
wenzelm@36936
    73
  val unfold_conv = if_true_conv is_bool_case (Conv.rewrs_conv thms)
boehmes@36898
    74
in
boehmes@36898
    75
fun rewrite_bool_cases ctxt =
boehmes@40161
    76
  map (apsnd ((Term.exists_subterm is_bool_case o Thm.prop_of) ??
boehmes@40161
    77
    Conv.fconv_rule (Conv.top_conv (K unfold_conv) ctxt)))
boehmes@36898
    78
end
boehmes@36898
    79
boehmes@36898
    80
boehmes@36898
    81
boehmes@36898
    82
(* normalization of numerals: rewriting of negative integer numerals into
boehmes@36898
    83
   positive numerals, Numeral0 into 0, Numeral1 into 1 *)
boehmes@36898
    84
boehmes@36898
    85
local
boehmes@36898
    86
  fun is_number_sort ctxt T =
boehmes@36898
    87
    Sign.of_sort (ProofContext.theory_of ctxt) (T, @{sort number_ring})
boehmes@36898
    88
boehmes@36898
    89
  fun is_strange_number ctxt (t as Const (@{const_name number_of}, _) $ _) =
boehmes@36898
    90
        (case try HOLogic.dest_number t of
boehmes@36898
    91
          SOME (T, i) => is_number_sort ctxt T andalso i < 2
boehmes@36898
    92
        | NONE => false)
boehmes@36898
    93
    | is_strange_number _ _ = false
boehmes@36898
    94
boehmes@36898
    95
  val pos_numeral_ss = HOL_ss
boehmes@36898
    96
    addsimps [@{thm Int.number_of_minus}, @{thm Int.number_of_Min}]
boehmes@36898
    97
    addsimps [@{thm Int.number_of_Pls}, @{thm Int.numeral_1_eq_1}]
boehmes@36898
    98
    addsimps @{thms Int.pred_bin_simps}
boehmes@36898
    99
    addsimps @{thms Int.normalize_bin_simps}
boehmes@36898
   100
    addsimps @{lemma
boehmes@36898
   101
      "Int.Min = - Int.Bit1 Int.Pls"
boehmes@36898
   102
      "Int.Bit0 (- Int.Pls) = - Int.Pls"
boehmes@36898
   103
      "Int.Bit0 (- k) = - Int.Bit0 k"
boehmes@36898
   104
      "Int.Bit1 (- k) = - Int.Bit1 (Int.pred k)"
boehmes@36898
   105
      by simp_all (simp add: pred_def)}
boehmes@36898
   106
boehmes@36898
   107
  fun pos_conv ctxt = if_conv (is_strange_number ctxt)
boehmes@36898
   108
    (Simplifier.rewrite (Simplifier.context ctxt pos_numeral_ss))
boehmes@36898
   109
    Conv.no_conv
boehmes@36898
   110
in
boehmes@36898
   111
fun normalize_numerals ctxt =
boehmes@40161
   112
  map (apsnd ((Term.exists_subterm (is_strange_number ctxt) o Thm.prop_of) ??
boehmes@40161
   113
    Conv.fconv_rule (Conv.top_sweep_conv pos_conv ctxt)))
boehmes@36898
   114
end
boehmes@36898
   115
boehmes@36898
   116
boehmes@36898
   117
boehmes@36898
   118
(* embedding of standard natural number operations into integer operations *)
boehmes@36898
   119
boehmes@36898
   120
local
boehmes@40161
   121
  val nat_embedding = map (pair ~1) @{lemma
boehmes@36898
   122
    "nat (int n) = n"
boehmes@36898
   123
    "i >= 0 --> int (nat i) = i"
boehmes@36898
   124
    "i < 0 --> int (nat i) = 0"
boehmes@36898
   125
    by simp_all}
boehmes@36898
   126
boehmes@36898
   127
  val nat_rewriting = @{lemma
boehmes@36898
   128
    "0 = nat 0"
boehmes@36898
   129
    "1 = nat 1"
boehmes@36898
   130
    "number_of i = nat (number_of i)"
boehmes@36898
   131
    "int (nat 0) = 0"
boehmes@36898
   132
    "int (nat 1) = 1"
boehmes@36898
   133
    "a < b = (int a < int b)"
boehmes@36898
   134
    "a <= b = (int a <= int b)"
boehmes@36898
   135
    "Suc a = nat (int a + 1)"
boehmes@36898
   136
    "a + b = nat (int a + int b)"
boehmes@36898
   137
    "a - b = nat (int a - int b)"
boehmes@36898
   138
    "a * b = nat (int a * int b)"
boehmes@36898
   139
    "a div b = nat (int a div int b)"
boehmes@36898
   140
    "a mod b = nat (int a mod int b)"
boehmes@36898
   141
    "min a b = nat (min (int a) (int b))"
boehmes@36898
   142
    "max a b = nat (max (int a) (int b))"
boehmes@36898
   143
    "int (nat (int a + int b)) = int a + int b"
boehmes@36898
   144
    "int (nat (int a * int b)) = int a * int b"
boehmes@36898
   145
    "int (nat (int a div int b)) = int a div int b"
boehmes@36898
   146
    "int (nat (int a mod int b)) = int a mod int b"
boehmes@36898
   147
    "int (nat (min (int a) (int b))) = min (int a) (int b)"
boehmes@36898
   148
    "int (nat (max (int a) (int b))) = max (int a) (int b)"
boehmes@36898
   149
    by (simp_all add: nat_mult_distrib nat_div_distrib nat_mod_distrib
boehmes@36898
   150
      int_mult[symmetric] zdiv_int[symmetric] zmod_int[symmetric])}
boehmes@36898
   151
boehmes@36898
   152
  fun on_positive num f x = 
boehmes@36898
   153
    (case try HOLogic.dest_number (Thm.term_of num) of
boehmes@36898
   154
      SOME (_, i) => if i >= 0 then SOME (f x) else NONE
boehmes@36898
   155
    | NONE => NONE)
boehmes@36898
   156
boehmes@36898
   157
  val cancel_int_nat_ss = HOL_ss
boehmes@36898
   158
    addsimps [@{thm Nat_Numeral.nat_number_of}]
boehmes@36898
   159
    addsimps [@{thm Nat_Numeral.int_nat_number_of}]
boehmes@36898
   160
    addsimps @{thms neg_simps}
boehmes@36898
   161
boehmes@36898
   162
  fun cancel_int_nat_simproc _ ss ct = 
boehmes@36898
   163
    let
boehmes@36898
   164
      val num = Thm.dest_arg (Thm.dest_arg ct)
boehmes@36898
   165
      val goal = Thm.mk_binop @{cterm "op == :: int => _"} ct num
boehmes@36898
   166
      val simpset = Simplifier.inherit_context ss cancel_int_nat_ss
boehmes@36898
   167
      fun tac _ = Simplifier.simp_tac simpset 1
boehmes@36898
   168
    in on_positive num (Goal.prove_internal [] goal) tac end
boehmes@36898
   169
boehmes@36898
   170
  val nat_ss = HOL_ss
boehmes@36898
   171
    addsimps nat_rewriting
boehmes@36898
   172
    addsimprocs [Simplifier.make_simproc {
boehmes@36898
   173
      name = "cancel_int_nat_num", lhss = [@{cpat "int (nat _)"}],
boehmes@36898
   174
      proc = cancel_int_nat_simproc, identifier = [] }]
boehmes@36898
   175
boehmes@36898
   176
  fun conv ctxt = Simplifier.rewrite (Simplifier.context ctxt nat_ss)
boehmes@36898
   177
boehmes@36898
   178
  val uses_nat_type = Term.exists_type (Term.exists_subtype (equal @{typ nat}))
boehmes@36898
   179
  val uses_nat_int =
boehmes@36898
   180
    Term.exists_subterm (member (op aconv) [@{term int}, @{term nat}])
boehmes@36898
   181
in
boehmes@36898
   182
fun nat_as_int ctxt =
boehmes@40161
   183
  map (apsnd ((uses_nat_type o Thm.prop_of) ?? Conv.fconv_rule (conv ctxt))) #>
boehmes@40161
   184
  exists (uses_nat_int o Thm.prop_of o snd) ?? append nat_embedding
boehmes@36898
   185
end
boehmes@36898
   186
boehmes@36898
   187
boehmes@36898
   188
boehmes@36898
   189
(* further normalizations: beta/eta, universal closure, atomize *)
boehmes@36898
   190
boehmes@36898
   191
val eta_expand_eq = @{lemma "f == (%x. f x)" by (rule reflexive)}
boehmes@36898
   192
boehmes@36898
   193
fun eta_expand_conv cv ctxt =
boehmes@36898
   194
  Conv.rewr_conv eta_expand_eq then_conv Conv.abs_conv (cv o snd) ctxt
boehmes@36898
   195
boehmes@36898
   196
local
boehmes@36898
   197
  val eta_conv = eta_expand_conv
boehmes@36898
   198
wenzelm@36936
   199
  fun keep_conv ctxt = Conv.binder_conv (norm_conv o snd) ctxt
boehmes@36898
   200
  and eta_binder_conv ctxt = Conv.arg_conv (eta_conv norm_conv ctxt)
boehmes@36898
   201
  and keep_let_conv ctxt = Conv.combination_conv
boehmes@36898
   202
    (Conv.arg_conv (norm_conv ctxt)) (Conv.abs_conv (norm_conv o snd) ctxt)
boehmes@36898
   203
  and unfold_let_conv ctxt = Conv.combination_conv
boehmes@36898
   204
    (Conv.arg_conv (norm_conv ctxt)) (eta_conv norm_conv ctxt)
boehmes@36898
   205
  and unfold_conv thm ctxt = Conv.rewr_conv thm then_conv keep_conv ctxt
boehmes@36898
   206
  and unfold_ex1_conv ctxt = unfold_conv @{thm Ex1_def} ctxt
boehmes@37786
   207
  and unfold_ball_conv ctxt = unfold_conv (mk_meta_eq @{thm Ball_def}) ctxt
boehmes@37786
   208
  and unfold_bex_conv ctxt = unfold_conv (mk_meta_eq @{thm Bex_def}) ctxt
boehmes@36898
   209
  and norm_conv ctxt ct =
boehmes@36898
   210
    (case Thm.term_of ct of
boehmes@36898
   211
      Const (@{const_name All}, _) $ Abs _ => keep_conv
boehmes@36898
   212
    | Const (@{const_name All}, _) $ _ => eta_binder_conv
boehmes@36898
   213
    | Const (@{const_name All}, _) => eta_conv eta_binder_conv
boehmes@36898
   214
    | Const (@{const_name Ex}, _) $ Abs _ => keep_conv
boehmes@36898
   215
    | Const (@{const_name Ex}, _) $ _ => eta_binder_conv
boehmes@36898
   216
    | Const (@{const_name Ex}, _) => eta_conv eta_binder_conv
boehmes@36898
   217
    | Const (@{const_name Let}, _) $ _ $ Abs _ => keep_let_conv
boehmes@36898
   218
    | Const (@{const_name Let}, _) $ _ $ _ => unfold_let_conv
boehmes@36898
   219
    | Const (@{const_name Let}, _) $ _ => eta_conv unfold_let_conv
boehmes@36898
   220
    | Const (@{const_name Let}, _) => eta_conv (eta_conv unfold_let_conv)
boehmes@36898
   221
    | Const (@{const_name Ex1}, _) $ _ => unfold_ex1_conv
boehmes@36898
   222
    | Const (@{const_name Ex1}, _) => eta_conv unfold_ex1_conv 
boehmes@36898
   223
    | Const (@{const_name Ball}, _) $ _ $ _ => unfold_ball_conv
boehmes@36898
   224
    | Const (@{const_name Ball}, _) $ _ => eta_conv unfold_ball_conv
boehmes@36898
   225
    | Const (@{const_name Ball}, _) => eta_conv (eta_conv unfold_ball_conv)
boehmes@36898
   226
    | Const (@{const_name Bex}, _) $ _ $ _ => unfold_bex_conv
boehmes@36898
   227
    | Const (@{const_name Bex}, _) $ _ => eta_conv unfold_bex_conv
boehmes@36898
   228
    | Const (@{const_name Bex}, _) => eta_conv (eta_conv unfold_bex_conv)
boehmes@36898
   229
    | Abs _ => Conv.abs_conv (norm_conv o snd)
boehmes@36898
   230
    | _ $ _ => Conv.comb_conv o norm_conv
boehmes@36898
   231
    | _ => K Conv.all_conv) ctxt ct
boehmes@36898
   232
boehmes@36898
   233
  fun is_normed t =
boehmes@36898
   234
    (case t of
boehmes@36898
   235
      Const (@{const_name All}, _) $ Abs (_, _, u) => is_normed u
boehmes@36898
   236
    | Const (@{const_name All}, _) $ _ => false
boehmes@36898
   237
    | Const (@{const_name All}, _) => false
boehmes@36898
   238
    | Const (@{const_name Ex}, _) $ Abs (_, _, u) => is_normed u
boehmes@36898
   239
    | Const (@{const_name Ex}, _) $ _ => false
boehmes@36898
   240
    | Const (@{const_name Ex}, _) => false
boehmes@36898
   241
    | Const (@{const_name Let}, _) $ u1 $ Abs (_, _, u2) =>
boehmes@36898
   242
        is_normed u1 andalso is_normed u2
boehmes@36898
   243
    | Const (@{const_name Let}, _) $ _ $ _ => false
boehmes@36898
   244
    | Const (@{const_name Let}, _) $ _ => false
boehmes@36898
   245
    | Const (@{const_name Let}, _) => false
boehmes@36898
   246
    | Const (@{const_name Ex1}, _) => false
boehmes@36898
   247
    | Const (@{const_name Ball}, _) => false
boehmes@36898
   248
    | Const (@{const_name Bex}, _) => false
boehmes@36898
   249
    | Abs (_, _, u) => is_normed u
boehmes@36898
   250
    | u1 $ u2 => is_normed u1 andalso is_normed u2
boehmes@36898
   251
    | _ => true)
boehmes@36898
   252
in
boehmes@36898
   253
fun norm_binder_conv ctxt = if_conv is_normed Conv.all_conv (norm_conv ctxt)
boehmes@36898
   254
end
boehmes@36898
   255
boehmes@36898
   256
fun norm_def ctxt thm =
boehmes@36898
   257
  (case Thm.prop_of thm of
haftmann@38864
   258
    @{term Trueprop} $ (Const (@{const_name HOL.eq}, _) $ _ $ Abs _) =>
boehmes@36898
   259
      norm_def ctxt (thm RS @{thm fun_cong})
boehmes@36898
   260
  | Const (@{const_name "=="}, _) $ _ $ Abs _ =>
boehmes@36898
   261
      norm_def ctxt (thm RS @{thm meta_eq_to_obj_eq})
boehmes@36898
   262
  | _ => thm)
boehmes@36898
   263
boehmes@36898
   264
fun atomize_conv ctxt ct =
boehmes@36898
   265
  (case Thm.term_of ct of
boehmes@36898
   266
    @{term "op ==>"} $ _ $ _ =>
boehmes@36898
   267
      Conv.binop_conv (atomize_conv ctxt) then_conv
boehmes@36898
   268
      Conv.rewr_conv @{thm atomize_imp}
boehmes@36898
   269
  | Const (@{const_name "=="}, _) $ _ $ _ =>
boehmes@36898
   270
      Conv.binop_conv (atomize_conv ctxt) then_conv
boehmes@36898
   271
      Conv.rewr_conv @{thm atomize_eq}
boehmes@36898
   272
  | Const (@{const_name all}, _) $ Abs _ =>
wenzelm@36936
   273
      Conv.binder_conv (atomize_conv o snd) ctxt then_conv
boehmes@36898
   274
      Conv.rewr_conv @{thm atomize_all}
boehmes@36898
   275
  | _ => Conv.all_conv) ct
boehmes@36898
   276
boehmes@36898
   277
fun normalize_rule ctxt =
boehmes@36898
   278
  Conv.fconv_rule (
boehmes@36898
   279
    (* reduce lambda abstractions, except at known binders: *)
boehmes@36898
   280
    Thm.beta_conversion true then_conv
boehmes@36898
   281
    Thm.eta_conversion then_conv
boehmes@36898
   282
    norm_binder_conv ctxt) #>
boehmes@36898
   283
  norm_def ctxt #>
boehmes@36898
   284
  Drule.forall_intr_vars #>
boehmes@36898
   285
  Conv.fconv_rule (atomize_conv ctxt)
boehmes@36898
   286
boehmes@36898
   287
boehmes@36898
   288
boehmes@36898
   289
(* lift lambda terms into additional rules *)
boehmes@36898
   290
boehmes@36898
   291
local
boehmes@36898
   292
  val meta_eq = @{cpat "op =="}
boehmes@36898
   293
  val meta_eqT = hd (Thm.dest_ctyp (Thm.ctyp_of_term meta_eq))
boehmes@36898
   294
  fun inst_meta cT = Thm.instantiate_cterm ([(meta_eqT, cT)], []) meta_eq
boehmes@36898
   295
  fun mk_meta_eq ct cu = Thm.mk_binop (inst_meta (Thm.ctyp_of_term ct)) ct cu
boehmes@36898
   296
boehmes@36898
   297
  fun cert ctxt = Thm.cterm_of (ProofContext.theory_of ctxt)
boehmes@36898
   298
boehmes@36898
   299
  fun used_vars cvs ct =
boehmes@36898
   300
    let
boehmes@36898
   301
      val lookup = AList.lookup (op aconv) (map (` Thm.term_of) cvs)
boehmes@36898
   302
      val add = (fn SOME ct => insert (op aconvc) ct | _ => I)
boehmes@36898
   303
    in Term.fold_aterms (add o lookup) (Thm.term_of ct) [] end
boehmes@36898
   304
boehmes@36898
   305
  fun apply cv thm = 
boehmes@36898
   306
    let val thm' = Thm.combination thm (Thm.reflexive cv)
boehmes@36898
   307
    in Thm.transitive thm' (Thm.beta_conversion false (Thm.rhs_of thm')) end
boehmes@36898
   308
  fun apply_def cvs eq = Thm.symmetric (fold apply cvs eq)
boehmes@36898
   309
boehmes@36898
   310
  fun replace_lambda cvs ct (cx as (ctxt, defs)) =
boehmes@36898
   311
    let
boehmes@36898
   312
      val cvs' = used_vars cvs ct
boehmes@36898
   313
      val ct' = fold_rev Thm.cabs cvs' ct
boehmes@36898
   314
    in
boehmes@36898
   315
      (case Termtab.lookup defs (Thm.term_of ct') of
boehmes@36898
   316
        SOME eq => (apply_def cvs' eq, cx)
boehmes@36898
   317
      | NONE =>
boehmes@36898
   318
          let
boehmes@36898
   319
            val {T, ...} = Thm.rep_cterm ct' and n = Name.uu
boehmes@36898
   320
            val (n', ctxt') = yield_singleton Variable.variant_fixes n ctxt
boehmes@36898
   321
            val cu = mk_meta_eq (cert ctxt (Free (n', T))) ct'
boehmes@36898
   322
            val (eq, ctxt'') = yield_singleton Assumption.add_assumes cu ctxt'
boehmes@36898
   323
            val defs' = Termtab.update (Thm.term_of ct', eq) defs
boehmes@36898
   324
          in (apply_def cvs' eq, (ctxt'', defs')) end)
boehmes@36898
   325
    end
boehmes@36898
   326
boehmes@36898
   327
  fun none ct cx = (Thm.reflexive ct, cx)
boehmes@36898
   328
  fun in_comb f g ct cx =
boehmes@36898
   329
    let val (cu1, cu2) = Thm.dest_comb ct
boehmes@36898
   330
    in cx |> f cu1 ||>> g cu2 |>> uncurry Thm.combination end
boehmes@36898
   331
  fun in_arg f = in_comb none f
boehmes@36898
   332
  fun in_abs f cvs ct (ctxt, defs) =
boehmes@36898
   333
    let
boehmes@36898
   334
      val (n, ctxt') = yield_singleton Variable.variant_fixes Name.uu ctxt
boehmes@36898
   335
      val (cv, cu) = Thm.dest_abs (SOME n) ct
boehmes@36898
   336
    in  (ctxt', defs) |> f (cv :: cvs) cu |>> Thm.abstract_rule n cv end
boehmes@36898
   337
boehmes@36898
   338
  fun traverse cvs ct =
boehmes@36898
   339
    (case Thm.term_of ct of
boehmes@36898
   340
      Const (@{const_name All}, _) $ Abs _ => in_arg (in_abs traverse cvs)
boehmes@36898
   341
    | Const (@{const_name Ex}, _) $ Abs _ => in_arg (in_abs traverse cvs)
boehmes@36898
   342
    | Const (@{const_name Let}, _) $ _ $ Abs _ =>
boehmes@36898
   343
        in_comb (in_arg (traverse cvs)) (in_abs traverse cvs)
boehmes@36898
   344
    | Abs _ => at_lambda cvs
boehmes@36898
   345
    | _ $ _ => in_comb (traverse cvs) (traverse cvs)
boehmes@36898
   346
    | _ => none) ct
boehmes@36898
   347
boehmes@36898
   348
  and at_lambda cvs ct =
boehmes@36898
   349
    in_abs traverse cvs ct #-> (fn thm =>
boehmes@36898
   350
    replace_lambda cvs (Thm.rhs_of thm) #>> Thm.transitive thm)
boehmes@36898
   351
boehmes@36898
   352
  fun has_free_lambdas t =
boehmes@36898
   353
    (case t of
boehmes@36898
   354
      Const (@{const_name All}, _) $ Abs (_, _, u) => has_free_lambdas u
boehmes@36898
   355
    | Const (@{const_name Ex}, _) $ Abs (_, _, u) => has_free_lambdas u
boehmes@36898
   356
    | Const (@{const_name Let}, _) $ u1 $ Abs (_, _, u2) =>
boehmes@36898
   357
        has_free_lambdas u1 orelse has_free_lambdas u2
boehmes@36898
   358
    | Abs _ => true
boehmes@36898
   359
    | u1 $ u2 => has_free_lambdas u1 orelse has_free_lambdas u2
boehmes@36898
   360
    | _ => false)
boehmes@36898
   361
boehmes@36898
   362
  fun lift_lm f thm cx =
boehmes@36898
   363
    if not (has_free_lambdas (Thm.prop_of thm)) then (thm, cx)
boehmes@36898
   364
    else cx |> f (Thm.cprop_of thm) |>> (fn thm' => Thm.equal_elim thm' thm)
boehmes@36898
   365
in
boehmes@40161
   366
fun lift_lambdas irules ctxt =
boehmes@36898
   367
  let
boehmes@36898
   368
    val cx = (ctxt, Termtab.empty)
boehmes@40161
   369
    val (idxs, thms) = split_list irules
boehmes@36898
   370
    val (thms', (ctxt', defs)) = fold_map (lift_lm (traverse [])) thms cx
boehmes@36898
   371
    val eqs = Termtab.fold (cons o normalize_rule ctxt' o snd) defs []
boehmes@40161
   372
  in (map (pair ~1) eqs @ (idxs ~~ thms'), ctxt') end
boehmes@36898
   373
end
boehmes@36898
   374
boehmes@36898
   375
boehmes@36898
   376
boehmes@36898
   377
(* make application explicit for functions with varying number of arguments *)
boehmes@36898
   378
boehmes@36898
   379
local
boehmes@36898
   380
  val const = prefix "c" and free = prefix "f"
boehmes@36898
   381
  fun min i (e as (_, j)) = if i <> j then (true, Int.min (i, j)) else e
boehmes@36898
   382
  fun add t i = Symtab.map_default (t, (false, i)) (min i)
boehmes@36898
   383
  fun traverse t =
boehmes@36898
   384
    (case Term.strip_comb t of
boehmes@36898
   385
      (Const (n, _), ts) => add (const n) (length ts) #> fold traverse ts 
boehmes@36898
   386
    | (Free (n, _), ts) => add (free n) (length ts) #> fold traverse ts
boehmes@36898
   387
    | (Abs (_, _, u), ts) => fold traverse (u :: ts)
boehmes@36898
   388
    | (_, ts) => fold traverse ts)
boehmes@40161
   389
  fun prune tab = Symtab.fold (fn (n, (true, i)) =>
boehmes@40161
   390
    Symtab.update (n, i) | _ => I) tab Symtab.empty
boehmes@36898
   391
boehmes@36898
   392
  fun binop_conv cv1 cv2 = Conv.combination_conv (Conv.arg_conv cv1) cv2
boehmes@36898
   393
  fun nary_conv conv1 conv2 ct =
boehmes@36898
   394
    (Conv.combination_conv (nary_conv conv1 conv2) conv2 else_conv conv1) ct
boehmes@36898
   395
  fun abs_conv conv tb = Conv.abs_conv (fn (cv, cx) =>
boehmes@36898
   396
    let val n = fst (Term.dest_Free (Thm.term_of cv))
boehmes@36898
   397
    in conv (Symtab.update (free n, 0) tb) cx end)
boehmes@37153
   398
  val fun_app_rule = @{lemma "f x == fun_app f x" by (simp add: fun_app_def)}
boehmes@36898
   399
in
boehmes@40161
   400
fun explicit_application ctxt irules =
boehmes@36898
   401
  let
boehmes@36898
   402
    fun sub_conv tb ctxt ct =
boehmes@36898
   403
      (case Term.strip_comb (Thm.term_of ct) of
boehmes@36898
   404
        (Const (n, _), ts) => app_conv tb (const n) (length ts) ctxt
boehmes@36898
   405
      | (Free (n, _), ts) => app_conv tb (free n) (length ts) ctxt
boehmes@36898
   406
      | (Abs _, _) => nary_conv (abs_conv sub_conv tb ctxt) (sub_conv tb ctxt)
boehmes@36898
   407
      | (_, _) => nary_conv Conv.all_conv (sub_conv tb ctxt)) ct
boehmes@36898
   408
    and app_conv tb n i ctxt =
boehmes@36898
   409
      (case Symtab.lookup tb n of
boehmes@36898
   410
        NONE => nary_conv Conv.all_conv (sub_conv tb ctxt)
boehmes@37153
   411
      | SOME j => fun_app_conv tb ctxt (i - j))
boehmes@37153
   412
    and fun_app_conv tb ctxt i ct = (
boehmes@36898
   413
      if i = 0 then nary_conv Conv.all_conv (sub_conv tb ctxt)
boehmes@36898
   414
      else
boehmes@37153
   415
        Conv.rewr_conv fun_app_rule then_conv
boehmes@37153
   416
        binop_conv (fun_app_conv tb ctxt (i-1)) (sub_conv tb ctxt)) ct
boehmes@36898
   417
boehmes@36898
   418
    fun needs_exp_app tab = Term.exists_subterm (fn
boehmes@36898
   419
        Bound _ $ _ => true
boehmes@36898
   420
      | Const (n, _) => Symtab.defined tab (const n)
boehmes@36898
   421
      | Free (n, _) => Symtab.defined tab (free n)
boehmes@36898
   422
      | _ => false)
boehmes@36898
   423
boehmes@36898
   424
    fun rewrite tab ctxt thm =
boehmes@36898
   425
      if not (needs_exp_app tab (Thm.prop_of thm)) then thm
boehmes@36898
   426
      else Conv.fconv_rule (sub_conv tab ctxt) thm
boehmes@36898
   427
boehmes@40161
   428
    val tab = prune (fold (traverse o Thm.prop_of o snd) irules Symtab.empty)
boehmes@40161
   429
  in map (apsnd (rewrite tab ctxt)) irules end
boehmes@36898
   430
end
boehmes@36898
   431
boehmes@36898
   432
boehmes@36898
   433
boehmes@39483
   434
(* add missing datatype selectors via hypothetical definitions *)
boehmes@39483
   435
boehmes@39483
   436
local
boehmes@39483
   437
  val add = (fn Type (n, _) => Symtab.update (n, ()) | _ => I)
boehmes@39483
   438
boehmes@39483
   439
  fun collect t =
boehmes@39483
   440
    (case Term.strip_comb t of
boehmes@39483
   441
      (Abs (_, T, t), _) => add T #> collect t
boehmes@39483
   442
    | (Const (_, T), ts) => collects T ts
boehmes@39483
   443
    | (Free (_, T), ts) => collects T ts
boehmes@39483
   444
    | _ => I)
boehmes@39483
   445
  and collects T ts =
boehmes@39483
   446
    let val ((Ts, Us), U) = Term.strip_type T |> apfst (chop (length ts))
boehmes@39483
   447
    in fold add Ts #> add (Us ---> U) #> fold collect ts end
boehmes@39483
   448
boehmes@39483
   449
  fun add_constructors thy n =
boehmes@39483
   450
    (case Datatype.get_info thy n of
boehmes@39483
   451
      NONE => I
boehmes@39483
   452
    | SOME {descr, ...} => fold (fn (_, (_, _, cs)) => fold (fn (n, ds) =>
boehmes@39483
   453
        fold (insert (op =) o pair n) (1 upto length ds)) cs) descr)
boehmes@39483
   454
boehmes@39483
   455
  fun add_selector (c as (n, i)) ctxt =
boehmes@39483
   456
    (case Datatype_Selectors.lookup_selector ctxt c of
boehmes@39483
   457
      SOME _ => ctxt
boehmes@39483
   458
    | NONE =>
boehmes@39483
   459
        let
boehmes@39483
   460
          val T = Sign.the_const_type (ProofContext.theory_of ctxt) n
boehmes@39483
   461
          val U = Term.body_type T --> nth (Term.binder_types T) (i-1)
boehmes@39483
   462
        in
boehmes@39483
   463
          ctxt
boehmes@39483
   464
          |> yield_singleton Variable.variant_fixes Name.uu
boehmes@39483
   465
          |>> pair ((n, T), i) o rpair U
boehmes@39483
   466
          |-> Context.proof_map o Datatype_Selectors.add_selector
boehmes@39483
   467
        end)
boehmes@39483
   468
in
boehmes@39483
   469
boehmes@40161
   470
fun datatype_selectors irules ctxt =
boehmes@39483
   471
  let
boehmes@40161
   472
    val ns = Symtab.keys (fold (collect o Thm.prop_of o snd) irules Symtab.empty)
boehmes@39483
   473
    val cs = fold (add_constructors (ProofContext.theory_of ctxt)) ns []
boehmes@40161
   474
  in (irules, fold add_selector cs ctxt) end
boehmes@39483
   475
    (* FIXME: also generate hypothetical definitions for the selectors *)
boehmes@39483
   476
boehmes@39483
   477
end
boehmes@39483
   478
boehmes@39483
   479
boehmes@39483
   480
boehmes@36898
   481
(* combined normalization *)
boehmes@36898
   482
boehmes@40161
   483
type extra_norm = (int * thm) list -> Proof.context ->
boehmes@40161
   484
  (int * thm) list * Proof.context
boehmes@36898
   485
boehmes@40161
   486
fun with_context f irules ctxt = (f ctxt irules, ctxt)
boehmes@36898
   487
boehmes@40161
   488
fun normalize extra_norm with_datatypes irules ctxt =
boehmes@40161
   489
  irules
boehmes@36898
   490
  |> trivial_distinct ctxt
boehmes@36898
   491
  |> rewrite_bool_cases ctxt
boehmes@36898
   492
  |> normalize_numerals ctxt
boehmes@36898
   493
  |> nat_as_int ctxt
boehmes@36898
   494
  |> rpair ctxt
boehmes@36898
   495
  |-> extra_norm
boehmes@40161
   496
  |-> with_context (fn cx => map (apsnd (normalize_rule cx)))
boehmes@36898
   497
  |-> SMT_Monomorph.monomorph
boehmes@36898
   498
  |-> lift_lambdas
boehmes@36898
   499
  |-> with_context explicit_application
boehmes@39483
   500
  |-> (if with_datatypes then datatype_selectors else pair)
boehmes@36898
   501
boehmes@36898
   502
end