Theory UpperPD

(*  Title:      HOL/HOLCF/UpperPD.thy
    Author:     Brian Huffman
*)

section ‹Upper powerdomain›

theory UpperPD
imports Compact_Basis
begin

subsection ‹Basis preorder›

definition
  upper_le :: "'a pd_basis  'a pd_basis  bool" (infix "≤♯" 50) where
  "upper_le = (λu v. yRep_pd_basis v. xRep_pd_basis u. x  y)"

lemma upper_le_refl [simp]: "t ≤♯ t"
unfolding upper_le_def by fast

lemma upper_le_trans: "t ≤♯ u; u ≤♯ v  t ≤♯ v"
unfolding upper_le_def
apply (rule ballI)
apply (drule (1) bspec, erule bexE)
apply (drule (1) bspec, erule bexE)
apply (erule rev_bexI)
apply (erule (1) below_trans)
done

interpretation upper_le: preorder upper_le
by (rule preorder.intro, rule upper_le_refl, rule upper_le_trans)

lemma upper_le_minimal [simp]: "PDUnit compact_bot ≤♯ t"
unfolding upper_le_def Rep_PDUnit by simp

lemma PDUnit_upper_mono: "x  y  PDUnit x ≤♯ PDUnit y"
unfolding upper_le_def Rep_PDUnit by simp

lemma PDPlus_upper_mono: "s ≤♯ t; u ≤♯ v  PDPlus s u ≤♯ PDPlus t v"
unfolding upper_le_def Rep_PDPlus by fast

lemma PDPlus_upper_le: "PDPlus t u ≤♯ t"
unfolding upper_le_def Rep_PDPlus by fast

lemma upper_le_PDUnit_PDUnit_iff [simp]:
  "(PDUnit a ≤♯ PDUnit b) = (a  b)"
unfolding upper_le_def Rep_PDUnit by fast

lemma upper_le_PDPlus_PDUnit_iff:
  "(PDPlus t u ≤♯ PDUnit a) = (t ≤♯ PDUnit a  u ≤♯ PDUnit a)"
unfolding upper_le_def Rep_PDPlus Rep_PDUnit by fast

lemma upper_le_PDPlus_iff: "(t ≤♯ PDPlus u v) = (t ≤♯ u  t ≤♯ v)"
unfolding upper_le_def Rep_PDPlus by fast

lemma upper_le_induct [induct set: upper_le]:
  assumes le: "t ≤♯ u"
  assumes 1: "a b. a  b  P (PDUnit a) (PDUnit b)"
  assumes 2: "t u a. P t (PDUnit a)  P (PDPlus t u) (PDUnit a)"
  assumes 3: "t u v. P t u; P t v  P t (PDPlus u v)"
  shows "P t u"
using le apply (induct u arbitrary: t rule: pd_basis_induct)
apply (erule rev_mp)
apply (induct_tac t rule: pd_basis_induct)
apply (simp add: 1)
apply (simp add: upper_le_PDPlus_PDUnit_iff)
apply (simp add: 2)
apply (subst PDPlus_commute)
apply (simp add: 2)
apply (simp add: upper_le_PDPlus_iff 3)
done


subsection ‹Type definition›

typedef 'a upper_pd  ("('(_')♯)") =
  "{S::'a pd_basis set. upper_le.ideal S}"
by (rule upper_le.ex_ideal)

instantiation upper_pd :: (bifinite) below
begin

definition
  "x  y  Rep_upper_pd x  Rep_upper_pd y"

instance ..
end

instance upper_pd :: (bifinite) po
using type_definition_upper_pd below_upper_pd_def
by (rule upper_le.typedef_ideal_po)

instance upper_pd :: (bifinite) cpo
using type_definition_upper_pd below_upper_pd_def
by (rule upper_le.typedef_ideal_cpo)

definition
  upper_principal :: "'a pd_basis  'a upper_pd" where
  "upper_principal t = Abs_upper_pd {u. u ≤♯ t}"

interpretation upper_pd:
  ideal_completion upper_le upper_principal Rep_upper_pd
using type_definition_upper_pd below_upper_pd_def
using upper_principal_def pd_basis_countable
by (rule upper_le.typedef_ideal_completion)

text ‹Upper powerdomain is pointed›

lemma upper_pd_minimal: "upper_principal (PDUnit compact_bot)  ys"
by (induct ys rule: upper_pd.principal_induct, simp, simp)

instance upper_pd :: (bifinite) pcpo
by intro_classes (fast intro: upper_pd_minimal)

lemma inst_upper_pd_pcpo: " = upper_principal (PDUnit compact_bot)"
by (rule upper_pd_minimal [THEN bottomI, symmetric])


subsection ‹Monadic unit and plus›

definition
  upper_unit :: "'a  'a upper_pd" where
  "upper_unit = compact_basis.extension (λa. upper_principal (PDUnit a))"

definition
  upper_plus :: "'a upper_pd  'a upper_pd  'a upper_pd" where
  "upper_plus = upper_pd.extension (λt. upper_pd.extension (λu.
      upper_principal (PDPlus t u)))"

