src/HOL/Partial_Function.thy
author krauss
Tue May 31 11:11:17 2011 +0200 (2011-05-31)
changeset 43081 1a39c9898ae6
parent 43080 73a1d6a7ef1d
child 43082 8d0c44de9773
permissions -rw-r--r--
admissibility on option type
krauss@40107
     1
(* Title:    HOL/Partial_Function.thy
krauss@40107
     2
   Author:   Alexander Krauss, TU Muenchen
krauss@40107
     3
*)
krauss@40107
     4
krauss@40107
     5
header {* Partial Function Definitions *}
krauss@40107
     6
krauss@40107
     7
theory Partial_Function
krauss@40107
     8
imports Complete_Partial_Order Option
krauss@40107
     9
uses 
krauss@40107
    10
  "Tools/Function/function_lib.ML" 
krauss@40107
    11
  "Tools/Function/partial_function.ML" 
krauss@40107
    12
begin
krauss@40107
    13
krauss@40107
    14
setup Partial_Function.setup
krauss@40107
    15
krauss@40107
    16
subsection {* Axiomatic setup *}
krauss@40107
    17
krauss@40107
    18
text {* This techical locale constains the requirements for function
krauss@42949
    19
  definitions with ccpo fixed points. *}
krauss@40107
    20
krauss@40107
    21
definition "fun_ord ord f g \<longleftrightarrow> (\<forall>x. ord (f x) (g x))"
krauss@40107
    22
definition "fun_lub L A = (\<lambda>x. L {y. \<exists>f\<in>A. y = f x})"
krauss@40107
    23
definition "img_ord f ord = (\<lambda>x y. ord (f x) (f y))"
krauss@40107
    24
definition "img_lub f g Lub = (\<lambda>A. g (Lub (f ` A)))"
krauss@40107
    25
krauss@43081
    26
lemma chain_fun: 
krauss@43081
    27
  assumes A: "chain (fun_ord ord) A"
krauss@43081
    28
  shows "chain ord {y. \<exists>f\<in>A. y = f a}" (is "chain ord ?C")
krauss@43081
    29
proof (rule chainI)
krauss@43081
    30
  fix x y assume "x \<in> ?C" "y \<in> ?C"
krauss@43081
    31
  then obtain f g where fg: "f \<in> A" "g \<in> A" 
krauss@43081
    32
    and [simp]: "x = f a" "y = g a" by blast
krauss@43081
    33
  from chainD[OF A fg]
krauss@43081
    34
  show "ord x y \<or> ord y x" unfolding fun_ord_def by auto
krauss@43081
    35
qed
krauss@43081
    36
krauss@40107
    37
lemma call_mono[partial_function_mono]: "monotone (fun_ord ord) ord (\<lambda>f. f t)"
krauss@40107
    38
by (rule monotoneI) (auto simp: fun_ord_def)
krauss@40107
    39
krauss@40288
    40
lemma let_mono[partial_function_mono]:
krauss@40288
    41
  "(\<And>x. monotone orda ordb (\<lambda>f. b f x))
krauss@40288
    42
  \<Longrightarrow> monotone orda ordb (\<lambda>f. Let t (b f))"
krauss@40288
    43
by (simp add: Let_def)
krauss@40288
    44
krauss@40107
    45
lemma if_mono[partial_function_mono]: "monotone orda ordb F 
krauss@40107
    46
\<Longrightarrow> monotone orda ordb G
krauss@40107
    47
\<Longrightarrow> monotone orda ordb (\<lambda>f. if c then F f else G f)"
krauss@40107
    48
unfolding monotone_def by simp
krauss@40107
    49
krauss@40107
    50
definition "mk_less R = (\<lambda>x y. R x y \<and> \<not> R y x)"
krauss@40107
    51
krauss@40107
    52
locale partial_function_definitions = 
krauss@40107
    53
  fixes leq :: "'a \<Rightarrow> 'a \<Rightarrow> bool"
krauss@40107
    54
  fixes lub :: "'a set \<Rightarrow> 'a"
krauss@40107
    55
  assumes leq_refl: "leq x x"
krauss@40107
    56
  assumes leq_trans: "leq x y \<Longrightarrow> leq y z \<Longrightarrow> leq x z"
krauss@40107
    57
  assumes leq_antisym: "leq x y \<Longrightarrow> leq y x \<Longrightarrow> x = y"
