cofe_solver.v 9.75 KB
Newer Older
1
From algebra Require Export cofe.
Robbert Krebbers's avatar
Robbert Krebbers committed
2

3
Record solution (F : cFunctor) := Solution {
4
  solution_car :> cofeT;
5 6
  solution_unfold : solution_car -n> F solution_car;
  solution_fold : F solution_car -n> solution_car;
7 8 9 10 11 12 13
  solution_fold_unfold X : solution_fold (solution_unfold X)  X;
  solution_unfold_fold X : solution_unfold (solution_fold X)  X
}.
Arguments solution_unfold {_} _.
Arguments solution_fold {_} _.

Module solver. Section solver.
14 15 16 17
Context (F : cFunctor) `{Finhab : Inhabited (F unitC)}.
Context (map_contractive :  {A1 A2 B1 B2},
  Contractive (@cFunctor_map F A1 A2 B1 B2)).
Notation map := (cFunctor_map F).
Robbert Krebbers's avatar
Robbert Krebbers committed
18 19

Fixpoint A (k : nat) : cofeT :=
20
  match k with 0 => unitC | S k => F (A k) end.
21 22 23 24 25 26 27 28 29 30 31
Fixpoint f (k : nat) : A k -n> A (S k) :=
  match k with 0 => CofeMor (λ _, inhabitant) | S k => map (g k,f k) end
with g (k : nat) : A (S k) -n> A k :=
  match k with 0 => CofeMor (λ _, ()) | S k => map (f k,g k) end.
Definition f_S k (x : A (S k)) : f (S k) x = map (g k,f k) x := eq_refl.
Definition g_S k (x : A (S (S k))) : g (S k) x = map (f k,g k) x := eq_refl.
Arguments A : simpl never.
Arguments f : simpl never.
Arguments g : simpl never.

Lemma gf {k} (x : A k) : g k (f k x)  x.
Robbert Krebbers's avatar
Robbert Krebbers committed
32 33
Proof.
  induction k as [|k IH]; simpl in *; [by destruct x|].
34
  rewrite -cFunctor_compose -{2}[x]cFunctor_id. by apply (contractive_proper map).
Robbert Krebbers's avatar
Robbert Krebbers committed
35
Qed.
36
Lemma fg {k} (x : A (S (S k))) : f (S k) (g (S k) x) {k} x.
Robbert Krebbers's avatar
Robbert Krebbers committed
37
Proof.
38
  induction k as [|k IH]; simpl.
39 40 41 42
  - rewrite f_S g_S -{2}[x]cFunctor_id -cFunctor_compose.
    apply (contractive_0 map).
  - rewrite f_S g_S -{2}[x]cFunctor_id -cFunctor_compose.
    by apply (contractive_S map).
Robbert Krebbers's avatar
Robbert Krebbers committed
43 44 45 46
Qed.

Record tower := {
  tower_car k :> A k;
47
  g_tower k : g k (tower_car (S k))  tower_car k
Robbert Krebbers's avatar
Robbert Krebbers committed
48 49
}.
Instance tower_equiv : Equiv tower := λ X Y,  k, X k  Y k.
50
Instance tower_dist : Dist tower := λ n X Y,  k, X k {n} Y k.
Robbert Krebbers's avatar
Robbert Krebbers committed
51 52 53 54 55 56
Program Definition tower_chain (c : chain tower) (k : nat) : chain (A k) :=
  {| chain_car i := c i k |}.
Next Obligation. intros c k n i ?; apply (chain_cauchy c n); lia. Qed.
Program Instance tower_compl : Compl tower := λ c,
  {| tower_car n := compl (tower_chain c n) |}.
Next Obligation.
57
  intros c k; apply equiv_dist=> n.
Robbert Krebbers's avatar
Robbert Krebbers committed
58
  by rewrite (conv_compl n (tower_chain c k))
59
    (conv_compl n (tower_chain c (S k))) /= (g_tower (c _) k).
Robbert Krebbers's avatar
Robbert Krebbers committed
60
Qed.
61
Definition tower_cofe_mixin : CofeMixin tower.
Robbert Krebbers's avatar
Robbert Krebbers committed
62 63
Proof.
  split.
64
  - intros X Y; split; [by intros HXY n k; apply equiv_dist|].
Robbert Krebbers's avatar
Robbert Krebbers committed
65
    intros HXY k; apply equiv_dist; intros n; apply HXY.
66
  - intros k; split.
Robbert Krebbers's avatar
Robbert Krebbers committed
67 68
    + by intros X n.
    + by intros X Y ? n.
69
    + by intros X Y Z ?? n; trans (Y n).
70
  - intros k X Y HXY n; apply dist_S.
Robbert Krebbers's avatar
Robbert Krebbers committed
71
    by rewrite -(g_tower X) (HXY (S n)) g_tower.
Robbert Krebbers's avatar
Robbert Krebbers committed
72
  - intros n c k; rewrite /= (conv_compl n (tower_chain c k)).
Robbert Krebbers's avatar
Robbert Krebbers committed
73 74
    apply (chain_cauchy c); lia.
Qed.
75
Definition T : cofeT := CofeT tower_cofe_mixin.
Robbert Krebbers's avatar
Robbert Krebbers committed
76 77

Fixpoint ff {k} (i : nat) : A k -n> A (i + k) :=
78
  match i with 0 => cid | S i => f (i + k)  ff i end.
Robbert Krebbers's avatar
Robbert Krebbers committed
79
Fixpoint gg {k} (i : nat) : A (i + k) -n> A k :=
80
  match i with 0 => cid | S i => gg i  g (i + k) end.
Robbert Krebbers's avatar
Robbert Krebbers committed
81
Lemma ggff {k i} (x : A k) : gg i (ff i x)  x.
Robbert Krebbers's avatar
Robbert Krebbers committed
82
Proof. induction i as [|i IH]; simpl; [done|by rewrite (gf (ff i x)) IH]. Qed.
83 84 85
Lemma f_tower k (X : tower) : f (S k) (X (S k)) {k} X (S (S k)).
Proof. intros. by rewrite -(fg (X (S (S k)))) -(g_tower X). Qed.
Lemma ff_tower k i (X : tower) : ff i (X (S k)) {k} X (i + S k).
Robbert Krebbers's avatar
Robbert Krebbers committed
86 87
Proof.
  intros; induction i as [|i IH]; simpl; [done|].
88
  by rewrite IH Nat.add_succ_r (dist_le _ _ _ _ (f_tower _ X)); last omega.
Robbert Krebbers's avatar
Robbert Krebbers committed
89 90
Qed.
Lemma gg_tower k i (X : tower) : gg i (X (i + k))  X k.
Robbert Krebbers's avatar
Robbert Krebbers committed
91
Proof. by induction i as [|i IH]; simpl; [done|rewrite g_tower IH]. Qed.
Robbert Krebbers's avatar
Robbert Krebbers committed
92 93 94 95 96 97 98 99 100 101 102 103 104

Instance tower_car_ne n k : Proper (dist n ==> dist n) (λ X, tower_car X k).
Proof. by intros X Y HX. Qed.
Definition project (k : nat) : T -n> A k := CofeMor (λ X : T, tower_car X k).

Definition coerce {i j} (H : i = j) : A i -n> A j :=
  eq_rect _ (λ i', A i -n> A i') cid _ H.
Lemma coerce_id {i} (H : i = i) (x : A i) : coerce H x = x.
Proof. unfold coerce. by rewrite (proof_irrel H (eq_refl i)). Qed.
Lemma coerce_proper {i j} (x y : A i) (H1 H2 : i = j) :
  x = y  coerce H1 x = coerce H2 y.
Proof. by destruct H1; rewrite !coerce_id. Qed.
Lemma g_coerce {k j} (H : S k = S j) (x : A (S k)) :
105
  g j (coerce H x) = coerce (Nat.succ_inj _ _ H) (g k x).
Robbert Krebbers's avatar
Robbert Krebbers committed
106 107
Proof. by assert (k = j) by lia; subst; rewrite !coerce_id. Qed.
Lemma coerce_f {k j} (H : S k = S j) (x : A k) :
108
  coerce H (f k x) = f j (coerce (Nat.succ_inj _ _ H) x).
Robbert Krebbers's avatar
Robbert Krebbers committed
109 110 111 112
Proof. by assert (k = j) by lia; subst; rewrite !coerce_id. Qed.
Lemma gg_gg {k i i1 i2 j} (H1 : k = i + j) (H2 : k = i2 + (i1 + j)) (x : A k) :
  gg i (coerce H1 x) = gg i1 (gg i2 (coerce H2 x)).
Proof.
113 114
  assert (i = i2 + i1) by lia; simplify_eq/=. revert j x H1.
  induction i2 as [|i2 IH]; intros j X H1; simplify_eq/=;
Robbert Krebbers's avatar
Robbert Krebbers committed
115
    [by rewrite coerce_id|by rewrite g_coerce IH].
Robbert Krebbers's avatar
Robbert Krebbers committed
116 117 118 119
Qed.
Lemma ff_ff {k i i1 i2 j} (H1 : i + k = j) (H2 : i1 + (i2 + k) = j) (x : A k) :
  coerce H1 (ff i x) = coerce H2 (ff i1 (ff i2 x)).
Proof.
120 121
  assert (i = i1 + i2) by lia; simplify_eq/=.
  induction i1 as [|i1 IH]; simplify_eq/=;
Robbert Krebbers's avatar
Robbert Krebbers committed
122
    [by rewrite coerce_id|by rewrite coerce_f IH].
Robbert Krebbers's avatar
Robbert Krebbers committed
123 124
Qed.

125
Definition embed_coerce {k} (i : nat) : A k -n> A i :=
Robbert Krebbers's avatar
Robbert Krebbers committed
126 127 128 129
  match le_lt_dec i k with
  | left H => gg (k-i)  coerce (eq_sym (Nat.sub_add _ _ H))
  | right H => coerce (Nat.sub_add k i (Nat.lt_le_incl _ _ H))  ff (i-k)
  end.
130 131
Lemma g_embed_coerce {k i} (x : A k) :
  g i (embed_coerce (S i) x)  embed_coerce i x.
Robbert Krebbers's avatar
Robbert Krebbers committed
132
Proof.
133
  unfold embed_coerce; destruct (le_lt_dec (S i) k), (le_lt_dec i k); simpl.
134 135 136
  - symmetry; by erewrite (@gg_gg _ _ 1 (k - S i)); simpl.
  - exfalso; lia.
  - assert (i = k) by lia; subst.
Robbert Krebbers's avatar
Robbert Krebbers committed
137
    rewrite (ff_ff _ (eq_refl (1 + (0 + k)))) /= gf.
Robbert Krebbers's avatar
Robbert Krebbers committed
138
    by rewrite (gg_gg _ (eq_refl (0 + (0 + k)))).
139
  - assert (H : 1 + ((i - k) + k) = S i) by lia.
Robbert Krebbers's avatar
Robbert Krebbers committed
140
    rewrite (ff_ff _ H) /= -{2}(gf (ff (i - k) x)) g_coerce.
Robbert Krebbers's avatar
Robbert Krebbers committed
141 142
    by erewrite coerce_proper by done.
Qed.
143 144 145 146 147
Program Definition embed (k : nat) (x : A k) : T :=
  {| tower_car n := embed_coerce n x |}.
Next Obligation. intros k x i. apply g_embed_coerce. Qed.
Instance: Params (@embed) 1.
Instance embed_ne k n : Proper (dist n ==> dist n) (embed k).
148
Proof. by intros x y Hxy i; rewrite /= Hxy. Qed.
149 150
Definition embed' (k : nat) : A k -n> T := CofeMor (embed k).
Lemma embed_f k (x : A k) : embed (S k) (f k x)  embed k x.
Robbert Krebbers's avatar
Robbert Krebbers committed
151
Proof.
152
  rewrite equiv_dist=> n i; rewrite /embed /= /embed_coerce.
Robbert Krebbers's avatar
Robbert Krebbers committed
153
  destruct (le_lt_dec i (S k)), (le_lt_dec i k); simpl.
154
  - assert (H : S k = S (k - i) + (0 + i)) by lia; rewrite (gg_gg _ H) /=.
Robbert Krebbers's avatar
Robbert Krebbers committed
155
    by erewrite g_coerce, gf, coerce_proper by done.
156
  - assert (S k = 0 + (0 + i)) as H by lia.
157
    rewrite (gg_gg _ H); simplify_eq/=.
Robbert Krebbers's avatar
Robbert Krebbers committed
158
    by rewrite (ff_ff _ (eq_refl (1 + (0 + k)))).
159 160
  - exfalso; lia.
  - assert (H : (i - S k) + (1 + k) = i) by lia; rewrite (ff_ff _ H) /=.
Robbert Krebbers's avatar
Robbert Krebbers committed
161 162
    by erewrite coerce_proper by done.
Qed.
163
Lemma embed_tower k (X : T) : embed (S k) (X (S k)) {k} X.
Robbert Krebbers's avatar
Robbert Krebbers committed
164
Proof.
165 166
  intros i; rewrite /= /embed_coerce.
  destruct (le_lt_dec i (S k)) as [H|H]; simpl.
167
  - rewrite -(gg_tower i (S k - i) X).
Robbert Krebbers's avatar
Robbert Krebbers committed
168
    apply (_ : Proper (_ ==> _) (gg _)); by destruct (eq_sym _).
169
  - rewrite (ff_tower k (i - S k) X). by destruct (Nat.sub_add _ _ _).
Robbert Krebbers's avatar
Robbert Krebbers committed
170 171
Qed.

172
Program Definition unfold_chain (X : T) : chain (F T) :=
173
  {| chain_car n := map (project n,embed' n) (X (S n)) |}.
Robbert Krebbers's avatar
Robbert Krebbers committed
174 175 176
Next Obligation.
  intros X n i Hi.
  assert ( k, i = k + n) as [k ?] by (exists (i - n); lia); subst; clear Hi.
177
  induction k as [|k IH]; simpl; first done.
178
  rewrite -IH -(dist_le _ _ _ _ (f_tower (k + n) _)); last lia.
179
  rewrite f_S -cFunctor_compose.
180
  by apply (contractive_ne map); split=> Y /=; rewrite ?g_tower ?embed_f.
Robbert Krebbers's avatar
Robbert Krebbers committed
181
Qed.
182
Definition unfold (X : T) : F T := compl (unfold_chain X).
183
Instance unfold_ne : Proper (dist n ==> dist n) unfold.
184
Proof.
Robbert Krebbers's avatar
Robbert Krebbers committed
185
  intros n X Y HXY. by rewrite /unfold (conv_compl n (unfold_chain X))
186
    (conv_compl n (unfold_chain Y)) /= (HXY (S n)).
187
Qed.
Robbert Krebbers's avatar
Robbert Krebbers committed
188

189
Program Definition fold (X : F T) : T :=
190
  {| tower_car n := g n (map (embed' n,project n) X) |}.
Robbert Krebbers's avatar
Robbert Krebbers committed
191
Next Obligation.
192
  intros X k. apply (_ : Proper (() ==> ()) (g k)).
193
  rewrite g_S -cFunctor_compose.
194
  apply (contractive_proper map); split=> Y; [apply embed_f|apply g_tower].
Robbert Krebbers's avatar
Robbert Krebbers committed
195
Qed.
196 197
Instance fold_ne : Proper (dist n ==> dist n) fold.
Proof. by intros n X Y HXY k; rewrite /fold /= HXY. Qed.
Robbert Krebbers's avatar
Robbert Krebbers committed
198

199
Theorem result : solution F.
Robbert Krebbers's avatar
Robbert Krebbers committed
200
Proof.
201
  apply (Solution F T (CofeMor unfold) (CofeMor fold)).
202
  - move=> X /=. rewrite equiv_dist=> n k; rewrite /unfold /fold /=.
203
    rewrite -g_tower -(gg_tower _ n); apply (_ : Proper (_ ==> _) (g _)).
204
    trans (map (ff n, gg n) (X (S (n + k)))).
Robbert Krebbers's avatar
Robbert Krebbers committed
205
    { rewrite /unfold (conv_compl n (unfold_chain X)).
206 207
      rewrite -(chain_cauchy (unfold_chain X) n (S (n + k))) /=; last lia.
      rewrite -(dist_le _ _ _ _ (f_tower (n + k) _)); last lia.
208
      rewrite f_S -!cFunctor_compose; apply (contractive_ne map); split=> Y.
209 210 211 212 213 214 215 216 217
      + rewrite /embed' /= /embed_coerce.
        destruct (le_lt_dec _ _); simpl; [exfalso; lia|].
        by rewrite (ff_ff _ (eq_refl (S n + (0 + k)))) /= gf.
      + rewrite /embed' /= /embed_coerce.
        destruct (le_lt_dec _ _); simpl; [|exfalso; lia].
        by rewrite (gg_gg _ (eq_refl (0 + (S n + k)))) /= gf. }
    assert ( i k (x : A (S i + k)) (H : S i + k = i + S k),
      map (ff i, gg i) x  gg i (coerce H x)) as map_ff_gg.
    { intros i; induction i as [|i IH]; intros k' x H; simpl.
218 219
      { by rewrite coerce_id cFunctor_id. }
      rewrite cFunctor_compose g_coerce; apply IH. }
220 221
    assert (H: S n + k = n + S k) by lia.
    rewrite (map_ff_gg _ _ _ H).
222
    apply (_ : Proper (_ ==> _) (gg _)); by destruct H.
223
  - intros X; rewrite equiv_dist=> n /=.
224
    rewrite /unfold /= (conv_compl' n (unfold_chain (fold X))) /=.
225
    rewrite g_S -!cFunctor_compose -{2}[X]cFunctor_id.
226
    apply (contractive_ne map); split => Y /=.
227
    + rewrite f_tower. apply dist_S. by rewrite embed_tower.
228
    + etrans; [apply embed_ne, equiv_dist, g_tower|apply embed_tower].
Robbert Krebbers's avatar
Robbert Krebbers committed
229
Qed.
230
End solver. End solver.