abbreviation
  upper_add :: "'a upper_pd  'a upper_pd  'a upper_pd"
    (infixl "∪♯" 65) where
  "xs ∪♯ ys == upper_plusxsys"

syntax
  "_upper_pd" :: "args  logic" ("{_}♯")

translations
  "{x,xs}♯" == "{x}♯ ∪♯ {xs}♯"
  "{x}♯" == "CONST upper_unitx"

lemma upper_unit_Rep_compact_basis [simp]:
  "{Rep_compact_basis a}♯ = upper_principal (PDUnit a)"
unfolding upper_unit_def
by (simp add: compact_basis.extension_principal PDUnit_upper_mono)

lemma upper_plus_principal [simp]:
  "upper_principal t ∪♯ upper_principal u = upper_principal (PDPlus t u)"
unfolding upper_plus_def
by (simp add: upper_pd.extension_principal
    upper_pd.extension_mono PDPlus_upper_mono)

interpretation upper_add: semilattice upper_add proof
  fix xs ys zs :: "'a upper_pd"
  show "(xs ∪♯ ys) ∪♯ zs = xs ∪♯ (ys ∪♯ zs)"
    apply (induct xs rule: upper_pd.principal_induct, simp)
    apply (induct ys rule: upper_pd.principal_induct, simp)
    apply (induct zs rule: upper_pd.principal_induct, simp)
    apply (simp add: PDPlus_assoc)
    done
  show "xs ∪♯ ys = ys ∪♯ xs"
    apply (induct xs rule: upper_pd.principal_induct, simp)
    apply (induct ys rule: upper_pd.principal_induct, simp)
    apply (simp add: PDPlus_commute)
    done
  show "xs ∪♯ xs = xs"
    apply (induct xs rule: upper_pd.principal_induct, simp)
    apply (simp add: PDPlus_absorb)
    done
qed

lemmas upper_plus_assoc = upper_add.assoc
lemmas upper_plus_commute = upper_add.commute
lemmas upper_plus_absorb = upper_add.idem
lemmas upper_plus_left_commute = upper_add.left_commute
lemmas upper_plus_left_absorb = upper_add.left_idem

text ‹Useful for simp add: upper_plus_ac›
lemmas upper_plus_ac =
  upper_plus_assoc upper_plus_commute upper_plus_left_commute

text ‹Useful for simp only: upper_plus_aci›
lemmas upper_plus_aci =
  upper_plus_ac upper_plus_absorb upper_plus_left_absorb

lemma upper_plus_below1: "xs ∪♯ ys  xs"
apply (induct xs rule: upper_pd.principal_induct, simp)
apply (induct ys rule: upper_pd.principal_induct, simp)
apply (simp add: PDPlus_upper_le)
done

lemma upper_plus_below2: "xs ∪♯ ys  ys"
by (subst upper_plus_commute, rule upper_plus_below1)

lemma upper_plus_greatest: "xs  ys; xs  zs  xs  ys ∪♯ zs"
apply (subst upper_plus_absorb [of xs, symmetric])
apply (erule (1) monofun_cfun [OF monofun_cfun_arg])
done

lemma upper_below_plus_iff [simp]:
  "xs  ys ∪♯ zs  xs  ys  xs  zs"
apply safe
apply (erule below_trans [OF _ upper_plus_below1])
apply (erule below_trans [OF _ upper_plus_below2])
apply (erule (1) upper_plus_greatest)
done