krauss@40107
    58
  assumes lub_upper: "chain leq A \<Longrightarrow> x \<in> A \<Longrightarrow> leq x (lub A)"
krauss@40107
    59
  assumes lub_least: "chain leq A \<Longrightarrow> (\<And>x. x \<in> A \<Longrightarrow> leq x z) \<Longrightarrow> leq (lub A) z"
krauss@40107
    60
krauss@40107
    61
lemma partial_function_lift:
krauss@40107
    62
  assumes "partial_function_definitions ord lb"
krauss@40107
    63
  shows "partial_function_definitions (fun_ord ord) (fun_lub lb)" (is "partial_function_definitions ?ordf ?lubf")
krauss@40107
    64
proof -
krauss@40107
    65
  interpret partial_function_definitions ord lb by fact
krauss@40107
    66
krauss@40107
    67
  show ?thesis
krauss@40107
    68
  proof
krauss@40107
    69
    fix x show "?ordf x x"
krauss@40107
    70
      unfolding fun_ord_def by (auto simp: leq_refl)
krauss@40107
    71
  next
krauss@40107
    72
    fix x y z assume "?ordf x y" "?ordf y z"
krauss@40107
    73
    thus "?ordf x z" unfolding fun_ord_def 
krauss@40107
    74
      by (force dest: leq_trans)
krauss@40107
    75
  next
krauss@40107
    76
    fix x y assume "?ordf x y" "?ordf y x"
krauss@40107
    77
    thus "x = y" unfolding fun_ord_def
krauss@43081
    78
      by (force intro!: dest: leq_antisym)
krauss@40107
    79
  next
krauss@40107
    80
    fix A f assume f: "f \<in> A" and A: "chain ?ordf A"
krauss@40107
    81
    thus "?ordf f (?lubf A)"
krauss@40107
    82
      unfolding fun_lub_def fun_ord_def
krauss@40107
    83
      by (blast intro: lub_upper chain_fun[OF A] f)
krauss@40107
    84
  next
krauss@40107
    85
    fix A :: "('b \<Rightarrow> 'a) set" and g :: "'b \<Rightarrow> 'a"
krauss@40107
    86
    assume A: "chain ?ordf A" and g: "\<And>f. f \<in> A \<Longrightarrow> ?ordf f g"
krauss@40107
    87
    show "?ordf (?lubf A) g" unfolding fun_lub_def fun_ord_def
krauss@40107
    88
      by (blast intro: lub_least chain_fun[OF A] dest: g[unfolded fun_ord_def])
krauss@40107
    89
   qed
krauss@40107
    90
qed
krauss@40107
    91
krauss@40107
    92
lemma ccpo: assumes "partial_function_definitions ord lb"
krauss@40107
    93
  shows "class.ccpo ord (mk_less ord) lb"
krauss@40107
    94
using assms unfolding partial_function_definitions_def mk_less_def
krauss@40107
    95
by unfold_locales blast+
krauss@40107
    96
krauss@40107
    97
lemma partial_function_image:
krauss@40107
    98
  assumes "partial_function_definitions ord Lub"
krauss@40107
    99
  assumes inj: "\<And>x y. f x = f y \<Longrightarrow> x = y"
krauss@40107
   100
  assumes inv: "\<And>x. f (g x) = x"
krauss@40107
   101
  shows "partial_function_definitions (img_ord f ord) (img_lub f g Lub)"
krauss@40107
   102
proof -
krauss@40107
   103
  let ?iord = "img_ord f ord"
krauss@40107
   104
  let ?ilub = "img_lub f g Lub"
krauss@40107
   105
krauss@40107
   106
  interpret partial_function_definitions ord Lub by fact
krauss@40107
   107
  show ?thesis
krauss@40107
   108
  proof
krauss@40107
   109
    fix A x assume "chain ?iord A" "x \<in> A"
krauss@40107
   110
    then have "chain ord (f ` A)" "f x \<in> f ` A"
krauss@40107
   111
      by (auto simp: img_ord_def intro: chainI dest: chainD)
krauss@40107
   112
    thus "?iord x (?ilub A)"
krauss@40107
   113
      unfolding inv img_lub_def img_ord_def by (rule lub_upper)
krauss@40107
   114
  next
krauss@40107
   115
    fix A x assume "chain ?iord A"
