src/HOLCF/Universal.thy
author huffman
Wed Oct 06 10:49:27 2010 -0700 (2010-10-06)
changeset 39974 b525988432e9
parent 36452 d37c6eed8117
child 39984 0300d5170622
permissions -rw-r--r--
major reorganization/simplification of HOLCF type classes:
removed profinite/bifinite classes and approx function;
universal domain uses approx_chain locale instead of bifinite class;
ideal_completion locale does not use 'take' functions, requires countable basis instead;
replaced type 'udom alg_defl' with type 'sfp';
replaced class 'rep' with class 'sfp';
renamed REP('a) to SFP('a);
huffman@27411
     1
(*  Title:      HOLCF/Universal.thy
huffman@27411
     2
    Author:     Brian Huffman
huffman@27411
     3
*)
huffman@27411
     4
huffman@35794
     5
header {* A universal bifinite domain *}
huffman@35794
     6
huffman@27411
     7
theory Universal
huffman@39974
     8
imports Completion Deflation Nat_Bijection
huffman@27411
     9
begin
huffman@27411
    10
huffman@39974
    11
subsection {* Basis for universal domain *}
huffman@39974
    12
huffman@39974
    13
subsubsection {* Basis datatype *}
huffman@27411
    14
huffman@27411
    15
types ubasis = nat
huffman@27411
    16
huffman@27411
    17
definition
huffman@27411
    18
  node :: "nat \<Rightarrow> ubasis \<Rightarrow> ubasis set \<Rightarrow> ubasis"
huffman@27411
    19
where
huffman@35701
    20
  "node i a S = Suc (prod_encode (i, prod_encode (a, set_encode S)))"
huffman@27411
    21
huffman@30505
    22
lemma node_not_0 [simp]: "node i a S \<noteq> 0"
huffman@27411
    23
unfolding node_def by simp
huffman@27411
    24
huffman@30505
    25
lemma node_gt_0 [simp]: "0 < node i a S"
huffman@27411
    26
unfolding node_def by simp
huffman@27411
    27
huffman@27411
    28
lemma node_inject [simp]:
huffman@30505
    29
  "\<lbrakk>finite S; finite T\<rbrakk>
huffman@30505
    30
    \<Longrightarrow> node i a S = node j b T \<longleftrightarrow> i = j \<and> a = b \<and> S = T"
huffman@35701
    31
unfolding node_def by (simp add: prod_encode_eq set_encode_eq)
huffman@27411
    32
huffman@30505
    33
lemma node_gt0: "i < node i a S"
huffman@27411
    34
unfolding node_def less_Suc_eq_le
huffman@35701
    35
by (rule le_prod_encode_1)
huffman@27411
    36
huffman@30505
    37
lemma node_gt1: "a < node i a S"
huffman@27411
    38
unfolding node_def less_Suc_eq_le
huffman@35701
    39
by (rule order_trans [OF le_prod_encode_1 le_prod_encode_2])
huffman@27411
    40
huffman@27411
    41
lemma nat_less_power2: "n < 2^n"
huffman@27411
    42
by (induct n) simp_all
huffman@27411
    43
huffman@30505
    44
lemma node_gt2: "\<lbrakk>finite S; b \<in> S\<rbrakk> \<Longrightarrow> b < node i a S"
huffman@35701
    45
unfolding node_def less_Suc_eq_le set_encode_def
huffman@35701
    46
apply (rule order_trans [OF _ le_prod_encode_2])
huffman@35701
    47
apply (rule order_trans [OF _ le_prod_encode_2])
huffman@30505
    48
apply (rule order_trans [where y="setsum (op ^ 2) {b}"])
huffman@27411
    49
apply (simp add: nat_less_power2 [THEN order_less_imp_le])
huffman@27411
    50
apply (erule setsum_mono2, simp, simp)
huffman@27411
    51
done
huffman@27411
    52
huffman@35701
    53
lemma eq_prod_encode_pairI:
huffman@35701
    54
  "\<lbrakk>fst (prod_decode x) = a; snd (prod_decode x) = b\<rbrakk> \<Longrightarrow> x = prod_encode (a, b)"
huffman@27411
    55
by (erule subst, erule subst, simp)
huffman@27411
    56
huffman@27411
    57
lemma node_cases:
huffman@27411
    58
  assumes 1: "x = 0 \<Longrightarrow> P"
huffman@30505
    59
  assumes 2: "\<And>i a S. \<lbrakk>finite S; x = node i a S\<rbrakk> \<Longrightarrow> P"
huffman@27411
    60
  shows "P"
huffman@27411
    61
 apply (cases x)
huffman@27411
    62
  apply (erule 1)
huffman@27411
    63
 apply (rule 2)
huffman@35701
    64
  apply (rule finite_set_decode)
huffman@27411
    65
 apply (simp add: node_def)
huffman@35701
    66
 apply (rule eq_prod_encode_pairI [OF refl])
huffman@35701
    67
 apply (rule eq_prod_encode_pairI [OF refl refl])
huffman@27411
    68
done
huffman@27411
    69
huffman@27411
    70
lemma node_induct:
huffman@27411
    71
  assumes 1: "P 0"
huffman@30505
    72
  assumes 2: "\<And>i a S. \<lbrakk>P a; finite S; \<forall>b\<in>S. P b\<rbrakk> \<Longrightarrow> P (node i a S)"
huffman@27411
    73
  shows "P x"
huffman@27411
    74
 apply (induct x rule: nat_less_induct)
huffman@27411
    75
 apply (case_tac n rule: node_cases)
huffman@27411
    76
  apply (simp add: 1)
huffman@27411
    77
 apply (simp add: 2 node_gt1 node_gt2)
huffman@27411
    78
done
huffman@27411
    79
huffman@39974
    80
subsubsection {* Basis ordering *}
huffman@27411
    81
huffman@27411
    82
inductive
huffman@27411
    83
  ubasis_le :: "nat \<Rightarrow> nat \<Rightarrow> bool"
huffman@27411
    84
where
huffman@30505
    85
  ubasis_le_refl: "ubasis_le a a"
huffman@27411
    86
| ubasis_le_trans:
huffman@30505
    87
    "\<lbrakk>ubasis_le a b; ubasis_le b c\<rbrakk> \<Longrightarrow> ubasis_le a c"
huffman@27411
    88
| ubasis_le_lower:
huffman@30505
    89
    "finite S \<Longrightarrow> ubasis_le a (node i a S)"
huffman@27411
    90
| ubasis_le_upper:
huffman@30505
    91
    "\<lbrakk>finite S; b \<in> S; ubasis_le a b\<rbrakk> \<Longrightarrow> ubasis_le (node i a S) b"
huffman@27411
    92
huffman@27411
    93
lemma ubasis_le_minimal: "ubasis_le 0 x"
huffman@27411
    94
apply (induct x rule: node_induct)
huffman@27411
    95
apply (rule ubasis_le_refl)
huffman@27411
    96
apply (erule ubasis_le_trans)
huffman@27411
    97
apply (erule ubasis_le_lower)
huffman@27411
    98
done
huffman@27411
    99
huffman@39974
   100
interpretation udom: preorder ubasis_le
huffman@39974
   101
apply default
huffman@39974
   102
apply (rule ubasis_le_refl)
huffman@39974
   103
apply (erule (1) ubasis_le_trans)
huffman@39974
   104
done
huffman@39974
   105
huffman@27411
   106
subsubsection {* Generic take function *}
huffman@27411
   107
huffman@27411
   108
function
huffman@27411
   109
  ubasis_until :: "(ubasis \<Rightarrow> bool) \<Rightarrow> ubasis \<Rightarrow> ubasis"
huffman@27411
   110
where
huffman@27411
   111
  "ubasis_until P 0 = 0"
huffman@30505
   112
| "finite S \<Longrightarrow> ubasis_until P (node i a S) =
huffman@30505
   113
    (if P (node i a S) then node i a S else ubasis_until P a)"
huffman@27411
   114
    apply clarify
huffman@27411
   115
    apply (rule_tac x=b in node_cases)
huffman@27411
   116
     apply simp
huffman@27411
   117
    apply simp
huffman@27411
   118
    apply fast
huffman@27411
   119
   apply simp
huffman@27411
   120
  apply simp
huffman@27411
   121
 apply simp
huffman@27411
   122
done
huffman@27411
   123
huffman@27411
   124
termination ubasis_until
huffman@27411
   125
apply (relation "measure snd")
huffman@27411
   126
