src/HOL/Tools/datatype_package/datatype_rep_proofs.ML
author haftmann
Sun Jun 21 08:38:58 2009 +0200 (2009-06-21)
changeset 31737 b3f63611784e
parent 31723 f5cafe803b55
permissions -rw-r--r--
simplified names of common datatype types
berghofe@5177
     1
(*  Title:      HOL/Tools/datatype_rep_proofs.ML
wenzelm@11539
     2
    Author:     Stefan Berghofer, TU Muenchen
berghofe@5177
     3
berghofe@5177
     4
Definitional introduction of datatypes
berghofe@5177
     5
Proof of characteristic theorems:
berghofe@5177
     6
berghofe@5177
     7
 - injectivity of constructors
berghofe@7228
     8
 - distinctness of constructors
berghofe@5177
     9
 - induction theorem
berghofe@5177
    10
*)
berghofe@5177
    11
berghofe@5177
    12
signature DATATYPE_REP_PROOFS =
berghofe@5177
    13
sig
haftmann@31737
    14
  include DATATYPE_COMMON
haftmann@27002
    15
  val distinctness_limit : int Config.T
haftmann@27002
    16
  val distinctness_limit_setup : theory -> theory
haftmann@31737
    17
  val representation_proofs : config -> info Symtab.table ->
haftmann@31668
    18
    string list -> descr list -> (string * sort) list ->
wenzelm@30345
    19
      (binding * mixfix) list -> (binding * mixfix) list list -> attribute
haftmann@18314
    20
        -> theory -> (thm list list * thm list list * thm list list *
haftmann@18314
    21
          DatatypeAux.simproc_dist list * thm) * theory
berghofe@5177
    22
end;
berghofe@5177
    23
berghofe@5177
    24
structure DatatypeRepProofs : DATATYPE_REP_PROOFS =
berghofe@5177
    25
struct
berghofe@5177
    26
berghofe@5177
    27
open DatatypeAux;
berghofe@5177
    28
haftmann@27002
    29
(*the kind of distinctiveness axioms depends on number of constructors*)
haftmann@27002
    30
val (distinctness_limit, distinctness_limit_setup) =
haftmann@27002
    31
  Attrib.config_int "datatype_distinctness_limit" 7;
haftmann@27002
    32
berghofe@5177
    33
val (_ $ (_ $ (_ $ (distinct_f $ _) $ _))) = hd (prems_of distinct_lemma);
berghofe@5177
    34
berghofe@21021
    35
val collect_simp = rewrite_rule [mk_meta_eq mem_Collect_eq];
berghofe@21021
    36
wenzelm@11435
    37
wenzelm@11435
    38
(** theory context references **)
berghofe@5177
    39
wenzelm@11435
    40
val f_myinv_f = thm "f_myinv_f";
wenzelm@11435
    41
val myinv_f_f = thm "myinv_f_f";
wenzelm@11435
    42
berghofe@5177
    43
haftmann@31737
    44
fun exh_thm_of (dt_info : info Symtab.table) tname =
wenzelm@17412
    45
  #exhaustion (the (Symtab.lookup dt_info tname));
berghofe@5177
    46
berghofe@5177
    47
(******************************************************************************)
berghofe@5177
    48
haftmann@31737
    49
fun representation_proofs (config : config) (dt_info : info Symtab.table)
wenzelm@8436
    50
      new_type_names descr sorts types_syntax constr_syntax case_names_induct thy =
berghofe@5177
    51
  let
wenzelm@19540
    52
    val Datatype_thy = ThyInfo.the_theory "Datatype" thy;
wenzelm@20820
    53
    val node_name = "Datatype.node";
wenzelm@20820
    54
    val In0_name = "Datatype.In0";
wenzelm@20820
    55
    val In1_name = "Datatype.In1";
wenzelm@20820
    56
    val Scons_name = "Datatype.Scons";
wenzelm@20820
    57
    val Leaf_name = "Datatype.Leaf";
wenzelm@20820
    58
    val Numb_name = "Datatype.Numb";
wenzelm@20820
    59
    val Lim_name = "Datatype.Lim";
berghofe@13641
    60
    val Suml_name = "Datatype.Suml";
berghofe@13641
    61
    val Sumr_name = "Datatype.Sumr";
berghofe@7015
    62
berghofe@13641
    63
    val [In0_inject, In1_inject, Scons_inject, Leaf_inject,
berghofe@13641
    64
         In0_eq, In1_eq, In0_not_In1, In1_not_In0,
wenzelm@26343
    65
         Lim_inject, Suml_inject, Sumr_inject] = map (PureThy.get_thm Datatype_thy)
wenzelm@26343
    66
          ["In0_inject", "In1_inject", "Scons_inject", "Leaf_inject",
wenzelm@26343
    67
           "In0_eq", "In1_eq", "In0_not_In1", "In1_not_In0",
wenzelm@26343
    68
           "Lim_inject", "Suml_inject", "Sumr_inject"];
berghofe@5177
    69
wenzelm@28362
    70
    val descr' = flat descr;
berghofe@5177
    71
berghofe@5661
    72
    val big_name = space_implode "_" new_type_names;