lemma upper_plus_below_unit_iff [simp]:
  "xs ∪♯ ys  {z}♯  xs  {z}♯  ys  {z}♯"
apply (induct xs rule: upper_pd.principal_induct, simp)
apply (induct ys rule: upper_pd.principal_induct, simp)
apply (induct z rule: compact_basis.principal_induct, simp)
apply (simp add: upper_le_PDPlus_PDUnit_iff)
done

lemma upper_unit_below_iff [simp]: "{x}♯  {y}♯  x  y"
apply (induct x rule: compact_basis.principal_induct, simp)
apply (induct y rule: compact_basis.principal_induct, simp)
apply simp
done

lemmas upper_pd_below_simps =
  upper_unit_below_iff
  upper_below_plus_iff
  upper_plus_below_unit_iff

lemma upper_unit_eq_iff [simp]: "{x}♯ = {y}♯  x = y"
unfolding po_eq_conv by simp

lemma upper_unit_strict [simp]: "{}♯ = "
using upper_unit_Rep_compact_basis [of compact_bot]
by (simp add: inst_upper_pd_pcpo)

lemma upper_plus_strict1 [simp]: " ∪♯ ys = "
by (rule bottomI, rule upper_plus_below1)

lemma upper_plus_strict2 [simp]: "xs ∪♯  = "
by (rule bottomI, rule upper_plus_below2)

lemma upper_unit_bottom_iff [simp]: "{x}♯ =   x = "
unfolding upper_unit_strict [symmetric] by (rule upper_unit_eq_iff)

lemma upper_plus_bottom_iff [simp]:
  "xs ∪♯ ys =   xs =   ys = "
apply (induct xs rule: upper_pd.principal_induct, simp)
apply (induct ys rule: upper_pd.principal_induct, simp)
apply (simp add: inst_upper_pd_pcpo upper_pd.principal_eq_iff
                 upper_le_PDPlus_PDUnit_iff)
done

lemma compact_upper_unit: "compact x  compact {x}♯"
by (auto dest!: compact_basis.compact_imp_principal)

lemma compact_upper_unit_iff [simp]: "compact {x}♯  compact x"
apply (safe elim!: compact_upper_unit)
apply (simp only: compact_def upper_unit_below_iff [symmetric])
apply (erule adm_subst [OF cont_Rep_cfun2])
done

lemma compact_upper_plus [simp]:
  "compact xs; compact ys  compact (xs ∪♯ ys)"
by (auto dest!: upper_pd.compact_imp_principal)


subsection ‹Induction rules›

lemma upper_pd_induct1:
  assumes P: "adm P"
  assumes unit: "x. P {x}♯"
  assumes insert: "x ys. P {x}♯; P ys  P ({x}♯ ∪♯ ys)"
  shows "P (xs::'a upper_pd)"
apply (induct xs rule: upper_pd.principal_induct, rule P)
apply (induct_tac a rule: pd_basis_induct1)
apply (simp only: upper_unit_Rep_compact_basis [symmetric])
apply (rule unit)
apply (simp only: upper_unit_Rep_compact_basis [symmetric]
                  upper_plus_principal [symmetric])
apply (erule insert [OF unit])
done

lemma upper_pd_induct
  [case_names adm upper_unit upper_plus, induct type: upper_pd]:
  assumes P: "adm P"
  assumes unit: "x. P {x}♯"
  assumes plus: "xs ys. P xs; P ys  P (xs ∪♯ ys)"
  shows "P (xs::'a upper_pd)"
apply (induct xs rule: upper_pd.principal_induct, rule P)
apply (induct_tac a rule: pd_basis_induct)
apply (simp only: upper_unit_Rep_compact_basis [symmetric] unit)
apply (simp only: upper_plus_principal [symmetric] plus)
done


subsection ‹Monadic bind›

definition
  upper_bind_basis ::
  "'a pd_basis  ('a  'b upper_pd)  'b upper_pd" where
  "upper_bind_basis = fold_pd
    (λa. Λ f. f(Rep_compact_basis a))
    (λx y. Λ f. xf ∪♯ yf)"

lemma ACI_upper_bind:
  "semilattice (λx y. Λ f. xf ∪♯ yf)"