apply (rule wf_measure)
huffman@27411
   127
apply (simp add: node_gt1)
huffman@27411
   128
done
huffman@27411
   129
huffman@27411
   130
lemma ubasis_until: "P 0 \<Longrightarrow> P (ubasis_until P x)"
huffman@27411
   131
by (induct x rule: node_induct) simp_all
huffman@27411
   132
huffman@27411
   133
lemma ubasis_until': "0 < ubasis_until P x \<Longrightarrow> P (ubasis_until P x)"
huffman@27411
   134
by (induct x rule: node_induct) auto
huffman@27411
   135
huffman@27411
   136
lemma ubasis_until_same: "P x \<Longrightarrow> ubasis_until P x = x"
huffman@27411
   137
by (induct x rule: node_induct) simp_all
huffman@27411
   138
huffman@27411
   139
lemma ubasis_until_idem:
huffman@27411
   140
  "P 0 \<Longrightarrow> ubasis_until P (ubasis_until P x) = ubasis_until P x"
huffman@27411
   141
by (rule ubasis_until_same [OF ubasis_until])
huffman@27411
   142
huffman@27411
   143
lemma ubasis_until_0:
huffman@27411
   144
  "\<forall>x. x \<noteq> 0 \<longrightarrow> \<not> P x \<Longrightarrow> ubasis_until P x = 0"
huffman@27411
   145
by (induct x rule: node_induct) simp_all
huffman@27411
   146
huffman@27411
   147
lemma ubasis_until_less: "ubasis_le (ubasis_until P x) x"
huffman@27411
   148
apply (induct x rule: node_induct)
huffman@27411
   149
apply (simp add: ubasis_le_refl)
huffman@27411
   150
apply (simp add: ubasis_le_refl)
huffman@27411
   151
apply (rule impI)
huffman@27411
   152
apply (erule ubasis_le_trans)
huffman@27411
   153
apply (erule ubasis_le_lower)
huffman@27411
   154
done
huffman@27411
   155
huffman@27411
   156
lemma ubasis_until_chain:
huffman@27411
   157
  assumes PQ: "\<And>x. P x \<Longrightarrow> Q x"
huffman@27411
   158
  shows "ubasis_le (ubasis_until P x) (ubasis_until Q x)"
huffman@27411
   159
apply (induct x rule: node_induct)
huffman@27411
   160
apply (simp add: ubasis_le_refl)
huffman@27411
   161
apply (simp add: ubasis_le_refl)
huffman@27411
   162
apply (simp add: PQ)
huffman@27411
   163
apply clarify
huffman@27411
   164
apply (rule ubasis_le_trans)
huffman@27411
   165
apply (rule ubasis_until_less)
huffman@27411
   166
apply (erule ubasis_le_lower)
huffman@27411
   167
done
huffman@27411
   168
huffman@27411
   169
lemma ubasis_until_mono:
huffman@30505
   170
  assumes "\<And>i a S b. \<lbrakk>finite S; P (node i a S); b \<in> S; ubasis_le a b\<rbrakk> \<Longrightarrow> P b"
huffman@30505
   171
  shows "ubasis_le a b \<Longrightarrow> ubasis_le (ubasis_until P a) (ubasis_until P b)"
huffman@30561
   172
proof (induct set: ubasis_le)
huffman@30561
   173
  case (ubasis_le_refl a) show ?case by (rule ubasis_le.ubasis_le_refl)
huffman@30561
   174
next
huffman@30561
   175
  case (ubasis_le_trans a b c) thus ?case by - (rule ubasis_le.ubasis_le_trans)
huffman@30561
   176
next
huffman@30561
   177
  case (ubasis_le_lower S a i) thus ?case
huffman@30561
   178
    apply (clarsimp simp add: ubasis_le_refl)
huffman@30561
   179
    apply (rule ubasis_le_trans [OF ubasis_until_less])
huffman@30561
   180
    apply (erule ubasis_le.ubasis_le_lower)
huffman@30561
   181
    done
huffman@30561
   182
next
huffman@30561
   183
  case (ubasis_le_upper S b a i) thus ?case
huffman@30561
   184
    apply clarsimp
huffman@30561
   185
    apply (subst ubasis_until_same)
huffman@30561
   186
     apply (erule (3) prems)
huffman@30561
   187
    apply (erule (2) ubasis_le.ubasis_le_upper)
huffman@30561
   188
    done
huffman@30561
   189
qed
huffman@27411
   190
huffman@27411
   191
lemma finite_range_ubasis_until:
huffman@27411
   192
  "finite {x. P x} \<Longrightarrow> finite (range (ubasis_until P))"
huffman@27411
   193
apply (rule finite_subset [where B="insert 0 {x. P x}"])
huffman@27411
   194
apply (clarsimp simp add: ubasis_until')
huffman@27411
   195
apply simp
huffman@27411
   196
done
huffman@27411
   197
huffman@27411
   198
huffman@27411
   199
subsection {* Defining the universal domain by ideal completion *}
huffman@27411
   200
huffman@27411
   201
typedef (open) udom = "{S. udom.ideal S}"
huffman@27411
   202
by (fast intro: udom.ideal_principal)
huffman@27411
   203
huffman@31076
   204
instantiation udom :: below
huffman@27411
   205
begin
huffman@27411
   206
huffman@27411
   207
definition
huffman@27411
   208
  "x \<sqsubseteq> y \<longleftrightarrow> Rep_udom x \<subseteq> Rep_udom y"
huffman@27411
   209
huffman@27411
   210
instance ..
huffman@27411
   211
end
huffman@27411
   212
huffman@27411
   213
instance udom :: po
huffman@39974
   214
using type_definition_udom below_udom_def
huffman@39974
   215
by (rule udom.typedef_ideal_po)
huffman@27411
   216
huffman@27411
   217
instance udom :: cpo
huffman@39974
   218
using type_definition_udom below_udom_def
huffman@39974
   219
by (rule udom.typedef_ideal_cpo)
huffman@27411
   220
huffman@27411
   221
lemma Rep_udom_lub:
huffman@27411
   222
  "chain Y \<Longrightarrow> Rep_udom (\<Squnion>i. Y i) = (\<Union>i. Rep_udom (Y i))"
huffman@39974
   223
using type_definition_udom below_udom_def
huffman@39974
   224
by (rule udom.typedef_ideal_rep_contlub)
huffman@27411
   225
huffman@27411
   226
lemma ideal_Rep_udom: "udom.ideal (Rep_udom xs)"
huffman@27411
   227
by (rule Rep_udom [unfolded mem_Collect_eq])
huffman@27411
   228
huffman@27411
   229
definition
huffman@27411
   230
  udom_principal :: "nat \<Rightarrow> udom" where
huffman@27411
   231
  "udom_principal t = Abs_udom {u. ubasis_le u t}"
huffman@27411
   232
huffman@27411
   233
lemma Rep_udom_principal:
huffman@27411
   234
  "Rep_udom (udom_principal t) = {u. ubasis_le u t}"
huffman@27411
   235
unfolding udom_principal_def
huffman@27411
   236
by (simp add: Abs_udom_inverse udom.ideal_principal)
huffman@27411
   237
wenzelm@30729
   238
interpretation udom:
huffman@39974
   239
  ideal_completion ubasis_le udom_principal Rep_udom
huffman@27411
   240
apply unfold_locales
huffman@27411
   241
apply (rule ideal_Rep_udom)
huffman@27411
   242
apply (erule Rep_udom_lub)
huffman@27411
   243
apply (rule Rep_udom_principal)
huffman@31076
   244
apply (simp only: below_udom_def)
huffman@39974
   245
apply (rule exI, rule inj_on_id)
huffman@27411
   246
done
huffman@27411
   247
huffman@27411
   248
text {* Universal domain is pointed *}
huffman@27411
   249
huffman@27411
   250
lemma udom_minimal: "udom_principal 0 \<sqsubseteq> x"
huffman@27411
   251
apply (induct x rule: udom.principal_induct)
huffman@27411
   252
apply (simp, simp add: ubasis_le_minimal)
huffman@27411
   253
done
huffman@27411
   254
huffman@27411
   255
instance udom :: pcpo
huffman@27411
   256
by intro_classes (fast intro: udom_minimal)
huffman@27411
   257
huffman@27411
   258
lemma inst_udom_pcpo: "\<bottom> = udom_principal 0"
huffman@27411
   259
by (rule udom_minimal [THEN UU_I, symmetric])
huffman@27411
   260
huffman@39974
   261
huffman@39974
   262
subsection {* Compact bases of domains *}
huffman@27411
   263
huffman@39974
   264
typedef (open) 'a compact_basis = "{x::'a::pcpo. compact x}"
huffman@39974
   265