krauss@40107
   116
      and 1: "\<And>z. z \<in> A \<Longrightarrow> ?iord z x"
krauss@40107
   117
    then have "chain ord (f ` A)"
krauss@40107
   118
      by (auto simp: img_ord_def intro: chainI dest: chainD)
krauss@40107
   119
    thus "?iord (?ilub A) x"
krauss@40107
   120
      unfolding inv img_lub_def img_ord_def
krauss@40107
   121
      by (rule lub_least) (auto dest: 1[unfolded img_ord_def])
krauss@40107
   122
  qed (auto simp: img_ord_def intro: leq_refl dest: leq_trans leq_antisym inj)
krauss@40107
   123
qed
krauss@40107
   124
krauss@40107
   125
context partial_function_definitions
krauss@40107
   126
begin
krauss@40107
   127
krauss@40107
   128
abbreviation "le_fun \<equiv> fun_ord leq"
krauss@40107
   129
abbreviation "lub_fun \<equiv> fun_lub lub"
krauss@40107
   130
abbreviation "fixp_fun == ccpo.fixp le_fun lub_fun"
krauss@40107
   131
abbreviation "mono_body \<equiv> monotone le_fun leq"
krauss@43081
   132
abbreviation "admissible \<equiv> ccpo.admissible le_fun lub_fun"
krauss@40107
   133
krauss@40107
   134
text {* Interpret manually, to avoid flooding everything with facts about
krauss@40107
   135
  orders *}
krauss@40107
   136
krauss@40107
   137
lemma ccpo: "class.ccpo le_fun (mk_less le_fun) lub_fun"
krauss@40107
   138
apply (rule ccpo)
krauss@40107
   139
apply (rule partial_function_lift)
krauss@40107
   140
apply (rule partial_function_definitions_axioms)
krauss@40107
   141
done
krauss@40107
   142
krauss@40107
   143
text {* The crucial fixed-point theorem *}
krauss@40107
   144
krauss@40107
   145
lemma mono_body_fixp: 
krauss@40107
   146
  "(\<And>x. mono_body (\<lambda>f. F f x)) \<Longrightarrow> fixp_fun F = F (fixp_fun F)"
krauss@40107
   147
by (rule ccpo.fixp_unfold[OF ccpo]) (auto simp: monotone_def fun_ord_def)
krauss@40107
   148
krauss@40107
   149
text {* Version with curry/uncurry combinators, to be used by package *}
krauss@40107
   150
krauss@40107
   151
lemma fixp_rule_uc:
krauss@40107
   152
  fixes F :: "'c \<Rightarrow> 'c" and
krauss@40107
   153
    U :: "'c \<Rightarrow> 'b \<Rightarrow> 'a" and
krauss@40107
   154
    C :: "('b \<Rightarrow> 'a) \<Rightarrow> 'c"
krauss@40107
   155
  assumes mono: "\<And>x. mono_body (\<lambda>f. U (F (C f)) x)"
krauss@40107
   156
  assumes eq: "f \<equiv> C (fixp_fun (\<lambda>f. U (F (C f))))"
krauss@40107
   157
  assumes inverse: "\<And>f. C (U f) = f"
krauss@40107
   158
  shows "f = F f"
krauss@40107
   159
proof -
krauss@40107
   160
  have "f = C (fixp_fun (\<lambda>f. U (F (C f))))" by (simp add: eq)
krauss@40107
   161
  also have "... = C (U (F (C (fixp_fun (\<lambda>f. U (F (C f)))))))"
krauss@40107
   162
    by (subst mono_body_fixp[of "%f. U (F (C f))", OF mono]) (rule refl)
krauss@40107
   163
  also have "... = F (C (fixp_fun (\<lambda>f. U (F (C f)))))" by (rule inverse)
krauss@40107
   164
  also have "... = F f" by (simp add: eq)
krauss@40107
   165
  finally show "f = F f" .
krauss@40107
   166
qed
krauss@40107
   167
krauss@40107
   168
text {* Rules for @{term mono_body}: *}
krauss@40107
   169
krauss@40107
   170
lemma const_mono[partial_function_mono]: "monotone ord leq (\<lambda>f. c)"
krauss@40107
   171
by (rule monotoneI) (rule leq_refl)
krauss@40107
   172
krauss@40107
   173