apply unfold_locales
apply (simp add: upper_plus_assoc)
apply (simp add: upper_plus_commute)
apply (simp add: eta_cfun)
done

lemma upper_bind_basis_simps [simp]:
  "upper_bind_basis (PDUnit a) =
    (Λ f. f(Rep_compact_basis a))"
  "upper_bind_basis (PDPlus t u) =
    (Λ f. upper_bind_basis tf ∪♯ upper_bind_basis uf)"
unfolding upper_bind_basis_def
apply -
apply (rule fold_pd_PDUnit [OF ACI_upper_bind])
apply (rule fold_pd_PDPlus [OF ACI_upper_bind])
done

lemma upper_bind_basis_mono:
  "t ≤♯ u  upper_bind_basis t  upper_bind_basis u"
unfolding cfun_below_iff
apply (erule upper_le_induct, safe)
apply (simp add: monofun_cfun)
apply (simp add: below_trans [OF upper_plus_below1])
apply simp
done

definition
  upper_bind :: "'a upper_pd  ('a  'b upper_pd)  'b upper_pd" where
  "upper_bind = upper_pd.extension upper_bind_basis"

syntax
  "_upper_bind" :: "[logic, logic, logic]  logic"
    ("(3⋃♯__./ _)" [0, 0, 10] 10)

translations
  "⋃♯xxs. e" == "CONST upper_bindxs(Λ x. e)"

lemma upper_bind_principal [simp]:
  "upper_bind(upper_principal t) = upper_bind_basis t"
unfolding upper_bind_def
apply (rule upper_pd.extension_principal)
apply (erule upper_bind_basis_mono)
done

lemma upper_bind_unit [simp]:
  "upper_bind{x}♯f = fx"
by (induct x rule: compact_basis.principal_induct, simp, simp)

lemma upper_bind_plus [simp]:
  "upper_bind(xs ∪♯ ys)f = upper_bindxsf ∪♯ upper_bindysf"
by (induct xs rule: upper_pd.principal_induct, simp,
    induct ys rule: upper_pd.principal_induct, simp, simp)

lemma upper_bind_strict [simp]: "upper_bindf = f"
unfolding upper_unit_strict [symmetric] by (rule upper_bind_unit)

lemma upper_bind_bind:
  "upper_bind(upper_bindxsf)g = upper_bindxs(Λ x. upper_bind(fx)g)"
by (induct xs, simp_all)


subsection ‹Map›

definition
  upper_map :: "('a  'b)  'a upper_pd  'b upper_pd" where
  "upper_map = (Λ f xs. upper_bindxs(Λ x. {fx}♯))"

lemma upper_map_unit [simp]:
  "upper_mapf{x}♯ = {fx}♯"
unfolding upper_map_def by simp

lemma upper_map_plus [simp]:
  "upper_mapf(xs ∪♯ ys) = upper_mapfxs ∪♯ upper_mapfys"
unfolding upper_map_def by simp

lemma upper_map_bottom [simp]: "upper_mapf = {f}♯"
unfolding upper_map_def by simp

lemma upper_map_ident: "upper_map(Λ x. x)xs = xs"
by (induct xs rule: upper_pd_induct, simp_all)

lemma upper_map_ID: "upper_mapID = ID"
by (simp add: cfun_eq_iff ID_def upper_map_ident)

lemma upper_map_map:
  "upper_mapf(upper_mapgxs) = upper_map(Λ x. f(gx))xs"
by (induct xs rule: upper_pd_induct, simp_all)

lemma upper_bind_map:
  "upper_bind(upper_mapfxs)g = upper_bindxs(Λ x. g(fx))"
by (simp add: upper_map_def upper_bind_bind)

lemma upper_map_bind:
  "upper_mapf(upper_bindxsg) = upper_bindxs(Λ x. upper_mapf(gx))"
by (simp add: upper_map_def upper_bind_bind)

lemma ep_pair_upper_map: "ep_pair e p  ep_pair (upper_mape) (upper_mapp)"
apply standard
apply (induct_tac x rule: upper_pd_induct, simp_all add: ep_pair.e_inverse)
apply (induct_tac y rule: upper_pd_induct)
apply (simp_all add: ep_pair.e_p_below monofun_cfun del: upper_below_plus_iff)
done