by auto
huffman@39974
   266
huffman@39974
   267
lemma compact_Rep_compact_basis: "compact (Rep_compact_basis a)"
huffman@39974
   268
by (rule Rep_compact_basis [unfolded mem_Collect_eq])
huffman@39974
   269
huffman@39974
   270
instantiation compact_basis :: (pcpo) below
huffman@27411
   271
begin
huffman@27411
   272
huffman@27411
   273
definition
huffman@39974
   274
  compact_le_def:
huffman@39974
   275
    "(op \<sqsubseteq>) \<equiv> (\<lambda>x y. Rep_compact_basis x \<sqsubseteq> Rep_compact_basis y)"
huffman@27411
   276
huffman@39974
   277
instance ..
huffman@27411
   278
end
huffman@27411
   279
huffman@39974
   280
instance compact_basis :: (pcpo) po
huffman@39974
   281
using type_definition_compact_basis compact_le_def
huffman@39974
   282
by (rule typedef_po)
huffman@39974
   283
huffman@39974
   284
definition
huffman@39974
   285
  approximants :: "'a \<Rightarrow> 'a compact_basis set" where
huffman@39974
   286
  "approximants = (\<lambda>x. {a. Rep_compact_basis a \<sqsubseteq> x})"
huffman@27411
   287
huffman@39974
   288
definition
huffman@39974
   289
  compact_bot :: "'a::pcpo compact_basis" where
huffman@39974
   290
  "compact_bot = Abs_compact_basis \<bottom>"
huffman@39974
   291
huffman@39974
   292
lemma Rep_compact_bot [simp]: "Rep_compact_basis compact_bot = \<bottom>"
huffman@39974
   293
unfolding compact_bot_def by (simp add: Abs_compact_basis_inverse)
huffman@39974
   294
huffman@39974
   295
lemma compact_bot_minimal [simp]: "compact_bot \<sqsubseteq> a"
huffman@39974
   296
unfolding compact_le_def Rep_compact_bot by simp
huffman@27411
   297
huffman@27411
   298
huffman@35900
   299
subsection {* Universality of \emph{udom} *}
huffman@27411
   300
huffman@39974
   301
text {* We use a locale to parameterize the construction over a chain
huffman@39974
   302
of approx functions on the type to be embedded. *}
huffman@39974
   303
huffman@39974
   304
locale approx_chain =
huffman@39974
   305
  fixes approx :: "nat \<Rightarrow> 'a::pcpo \<rightarrow> 'a"
huffman@39974
   306
  assumes chain_approx [simp]: "chain (\<lambda>i. approx i)"
huffman@39974
   307
  assumes lub_approx [simp]: "(\<Squnion>i. approx i) = ID"
huffman@39974
   308
  assumes finite_deflation_approx: "\<And>i. finite_deflation (approx i)"
huffman@39974
   309
begin
huffman@27411
   310
huffman@27411
   311
subsubsection {* Choosing a maximal element from a finite set *}
huffman@27411
   312
huffman@27411
   313
lemma finite_has_maximal:
huffman@39974
   314
  fixes A :: "'a compact_basis set"
huffman@27411
   315
  shows "\<lbrakk>finite A; A \<noteq> {}\<rbrakk> \<Longrightarrow> \<exists>x\<in>A. \<forall>y\<in>A. x \<sqsubseteq> y \<longrightarrow> x = y"
huffman@27411
   316
proof (induct rule: finite_ne_induct)
huffman@27411
   317
  case (singleton x)
huffman@27411
   318
    show ?case by simp
huffman@27411
   319
next
huffman@27411
   320
  case (insert a A)
huffman@27411
   321
  from `\<exists>x\<in>A. \<forall>y\<in>A. x \<sqsubseteq> y \<longrightarrow> x = y`
huffman@27411
   322
  obtain x where x: "x \<in> A"
huffman@27411
   323
           and x_eq: "\<And>y. \<lbrakk>y \<in> A; x \<sqsubseteq> y\<rbrakk> \<Longrightarrow> x = y" by fast
huffman@27411
   324
  show ?case
huffman@27411
   325
  proof (intro bexI ballI impI)
huffman@27411
   326
    fix y
huffman@27411
   327
    assume "y \<in> insert a A" and "(if x \<sqsubseteq> a then a else x) \<sqsubseteq> y"
huffman@27411
   328
    thus "(if x \<sqsubseteq> a then a else x) = y"
huffman@27411
   329
      apply auto
huffman@31076
   330
      apply (frule (1) below_trans)
huffman@27411
   331
      apply (frule (1) x_eq)
huffman@31076
   332
      apply (rule below_antisym, assumption)
huffman@27411
   333
      apply simp
huffman@27411
   334
      apply (erule (1) x_eq)
huffman@27411
   335
      done
huffman@27411
   336
  next
huffman@27411
   337
    show "(if x \<sqsubseteq> a then a else x) \<in> insert a A"
huffman@27411
   338
      by (simp add: x)
huffman@27411
   339
  qed
huffman@27411
   340
qed
huffman@27411
   341
huffman@27411
   342
definition
huffman@27411
   343
  choose :: "'a compact_basis set \<Rightarrow> 'a compact_basis"
huffman@27411
   344
where
huffman@27411
   345
  "choose A = (SOME x. x \<in> {x\<in>A. \<forall>y\<in>A. x \<sqsubseteq> y \<longrightarrow> x = y})"
huffman@27411
   346
huffman@27411
   347
lemma choose_lemma:
huffman@27411
   348
  "\<lbrakk>finite A; A \<noteq> {}\<rbrakk> \<Longrightarrow> choose A \<in> {x\<in>A. \<forall>y\<in>A. x \<sqsubseteq> y \<longrightarrow> x = y}"
huffman@27411
   349
unfolding choose_def
huffman@27411
   350
apply (rule someI_ex)
huffman@27411
   351
apply (frule (1) finite_has_maximal, fast)
huffman@27411
   352
done
huffman@27411
   353
huffman@27411
   354
lemma maximal_choose:
huffman@27411
   355
  "\<lbrakk>finite A; y \<in> A; choose A \<sqsubseteq> y\<rbrakk> \<Longrightarrow> choose A = y"
huffman@27411
   356
apply (cases "A = {}", simp)
huffman@27411
   357
apply (frule (1) choose_lemma, simp)
huffman@27411
   358
done
huffman@27411
   359
huffman@27411
   360
lemma choose_in: "\<lbrakk>finite A; A \<noteq> {}\<rbrakk> \<Longrightarrow> choose A \<in> A"
huffman@27411
   361
by (frule (1) choose_lemma, simp)
huffman@27411
   362
huffman@27411
   363
function
huffman@27411
   364
  choose_pos :: "'a compact_basis set \<Rightarrow> 'a compact_basis \<Rightarrow> nat"
huffman@27411
   365
where
huffman@27411
   366
  "choose_pos A x =
huffman@27411
   367
    (if finite A \<and> x \<in> A \<and> x \<noteq> choose A
huffman@27411
   368
      then Suc (choose_pos (A - {choose A}) x) else 0)"
huffman@27411
   369
by auto
huffman@27411
   370
huffman@27411
   371
termination choose_pos
huffman@27411
   372
apply (relation "measure (card \<circ> fst)", simp)
huffman@27411
   373
apply clarsimp
huffman@27411
   374
apply (rule card_Diff1_less)
huffman@27411
   375
apply assumption
huffman@27411
   376
apply (erule choose_in)
huffman@27411
   377
apply clarsimp
huffman@27411
   378
done
huffman@27411
   379
huffman@27411
   380
declare choose_pos.simps [simp del]
huffman@27411
   381
huffman@27411
   382
lemma choose_pos_choose: "finite A \<Longrightarrow> choose_pos A (choose A) = 0"
huffman@27411
   383
by (simp add: choose_pos.simps)
huffman@27411
   384
huffman@27411
   385
lemma inj_on_choose_pos [OF refl]:
huffman@27411
   386
  "\<lbrakk>card A = n; finite A\<rbrakk> \<Longrightarrow> inj_on (choose_pos A) A"
huffman@27411
   387
 apply (induct n arbitrary: A)
huffman@27411
   388
  apply simp
huffman@27411
   389
 apply (case_tac "A = {}", simp)