end
krauss@40107
   174
krauss@40107
   175
krauss@40107
   176
subsection {* Flat interpretation: tailrec and option *}
krauss@40107
   177
krauss@40107
   178
definition 
krauss@40107
   179
  "flat_ord b x y \<longleftrightarrow> x = b \<or> x = y"
krauss@40107
   180
krauss@40107
   181
definition 
krauss@40107
   182
  "flat_lub b A = (if A \<subseteq> {b} then b else (THE x. x \<in> A - {b}))"
krauss@40107
   183
krauss@40107
   184
lemma flat_interpretation:
krauss@40107
   185
  "partial_function_definitions (flat_ord b) (flat_lub b)"
krauss@40107
   186
proof
krauss@40107
   187
  fix A x assume 1: "chain (flat_ord b) A" "x \<in> A"
krauss@40107
   188
  show "flat_ord b x (flat_lub b A)"
krauss@40107
   189
  proof cases
krauss@40107
   190
    assume "x = b"
krauss@40107
   191
    thus ?thesis by (simp add: flat_ord_def)
krauss@40107
   192
  next
krauss@40107
   193
    assume "x \<noteq> b"
krauss@40107
   194
    with 1 have "A - {b} = {x}"
krauss@40107
   195
      by (auto elim: chainE simp: flat_ord_def)
krauss@40107
   196
    then have "flat_lub b A = x"
krauss@40107
   197
      by (auto simp: flat_lub_def)
krauss@40107
   198
    thus ?thesis by (auto simp: flat_ord_def)
krauss@40107
   199
  qed
krauss@40107
   200
next
krauss@40107
   201
  fix A z assume A: "chain (flat_ord b) A"
krauss@40107
   202
    and z: "\<And>x. x \<in> A \<Longrightarrow> flat_ord b x z"
krauss@40107
   203
  show "flat_ord b (flat_lub b A) z"
krauss@40107
   204
  proof cases
krauss@40107
   205
    assume "A \<subseteq> {b}"
krauss@40107
   206
    thus ?thesis
krauss@40107
   207
      by (auto simp: flat_lub_def flat_ord_def)
krauss@40107
   208
  next
krauss@40107
   209
    assume nb: "\<not> A \<subseteq> {b}"
krauss@40107
   210
    then obtain y where y: "y \<in> A" "y \<noteq> b" by auto
krauss@40107
   211
    with A have "A - {b} = {y}"
krauss@40107
   212
      by (auto elim: chainE simp: flat_ord_def)
krauss@40107
   213
    with nb have "flat_lub b A = y"
krauss@40107
   214
      by (auto simp: flat_lub_def)
krauss@40107
   215
    with z y show ?thesis by auto    
krauss@40107
   216
  qed
krauss@40107
   217
qed (auto simp: flat_ord_def)
krauss@40107
   218
krauss@40107
   219
interpretation tailrec!:
krauss@40107
   220
  partial_function_definitions "flat_ord undefined" "flat_lub undefined"
krauss@40107
   221
by (rule flat_interpretation)
krauss@40107
   222
krauss@40107
   223
interpretation option!:
krauss@40107
   224
  partial_function_definitions "flat_ord None" "flat_lub None"
krauss@40107
   225
by (rule flat_interpretation)
krauss@40107
   226
krauss@42949
   227
declaration {* Partial_Function.init "tailrec" @{term tailrec.fixp_fun}
krauss@43080
   228
  @{term tailrec.mono_body} @{thm tailrec.fixp_rule_uc} NONE *}
krauss@42949
   229
krauss@42949
   230
declaration {* Partial_Function.init "option" @{term option.fixp_fun}
krauss@43080
   231
  @{term option.mono_body} @{thm option.fixp_rule_uc} NONE *}
krauss@42949
   232
krauss@42949
   233
krauss@40107
   234
abbreviation "option_ord \<equiv> flat_ord None"
krauss@40107
   235
abbreviation "mono_option \<equiv> monotone (fun_ord option_ord) option_ord"
krauss@40107
   236
krauss@40107
   237
lemma bind_mono[partial_function_mono]:
krauss@40107
   238
assumes mf: "mono_option B" and mg: "\<And>y. mono_option (\<lambda>f. C y f)"
krauss@40107
   239
shows "mono_option (\<lambda>f. Option.bind (B f) (\<lambda>y. C y f))"
krauss@40107
   240