haftmann@31668
    73
    val thy1 = add_path (#flat_names config) big_name thy;
berghofe@5661
    74
    val big_rec_name = big_name ^ "_rep_set";
berghofe@21021
    75
    val rep_set_names' =
berghofe@5177
    76
      (if length descr' = 1 then [big_rec_name] else
berghofe@5177
    77
        (map ((curry (op ^) (big_rec_name ^ "_")) o string_of_int)
berghofe@5177
    78
          (1 upto (length descr'))));
haftmann@28965
    79
    val rep_set_names = map (Sign.full_bname thy1) rep_set_names';
berghofe@5177
    80
berghofe@5661
    81
    val tyvars = map (fn (_, (_, Ts, _)) => map dest_DtTFree Ts) (hd descr);
berghofe@5661
    82
    val leafTs' = get_nonrec_types descr' sorts;
berghofe@7015
    83
    val branchTs = get_branching_types descr' sorts;
berghofe@7015
    84
    val branchT = if null branchTs then HOLogic.unitT
wenzelm@23419
    85
      else BalancedTree.make (fn (T, U) => Type ("+", [T, U])) branchTs;
berghofe@13641
    86
    val arities = get_arities descr' \ 0;
wenzelm@30190
    87
    val unneeded_vars = hd tyvars \\ List.foldr OldTerm.add_typ_tfree_names [] (leafTs' @ branchTs);
haftmann@17521
    88
    val leafTs = leafTs' @ (map (fn n => TFree (n, (the o AList.lookup (op =) sorts) n)) unneeded_vars);
berghofe@5177
    89
    val recTs = get_rec_types descr' sorts;
skalberg@15570
    90
    val newTs = Library.take (length (hd descr), recTs);
skalberg@15570
    91
    val oldTs = Library.drop (length (hd descr), recTs);
berghofe@5177
    92
    val sumT = if null leafTs then HOLogic.unitT
wenzelm@23419
    93
      else BalancedTree.make (fn (T, U) => Type ("+", [T, U])) leafTs;
berghofe@7015
    94
    val Univ_elT = HOLogic.mk_setT (Type (node_name, [sumT, branchT]));
berghofe@5177
    95
    val UnivT = HOLogic.mk_setT Univ_elT;
berghofe@21021
    96
    val UnivT' = Univ_elT --> HOLogic.boolT;
berghofe@21021
    97
    val Collect = Const ("Collect", UnivT' --> UnivT);
berghofe@5177
    98
berghofe@5177
    99
    val In0 = Const (In0_name, Univ_elT --> Univ_elT);
berghofe@5177
   100
    val In1 = Const (In1_name, Univ_elT --> Univ_elT);
berghofe@5177
   101
    val Leaf = Const (Leaf_name, sumT --> Univ_elT);
berghofe@7015
   102
    val Lim = Const (Lim_name, (branchT --> Univ_elT) --> Univ_elT);
berghofe@5177
   103
berghofe@5177
   104
    (* make injections needed for embedding types in leaves *)
berghofe@5177
   105
berghofe@5177
   106
    fun mk_inj T' x =
berghofe@5177
   107
      let
berghofe@5177
   108
        fun mk_inj' T n i =
berghofe@5177
   109
          if n = 1 then x else
berghofe@5177
   110
          let val n2 = n div 2;
berghofe@5177
   111
              val Type (_, [T1, T2]) = T
berghofe@5177
   112
          in
berghofe@5177
   113
            if i <= n2 then
paulson@15391
   114
              Const ("Sum_Type.Inl", T1 --> T) $ (mk_inj' T1 n2 i)
berghofe@5177
   115
            else
paulson@15391
   116
              Const ("Sum_Type.Inr", T2 --> T) $ (mk_inj' T2 (n - n2) (i - n2))
berghofe@5177
   117
          end
berghofe@5177
   118
      in mk_inj' sumT (length leafTs) (1 + find_index_eq T' leafTs)
berghofe@5177
   119
      end;
berghofe@5177
   120
berghofe@5177
   121
    (* make injections for constructors *)
berghofe@5177
   122
wenzelm@23419
   123
    fun mk_univ_inj ts = BalancedTree.access
wenzelm@23419
   124
      {left = fn t => In0 $ t,
wenzelm@23419
   125
        right = fn t => In1 $ t,
wenzelm@23419
   126
        init =
haftmann@28524
   127
          if ts = [] then Const (@{const_name undefined}, Univ_elT)
wenzelm@23419
   128
          else foldr1 (HOLogic.mk_binop Scons_name) ts};
berghofe@5177
   129
berghofe@7015
   130
    (* function spaces *)
berghofe@7015
   131
berghofe@7015
   132
    fun mk_fun_inj T' x =
berghofe@7015
   133
      let
berghofe@7015
   134
        fun mk_inj T n i =
berghofe@7015
   135
          if n = 1 then x else
berghofe@7015
   136
          let
berghofe@7015
   137
            val n2 = n div 2;
berghofe@7015
   138
            val Type (_, [T1, T2]) = T;
berghofe@13641
   139
            fun mkT U = (U --> Univ_elT) --> T --> Univ_elT
berghofe@7015
   140
          in
berghofe@13641
   141
            if i <= n2 then Const (Suml_name, mkT T1) $ mk_inj T1 n2 i
berghofe@13641
   142
            else Const (Sumr_name, mkT T2) $ mk_inj T2 (n - n2) (i - n2)
berghofe@7015
   143
          end
berghofe@7015
   144
      in mk_inj branchT (length branchTs) (1 + find_index_eq T' branchTs)
berghofe@7015
   145
      end;
berghofe@7015
   146
wenzelm@30190
   147
    val mk_lim = List.foldr (fn (T, t) => Lim $ mk_fun_inj T (Abs ("x", T, t)));
berghofe@13641
   148
berghofe@5177
   149
    (************** generate introduction rules for representing set **********)
berghofe@5177
   150
haftmann@31668
   151
    val _ = message config "Constructing representing sets ...";
berghofe@5177
   152
berghofe@5177
   153
    (* make introduction rule for a single constructor *)
berghofe@5177
   154
berghofe@5177
   155
    fun make_intr s n (i, (_, cargs)) =
berghofe@5177
   156
      let
berghofe@13641
   157
        fun mk_prem (dt, (j, prems, ts)) = (case strip_dtyp dt of
berghofe@13641
   158
            (dts, DtRec k) =>
berghofe@13641
   159
              let
berghofe@13641
   160
                val Ts = map (typ_of_dtyp descr' sorts) dts;
berghofe@13641
   161
                val free_t =
berghofe@13641
   162
                  app_bnds (mk_Free "x" (Ts ---> Univ_elT) j) (length Ts)
berghofe@13641
   163
              in (j + 1, list_all (map (pair "x") Ts,
berghofe@21021
   164
                  HOLogic.mk_Trueprop
berghofe@21021
   165
                    (Free (List.nth (rep_set_names', k), UnivT') $ free_t)) :: prems,
skalberg@15574
   166
                mk_lim free_t Ts :: ts)
berghofe@5177
   167
              end
berghofe@13641
   168
          | _ =>
berghofe@5177
   169
              let val T = typ_of_dtyp descr' sorts dt
berghofe@5177
   170
              in (j + 1, prems, (Leaf $ mk_inj T (mk_Free "x" T j))::ts)
berghofe@13641
   171
              end);
berghofe@5177
   172
wenzelm@30190
   173
        val (_, prems, ts) = List.foldr mk_prem (1, [], []) cargs;
berghofe@21021
   174
        val concl = HOLogic.mk_Trueprop
berghofe@21021
   175
          (Free (s, UnivT') $ mk_univ_inj ts n i)
berghofe@13641
   176
      in Logic.list_implies (prems, concl)
berghofe@5177
   177
      end;
berghofe@5177
   178
wenzelm@28362
   179
    val intr_ts = maps (fn ((_, (_, _, constrs)), rep_set_name) =>
berghofe@5177
   180
      map (make_intr rep_set_name (length constrs))
wenzelm@28362
   181
        ((1 upto (length constrs)) ~~ constrs)) (descr' ~~ rep_set_names');
berghofe@5177
   182
wenzelm@21365
   183
    val ({raw_induct = rep_induct, intrs = rep_intrs, ...}, thy2) =
haftmann@31723
   184
        Inductive.add_inductive_global (serial_string ())
haftmann@31668
   185
          {quiet_mode = #quiet config, verbose = false, kind = Thm.internalK,
haftmann@28965
   186
           alt_name = Binding.name big_rec_name, coind = false, no_elim = true, no_ind = false,
wenzelm@29389
   187
           skip_mono = true, fork_mono = false}
haftmann@28965
   188
          (map (fn s => ((Binding.name s, UnivT'), NoSyn)) rep_set_names') []
haftmann@28965
   189
          (map (fn x => (Attrib.empty_binding, x)) intr_ts) [] thy1;
berghofe@5177
   190
berghofe@5177
   191
    (********************************* typedef ********************************)
berghofe@5177
   192
berghofe@21021
   193
    val (typedefs, thy3) = thy2 |>
haftmann@31668
   194
      parent_path (#flat_names config) |>
berghofe@21021
   195
      fold_map (fn ((((name, mx), tvs), c), name') =>
haftmann@31723
   196
          Typedef.add_typedef false (SOME (Binding.name name')) (name, tvs, mx)
berghofe@21021
   197
            (Collect $ Const (c, UnivT')) NONE
berghofe@21021
   198
            (rtac exI 1 THEN rtac CollectI 1 THEN
berghofe@21021
   199
              QUIET_BREADTH_FIRST (has_fewer_prems 1)
wenzelm@26475
   200
              (resolve_tac rep_intrs 1)))
berghofe@21021
   201
                (types_syntax ~~ tyvars ~~
berghofe@21021
   202
                  (Library.take (length newTs, rep_set_names)) ~~ new_type_names) ||>
haftmann@31668
   203
      add_path (#flat_names config) big_name;
berghofe@5177
   204
berghofe@5177
   205
    (*********************** definition of constructors ***********************)
berghofe@5177
   206
berghofe@5177
   207
    val big_rep_name = (space_implode "_" new_type_names) ^ "_Rep_";
berghofe@5177
   208
    val rep_names = map (curry op ^ "Rep_") new_type_names;
berghofe@5177
   209
    val rep_names' = map (fn i => big_rep_name ^ (string_of_int i))
wenzelm@28362
   210
      (1 upto (length (flat (tl descr))));
wenzelm@22578
   211
    val all_rep_names = map (Sign.intern_const thy3) rep_names @
haftmann@28965
   212
      map (Sign.full_bname thy3) rep_names';
berghofe@5177
   213
berghofe@5177
   214
    (* isomorphism declarations *)
berghofe@5177
   215
wenzelm@30345
   216
    val iso_decls = map (fn (T, s) => (Binding.name s, T --> Univ_elT, NoSyn))
berghofe@5177
   217
      (oldTs ~~ rep_names');
berghofe@5177
   218
berghofe@5177
   219
    (* constructor definitions *)
berghofe@5177
   220
berghofe@5177
   221
    fun make_constr_def tname T n ((thy, defs, eqns, i), ((cname, cargs), (cname', mx))) =
berghofe@5177
   222
      let
berghofe@5177
   223
        fun constr_arg (dt, (j, l_args, r_args)) =
berghofe@5177
   224
          let val T = typ_of_dtyp descr' sorts dt;
berghofe@5177
   225
              val free_t = mk_Free "x" T j
berghofe@13641
   226
          in (case (strip_dtyp dt, strip_type T) of
skalberg@15574
   227
              ((_, DtRec m), (Us, U)) => (j + 1, free_t :: l_args, mk_lim
skalberg@15574
   228
                (Const (List.nth (all_rep_names, m), U --> Univ_elT) $
skalberg@15574
   229
                   app_bnds free_t (length Us)) Us :: r_args)
berghofe@5177
   230
            | _ => (j + 1, free_t::l_args, (Leaf $ mk_inj T free_t)::r_args))
berghofe@5177
   231
          end;
berghofe@5177
   232
wenzelm@30190
   233
        val (_, l_args, r_args) = List.foldr constr_arg (1, [], []) cargs;
berghofe@5177
   234
        val constrT = (map (typ_of_dtyp descr' sorts) cargs) ---> T;
wenzelm@22578
   235
        val abs_name = Sign.intern_const thy ("Abs_" ^ tname);
wenzelm@22578
   236
        val rep_name = Sign.intern_const thy ("Rep_" ^ tname);
berghofe@5177
   237
        val lhs = list_comb (Const (cname, constrT), l_args);
berghofe@5177
   238
        val rhs = mk_univ_inj r_args n i;
wenzelm@27330
   239
        val def = Logic.mk_equals (lhs, Const (abs_name, Univ_elT --> T) $ rhs);
wenzelm@30364
   240
        val def_name = Long_Name.base_name cname ^ "_def";
berghofe@5177
   241
        val eqn = HOLogic.mk_Trueprop (HOLogic.mk_eq
berghofe@5177
   242
          (Const (rep_name, T --> Univ_elT) $ lhs, rhs));
haftmann@18377
   243
        val ([def_thm], thy') =
haftmann@18377
   244
          thy
wenzelm@24712
   245
          |> Sign.add_consts_i [(cname', constrT, mx)]
haftmann@29579
   246
          |> (PureThy.add_defs false o map Thm.no_attributes) [(Binding.name def_name, def)];
berghofe@5177
   247
wenzelm@8436
   248
      in (thy', defs @ [def_thm], eqns @ [eqn], i + 1) end;
berghofe@5177
   249
berghofe@5177
   250
    (* constructor definitions for datatype *)
berghofe@5177
   251
berghofe@5177
   252
    fun dt_constr_defs ((thy, defs, eqns, rep_congs, dist_lemmas),
berghofe@5177
   253
        ((((_, (_, _, constrs)), tname), T), constr_syntax)) =
berghofe@5177
   254
      let
berghofe@5177
   255
        val _ $ (_ $ (cong_f $ _) $ _) = concl_of arg_cong;
wenzelm@22578
   256
        val rep_const = cterm_of thy
wenzelm@22578
   257
          (Const (Sign.intern_const thy ("Rep_" ^ tname), T --> Univ_elT));
wenzelm@22578
   258
        val cong' = standard (cterm_instantiate [(cterm_of thy cong_f, rep_const)] arg_cong);
wenzelm@22578
   259
        val dist = standard (cterm_instantiate [(cterm_of thy distinct_f, rep_const)] distinct_lemma);
skalberg@15570
   260
        val (thy', defs', eqns', _) = Library.foldl ((make_constr_def tname T) (length constrs))
haftmann@31668
   261
          ((add_path (#flat_names config) tname thy, defs, [], 1), constrs ~~ constr_syntax)
berghofe@5177
   262
      in
haftmann@31668
   263
        (parent_path (#flat_names config) thy', defs', eqns @ [eqns'],
berghofe@5177
   264
          rep_congs @ [cong'], dist_lemmas @ [dist])
berghofe@5177
   265
      end;
berghofe@5177
   266
skalberg@15570
   267
    val (thy4, constr_defs, constr_rep_eqns, rep_congs, dist_lemmas) = Library.foldl dt_constr_defs
haftmann@31668
   268
      ((thy3 |> Sign.add_consts_i iso_decls |> parent_path (#flat_names config), [], [], [], []),
berghofe@5177
   269
        hd descr ~~ new_type_names ~~ newTs ~~ constr_syntax);
berghofe@5177
   270
berghofe@5177
   271
    (*********** isomorphisms for new types (introduced by typedef) ***********)
berghofe@5177
   272
haftmann@31668
   273
    val _ = message config "Proving isomorphism properties ...";
berghofe@5177
   274
berghofe@21021
   275
    val newT_iso_axms = map (fn (_, td) =>
berghofe@21021
   276
      (collect_simp (#Abs_inverse td), #Rep_inverse td,
berghofe@21021
   277
       collect_simp (#Rep td))) typedefs;
berghofe@5177
   278
berghofe@21021
   279
    val newT_iso_inj_thms = map (fn (_, td) =>
berghofe@21021
   280
      (collect_simp (#Abs_inject td) RS iffD1, #Rep_inject td RS iffD1)) typedefs;
berghofe@5177
   281
berghofe@5177
   282
    (********* isomorphisms between existing types and "unfolded" types *******)
berghofe@5177
   283
berghofe@5177
   284
    (*---------------------------------------------------------------------*)
berghofe@5177
   285
    (* isomorphisms are defined using primrec-combinators:                 *)
berghofe@5177
   286
    (* generate appropriate functions for instantiating primrec-combinator *)
berghofe@5177
   287
    (*                                                                     *)
berghofe@13641
   288
    (*   e.g.  dt_Rep_i = list_rec ... (%h t y. In1 (Scons (Leaf h) y))    *)
berghofe@5177
   289
    (*                                                                     *)
berghofe@5177
   290
    (* also generate characteristic equations for isomorphisms             *)
berghofe@5177
   291
    (*                                                                     *)
berghofe@13641
   292
    (*   e.g.  dt_Rep_i (cons h t) = In1 (Scons (dt_Rep_j h) (dt_Rep_i t)) *)
berghofe@5177
   293
    (*---------------------------------------------------------------------*)
berghofe@5177
   294
berghofe@5177
   295
    fun make_iso_def k ks n ((fs, eqns, i), (cname, cargs)) =
berghofe@5177
   296
      let
berghofe@5177
   297
        val argTs = map (typ_of_dtyp descr' sorts) cargs;
skalberg@15570
   298
        val T = List.nth (recTs, k);
skalberg@15570
   299
        val rep_name = List.nth (all_rep_names, k);
berghofe@5177
   300
        val rep_const = Const (rep_name, T --> Univ_elT);
berghofe@5177
   301
        val constr = Const (cname, argTs ---> T);
berghofe@5177
   302
berghofe@7015
   303
        fun process_arg ks' ((i2, i2', ts, Ts), dt) =
berghofe@13641
   304
          let
berghofe@13641
   305
            val T' = typ_of_dtyp descr' sorts dt;
berghofe@13641
   306
            val (Us, U) = strip_type T'
berghofe@13641
   307
          in (case strip_dtyp dt of
berghofe@13641
   308
              (_, DtRec j) => if j mem ks' then
skalberg@15574
   309
                  (i2 + 1, i2' + 1, ts @ [mk_lim (app_bnds
skalberg@15574
   310
                     (mk_Free "y" (Us ---> Univ_elT) i2') (length Us)) Us],
berghofe@13641
   311
                   Ts @ [Us ---> Univ_elT])
berghofe@5177
   312
                else
skalberg@15574
   313
                  (i2 + 1, i2', ts @ [mk_lim
skalberg@15574
   314
                     (Const (List.nth (all_rep_names, j), U --> Univ_elT) $
skalberg@15574
   315
                        app_bnds (mk_Free "x" T' i2) (length Us)) Us], Ts)
berghofe@7015
   316
            | _ => (i2 + 1, i2', ts @ [Leaf $ mk_inj T' (mk_Free "x" T' i2)], Ts))
berghofe@5177
   317
          end;
berghofe@5177
   318
skalberg@15570
   319
        val (i2, i2', ts, Ts) = Library.foldl (process_arg ks) ((1, 1, [], []), cargs);
berghofe@5177
   320
        val xs = map (uncurry (mk_Free "x")) (argTs ~~ (1 upto (i2 - 1)));
berghofe@7015
   321
        val ys = map (uncurry (mk_Free "y")) (Ts ~~ (1 upto (i2' - 1)));
berghofe@5177
   322
        val f = list_abs_free (map dest_Free (xs @ ys), mk_univ_inj ts n i);
berghofe@5177
   323
skalberg@15570
   324
        val (_, _, ts', _) = Library.foldl (process_arg []) ((1, 1, [], []), cargs);
berghofe@5177
   325
        val eqn = HOLogic.mk_Trueprop (HOLogic.mk_eq
berghofe@5177
   326
          (rep_const $ list_comb (constr, xs), mk_univ_inj ts' n i))
berghofe@5177
   327
berghofe@5177
   328
      in (fs @ [f], eqns @ [eqn], i + 1) end;
berghofe@5177
   329
berghofe@5177
   330
    (* define isomorphisms for all mutually recursive datatypes in list ds *)
berghofe@5177
   331
berghofe@5177
   332
    fun make_iso_defs (ds, (thy, char_thms)) =
berghofe@5177
   333
      let
berghofe@5177
   334
        val ks = map fst ds;
berghofe@5177
   335
        val (_, (tname, _, _)) = hd ds;
wenzelm@17412
   336
        val {rec_rewrites, rec_names, ...} = the (Symtab.lookup dt_info tname);
berghofe@5177
   337
berghofe@5177
   338
        fun process_dt ((fs, eqns, isos), (k, (tname, _, constrs))) =
berghofe@5177
   339
          let
skalberg@15570
   340
            val (fs', eqns', _) = Library.foldl (make_iso_def k ks (length constrs))
berghofe@5177
   341
              ((fs, eqns, 1), constrs);
skalberg@15570
   342
            val iso = (List.nth (recTs, k), List.nth (all_rep_names, k))
berghofe@5177
   343
          in (fs', eqns', isos @ [iso]) end;
berghofe@5177
   344
        
skalberg@15570
   345
        val (fs, eqns, isos) = Library.foldl process_dt (([], [], []), ds);
berghofe@5177
   346
        val fTs = map fastype_of fs;
wenzelm@30364
   347
        val defs = map (fn (rec_name, (T, iso_name)) => (Binding.name (Long_Name.base_name iso_name ^ "_def"),
wenzelm@27330
   348
          Logic.mk_equals (Const (iso_name, T --> Univ_elT),
wenzelm@27330
   349
            list_comb (Const (rec_name, fTs @ [T] ---> Univ_elT), fs)))) (rec_names ~~ isos);
wenzelm@28362
   350
        val (def_thms, thy') =
wenzelm@28362
   351
          apsnd Theory.checkpoint ((PureThy.add_defs false o map Thm.no_attributes) defs thy);
berghofe@5177
   352
berghofe@5177
   353
        (* prove characteristic equations *)
berghofe@5177
   354
oheimb@5553
   355
        val rewrites = def_thms @ (map mk_meta_eq rec_rewrites);
berghofe@26531
   356
        val char_thms' = map (fn eqn => SkipProof.prove_global thy' [] [] eqn
wenzelm@20046
   357
          (fn _ => EVERY [rewrite_goals_tac rewrites, rtac refl 1])) eqns;
berghofe@5177
   358
berghofe@5177
   359
      in (thy', char_thms' @ char_thms) end;
berghofe@5177
   360
wenzelm@30190
   361
    val (thy5, iso_char_thms) = apfst Theory.checkpoint (List.foldr make_iso_defs
haftmann@31668
   362
      (add_path (#flat_names config) big_name thy4, []) (tl descr));
berghofe@5177
   363
berghofe@5177
   364
    (* prove isomorphism properties *)
berghofe@5177
   365
wenzelm@28362
   366
    fun mk_funs_inv thy thm =
berghofe@7015
   367
      let
wenzelm@26626
   368
        val prop = Thm.prop_of thm;
berghofe@21021
   369
        val _ $ (_ $ ((S as Const (_, Type (_, [U, _]))) $ _ )) $
wenzelm@16287
   370
          (_ $ (_ $ (r $ (a $ _)) $ _)) = Type.freeze prop;
wenzelm@29270
   371
        val used = OldTerm.add_term_tfree_names (a, []);
berghofe@13641
   372
berghofe@13641
   373
        fun mk_thm i =
berghofe@13641
   374
          let
berghofe@13641
   375
            val Ts = map (TFree o rpair HOLogic.typeS)
wenzelm@20071
   376
              (Name.variant_list used (replicate i "'t"));
berghofe@13641
   377
            val f = Free ("f", Ts ---> U)
berghofe@26531
   378
          in SkipProof.prove_global thy [] [] (Logic.mk_implies
berghofe@13641
   379
            (HOLogic.mk_Trueprop (HOLogic.list_all
berghofe@21021
   380
               (map (pair "x") Ts, S $ app_bnds f i)),
berghofe@13641
   381
             HOLogic.mk_Trueprop (HOLogic.mk_eq (list_abs (map (pair "x") Ts,
wenzelm@17985
   382
               r $ (a $ app_bnds f i)), f))))
berghofe@26806
   383
            (fn _ => EVERY [REPEAT_DETERM_N i (rtac ext 1),
berghofe@26806
   384
               REPEAT (etac allE 1), rtac thm 1, atac 1])
berghofe@13641
   385
          end
berghofe@13641
   386
      in map (fn r => r RS subst) (thm :: map mk_thm arities) end;
berghofe@7015
   387
berghofe@5177
   388
    (* prove  inj dt_Rep_i  and  dt_Rep_i x : dt_rep_set_i *)
berghofe@5177
   389
berghofe@26806
   390
    val fun_congs = map (fn T => make_elim (Drule.instantiate'
berghofe@26806
   391
      [SOME (ctyp_of thy5 T)] [] fun_cong)) branchTs;
berghofe@26806
   392
berghofe@5177
   393
    fun prove_iso_thms (ds, (inj_thms, elem_thms)) =
berghofe@5177
   394
      let
berghofe@5177
   395
        val (_, (tname, _, _)) = hd ds;
wenzelm@17412
   396
        val {induction, ...} = the (Symtab.lookup dt_info tname);
berghofe@5177
   397
berghofe@5177
   398
        fun mk_ind_concl (i, _) =
berghofe@5177
   399
          let
skalberg@15570
   400
            val T = List.nth (recTs, i);
skalberg@15570
   401
            val Rep_t = Const (List.nth (all_rep_names, i), T --> Univ_elT);
skalberg@15570
   402
            val rep_set_name = List.nth (rep_set_names, i)
berghofe@5177
   403
          in (HOLogic.all_const T $ Abs ("y", T, HOLogic.imp $
berghofe@5177
   404
                HOLogic.mk_eq (Rep_t $ mk_Free "x" T i, Rep_t $ Bound 0) $
berghofe@5177
   405
                  HOLogic.mk_eq (mk_Free "x" T i, Bound 0)),
berghofe@21021
   406
              Const (rep_set_name, UnivT') $ (Rep_t $ mk_Free "x" T i))
berghofe@5177
   407
          end;
berghofe@5177
   408
berghofe@5177
   409
        val (ind_concl1, ind_concl2) = ListPair.unzip (map mk_ind_concl ds);
berghofe@5177
   410
oheimb@5553
   411
        val rewrites = map mk_meta_eq iso_char_thms;
berghofe@21021
   412
        val inj_thms' = map snd newT_iso_inj_thms @
haftmann@26359
   413
          map (fn r => r RS @{thm injD}) inj_thms;
berghofe@5177
   414
berghofe@26531
   415
        val inj_thm = SkipProof.prove_global thy5 [] []
wenzelm@17985
   416
          (HOLogic.mk_Trueprop (mk_conj ind_concl1)) (fn _ => EVERY
berghofe@25678
   417
            [(indtac induction [] THEN_ALL_NEW ObjectLogic.atomize_prems_tac) 1,
berghofe@5177
   418
             REPEAT (EVERY
berghofe@5177
   419
               [rtac allI 1, rtac impI 1,
berghofe@5177
   420
                exh_tac (exh_thm_of dt_info) 1,
berghofe@5177
   421
                REPEAT (EVERY
berghofe@5177
   422
                  [hyp_subst_tac 1,
berghofe@5177
   423
                   rewrite_goals_tac rewrites,
berghofe@5177
   424
                   REPEAT (dresolve_tac [In0_inject, In1_inject] 1),
berghofe@5177
   425
                   (eresolve_tac [In0_not_In1 RS notE, In1_not_In0 RS notE] 1)
berghofe@5177
   426
                   ORELSE (EVERY
berghofe@13641
   427
                     [REPEAT (eresolve_tac (Scons_inject ::
berghofe@13641
   428
                        map make_elim [Leaf_inject, Inl_inject, Inr_inject]) 1),
berghofe@13641
   429
                      REPEAT (cong_tac 1), rtac refl 1,
berghofe@13641
   430
                      REPEAT (atac 1 ORELSE (EVERY
berghofe@13641
   431
                        [REPEAT (rtac ext 1),
berghofe@13641
   432
                         REPEAT (eresolve_tac (mp :: allE ::
berghofe@13641
   433
                           map make_elim (Suml_inject :: Sumr_inject ::
berghofe@26806
   434
                             Lim_inject :: inj_thms') @ fun_congs) 1),
wenzelm@20046
   435
                         atac 1]))])])])]);
berghofe@5177
   436
haftmann@26359
   437
        val inj_thms'' = map (fn r => r RS @{thm datatype_injI})
paulson@6171
   438
                             (split_conj_thm inj_thm);
berghofe@5177
   439
paulson@6171
   440
        val elem_thm = 
berghofe@26531
   441
            SkipProof.prove_global thy5 [] [] (HOLogic.mk_Trueprop (mk_conj ind_concl2))
wenzelm@20046
   442
              (fn _ =>
berghofe@25678
   443
               EVERY [(indtac induction [] THEN_ALL_NEW ObjectLogic.atomize_prems_tac) 1,
wenzelm@20046
   444
                rewrite_goals_tac rewrites,
wenzelm@20046
   445
                REPEAT ((resolve_tac rep_intrs THEN_ALL_NEW
wenzelm@20046
   446
                  ((REPEAT o etac allE) THEN' ares_tac elem_thms)) 1)]);
berghofe@5177
   447
berghofe@11471
   448
      in (inj_thms'' @ inj_thms, elem_thms @ (split_conj_thm elem_thm))
berghofe@11471
   449
      end;
berghofe@11471
   450
wenzelm@30190
   451
    val (iso_inj_thms_unfolded, iso_elem_thms) = List.foldr prove_iso_thms
skalberg@15574
   452
      ([], map #3 newT_iso_axms) (tl descr);
berghofe@21021
   453
    val iso_inj_thms = map snd newT_iso_inj_thms @
haftmann@26359
   454
      map (fn r => r RS @{thm injD}) iso_inj_thms_unfolded;
berghofe@11471
   455
berghofe@21021
   456
    (* prove  dt_rep_set_i x --> x : range dt_Rep_i *)
berghofe@11471
   457
berghofe@11471
   458
    fun mk_iso_t (((set_name, iso_name), i), T) =
berghofe@11471
   459
      let val isoT = T --> Univ_elT
berghofe@11471
   460
      in HOLogic.imp $ 
berghofe@21021
   461
        (Const (set_name, UnivT') $ mk_Free "x" Univ_elT i) $
haftmann@31643
   462
          (if i < length newTs then HOLogic.true_const
berghofe@11471
   463
           else HOLogic.mk_mem (mk_Free "x" Univ_elT i,
haftmann@31643
   464
             Const (@{const_name image}, isoT --> HOLogic.mk_setT T --> UnivT) $
haftmann@30304
   465
               Const (iso_name, isoT) $ Const (@{const_name UNIV}, HOLogic.mk_setT T)))
berghofe@5177
   466
      end;
berghofe@5177
   467
berghofe@11471
   468
    val iso_t = HOLogic.mk_Trueprop (mk_conj (map mk_iso_t
berghofe@11471
   469
      (rep_set_names ~~ all_rep_names ~~ (0 upto (length descr' - 1)) ~~ recTs)));
berghofe@11471
   470
berghofe@11471
   471
    (* all the theorems are proved by one single simultaneous induction *)
berghofe@11471
   472
haftmann@26359
   473
    val range_eqs = map (fn r => mk_meta_eq (r RS @{thm range_ex1_eq}))
berghofe@13641
   474
      iso_inj_thms_unfolded;
berghofe@13641
   475
berghofe@11471
   476
    val iso_thms = if length descr = 1 then [] else
skalberg@15570
   477
      Library.drop (length newTs, split_conj_thm
berghofe@26531
   478
        (SkipProof.prove_global thy5 [] [] iso_t (fn _ => EVERY
berghofe@25678
   479
           [(indtac rep_induct [] THEN_ALL_NEW ObjectLogic.atomize_prems_tac) 1,
berghofe@11471
   480
            REPEAT (rtac TrueI 1),
berghofe@13641
   481
            rewrite_goals_tac (mk_meta_eq choice_eq ::
haftmann@26359
   482
              symmetric (mk_meta_eq @{thm expand_fun_eq}) :: range_eqs),
berghofe@13641
   483
            rewrite_goals_tac (map symmetric range_eqs),
berghofe@11471
   484
            REPEAT (EVERY
berghofe@13641
   485
              [REPEAT (eresolve_tac ([rangeE, ex1_implies_ex RS exE] @
wenzelm@28362
   486
                 maps (mk_funs_inv thy5 o #1) newT_iso_axms) 1),
berghofe@11471
   487
               TRY (hyp_subst_tac 1),
berghofe@11471
   488
               rtac (sym RS range_eqI) 1,
wenzelm@20046
   489
               resolve_tac iso_char_thms 1])])));
wenzelm@11435
   490
wenzelm@11435
   491
    val Abs_inverse_thms' =
wenzelm@11435
   492
      map #1 newT_iso_axms @
haftmann@18330
   493
      map2 (fn r_inj => fn r => f_myinv_f OF [r_inj, r RS mp])
haftmann@18330
   494
        iso_inj_thms_unfolded iso_thms;
wenzelm@11435
   495
wenzelm@28362
   496
    val Abs_inverse_thms = maps (mk_funs_inv thy5) Abs_inverse_thms';
berghofe@5177
   497
berghofe@5177
   498
    (******************* freeness theorems for constructors *******************)
berghofe@5177
   499
haftmann@31668
   500
    val _ = message config "Proving freeness of constructors ...";
berghofe@5177
   501
berghofe@5177
   502
    (* prove theorem  Rep_i (Constr_j ...) = Inj_j ...  *)
berghofe@5177
   503
    
berghofe@5177
   504
    fun prove_constr_rep_thm eqn =
berghofe@5177
   505
      let
berghofe@21021
   506
        val inj_thms = map fst newT_iso_inj_thms;
haftmann@26359
   507
        val rewrites = @{thm o_def} :: constr_defs @ (map (mk_meta_eq o #2) newT_iso_axms)
berghofe@26531
   508
      in SkipProof.prove_global thy5 [] [] eqn (fn _ => EVERY
berghofe@5177
   509
        [resolve_tac inj_thms 1,
berghofe@5177
   510
         rewrite_goals_tac rewrites,
berghofe@21021
   511
         rtac refl 3,
berghofe@5177
   512
         resolve_tac rep_intrs 2,
wenzelm@20046
   513
         REPEAT (resolve_tac iso_elem_thms 1)])
berghofe@5177
   514
      end;
berghofe@5177
   515
berghofe@5177
   516
    (*--------------------------------------------------------------*)
berghofe@5177
   517
    (* constr_rep_thms and rep_congs are used to prove distinctness *)
berghofe@7015
   518
    (* of constructors.                                             *)
berghofe@5177
   519
    (*--------------------------------------------------------------*)
berghofe@5177
   520
berghofe@5177
   521
    val constr_rep_thms = map (map prove_constr_rep_thm) constr_rep_eqns;
berghofe@5177
   522
berghofe@5177
   523
    val dist_rewrites = map (fn (rep_thms, dist_lemma) =>
berghofe@5177
   524
      dist_lemma::(rep_thms @ [In0_eq, In1_eq, In0_not_In1, In1_not_In0]))
berghofe@5177
   525
        (constr_rep_thms ~~ dist_lemmas);
berghofe@5177
   526
haftmann@26969
   527
    fun prove_distinct_thms _ _ (_, []) = []
haftmann@26969
   528
      | prove_distinct_thms lim dist_rewrites' (k, ts as _ :: _) =
haftmann@26969
   529
          if k >= lim then [] else let
haftmann@27002
   530
            (*number of constructors < distinctness_limit : C_i ... ~= C_j ...*)
haftmann@26969
   531
            fun prove [] = []
haftmann@27300
   532
              | prove (t :: ts) =
haftmann@26969
   533
                  let
haftmann@26969
   534
                    val dist_thm = SkipProof.prove_global thy5 [] [] t (fn _ =>
haftmann@26969
   535
                      EVERY [simp_tac (HOL_ss addsimps dist_rewrites') 1])
haftmann@26969
   536
                  in dist_thm :: standard (dist_thm RS not_sym) :: prove ts end;
haftmann@26969
   537
          in prove ts end;
berghofe@7015
   538
haftmann@27300
   539
    val distinct_thms = DatatypeProp.make_distincts descr sorts
haftmann@27300
   540
      |> map2 (prove_distinct_thms
haftmann@27300
   541
           (Config.get_thy thy5 distinctness_limit)) dist_rewrites;
berghofe@7015
   542
berghofe@7015
   543
    val simproc_dists = map (fn ((((_, (_, _, constrs)), rep_thms), congr), dists) =>
haftmann@27002
   544
      if length constrs < Config.get_thy thy5 distinctness_limit
wenzelm@24098
   545
      then FewConstrs dists
berghofe@7015
   546
      else ManyConstrs (congr, HOL_basic_ss addsimps rep_thms)) (hd descr ~~
berghofe@7015
   547
        constr_rep_thms ~~ rep_congs ~~ distinct_thms);
berghofe@7015
   548
berghofe@5177
   549
    (* prove injectivity of constructors *)
berghofe@5177
   550
berghofe@5177
   551
    fun prove_constr_inj_thm rep_thms t =
berghofe@13641
   552
      let val inj_thms = Scons_inject :: (map make_elim
berghofe@21021
   553
        (iso_inj_thms @
berghofe@13641
   554
          [In0_inject, In1_inject, Leaf_inject, Inl_inject, Inr_inject,
berghofe@13641
   555
           Lim_inject, Suml_inject, Sumr_inject]))
berghofe@26531
   556
      in SkipProof.prove_global thy5 [] [] t (fn _ => EVERY
berghofe@5177
   557
        [rtac iffI 1,
berghofe@5177
   558
         REPEAT (etac conjE 2), hyp_subst_tac 2, rtac refl 2,
berghofe@5177
   559
         dresolve_tac rep_congs 1, dtac box_equals 1,
berghofe@13641
   560
         REPEAT (resolve_tac rep_thms 1),
berghofe@5177
   561
         REPEAT (eresolve_tac inj_thms 1),
berghofe@13641
   562
         REPEAT (ares_tac [conjI] 1 ORELSE (EVERY [REPEAT (rtac ext 1),
berghofe@13641
   563
           REPEAT (eresolve_tac (make_elim fun_cong :: inj_thms) 1),
wenzelm@20046
   564
           atac 1]))])
berghofe@5177
   565
      end;
berghofe@5177
   566
berghofe@5177
   567
    val constr_inject = map (fn (ts, thms) => map (prove_constr_inj_thm thms) ts)
berghofe@5177
   568
      ((DatatypeProp.make_injs descr sorts) ~~ constr_rep_thms);
berghofe@5177
   569
haftmann@18314
   570
    val ((constr_inject', distinct_thms'), thy6) =
haftmann@18314
   571
      thy5
haftmann@31668
   572
      |> parent_path (#flat_names config)
haftmann@18314
   573
      |> store_thmss "inject" new_type_names constr_inject
haftmann@18314
   574
      ||>> store_thmss "distinct" new_type_names distinct_thms;
berghofe@5177
   575
berghofe@5177
   576
    (*************************** induction theorem ****************************)
berghofe@5177
   577
haftmann@31668
   578
    val _ = message config "Proving induction rule for datatypes ...";
berghofe@5177
   579
berghofe@5177
   580
    val Rep_inverse_thms = (map (fn (_, iso, _) => iso RS subst) newT_iso_axms) @
berghofe@11471
   581
      (map (fn r => r RS myinv_f_f RS subst) iso_inj_thms_unfolded);
berghofe@11471
   582
    val Rep_inverse_thms' = map (fn r => r RS myinv_f_f) iso_inj_thms_unfolded;
berghofe@5177
   583
berghofe@5177
   584
    fun mk_indrule_lemma ((prems, concls), ((i, _), T)) =
berghofe@5177
   585
      let
skalberg@15570
   586
        val Rep_t = Const (List.nth (all_rep_names, i), T --> Univ_elT) $
berghofe@5177
   587
          mk_Free "x" T i;
berghofe@5177
   588
berghofe@5177
   589
        val Abs_t = if i < length newTs then
wenzelm@22578
   590
            Const (Sign.intern_const thy6
skalberg@15570
   591
              ("Abs_" ^ (List.nth (new_type_names, i))), Univ_elT --> T)
wenzelm@11435
   592
          else Const ("Inductive.myinv", [T --> Univ_elT, Univ_elT] ---> T) $
skalberg@15570
   593
            Const (List.nth (all_rep_names, i), T --> Univ_elT)
berghofe@5177
   594
berghofe@21021
   595
      in (prems @ [HOLogic.imp $
berghofe@21021
   596
            (Const (List.nth (rep_set_names, i), UnivT') $ Rep_t) $
berghofe@5177
   597
              (mk_Free "P" (T --> HOLogic.boolT) (i + 1) $ (Abs_t $ Rep_t))],
berghofe@5177
   598
          concls @ [mk_Free "P" (T --> HOLogic.boolT) (i + 1) $ mk_Free "x" T i])
berghofe@5177
   599
      end;
berghofe@5177
   600
berghofe@5177
   601
    val (indrule_lemma_prems, indrule_lemma_concls) =
skalberg@15570
   602
      Library.foldl mk_indrule_lemma (([], []), (descr' ~~ recTs));
berghofe@5177
   603
wenzelm@22578
   604
    val cert = cterm_of thy6;
berghofe@5177
   605
berghofe@26531
   606
    val indrule_lemma = SkipProof.prove_global thy6 [] []
berghofe@5177
   607
      (Logic.mk_implies
berghofe@5177
   608
        (HOLogic.mk_Trueprop (mk_conj indrule_lemma_prems),
wenzelm@17985
   609
         HOLogic.mk_Trueprop (mk_conj indrule_lemma_concls))) (fn _ => EVERY
wenzelm@17985
   610
           [REPEAT (etac conjE 1),
berghofe@5177
   611
            REPEAT (EVERY
berghofe@5177
   612
              [TRY (rtac conjI 1), resolve_tac Rep_inverse_thms 1,
wenzelm@20046
   613
               etac mp 1, resolve_tac iso_elem_thms 1])]);
berghofe@5177
   614
wenzelm@8305
   615
    val Ps = map head_of (HOLogic.dest_conj (HOLogic.dest_Trueprop (concl_of indrule_lemma)));
berghofe@5177
   616
    val frees = if length Ps = 1 then [Free ("P", snd (dest_Var (hd Ps)))] else
berghofe@5177
   617
      map (Free o apfst fst o dest_Var) Ps;
berghofe@5177
   618
    val indrule_lemma' = cterm_instantiate (map cert Ps ~~ map cert frees) indrule_lemma;
berghofe@5177
   619
wenzelm@17985
   620
    val dt_induct_prop = DatatypeProp.make_ind descr sorts;
berghofe@26531
   621
    val dt_induct = SkipProof.prove_global thy6 []
wenzelm@17985
   622
      (Logic.strip_imp_prems dt_induct_prop) (Logic.strip_imp_concl dt_induct_prop)
wenzelm@26711
   623
      (fn {prems, ...} => EVERY
berghofe@13641
   624
        [rtac indrule_lemma' 1,
berghofe@25678
   625
         (indtac rep_induct [] THEN_ALL_NEW ObjectLogic.atomize_prems_tac) 1,
berghofe@5177
   626
         EVERY (map (fn (prem, r) => (EVERY
berghofe@13641
   627
           [REPEAT (eresolve_tac Abs_inverse_thms 1),
berghofe@5177
   628
            simp_tac (HOL_basic_ss addsimps ((symmetric r)::Rep_inverse_thms')) 1,
berghofe@13641
   629
            DEPTH_SOLVE_1 (ares_tac [prem] 1 ORELSE etac allE 1)]))
wenzelm@20046
   630
                (prems ~~ (constr_defs @ (map mk_meta_eq iso_char_thms))))]);
berghofe@5177
   631
haftmann@18377
   632
    val ([dt_induct'], thy7) =
haftmann@18377
   633
      thy6
wenzelm@24712
   634
      |> Sign.add_path big_name
haftmann@29579
   635
      |> PureThy.add_thms [((Binding.name "induct", dt_induct), [case_names_induct])]
wenzelm@28362
   636
      ||> Sign.parent_path
wenzelm@28362
   637
      ||> Theory.checkpoint;
berghofe@5177
   638
haftmann@18314
   639
  in
haftmann@18314
   640
    ((constr_inject', distinct_thms', dist_rewrites, simproc_dists, dt_induct'), thy7)
berghofe@5177
   641
  end;
berghofe@5177
   642
berghofe@5177
   643
end;