huffman@27411
   390
 apply (frule (1) choose_in)
huffman@27411
   391
 apply (rule inj_onI)
huffman@27411
   392
 apply (drule_tac x="A - {choose A}" in meta_spec, simp)
huffman@27411
   393
 apply (simp add: choose_pos.simps)
huffman@27411
   394
 apply (simp split: split_if_asm)
huffman@27411
   395
 apply (erule (1) inj_onD, simp, simp)
huffman@27411
   396
done
huffman@27411
   397
huffman@27411
   398
lemma choose_pos_bounded [OF refl]:
huffman@27411
   399
  "\<lbrakk>card A = n; finite A; x \<in> A\<rbrakk> \<Longrightarrow> choose_pos A x < n"
huffman@27411
   400
apply (induct n arbitrary: A)
huffman@27411
   401
apply simp
huffman@27411
   402
 apply (case_tac "A = {}", simp)
huffman@27411
   403
 apply (frule (1) choose_in)
huffman@27411
   404
apply (subst choose_pos.simps)
huffman@27411
   405
apply simp
huffman@27411
   406
done
huffman@27411
   407
huffman@27411
   408
lemma choose_pos_lessD:
huffman@27411
   409
  "\<lbrakk>choose_pos A x < choose_pos A y; finite A; x \<in> A; y \<in> A\<rbrakk> \<Longrightarrow> \<not> x \<sqsubseteq> y"
huffman@27411
   410
 apply (induct A x arbitrary: y rule: choose_pos.induct)
huffman@27411
   411
 apply simp
huffman@27411
   412
 apply (case_tac "x = choose A")
huffman@27411
   413
  apply simp
huffman@27411
   414
  apply (rule notI)
huffman@27411
   415
  apply (frule (2) maximal_choose)
huffman@27411
   416
  apply simp
huffman@27411
   417
 apply (case_tac "y = choose A")
huffman@27411
   418
  apply (simp add: choose_pos_choose)
huffman@27411
   419
 apply (drule_tac x=y in meta_spec)
huffman@27411
   420
 apply simp
huffman@27411
   421
 apply (erule meta_mp)
huffman@27411
   422
 apply (simp add: choose_pos.simps)
huffman@27411
   423
done
huffman@27411
   424
huffman@39974
   425
subsubsection {* Properties of approx function *}
huffman@39974
   426
huffman@39974
   427
lemma deflation_approx: "deflation (approx i)"
huffman@39974
   428
using finite_deflation_approx by (rule finite_deflation_imp_deflation)
huffman@39974
   429
huffman@39974
   430
lemma approx_idem: "approx i\<cdot>(approx i\<cdot>x) = approx i\<cdot>x"
huffman@39974
   431
using deflation_approx by (rule deflation.idem)
huffman@39974
   432
huffman@39974
   433
lemma approx_below: "approx i\<cdot>x \<sqsubseteq> x"
huffman@39974
   434
using deflation_approx by (rule deflation.below)
huffman@39974
   435
huffman@39974
   436
lemma finite_range_approx: "finite (range (\<lambda>x. approx i\<cdot>x))"
huffman@39974
   437
apply (rule finite_deflation.finite_range)
huffman@39974
   438
apply (rule finite_deflation_approx)
huffman@39974
   439
done
huffman@39974
   440
huffman@39974
   441
lemma compact_approx: "compact (approx n\<cdot>x)"
huffman@39974
   442
apply (rule finite_deflation.compact)
huffman@39974
   443
apply (rule finite_deflation_approx)
huffman@39974
   444
done
huffman@39974
   445
huffman@39974
   446
lemma compact_eq_approx: "compact x \<Longrightarrow> \<exists>i. approx i\<cdot>x = x"
huffman@39974
   447
by (rule admD2, simp_all)
huffman@39974
   448
huffman@39974
   449
subsubsection {* Compact basis take function *}
huffman@27411
   450
huffman@27411
   451
primrec
huffman@39974
   452
  cb_take :: "nat \<Rightarrow> 'a compact_basis \<Rightarrow> 'a compact_basis" where
huffman@27411
   453
  "cb_take 0 = (\<lambda>x. compact_bot)"
huffman@39974
   454
| "cb_take (Suc n) = (\<lambda>a. Abs_compact_basis (approx n\<cdot>(Rep_compact_basis a)))"
huffman@39974
   455
huffman@39974
   456
declare cb_take.simps [simp del]
huffman@39974
   457
huffman@39974
   458
lemma cb_take_zero [simp]: "cb_take 0 a = compact_bot"
huffman@39974
   459
by (simp only: cb_take.simps)
huffman@39974
   460
huffman@39974
   461
lemma Rep_cb_take:
huffman@39974
   462
  "Rep_compact_basis (cb_take (Suc n) a) = approx n\<cdot>(Rep_compact_basis a)"
huffman@39974
   463
by (simp add: Abs_compact_basis_inverse cb_take.simps(2) compact_approx)
huffman@39974
   464
huffman@39974
   465
lemmas approx_Rep_compact_basis = Rep_cb_take [symmetric]
huffman@27411
   466
huffman@27411
   467
lemma cb_take_covers: "\<exists>n. cb_take n x = x"
huffman@39974
   468
apply (subgoal_tac "\<exists>n. cb_take (Suc n) x = x", fast)
huffman@39974
   469
apply (simp add: Rep_compact_basis_inject [symmetric])
huffman@39974
   470
apply (simp add: Rep_cb_take)
huffman@39974
   471
apply (rule compact_eq_approx)
huffman@39974
   472
apply (rule compact_Rep_compact_basis)
huffman@27411
   473
done
huffman@27411
   474
huffman@27411
   475
lemma cb_take_less: "cb_take n x \<sqsubseteq> x"
huffman@39974
   476
unfolding compact_le_def
huffman@39974
   477
by (cases n, simp, simp add: Rep_cb_take approx_below)
huffman@27411
   478
huffman@27411
   479
lemma cb_take_idem: "cb_take n (cb_take n x) = cb_take n x"
huffman@39974
   480
unfolding Rep_compact_basis_inject [symmetric]
huffman@39974
   481
by (cases n, simp, simp add: Rep_cb_take approx_idem)
huffman@27411
   482
huffman@27411
   483
lemma cb_take_mono: "x \<sqsubseteq> y \<Longrightarrow> cb_take n x \<sqsubseteq> cb_take n y"
huffman@39974
   484
unfolding compact_le_def
huffman@39974
   485
by (cases n, simp, simp add: Rep_cb_take monofun_cfun_arg)
huffman@27411
   486
huffman@27411
   487
lemma cb_take_chain_le: "m \<le> n \<Longrightarrow> cb_take m x \<sqsubseteq> cb_take n x"
huffman@39974
   488
unfolding compact_le_def
huffman@39974
   489
apply (cases m, simp, cases n, simp)
huffman@39974
   490
apply (simp add: Rep_cb_take, rule chain_mono, simp, simp)
huffman@27411
   491
done
huffman@27411
   492
huffman@27411
   493
lemma finite_range_cb_take: "finite (range (cb_take n))"
huffman@27411
   494
apply (cases n)
huffman@39974
   495
apply (subgoal_tac "range (cb_take 0) = {compact_bot}", simp, force)
huffman@39974
   496
apply (rule finite_imageD [where f="Rep_compact_basis"])
huffman@39974
   497
apply (rule finite_subset [where B="range (\<lambda>x. approx (n - 1)\<cdot>x)"])
huffman@39974
   498
apply (clarsimp simp add: Rep_cb_take)
huffman@39974
   499
apply (rule finite_range_approx)
huffman@39974
   500
apply (rule inj_onI, simp add: Rep_compact_basis_inject)
huffman@27411
   501
done
huffman@27411
   502
huffman@39974
   503
subsubsection {* Rank of basis elements *}
huffman@39974
   504
huffman@27411
   505
definition
huffman@27411
   506
  rank :: "'a compact_basis \<Rightarrow> nat"
huffman@27411
   507
where
huffman@27411
   508
  "rank x = (LEAST n. cb_take n x = x)"
huffman@27411
   509
huffman@27411
   510
lemma compact_approx_rank: "cb_take (rank x) x = x"
huffman@27411
   511
unfolding rank_def
huffman@27411
   512
apply (rule LeastI_ex)
huffman@27411
   513
apply (rule cb_take_covers)
huffman@27411
   514
done
huffman@27411
   515
huffman@27411
   516
lemma rank_leD: "rank x \<le> n \<Longrightarrow> cb_take n x = x"
huffman@31076
   517