proof (rule monotoneI)
krauss@40107
   241
  fix f g :: "'a \<Rightarrow> 'b option" assume fg: "fun_ord option_ord f g"
krauss@40107
   242
  with mf
krauss@40107
   243
  have "option_ord (B f) (B g)" by (rule monotoneD[of _ _ _ f g])
krauss@40107
   244
  then have "option_ord (Option.bind (B f) (\<lambda>y. C y f)) (Option.bind (B g) (\<lambda>y. C y f))"
krauss@40107
   245
    unfolding flat_ord_def by auto    
krauss@40107
   246
  also from mg
krauss@40107
   247
  have "\<And>y'. option_ord (C y' f) (C y' g)"
krauss@40107
   248
    by (rule monotoneD) (rule fg)
krauss@40107
   249
  then have "option_ord (Option.bind (B g) (\<lambda>y'. C y' f)) (Option.bind (B g) (\<lambda>y'. C y' g))"
krauss@40107
   250
    unfolding flat_ord_def by (cases "B g") auto
krauss@40107
   251
  finally (option.leq_trans)
krauss@40107
   252
  show "option_ord (Option.bind (B f) (\<lambda>y. C y f)) (Option.bind (B g) (\<lambda>y'. C y' g))" .
krauss@40107
   253
qed
krauss@40107
   254
krauss@43081
   255
lemma flat_lub_in_chain:
krauss@43081
   256
  assumes ch: "chain (flat_ord b) A "
krauss@43081
   257
  assumes lub: "flat_lub b A = a"
krauss@43081
   258
  shows "a = b \<or> a \<in> A"
krauss@43081
   259
proof (cases "A \<subseteq> {b}")
krauss@43081
   260
  case True
krauss@43081
   261
  then have "flat_lub b A = b" unfolding flat_lub_def by simp
krauss@43081
   262
  with lub show ?thesis by simp
krauss@43081
   263
next
krauss@43081
   264
  case False
krauss@43081
   265
  then obtain c where "c \<in> A" and "c \<noteq> b" by auto
krauss@43081
   266
  { fix z assume "z \<in> A"
krauss@43081
   267
    from chainD[OF ch `c \<in> A` this] have "z = c \<or> z = b"
krauss@43081
   268
      unfolding flat_ord_def using `c \<noteq> b` by auto }
krauss@43081
   269
  with False have "A - {b} = {c}" by auto
krauss@43081
   270
  with False have "flat_lub b A = c" by (auto simp: flat_lub_def)
krauss@43081
   271
  with `c \<in> A` lub show ?thesis by simp
krauss@43081
   272
qed
krauss@43081
   273
krauss@43081
   274
lemma option_admissible: "option.admissible (%(f::'a \<Rightarrow> 'b option).
krauss@43081
   275
  (\<forall>x y. f x = Some y \<longrightarrow> P x y))"
krauss@43081
   276
proof (rule ccpo.admissibleI[OF option.ccpo])
krauss@43081
   277
  fix A :: "('a \<Rightarrow> 'b option) set"
krauss@43081
   278
  assume ch: "chain option.le_fun A"
krauss@43081
   279
    and IH: "\<forall>f\<in>A. \<forall>x y. f x = Some y \<longrightarrow> P x y"
krauss@43081
   280
  from ch have ch': "\<And>x. chain option_ord {y. \<exists>f\<in>A. y = f x}" by (rule chain_fun)
krauss@43081
   281
  show "\<forall>x y. option.lub_fun A x = Some y \<longrightarrow> P x y"
krauss@43081
   282
  proof (intro allI impI)
krauss@43081
   283
    fix x y assume "option.lub_fun A x = Some y"
krauss@43081
   284
    from flat_lub_in_chain[OF ch' this[unfolded fun_lub_def]]
krauss@43081
   285
    have "Some y \<in> {y. \<exists>f\<in>A. y = f x}" by simp
krauss@43081
   286
    then have "\<exists>f\<in>A. f x = Some y" by auto
krauss@43081
   287
    with IH show "P x y" by auto
krauss@43081
   288
  qed
krauss@43081
   289
qed
krauss@43081
   290
krauss@43081
   291
krauss@40252
   292
hide_const (open) chain
krauss@40107
   293
krauss@40107
   294
end
krauss@40107
   295