lemma deflation_upper_map: "deflation d  deflation (upper_mapd)"
apply standard
apply (induct_tac x rule: upper_pd_induct, simp_all add: deflation.idem)
apply (induct_tac x rule: upper_pd_induct)
apply (simp_all add: deflation.below monofun_cfun del: upper_below_plus_iff)
done

(* FIXME: long proof! *)
lemma finite_deflation_upper_map:
  assumes "finite_deflation d" shows "finite_deflation (upper_mapd)"
proof (rule finite_deflation_intro)
  interpret d: finite_deflation d by fact
  from d.deflation_axioms show "deflation (upper_mapd)"
    by (rule deflation_upper_map)
  have "finite (range (λx. dx))" by (rule d.finite_range)
  hence "finite (Rep_compact_basis -` range (λx. dx))"
    by (rule finite_vimageI, simp add: inj_on_def Rep_compact_basis_inject)
  hence "finite (Pow (Rep_compact_basis -` range (λx. dx)))" by simp
  hence "finite (Rep_pd_basis -` (Pow (Rep_compact_basis -` range (λx. dx))))"
    by (rule finite_vimageI, simp add: inj_on_def Rep_pd_basis_inject)
  hence *: "finite (upper_principal ` Rep_pd_basis -` (Pow (Rep_compact_basis -` range (λx. dx))))" by simp
  hence "finite (range (λxs. upper_mapdxs))"
    apply (rule rev_finite_subset)
    apply clarsimp
    apply (induct_tac xs rule: upper_pd.principal_induct)
    apply (simp add: adm_mem_finite *)
    apply (rename_tac t, induct_tac t rule: pd_basis_induct)
    apply (simp only: upper_unit_Rep_compact_basis [symmetric] upper_map_unit)
    apply simp
    apply (subgoal_tac "b. d(Rep_compact_basis a) = Rep_compact_basis b")
    apply clarsimp
    apply (rule imageI)
    apply (rule vimageI2)
    apply (simp add: Rep_PDUnit)
    apply (rule range_eqI)
    apply (erule sym)
    apply (rule exI)
    apply (rule Abs_compact_basis_inverse [symmetric])
    apply (simp add: d.compact)
    apply (simp only: upper_plus_principal [symmetric] upper_map_plus)
    apply clarsimp
    apply (rule imageI)
    apply (rule vimageI2)
    apply (simp add: Rep_PDPlus)
    done
  thus "finite {xs. upper_mapdxs = xs}"
    by (rule finite_range_imp_finite_fixes)
qed

subsection ‹Upper powerdomain is bifinite›

lemma approx_chain_upper_map:
  assumes "approx_chain a"
  shows "approx_chain (λi. upper_map(a i))"
  using assms unfolding approx_chain_def
  by (simp add: lub_APP upper_map_ID finite_deflation_upper_map)

instance upper_pd :: (bifinite) bifinite
proof
  show "(a::nat  'a upper_pd  'a upper_pd). approx_chain a"
    using bifinite [where 'a='a]
    by (fast intro!: approx_chain_upper_map)
qed

subsection ‹Join›

definition
  upper_join :: "'a upper_pd upper_pd  'a upper_pd" where
  "upper_join = (Λ xss. upper_bindxss(Λ xs. xs))"

lemma upper_join_unit [simp]:
  "upper_join{xs}♯ = xs"
unfolding upper_join_def by simp

lemma upper_join_plus [simp]:
  "upper_join(xss ∪♯ yss) = upper_joinxss ∪♯ upper_joinyss"
unfolding upper_join_def by simp

lemma upper_join_bottom [simp]: "upper_join = "
unfolding upper_join_def by simp

lemma upper_join_map_unit:
  "upper_join(upper_mapupper_unitxs) = xs"
by (induct xs rule: upper_pd_induct, simp_all)

lemma upper_join_map_join:
  "upper_join(upper_mapupper_joinxsss) = upper_join(upper_joinxsss)"
by (induct xsss rule: upper_pd_induct, simp_all)

lemma upper_join_map_map:
  "upper_join(upper_map(upper_mapf)xss) =
   upper_mapf(upper_joinxss)"
by (induct xss rule: upper_pd_induct, simp_all)

end