apply (rule below_antisym [OF cb_take_less])
huffman@27411
   518
apply (subst compact_approx_rank [symmetric])
huffman@27411
   519
apply (erule cb_take_chain_le)
huffman@27411
   520
done
huffman@27411
   521
huffman@27411
   522
lemma rank_leI: "cb_take n x = x \<Longrightarrow> rank x \<le> n"
huffman@27411
   523
unfolding rank_def by (rule Least_le)
huffman@27411
   524
huffman@27411
   525
lemma rank_le_iff: "rank x \<le> n \<longleftrightarrow> cb_take n x = x"
huffman@27411
   526
by (rule iffI [OF rank_leD rank_leI])
huffman@27411
   527
huffman@30505
   528
lemma rank_compact_bot [simp]: "rank compact_bot = 0"
huffman@30505
   529
using rank_leI [of 0 compact_bot] by simp
huffman@30505
   530
huffman@30505
   531
lemma rank_eq_0_iff [simp]: "rank x = 0 \<longleftrightarrow> x = compact_bot"
huffman@30505
   532
using rank_le_iff [of x 0] by auto
huffman@30505
   533
huffman@27411
   534
definition
huffman@27411
   535
  rank_le :: "'a compact_basis \<Rightarrow> 'a compact_basis set"
huffman@27411
   536
where
huffman@27411
   537
  "rank_le x = {y. rank y \<le> rank x}"
huffman@27411
   538
huffman@27411
   539
definition
huffman@27411
   540
  rank_lt :: "'a compact_basis \<Rightarrow> 'a compact_basis set"
huffman@27411
   541
where
huffman@27411
   542
  "rank_lt x = {y. rank y < rank x}"
huffman@27411
   543
huffman@27411
   544
definition
huffman@27411
   545
  rank_eq :: "'a compact_basis \<Rightarrow> 'a compact_basis set"
huffman@27411
   546
where
huffman@27411
   547
  "rank_eq x = {y. rank y = rank x}"
huffman@27411
   548
huffman@27411
   549
lemma rank_eq_cong: "rank x = rank y \<Longrightarrow> rank_eq x = rank_eq y"
huffman@27411
   550
unfolding rank_eq_def by simp
huffman@27411
   551
huffman@27411
   552
lemma rank_lt_cong: "rank x = rank y \<Longrightarrow> rank_lt x = rank_lt y"
huffman@27411
   553
unfolding rank_lt_def by simp
huffman@27411
   554
huffman@27411
   555
lemma rank_eq_subset: "rank_eq x \<subseteq> rank_le x"
huffman@27411
   556
unfolding rank_eq_def rank_le_def by auto
huffman@27411
   557
huffman@27411
   558
lemma rank_lt_subset: "rank_lt x \<subseteq> rank_le x"
huffman@27411
   559
unfolding rank_lt_def rank_le_def by auto
huffman@27411
   560
huffman@27411
   561
lemma finite_rank_le: "finite (rank_le x)"
huffman@27411
   562
unfolding rank_le_def
huffman@27411
   563
apply (rule finite_subset [where B="range (cb_take (rank x))"])
huffman@27411
   564
apply clarify
huffman@27411
   565
apply (rule range_eqI)
huffman@27411
   566
apply (erule rank_leD [symmetric])
huffman@27411
   567
apply (rule finite_range_cb_take)
huffman@27411
   568
done
huffman@27411
   569
huffman@27411
   570
lemma finite_rank_eq: "finite (rank_eq x)"
huffman@27411
   571
by (rule finite_subset [OF rank_eq_subset finite_rank_le])
huffman@27411
   572
huffman@27411
   573
lemma finite_rank_lt: "finite (rank_lt x)"
huffman@27411
   574
by (rule finite_subset [OF rank_lt_subset finite_rank_le])
huffman@27411
   575
huffman@27411
   576
lemma rank_lt_Int_rank_eq: "rank_lt x \<inter> rank_eq x = {}"
huffman@27411
   577
unfolding rank_lt_def rank_eq_def rank_le_def by auto
huffman@27411
   578
huffman@27411
   579
lemma rank_lt_Un_rank_eq: "rank_lt x \<union> rank_eq x = rank_le x"
huffman@27411
   580
unfolding rank_lt_def rank_eq_def rank_le_def by auto
huffman@27411
   581
huffman@30505
   582
subsubsection {* Sequencing basis elements *}
huffman@27411
   583
huffman@27411
   584
definition
huffman@30505
   585
  place :: "'a compact_basis \<Rightarrow> nat"
huffman@27411
   586
where
huffman@30505
   587
  "place x = card (rank_lt x) + choose_pos (rank_eq x) x"
huffman@27411
   588
huffman@30505
   589
lemma place_bounded: "place x < card (rank_le x)"
huffman@30505
   590
unfolding place_def
huffman@27411
   591
 apply (rule ord_less_eq_trans)
huffman@27411
   592
  apply (rule add_strict_left_mono)
huffman@27411
   593
  apply (rule choose_pos_bounded)
huffman@27411
   594
   apply (rule finite_rank_eq)
huffman@27411
   595
  apply (simp add: rank_eq_def)
huffman@27411
   596
 apply (subst card_Un_disjoint [symmetric])
huffman@27411
   597
    apply (rule finite_rank_lt)
huffman@27411
   598
   apply (rule finite_rank_eq)
huffman@27411
   599
  apply (rule rank_lt_Int_rank_eq)
huffman@27411
   600
 apply (simp add: rank_lt_Un_rank_eq)
huffman@27411
   601
done
huffman@27411
   602
huffman@30505
   603
lemma place_ge: "card (rank_lt x) \<le> place x"
huffman@30505
   604
unfolding place_def by simp
huffman@27411
   605
huffman@30505
   606
lemma place_rank_mono:
huffman@27411
   607
  fixes x y :: "'a compact_basis"
huffman@30505
   608
  shows "rank x < rank y \<Longrightarrow> place x < place y"
huffman@30505
   609
apply (rule less_le_trans [OF place_bounded])
huffman@30505
   610
apply (rule order_trans [OF _ place_ge])
huffman@27411
   611
apply (rule card_mono)
huffman@27411
   612
apply (rule finite_rank_lt)
huffman@27411
   613
apply (simp add: rank_le_def rank_lt_def subset_eq)
huffman@27411
   614
done
huffman@27411
   615
huffman@30505
   616
lemma place_eqD: "place x = place y \<Longrightarrow> x = y"
huffman@27411
   617
 apply (rule linorder_cases [where x="rank x" and y="rank y"])
huffman@30505
   618
   apply (drule place_rank_mono, simp)
huffman@30505
   619
  apply (simp add: place_def)
huffman@27411
   620
  apply (rule inj_on_choose_pos [where A="rank_eq x", THEN inj_onD])
huffman@27411
   621
     apply (rule finite_rank_eq)
huffman@27411
   622
    apply (simp cong: rank_lt_cong rank_eq_cong)
huffman@27411
   623
   apply (simp add: rank_eq_def)
huffman@27411
   624
  apply (simp add: rank_eq_def)
huffman@30505
   625
 apply (drule place_rank_mono, simp)
huffman@27411
   626
done
huffman@27411
   627
huffman@30505
   628
lemma inj_place: "inj place"
huffman@30505
   629
by (rule inj_onI, erule place_eqD)
huffman@27411
   630
huffman@27411
   631
subsubsection {* Embedding and projection on basis elements *}
huffman@27411
   632
huffman@30505
   633
definition
huffman@30505
   634
  sub :: "'a compact_basis \<Rightarrow> 'a compact_basis"
huffman@30505
   635
where
huffman@30505
   636
  "sub x = (case rank x of 0 \<Rightarrow> compact_bot | Suc k \<Rightarrow> cb_take k x)"
huffman@30505
   637
huffman@30505
   638
lemma rank_sub_less: "x \<noteq> compact_bot \<Longrightarrow> rank (sub x) < rank x"
huffman@30505
   639
unfolding sub_def
huffman@30505
   640
apply (cases "rank x", simp)
huffman@30505
   641
apply (simp add: less_Suc_eq_le)
huffman@30505
   642
apply (rule rank_leI)
huffman@30505
   643
apply (rule cb_take_idem)
huffman@30505
   644
done
huffman@30505
   645
huffman@30505
   646
lemma place_sub_less: "x \<noteq> compact_bot \<Longrightarrow> place (sub x) < place x"
huffman@30505
   647
apply (rule place_rank_mono)
huffman@30505
   648
apply (erule rank_sub_less)
huffman@30505
   649
done
huffman@30505
   650
huffman@30505
   651
lemma sub_below: "sub x \<sqsubseteq> x"
huffman@30505
   652
unfolding sub_def by (cases "rank x", simp_all add: cb_take_less)
huffman@30505
   653
huffman@30505
   654
lemma rank_less_imp_below_sub: "\<lbrakk>x \<sqsubseteq> y; rank x < rank y\<rbrakk> \<Longrightarrow> x \<sqsubseteq> sub y"
huffman@30505
   655
unfolding sub_def
huffman@30505
   656
apply (cases "rank y", simp)
huffman@30505
   657
apply (simp add: less_Suc_eq_le)
huffman@30505
   658
apply (subgoal_tac "cb_take nat x \<sqsubseteq> cb_take nat y")
huffman@30505
   659
apply (simp add: rank_leD)
huffman@30505
   660
apply (erule cb_take_mono)
huffman@30505
   661
done
huffman@30505
   662
huffman@27411
   663
function
huffman@27411
   664
  basis_emb :: "'a compact_basis \<Rightarrow> ubasis"
huffman@27411
   665
where
huffman@27411
   666
  "basis_emb x = (if x = compact_bot then 0 else
huffman@30505
   667
    node (place x) (basis_emb (sub x))
huffman@30505
   668
      (basis_emb ` {y. place y < place x \<and> x \<sqsubseteq> y}))"
huffman@27411
   669
by auto
huffman@27411
   670
huffman@27411
   671
termination basis_emb
huffman@30505
   672
apply (relation "measure place", simp)
huffman@30505
   673
apply (simp add: place_sub_less)
huffman@27411
   674
apply simp
huffman@27411
   675
done
huffman@27411
   676
huffman@27411
   677
declare basis_emb.simps [simp del]
huffman@27411
   678
huffman@27411
   679
lemma basis_emb_compact_bot [simp]: "basis_emb compact_bot = 0"
huffman@27411
   680
by (simp add: basis_emb.simps)
huffman@27411
   681
huffman@30505
   682
lemma fin1: "finite {y. place y < place x \<and> x \<sqsubseteq> y}"
huffman@27411
   683
apply (subst Collect_conj_eq)
huffman@27411
   684
apply (rule finite_Int)
huffman@27411
   685
apply (rule disjI1)
huffman@30505
   686
apply (subgoal_tac "finite (place -` {n. n < place x})", simp)
huffman@30505
   687
apply (rule finite_vimageI [OF _ inj_place])
huffman@27411
   688
apply (simp add: lessThan_def [symmetric])
huffman@27411
   689
done
huffman@27411
   690
huffman@30505
   691
lemma fin2: "finite (basis_emb ` {y. place y < place x \<and> x \<sqsubseteq> y})"
huffman@27411
   692
by (rule finite_imageI [OF fin1])
huffman@27411
   693
huffman@30505
   694
lemma rank_place_mono:
huffman@30505
   695
  "\<lbrakk>place x < place y; x \<sqsubseteq> y\<rbrakk> \<Longrightarrow> rank x < rank y"
huffman@30505
   696
apply (rule linorder_cases, assumption)
huffman@30505
   697
apply (simp add: place_def cong: rank_lt_cong rank_eq_cong)
huffman@30505
   698
apply (drule choose_pos_lessD)
huffman@30505
   699
apply (rule finite_rank_eq)
huffman@30505
   700
apply (simp add: rank_eq_def)
huffman@30505
   701
apply (simp add: rank_eq_def)
huffman@30505
   702
apply simp
huffman@30505
   703
apply (drule place_rank_mono, simp)
huffman@30505
   704
done
huffman@30505
   705
huffman@30505
   706
lemma basis_emb_mono:
huffman@30505
   707
  "x \<sqsubseteq> y \<Longrightarrow> ubasis_le (basis_emb x) (basis_emb y)"
berghofe@34915
   708
proof (induct "max (place x) (place y)" arbitrary: x y rule: less_induct)
berghofe@34915
   709
  case less
huffman@30505
   710
  show ?case proof (rule linorder_cases)
huffman@30505
   711
    assume "place x < place y"
huffman@30505
   712
    then have "rank x < rank y"
huffman@30505
   713
      using `x \<sqsubseteq> y` by (rule rank_place_mono)
huffman@30505
   714
    with `place x < place y` show ?case
huffman@30505
   715
      apply (case_tac "y = compact_bot", simp)
huffman@30505
   716
      apply (simp add: basis_emb.simps [of y])
huffman@30505
   717
      apply (rule ubasis_le_trans [OF _ ubasis_le_lower [OF fin2]])
berghofe@34915
   718
      apply (rule less)
huffman@30505
   719
       apply (simp add: less_max_iff_disj)
huffman@30505
   720
       apply (erule place_sub_less)
huffman@30505
   721
      apply (erule rank_less_imp_below_sub [OF `x \<sqsubseteq> y`])
huffman@27411
   722
      done
huffman@30505
   723
  next
huffman@30505
   724
    assume "place x = place y"
huffman@30505
   725
    hence "x = y" by (rule place_eqD)
huffman@30505
   726
    thus ?case by (simp add: ubasis_le_refl)
huffman@30505
   727
  next
huffman@30505
   728
    assume "place x > place y"
huffman@30505
   729
    with `x \<sqsubseteq> y` show ?case
huffman@30505
   730
      apply (case_tac "x = compact_bot", simp add: ubasis_le_minimal)
huffman@30505
   731
      apply (simp add: basis_emb.simps [of x])
huffman@30505
   732
      apply (rule ubasis_le_upper [OF fin2], simp)
berghofe@34915
   733
      apply (rule less)
huffman@30505
   734
       apply (simp add: less_max_iff_disj)
huffman@30505
   735
       apply (erule place_sub_less)
huffman@31076
   736
      apply (erule rev_below_trans)
huffman@30505
   737
      apply (rule sub_below)
huffman@30505
   738
      done
huffman@27411
   739
  qed
huffman@27411
   740
qed
huffman@27411
   741
huffman@27411
   742
lemma inj_basis_emb: "inj basis_emb"
huffman@27411
   743
 apply (rule inj_onI)
huffman@27411
   744
 apply (case_tac "x = compact_bot")
huffman@27411
   745
  apply (case_tac [!] "y = compact_bot")
huffman@27411
   746
    apply simp
huffman@27411
   747
   apply (simp add: basis_emb.simps)
huffman@27411
   748
  apply (simp add: basis_emb.simps)
huffman@27411
   749
 apply (simp add: basis_emb.simps)
huffman@30505
   750
 apply (simp add: fin2 inj_eq [OF inj_place])
huffman@27411
   751
done
huffman@27411
   752
huffman@27411
   753
definition
huffman@30505
   754
  basis_prj :: "ubasis \<Rightarrow> 'a compact_basis"
huffman@27411
   755
where
huffman@27411
   756
  "basis_prj x = inv basis_emb
huffman@30505
   757
    (ubasis_until (\<lambda>x. x \<in> range (basis_emb :: 'a compact_basis \<Rightarrow> ubasis)) x)"
huffman@27411
   758
huffman@27411
   759
lemma basis_prj_basis_emb: "\<And>x. basis_prj (basis_emb x) = x"
huffman@27411
   760
unfolding basis_prj_def
huffman@27411
   761
 apply (subst ubasis_until_same)
huffman@27411
   762
  apply (rule rangeI)
huffman@27411
   763
 apply (rule inv_f_f)
huffman@27411
   764
 apply (rule inj_basis_emb)
huffman@27411
   765
done
huffman@27411
   766
huffman@27411
   767
lemma basis_prj_node:
huffman@30505
   768
  "\<lbrakk>finite S; node i a S \<notin> range (basis_emb :: 'a compact_basis \<Rightarrow> nat)\<rbrakk>
huffman@30505
   769
    \<Longrightarrow> basis_prj (node i a S) = (basis_prj a :: 'a compact_basis)"
huffman@27411
   770
unfolding basis_prj_def by simp
huffman@27411
   771
huffman@27411
   772
lemma basis_prj_0: "basis_prj 0 = compact_bot"
huffman@27411
   773
apply (subst basis_emb_compact_bot [symmetric])
huffman@27411
   774
apply (rule basis_prj_basis_emb)
huffman@27411
   775
done
huffman@27411
   776
huffman@30505
   777
lemma node_eq_basis_emb_iff:
huffman@30505
   778
  "finite S \<Longrightarrow> node i a S = basis_emb x \<longleftrightarrow>
huffman@30505
   779
    x \<noteq> compact_bot \<and> i = place x \<and> a = basis_emb (sub x) \<and>
huffman@30505
   780
        S = basis_emb ` {y. place y < place x \<and> x \<sqsubseteq> y}"
huffman@30505
   781
apply (cases "x = compact_bot", simp)
huffman@30505
   782
apply (simp add: basis_emb.simps [of x])
huffman@30505
   783
apply (simp add: fin2)
huffman@27411
   784
done
huffman@27411
   785
huffman@30505
   786
lemma basis_prj_mono: "ubasis_le a b \<Longrightarrow> basis_prj a \<sqsubseteq> basis_prj b"
huffman@30505
   787
proof (induct a b rule: ubasis_le.induct)
huffman@31076
   788
  case (ubasis_le_refl a) show ?case by (rule below_refl)
huffman@30505
   789
next
huffman@31076
   790
  case (ubasis_le_trans a b c) thus ?case by - (rule below_trans)
huffman@30505
   791
next
huffman@30505
   792
  case (ubasis_le_lower S a i) thus ?case
huffman@30561
   793
    apply (cases "node i a S \<in> range (basis_emb :: 'a compact_basis \<Rightarrow> nat)")
huffman@30505
   794
     apply (erule rangeE, rename_tac x)
huffman@30505
   795
     apply (simp add: basis_prj_basis_emb)
huffman@30505
   796
     apply (simp add: node_eq_basis_emb_iff)
huffman@30505
   797
     apply (simp add: basis_prj_basis_emb)
huffman@30505
   798
     apply (rule sub_below)
huffman@30505
   799
    apply (simp add: basis_prj_node)
huffman@30505
   800
    done
huffman@30505
   801
next
huffman@30505
   802
  case (ubasis_le_upper S b a i) thus ?case
huffman@30561
   803
    apply (cases "node i a S \<in> range (basis_emb :: 'a compact_basis \<Rightarrow> nat)")
huffman@30505
   804
     apply (erule rangeE, rename_tac x)
huffman@30505
   805
     apply (simp add: basis_prj_basis_emb)
huffman@30505
   806
     apply (clarsimp simp add: node_eq_basis_emb_iff)
huffman@30505
   807
     apply (simp add: basis_prj_basis_emb)
huffman@30505
   808
    apply (simp add: basis_prj_node)
huffman@30505
   809
    done
huffman@30505
   810
qed
huffman@30505
   811
huffman@27411
   812
lemma basis_emb_prj_less: "ubasis_le (basis_emb (basis_prj x)) x"
huffman@27411
   813
unfolding basis_prj_def
wenzelm@33071
   814
 apply (subst f_inv_into_f [where f=basis_emb])
huffman@27411
   815
  apply (rule ubasis_until)
huffman@27411
   816
  apply (rule range_eqI [where x=compact_bot])
huffman@27411
   817
  apply simp
huffman@27411
   818
 apply (rule ubasis_until_less)
huffman@27411
   819
done
huffman@27411
   820
huffman@39974
   821
end
huffman@39974
   822
huffman@39974
   823
sublocale approx_chain \<subseteq> compact_basis!:
huffman@39974
   824
  ideal_completion below Rep_compact_basis
huffman@39974
   825
    "approximants :: 'a \<Rightarrow> 'a compact_basis set"
huffman@39974
   826
proof
huffman@39974
   827
  fix w :: "'a"
huffman@39974
   828
  show "below.ideal (approximants w)"
huffman@39974
   829
  proof (rule below.idealI)
huffman@39974
   830
    show "\<exists>x. x \<in> approximants w"
huffman@39974
   831
      unfolding approximants_def
huffman@39974
   832
      apply (rule_tac x="Abs_compact_basis (approx 0\<cdot>w)" in exI)
huffman@39974
   833
      apply (simp add: Abs_compact_basis_inverse approx_below compact_approx)
huffman@39974
   834
      done
huffman@39974
   835
  next
huffman@39974
   836
    fix x y :: "'a compact_basis"
huffman@39974
   837
    assume "x \<in> approximants w" "y \<in> approximants w"
huffman@39974
   838
    thus "\<exists>z \<in> approximants w. x \<sqsubseteq> z \<and> y \<sqsubseteq> z"
huffman@39974
   839
      unfolding approximants_def
huffman@39974
   840
      apply simp
huffman@39974
   841
      apply (cut_tac a=x in compact_Rep_compact_basis)
huffman@39974
   842
      apply (cut_tac a=y in compact_Rep_compact_basis)
huffman@39974
   843
      apply (drule compact_eq_approx)
huffman@39974
   844
      apply (drule compact_eq_approx)
huffman@39974
   845
      apply (clarify, rename_tac i j)
huffman@39974
   846
      apply (rule_tac x="Abs_compact_basis (approx (max i j)\<cdot>w)" in exI)
huffman@39974
   847
      apply (simp add: compact_le_def)
huffman@39974
   848
      apply (simp add: Abs_compact_basis_inverse approx_below compact_approx)
huffman@39974
   849
      apply (erule subst, erule subst)
huffman@39974
   850
      apply (simp add: monofun_cfun chain_mono [OF chain_approx])
huffman@39974
   851
      done
huffman@39974
   852
  next
huffman@39974
   853
    fix x y :: "'a compact_basis"
huffman@39974
   854
    assume "x \<sqsubseteq> y" "y \<in> approximants w" thus "x \<in> approximants w"
huffman@39974
   855
      unfolding approximants_def
huffman@39974
   856
      apply simp
huffman@39974
   857
      apply (simp add: compact_le_def)
huffman@39974
   858
      apply (erule (1) below_trans)
huffman@39974
   859
      done
huffman@39974
   860
  qed
huffman@39974
   861
next
huffman@39974
   862
  fix Y :: "nat \<Rightarrow> 'a"
huffman@39974
   863
  assume Y: "chain Y"
huffman@39974
   864
  show "approximants (\<Squnion>i. Y i) = (\<Union>i. approximants (Y i))"
huffman@39974
   865
    unfolding approximants_def
huffman@39974
   866
    apply safe
huffman@39974
   867
    apply (simp add: compactD2 [OF compact_Rep_compact_basis Y])
huffman@39974
   868
    apply (erule below_trans, rule is_ub_thelub [OF Y])
huffman@39974
   869
    done
huffman@39974
   870
next
huffman@39974
   871
  fix a :: "'a compact_basis"
huffman@39974
   872
  show "approximants (Rep_compact_basis a) = {b. b \<sqsubseteq> a}"
huffman@39974
   873
    unfolding approximants_def compact_le_def ..
huffman@39974
   874
next
huffman@39974
   875
  fix x y :: "'a"
huffman@39974
   876
  assume "approximants x \<subseteq> approximants y" thus "x \<sqsubseteq> y"
huffman@39974
   877
    apply (subgoal_tac "(\<Squnion>i. approx i\<cdot>x) \<sqsubseteq> y")
huffman@39974
   878
    apply (simp add: lub_distribs)
huffman@39974
   879
    apply (rule admD, simp, simp)
huffman@39974
   880
    apply (drule_tac c="Abs_compact_basis (approx i\<cdot>x)" in subsetD)
huffman@39974
   881
    apply (simp add: approximants_def Abs_compact_basis_inverse
huffman@39974
   882
                     approx_below compact_approx)
huffman@39974
   883
    apply (simp add: approximants_def Abs_compact_basis_inverse compact_approx)
huffman@39974
   884
    done
huffman@39974
   885
next
huffman@39974
   886
  show "\<exists>f::'a compact_basis \<Rightarrow> nat. inj f"
huffman@39974
   887
    by (rule exI, rule inj_place)
huffman@39974
   888
qed
huffman@27411
   889
huffman@35900
   890
subsubsection {* EP-pair from any bifinite domain into \emph{udom} *}
huffman@27411
   891
huffman@39974
   892
context approx_chain begin
huffman@39974
   893
huffman@27411
   894
definition
huffman@39974
   895
  udom_emb :: "'a \<rightarrow> udom"
huffman@27411
   896
where
huffman@27411
   897
  "udom_emb = compact_basis.basis_fun (\<lambda>x. udom_principal (basis_emb x))"
huffman@27411
   898
huffman@27411
   899
definition
huffman@39974
   900
  udom_prj :: "udom \<rightarrow> 'a"
huffman@27411
   901
where
huffman@27411
   902
  "udom_prj = udom.basis_fun (\<lambda>x. Rep_compact_basis (basis_prj x))"
huffman@27411
   903
huffman@27411
   904
lemma udom_emb_principal:
huffman@27411
   905
  "udom_emb\<cdot>(Rep_compact_basis x) = udom_principal (basis_emb x)"
huffman@27411
   906
unfolding udom_emb_def
huffman@27411
   907
apply (rule compact_basis.basis_fun_principal)
huffman@27411
   908
apply (rule udom.principal_mono)
huffman@27411
   909
apply (erule basis_emb_mono)
huffman@27411
   910
done
huffman@27411
   911
huffman@27411
   912
lemma udom_prj_principal:
huffman@27411
   913
  "udom_prj\<cdot>(udom_principal x) = Rep_compact_basis (basis_prj x)"
huffman@27411
   914
unfolding udom_prj_def
huffman@27411
   915
apply (rule udom.basis_fun_principal)
huffman@27411
   916
apply (rule compact_basis.principal_mono)
huffman@27411
   917
apply (erule basis_prj_mono)
huffman@27411
   918
done
huffman@27411
   919
huffman@27411
   920
lemma ep_pair_udom: "ep_pair udom_emb udom_prj"
huffman@27411
   921
 apply default
huffman@27411
   922
  apply (rule compact_basis.principal_induct, simp)
huffman@27411
   923
  apply (simp add: udom_emb_principal udom_prj_principal)
huffman@27411
   924
  apply (simp add: basis_prj_basis_emb)
huffman@27411
   925
 apply (rule udom.principal_induct, simp)
huffman@27411
   926
 apply (simp add: udom_emb_principal udom_prj_principal)
huffman@27411
   927
 apply (rule basis_emb_prj_less)
huffman@27411
   928
done
huffman@27411
   929
huffman@27411
   930
end
huffman@39974
   931
huffman@39974
   932
abbreviation "udom_emb \<equiv> approx_chain.udom_emb"
huffman@39974
   933
abbreviation "udom_prj \<equiv> approx_chain.udom_prj"
huffman@39974
   934
huffman@39974
   935
lemmas ep_pair_udom = approx_chain.ep_pair_udom
huffman@39974
   936
huffman@39974
   937
subsection {* Chain of approx functions for type \emph{udom} *}
huffman@39974
   938
huffman@39974
   939
definition
huffman@39974
   940
  udom_approx :: "nat \<Rightarrow> udom \<rightarrow> udom"
huffman@39974
   941
where
huffman@39974
   942
  "udom_approx i =
huffman@39974
   943
    udom.basis_fun (\<lambda>x. udom_principal (ubasis_until (\<lambda>y. y \<le> i) x))"
huffman@39974
   944
huffman@39974
   945
lemma udom_approx_mono:
huffman@39974
   946
  "ubasis_le a b \<Longrightarrow>
huffman@39974
   947
    udom_principal (ubasis_until (\<lambda>y. y \<le> i) a) \<sqsubseteq>
huffman@39974
   948
    udom_principal (ubasis_until (\<lambda>y. y \<le> i) b)"
huffman@39974
   949
apply (rule udom.principal_mono)
huffman@39974
   950
apply (rule ubasis_until_mono)
huffman@39974
   951
apply (frule (2) order_less_le_trans [OF node_gt2])
huffman@39974
   952
apply (erule order_less_imp_le)
huffman@39974
   953
apply assumption
huffman@39974
   954
done
huffman@39974
   955
huffman@39974
   956
lemma adm_mem_finite: "\<lbrakk>cont f; finite S\<rbrakk> \<Longrightarrow> adm (\<lambda>x. f x \<in> S)"
huffman@39974
   957
by (erule adm_subst, induct set: finite, simp_all)
huffman@39974
   958
huffman@39974
   959
lemma udom_approx_principal:
huffman@39974
   960
  "udom_approx i\<cdot>(udom_principal x) =
huffman@39974
   961
    udom_principal (ubasis_until (\<lambda>y. y \<le> i) x)"
huffman@39974
   962
unfolding udom_approx_def
huffman@39974
   963
apply (rule udom.basis_fun_principal)
huffman@39974
   964
apply (erule udom_approx_mono)
huffman@39974
   965
done
huffman@39974
   966
huffman@39974
   967
lemma finite_deflation_udom_approx: "finite_deflation (udom_approx i)"
huffman@39974
   968
proof
huffman@39974
   969
  fix x show "udom_approx i\<cdot>(udom_approx i\<cdot>x) = udom_approx i\<cdot>x"
huffman@39974
   970
    by (induct x rule: udom.principal_induct, simp)
huffman@39974
   971
       (simp add: udom_approx_principal ubasis_until_idem)
huffman@39974
   972
next
huffman@39974
   973
  fix x show "udom_approx i\<cdot>x \<sqsubseteq> x"
huffman@39974
   974
    by (induct x rule: udom.principal_induct, simp)
huffman@39974
   975
       (simp add: udom_approx_principal ubasis_until_less)
huffman@39974
   976
next
huffman@39974
   977
  have *: "finite (range (\<lambda>x. udom_principal (ubasis_until (\<lambda>y. y \<le> i) x)))"
huffman@39974
   978
    apply (subst range_composition [where f=udom_principal])
huffman@39974
   979
    apply (simp add: finite_range_ubasis_until)
huffman@39974
   980
    done
huffman@39974
   981
  show "finite {x. udom_approx i\<cdot>x = x}"
huffman@39974
   982
    apply (rule finite_range_imp_finite_fixes)
huffman@39974
   983
    apply (rule rev_finite_subset [OF *])
huffman@39974
   984
    apply (clarsimp, rename_tac x)
huffman@39974
   985
    apply (induct_tac x rule: udom.principal_induct)
huffman@39974
   986
    apply (simp add: adm_mem_finite *)
huffman@39974
   987
    apply (simp add: udom_approx_principal)
huffman@39974
   988
    done
huffman@39974
   989
qed
huffman@39974
   990
huffman@39974
   991
interpretation udom_approx: finite_deflation "udom_approx i"
huffman@39974
   992
by (rule finite_deflation_udom_approx)
huffman@39974
   993
huffman@39974
   994
lemma chain_udom_approx [simp]: "chain (\<lambda>i. udom_approx i)"
huffman@39974
   995
unfolding udom_approx_def
huffman@39974
   996
apply (rule chainI)
huffman@39974
   997
apply (rule udom.basis_fun_mono)
huffman@39974
   998
apply (erule udom_approx_mono)
huffman@39974
   999
apply (erule udom_approx_mono)
huffman@39974
  1000
apply (rule udom.principal_mono)
huffman@39974
  1001
apply (rule ubasis_until_chain, simp)
huffman@39974
  1002
done
huffman@39974
  1003
huffman@39974
  1004
lemma lub_udom_approx [simp]: "(\<Squnion>i. udom_approx i) = ID"
huffman@39974
  1005
apply (rule ext_cfun, simp add: contlub_cfun_fun)
huffman@39974
  1006
apply (rule below_antisym)
huffman@39974
  1007
apply (rule is_lub_thelub)
huffman@39974
  1008
apply (simp)
huffman@39974
  1009
apply (rule ub_rangeI)
huffman@39974
  1010
apply (rule udom_approx.below)
huffman@39974
  1011
apply (rule_tac x=x in udom.principal_induct)
huffman@39974
  1012
apply (simp add: lub_distribs)
huffman@39974
  1013
apply (rule rev_below_trans)
huffman@39974
  1014
apply (rule_tac x=a in is_ub_thelub)
huffman@39974
  1015
apply simp
huffman@39974
  1016
apply (simp add: udom_approx_principal)
huffman@39974
  1017
apply (simp add: ubasis_until_same ubasis_le_refl)
huffman@39974
  1018
done
huffman@39974
  1019
 
huffman@39974
  1020
lemma udom_approx: "approx_chain udom_approx"
huffman@39974
  1021
proof
huffman@39974
  1022
  show "chain (\<lambda>i. udom_approx i)"
huffman@39974
  1023
    by (rule chain_udom_approx)
huffman@39974
  1024
  show "(\<Squnion>i. udom_approx i) = ID"
huffman@39974
  1025
    by (rule lub_udom_approx)
huffman@39974
  1026
qed
huffman@39974
  1027
huffman@39974
  1028
hide_const (open) node
huffman@39974
  1029
huffman@39974
  1030
end