Hoare: Hoare Logic
Over the past two weeks, we've begun applying the mathematical
tools developed in the first part of the course to studying
the theory of a small programming language. (This language is
often called Imp, since it is a language of pure imperative
programs.)
- We defined a type of abstract syntax trees for Imp,
together with an evaluation relation (a partial function
on states) that specifies the operational semantics of
programs.
If the course had ended at this point, we would still have gotten to something extremely useful: a set of tools for defining and discussing programming languages and language features that are mathematically precise, easy to work with, and very flexible.
The language we defined, though small, captures some of the key features of full-blown languages like C, C++, and Java, including the fundamental notion of mutable state and some common control structures.
- We proved a number of metatheoretic properties ("meta" in
the sense that they are properties of the language as a
whole, rather than properties of particular programs in the
language), including:
- determinacy of evaluation
- equivalence of some different ways of writing down the
definition
- guaranteed termination of certain classes of programs
- correctness (in the sense of preserving behavioral
equivalence) of a number of useful program
transformations
We'll return to this theme later in the course when we discuss types and type soundness.
- determinacy of evaluation
- We saw a couple of examples of program verification -- using the precise definition of Imp to prove formally that certain particular programs (factorial and slow subtraction) satisfied particular specifications of their behavior.
This week, we're going to take this idea further. We'll
develop a reasoning system called Floyd-Hoare Logic
(often shortened to just Hoare Logic), in which each of
the syntactic constructs of Imp is equipped with a single,
generic "proof rule" that can be used to reason about
programs involving this construct.
Hoare Logic has a long history, dating back to the 1960s, and it has been the subject of intensive research right up to the present day. It lies at the core of a huge variety of tools that are now being used to specify and verify real software systems.
Hoare Logic has a long history, dating back to the 1960s, and it has been the subject of intensive research right up to the present day. It lies at the core of a huge variety of tools that are now being used to specify and verify real software systems.
If we're going to talk about specifications of programs, the
first thing we'll want is a way of making assertions about
properties that hold at particular points in time -- i.e.,
properties that may or may not be true of a given state.
Paraphrase the following assertions in English.
fun st => st X = 3
fun st => st X = x
fun st => st X <= st Y
fun st => st X = 3 \/ st X <= st Y
fun st => (mult (st Z) (st Z)) <= x
/\ ~ ((mult (S (st Z)) (S (st Z))) <= x)
fun st => True
fun st => False
(Remember that one-star exercises do not need to be handed in.)
fun st => st X = 3
fun st => st X = x
fun st => st X <= st Y
fun st => st X = 3 \/ st X <= st Y
fun st => (mult (st Z) (st Z)) <= x
/\ ~ ((mult (S (st Z)) (S (st Z))) <= x)
fun st => True
fun st => False
(Remember that one-star exercises do not need to be handed in.)
☐
This way of writing assertions is correct -- it precisely
captures what we mean, and it is exactly what we will use in
Coq proofs -- but it is not very nice to look at: every
single assertion that we ever write is going to begin with
fun st => , and everyplace we refer to a variable in the
current state it is written st X.
Moreover, this is the *only* way we use states and the only way we refer to the values of variables in the current state. (We never need to talk about two states at the same time, etc.) So when we are writing down assertions in informal contexts, we can make some simplifications: drop the initial fun st => and write just X instead of st X. For example, instead of
fun st => (mult (st Z) (st Z)) <= x
/\ ~ ((mult (S (st Z)) (S (st Z))) <= x)
we'll write just
(mult Z Z) <= x /\ ~ ((mult (S Z) (S Z)) <= x).
Moreover, this is the *only* way we use states and the only way we refer to the values of variables in the current state. (We never need to talk about two states at the same time, etc.) So when we are writing down assertions in informal contexts, we can make some simplifications: drop the initial fun st => and write just X instead of st X. For example, instead of
fun st => (mult (st Z) (st Z)) <= x
/\ ~ ((mult (S (st Z)) (S (st Z))) <= x)
we'll write just
(mult Z Z) <= x /\ ~ ((mult (S Z) (S Z)) <= x).
Next, we need a way of specifying -- making a general claim
about -- the behavior of a command. Since we've defined
assertions as a way of making general claims about the
properties of states, and since the behavior of a command is
to transform one state to another, a general claim about a
command takes the following form:
Since we'll be working a lot with Hoare triples, it's useful to have a compact notation:
{{ P }} c {{ Q }}.
(Traditionally, Hoare triples are written {P} c {Q}, but single braces are already used in Coq.)
- "If c is started in a state satisfying assertion P, then when (and if) c terminates, the final state is guaranteed to satisfy the assertion Q."
Since we'll be working a lot with Hoare triples, it's useful to have a compact notation:
{{ P }} c {{ Q }}.
(Traditionally, Hoare triples are written {P} c {Q}, but single braces are already used in Coq.)
Definition hoare_triple (P:Assertion) (c:com) (Q:Assertion) : Prop :=
forall st st',
c / st ~~> st'
-> P st
-> Q st'.
Notation "{{ P }} c {{ Q }}" := (hoare_triple P c Q) (at level 90).
Paraphrase the following Hoare triples in English. (Note
that we're using the informal convention for writing down
assertions.)
{{True}} c {{X = 5}}
{{X = x}} c {{X = plus x 5)}}
{{X <= Y}} c {{Y <= X}}
{{True}} c {{False}}
{{X = x}}
c
{{Y = real_fact x}}.
{{True}}
c
{{(mult Z Z) <= x /\ ~ ((mult (S Z) (S Z)) <= x)}}
{{True}} c {{X = 5}}
{{X = x}} c {{X = plus x 5)}}
{{X <= Y}} c {{Y <= X}}
{{True}} c {{False}}
{{X = x}}
c
{{Y = real_fact x}}.
{{True}}
c
{{(mult Z Z) <= x /\ ~ ((mult (S Z) (S Z)) <= x)}}
☐
Which of the following Hoare triples are valid (i.e., the
claimed relation between P, c, and Q is true)?
{{True}} X ::= A5 {{X = 5}}
{{X = 2}} X ::= !X +++ 1 {{X = 3}}
{{True}} X ::= A5; Y ::= A0 {{X = 5}}
{{X = 2 /\ X = 3}} X ::= A5 {{X = 0}}
{{True}} SKIP {{False}}
{{False}} SKIP {{True}}
{{True}} WHILE BTrue DO SKIP LOOP {{False}}
{{X = 0}} WHILE !X === A0 DO X ::= !X +++ 1 LOOP {{X = 1}}
{{X = 1}}
WHILE BNot (!X === A0) DO X ::= !X +++ 1 LOOP
{{X = 100}}
{{True}} X ::= A5 {{X = 5}}
{{X = 2}} X ::= !X +++ 1 {{X = 3}}
{{True}} X ::= A5; Y ::= A0 {{X = 5}}
{{X = 2 /\ X = 3}} X ::= A5 {{X = 0}}
{{True}} SKIP {{False}}
{{False}} SKIP {{True}}
{{True}} WHILE BTrue DO SKIP LOOP {{False}}
{{X = 0}} WHILE !X === A0 DO X ::= !X +++ 1 LOOP {{X = 1}}
{{X = 1}}
WHILE BNot (!X === A0) DO X ::= !X +++ 1 LOOP
{{X = 100}}
☐
To get us warmed up, here are two simple lemmas about Hoare
triples
Theorem hoare_post_true : forall (P Q : Assertion) c,
(forall st, Q st) ->
{{P}} c {{Q}}.
Proof.
intros P Q c H. unfold hoare_triple.
intros st st' Heval HP.
apply H. Qed.
Theorem hoare_pre_false : forall (P Q : Assertion) c,
(forall st, ~(P st)) ->
{{P}} c {{Q}}.
Proof.
intros P Q c H. unfold hoare_triple.
intros st st' Heval HP.
unfold not in H. apply H in HP.
destruct HP. Qed.
Some Hoare triples are more interesting than others. For
example,
{{ False }} X ::= !Y + A1 {{ X <= 5 }}
is perfectly valid, but it tells us nothing useful -- in particular, it doesn't tell us how we can use the command X ::= !Y + A1 to achieve the postcondition X <= 5.
By contrast,
{{ Y <= 4 /\ Z = 0 }} X ::= !Y + A1 {{ X <= 5 }}
is useful: it tells us that, if we can somehow create a situation in which we know that Y <= 4 /\ Z = 0, then running this command will produce a state satisfying this postcondition. However, this triple is still not as useful as it could be, because the Z = 0 clause in the precondition actually has nothing to do with the postcondition X <= 5. The *most* useful triple (with the same command and postcondition) is this one:
{{ Y <= 4 }} X ::= !Y + A1 {{ X <= 5 }}
We say that Y <= 4 is the weakest precondition of the command X ::= !Y + A1 for the postcondition X <= 5.
In general, we say that "P is the weakest precondition of c for Q" if
{{ False }} X ::= !Y + A1 {{ X <= 5 }}
is perfectly valid, but it tells us nothing useful -- in particular, it doesn't tell us how we can use the command X ::= !Y + A1 to achieve the postcondition X <= 5.
By contrast,
{{ Y <= 4 /\ Z = 0 }} X ::= !Y + A1 {{ X <= 5 }}
is useful: it tells us that, if we can somehow create a situation in which we know that Y <= 4 /\ Z = 0, then running this command will produce a state satisfying this postcondition. However, this triple is still not as useful as it could be, because the Z = 0 clause in the precondition actually has nothing to do with the postcondition X <= 5. The *most* useful triple (with the same command and postcondition) is this one:
{{ Y <= 4 }} X ::= !Y + A1 {{ X <= 5 }}
We say that Y <= 4 is the weakest precondition of the command X ::= !Y + A1 for the postcondition X <= 5.
In general, we say that "P is the weakest precondition of c for Q" if
- {{P}} c {{Q}}, and
- if P' is an assertion such that {{P'}} c {{Q}}, then, for all states st, P' st implies P st.
What are the weakest preconditions of the following commands
for the following postconditions?
{{ ? }} SKIP {{ X = 5 }}
{{ ? }} X ::= !Y +++ !Z {{ X = 5 }}
{{ ? }} X ::= !Y {{ X = Y }}
{{ ? }}
IFB !X === A0
THEN Y ::= !Z +++ A1
ELSE Y ::= !W +++ A2
{{ Y = 5 }}
{{ ? }}
X ::= 5
{{ X = 0 }}
{{ ? }}
WHILE BTrue DO X LOOP ::= 0
{{ X = 0 }}
{{ ? }} SKIP {{ X = 5 }}
{{ ? }} X ::= !Y +++ !Z {{ X = 5 }}
{{ ? }} X ::= !Y {{ X = Y }}
{{ ? }}
IFB !X === A0
THEN Y ::= !Z +++ A1
ELSE Y ::= !W +++ A2
{{ Y = 5 }}
{{ ? }}
X ::= 5
{{ X = 0 }}
{{ ? }}
WHILE BTrue DO X LOOP ::= 0
{{ X = 0 }}
☐
The rule for reasoning about assignment is the most basic of
the Hoare rules... and probably the trickiest! Here's how it
works.
Consider this (valid) Hoare triple:
{{ Y = 1 }} X ::= !Y {{ X = 1 }}
In English: if we start out in a state where the value of Y is 1 and we assign !Y to X, then we'll finish in a state where X is 1. That is, the property of being equal to 1 gets transferred from Y to X.
Similarly, in
{{ Y + Z = 1 }} X ::= !Y +++ !Z {{ X = 1 }}
the same property (being equal to one) gets transferred to X from the expression !Y +++ !Z on the right-hand side of the assignment.
More generally, if a is *any* arithmetic expression, then
{{ a = 1 }} X ::= a {{ X = 1 }}
is a valid Hoare triple.
Even more generally, if P is *any* property of numbers and a is any arithmetic expression, then
{{ P(a) }} X ::= a {{ P(X) }}
is a valid Hoare triple.
Rephrasing this a bit leads us to the general Hoare rule for assignment:
{{ P where a is substituted for X }} X ::= a {{ P }}
This rule looks backwards to everyone at first -- what it's saying is that, if P holds in an environment where X is replaced by the value of a, then P still holds after executing X ::= a.
For example, these are valid applications of the assignment rule:
{{ X + 1 <= 5 }} X ::= !X +++ A1 {{ X <= 5 }}
{{ 3 = 3 }} X ::= A3 {{ X = 3 }}
{{ 0 <= 3 /\ 3 <= 5 }} X ::= A3 {{ 0 <= X /\ X <= 5 }}
Consider this (valid) Hoare triple:
{{ Y = 1 }} X ::= !Y {{ X = 1 }}
In English: if we start out in a state where the value of Y is 1 and we assign !Y to X, then we'll finish in a state where X is 1. That is, the property of being equal to 1 gets transferred from Y to X.
Similarly, in
{{ Y + Z = 1 }} X ::= !Y +++ !Z {{ X = 1 }}
the same property (being equal to one) gets transferred to X from the expression !Y +++ !Z on the right-hand side of the assignment.
More generally, if a is *any* arithmetic expression, then
{{ a = 1 }} X ::= a {{ X = 1 }}
is a valid Hoare triple.
Even more generally, if P is *any* property of numbers and a is any arithmetic expression, then
{{ P(a) }} X ::= a {{ P(X) }}
is a valid Hoare triple.
Rephrasing this a bit leads us to the general Hoare rule for assignment:
{{ P where a is substituted for X }} X ::= a {{ P }}
This rule looks backwards to everyone at first -- what it's saying is that, if P holds in an environment where X is replaced by the value of a, then P still holds after executing X ::= a.
For example, these are valid applications of the assignment rule:
{{ X + 1 <= 5 }} X ::= !X +++ A1 {{ X <= 5 }}
{{ 3 = 3 }} X ::= A3 {{ X = 3 }}
{{ 0 <= 3 /\ 3 <= 5 }} X ::= A3 {{ 0 <= X /\ X <= 5 }}
We could try to formalize the assignment rule directly in Coq
by treating P as a family of assertions indexed by
arithmetic expressions -- something like this:
Theorem hoare_asgn_firsttry :
forall (P : aexp -> Assertion) V a,
{{fun st => P a st}} (V ::= a) {{fun st => P (AId V) st}}.
But this formulation is not very nice, for two reasons. First, it is not clear how we'd prove it is valid (we'd need to somehow reason about all possible propositions). And second, even if we could prove it, it would be awkward to use.
A much smoother way of formalizing the rule arises from the following obervation:
Theorem hoare_asgn_firsttry :
forall (P : aexp -> Assertion) V a,
{{fun st => P a st}} (V ::= a) {{fun st => P (AId V) st}}.
But this formulation is not very nice, for two reasons. First, it is not clear how we'd prove it is valid (we'd need to somehow reason about all possible propositions). And second, even if we could prove it, it would be awkward to use.
A much smoother way of formalizing the rule arises from the following obervation:
- For all states st, "P where a is substituted for X" holds in the state st if and only if P holds in the state override st V (aeval st a).
Substitution:
The proof rule for assignment:
Theorem hoare_asgn : forall P V a,
{{assn_sub P V a}} (V ::= a) {{P}}.
Proof.
unfold hoare_triple.
intros P V a st st' HE HP.
inversion HE. subst.
unfold assn_sub in HP. assumption. Qed.
Here's a first formal proof using this rule:
Example assn_sub_example :
{{fun st => 3 = 3}} (X ::= A3) {{fun st => st X = 3}}.
Proof.
assert ((fun st => 3 = 3) =
(assn_sub (fun st => st X = 3) X A3)).
Case "Proof of assertion".
unfold assn_sub. reflexivity.
rewrite -> H. apply hoare_asgn. Qed.
Unfortunately, the hoare_asgn rule doesn't literally apply
to the initial goal: it only works with triples whose
precondition has precisely the form assn_sub P V a for some
P, V, and a. So we start with asserting a little lemma
to get the goal into this form.
Doing this kind of fiddling with the goal state every time we
want to use hoare_asgn would gets tiresome pretty quickly.
Fortunately, there are easier alternatives. One simple one is
to state a slightly more general theorem that introduces an
explicit equality hypothesis:
Theorem hoare_asgn_eq : forall P P' V a,
P' = assn_sub P V a
-> {{P'}} (V ::= a) {{P}}.
Proof.
intros P P' V a H. rewrite H. apply hoare_asgn. Qed.
With this version of hoare_asgn, we can do the proof much
more smoothly.
Example assn_sub_example' :
{{fun st => 3 = 3}} (X ::= A3) {{fun st => st X = 3}}.
Proof.
apply hoare_asgn_eq. reflexivity. Qed.
Translate these informal Hoare triples:
{{ X + 1 <= 5 }} X ::= !X +++ A1 {{ X <= 5 }}
{{ 0 <= 3 /\ 3 <= 5 }} X ::= A3 {{ 0 <= X /\ X <= 5 }}
into formal statements and use hoare_asgn_eq to prove them.
{{ X + 1 <= 5 }} X ::= !X +++ A1 {{ X <= 5 }}
{{ 0 <= 3 /\ 3 <= 5 }} X ::= A3 {{ 0 <= X /\ X <= 5 }}
into formal statements and use hoare_asgn_eq to prove them.
(* FILL IN HERE *)
☐
If the assignment rule still seems "backward", it may help to
think a little about alternative "forward" rules. Here is a
seemingly natural one:
{{ True }} X ::= a {{ X = a }}
Explain what is wrong with this rule.
(* FILL IN HERE *)
{{ True }} X ::= a {{ X = a }}
Explain what is wrong with this rule.
(* FILL IN HERE *)
☐
(* Check hoare_asgn_eq. *)
The above discussion about the awkwardness of applying the
assignment rule illustrates a more general point: sometimes
the preconditions and postconditions we get from the Hoare
rules won't quite be the ones we want -- they may (as in the
above example) be logically equivalent but have a different
syntactic form that fails to unify with the goal we are
trying to prove, or they actually may be logically weaker or
stronger than the goal.
For instance, while
{{3 = 3}} (X ::= A3) {{X = 3}},
is a valid Hoare triple, what we probably have in mind when we think about the effect of this assignment is something more like this:
{{True}} (X ::= A3) {{X = 3}}.
This triple is also valid, but we can't derive it from hoare_asgn (or hoare_asgn_eq) because True and 3 = 3 are not equal, even after simplification.
In general, if we can derive {{P}} c {{Q}}, it is valid to change P to P' as long as P' is still enough to show P, and change Q to Q' as long as Q is enough to show Q'.
This observation is captured by the following Rule of Consequence.
For instance, while
{{3 = 3}} (X ::= A3) {{X = 3}},
is a valid Hoare triple, what we probably have in mind when we think about the effect of this assignment is something more like this:
{{True}} (X ::= A3) {{X = 3}}.
This triple is also valid, but we can't derive it from hoare_asgn (or hoare_asgn_eq) because True and 3 = 3 are not equal, even after simplification.
In general, if we can derive {{P}} c {{Q}}, it is valid to change P to P' as long as P' is still enough to show P, and change Q to Q' as long as Q is enough to show Q'.
This observation is captured by the following Rule of Consequence.
{{P'}} c {{Q'}} P implies P' (in every state) Q' implies Q (in every state) ----------------------------- {{P}} c {{Q}}For convenience, here are two more consequence rules -- one for situations where we want to just strengthen the precondition, and one for when we want to just loosen the postcondition.
{{P'}} c {{Q}} P implies P' (in every state) ----------------------------- {{P}} c {{Q}} {{P}} c {{Q'}} Q' implies Q (in every state) ----------------------------- {{P}} c {{Q}}
Theorem hoare_consequence : forall (P P' Q Q' : Assertion) c,
{{P'}} c {{Q'}} ->
(forall st, P st -> P' st) ->
(forall st, Q' st -> Q st) ->
{{P}} c {{Q}}.
Proof.
unfold hoare_triple.
intros P P' Q Q' c Hht HPP' HQ'Q.
intros st st' Hc HP.
apply HQ'Q. apply (Hht st st'). try assumption.
apply HPP'. assumption. Qed.
Theorem hoare_consequence_pre : forall (P P' Q : Assertion) c,
{{P'}} c {{Q}} ->
(forall st, P st -> P' st) ->
{{P}} c {{Q}}.
Proof.
intros P P' Q c Hhoare Himp.
apply hoare_consequence with (P' := P') (Q' := Q);
try assumption.
intros st H. apply H. Qed.
Theorem hoare_consequence_post : forall (P Q Q' : Assertion) c,
{{P}} c {{Q'}} ->
(forall st, Q' st -> Q st) ->
{{P}} c {{Q}}.
Proof.
intros P Q Q' c Hhoare Himp.
apply hoare_consequence with (P' := P) (Q' := Q');
try assumption.
intros st H. apply H. Qed.
For example, we might use (the "_pre" version of) the
consequence rule like this:
{{ True }} =>
{{ 1 = 1 }}
X ::= A1
{{ X = 1 }}
Or, formally...
{{ True }} =>
{{ 1 = 1 }}
X ::= A1
{{ X = 1 }}
Or, formally...
Example hoare_asgn_example1 :
{{fun st => True}} (X ::= A1) {{fun st => st X = 1}}.
Proof.
apply hoare_consequence_pre
with (P' := assn_sub (fun st => st X = 1) X A1).
apply hoare_asgn.
intros st H. unfold assn_sub. apply override_eq. Qed.
This is a good moment to introduce another convenient feature
of Coq. Having to write P' explicitly in the example above
is a bit annoying because the very next thing we are going to
do -- applying the hoare_asgn rule -- is going to determine
exactly what it should be. We can use eapply instead of
apply to tell Coq, essentially, "The missing part is going
to be filled in later."
Example hoare_asgn_example1' :
{{fun st => True}} (X ::= A1) {{fun st => st X = 1}}.
Proof.
eapply hoare_consequence_pre.
apply hoare_asgn.
intros st H. unfold assn_sub.
apply override_eq. Qed.
In general, eapply H tactic works just like apply H
except that, instead of failing if unifying the goal with the
conclusion of H does not determine how to instantiate all
of the variables appearing in the premises of H, eapply H
will replace these variables with existential variables
(written ?nnn) as placeholders for expressions that will be
determined (by further unification) later in the proof.
There is also an eassumption tactic that works similarly.
There is also an eassumption tactic that works similarly.
Example hoare_asgn_example2 : forall a n,
{{fun st => aeval st a = n}} (X ::= a)
{{fun st => st X = n}}.
Proof.
intros a n.
eapply hoare_consequence_pre. apply hoare_asgn.
intros st H. unfold assn_sub. subst. apply override_eq. Qed.
Informally:
{{ True }} =>
{{ a = a }}
X ::= a
{{ X = a }}
{{ True }} =>
{{ a = a }}
X ::= a
{{ X = a }}
Now we come to some easier rules...
Since SKIP doesn't change the state, it preserves any
property P:
{{ P }} SKIP {{ P }}
{{ P }} SKIP {{ P }}
Theorem hoare_skip : forall P,
{{P}} SKIP {{P}}.
Proof.
unfold hoare_triple. intros P st st' H HP. inversion H. subst.
assumption. Qed.
{{P}} SKIP {{P}}.
Proof.
unfold hoare_triple. intros P st st' H HP. inversion H. subst.
assumption. Qed.
More interestingly, if the command c1 takes any state where
P holds to a state where Q holds, and if c2 takes any
state where Q holds to one where R holds, then doing c1
followed by c2 will take any state where P holds to one
where R holds:
{{ P }} c1 {{ Q }}
{{ Q }} c2 {{ R }}
---------------------
{{ P }} c1;c2 {{ R }}
{{ P }} c1 {{ Q }}
{{ Q }} c2 {{ R }}
---------------------
{{ P }} c1;c2 {{ R }}
Theorem hoare_seq : forall P Q R c1 c2,
{{Q}} c2 {{R}}
-> {{P}} c1 {{Q}}
-> {{P}} c1;c2 {{R}}.
Proof.
unfold hoare_triple.
intros P Q R c1 c2 H1 H2 st st' H12 Pre.
inversion H12; subst.
apply (H1 st'0 st'); try assumption.
apply (H2 st st'0); try assumption. Qed.
{{Q}} c2 {{R}}
-> {{P}} c1 {{Q}}
-> {{P}} c1;c2 {{R}}.
Proof.
unfold hoare_triple.
intros P Q R c1 c2 H1 H2 st st' H12 Pre.
inversion H12; subst.
apply (H1 st'0 st'); try assumption.
apply (H2 st st'0); try assumption. Qed.
Note that, in the formal rule hoare_seq, the premises are
given in "backwards" order (c2 before c1). This matches
the natural flow of information in many of the situations
where we'll use the rule.
Informally, a nice way of recording a proof using this rule
is as a "decorated program" where the intermediate assertion
Q is written between c1 and c2:
{{ a = n }}
X ::= a;
{{ X = n }} <---- Q
SKIP
{{ X = n }}
{{ a = n }}
X ::= a;
{{ X = n }} <---- Q
SKIP
{{ X = n }}
Example hoare_asgn_example3 : forall a n,
{{fun st => aeval st a = n}}
(X ::= a; SKIP)
{{fun st => st X = n}}.
Proof.
intros a n. eapply hoare_seq.
Case "right part of seq".
apply hoare_skip.
Case "left part of seq".
eapply hoare_consequence_pre. apply hoare_asgn.
intros st H. unfold assn_sub. subst. apply override_eq. Qed.
{{fun st => aeval st a = n}}
(X ::= a; SKIP)
{{fun st => st X = n}}.
Proof.
intros a n. eapply hoare_seq.
Case "right part of seq".
apply hoare_skip.
Case "left part of seq".
eapply hoare_consequence_pre. apply hoare_asgn.
intros st H. unfold assn_sub. subst. apply override_eq. Qed.
Translate this decorated program into a formal proof:
{{ True }} =>
{{ 1 = 1 }}
X ::= A1;
{{ X = 1 }} =>
{{ X = 1 /\ 2 = 2 }}
Y ::= A2
{{ X = 1 /\ Y = 2 }}
{{ True }} =>
{{ 1 = 1 }}
X ::= A1;
{{ X = 1 }} =>
{{ X = 1 /\ 2 = 2 }}
Y ::= A2
{{ X = 1 /\ Y = 2 }}
Example hoare_asgn_example4 :
{{fun st => True}} (X ::= A1; Y ::= A2)
{{fun st => st X = 1 /\ st Y = 2}}.
Proof.
(* FILL IN HERE *) Admitted.
{{fun st => True}} (X ::= A1; Y ::= A2)
{{fun st => st X = 1 /\ st Y = 2}}.
Proof.
(* FILL IN HERE *) Admitted.
☐
Write an IMP program c that swaps the values of X and Y
and show (in Coq) that it satisfies the following
specification:
{{X <= Y}} c {{Y <= X}}
{{X <= Y}} c {{Y <= X}}
(* FILL IN HERE *)
☐
Explain why the following proposition can't be proven:
forall (a : aexp) (n : nat),
{{fun st => aeval st a = n}} (X ::= A3; Y ::= a)
{{fun st => st Y = n}}.
forall (a : aexp) (n : nat),
{{fun st => aeval st a = n}} (X ::= A3; Y ::= a)
{{fun st => st Y = n}}.
(* FILL IN HERE *)
☐
What sort of rule do we want for reasoning about conditional
commands? Certainly, if the same assertion Q holds after
executing either branch, then it holds after the whole
conditional. So we might be tempted to write:
{{P}} c1 {{Q}}
{{P}} c2 {{Q}}
--------------------------------
{{P}} IFB b THEN c1 ELSE c2 {{Q}}
However, this is rather weak. For example, using this rule, we cannot show that:
{{True}}
IFB !X === A0
THEN Y ::= A2
ELSE Y ::= !X +++ 1
{{ X <= Y }}
since the rule tells us nothing about the state in which the assignments take place in the then- and else- branches.
But, actually, we can say something more precise. In the "then" branch, we know that the boolean expression b evaluates to true, and in the "else" branch, we know it evaluates to false. Making this information available in the premises of the lemma gives us more information to work with when reasoning about the behavior of c1 and c2 (i.e., the reasons why they establish the postcondtion Q).
{{P /\ b}} c1 {{Q}}
{{P /\ ~b}} c2 {{Q}}
--------------------------------
{{P}} IFB b THEN c1 ELSE c2 {{Q}}
{{P}} c1 {{Q}}
{{P}} c2 {{Q}}
--------------------------------
{{P}} IFB b THEN c1 ELSE c2 {{Q}}
However, this is rather weak. For example, using this rule, we cannot show that:
{{True}}
IFB !X === A0
THEN Y ::= A2
ELSE Y ::= !X +++ 1
{{ X <= Y }}
since the rule tells us nothing about the state in which the assignments take place in the then- and else- branches.
But, actually, we can say something more precise. In the "then" branch, we know that the boolean expression b evaluates to true, and in the "else" branch, we know it evaluates to false. Making this information available in the premises of the lemma gives us more information to work with when reasoning about the behavior of c1 and c2 (i.e., the reasons why they establish the postcondtion Q).
{{P /\ b}} c1 {{Q}}
{{P /\ ~b}} c2 {{Q}}
--------------------------------
{{P}} IFB b THEN c1 ELSE c2 {{Q}}
To interpret this rule formally, we need to do a little work.
Strictly speaking, what we've written -- P /\ b -- is the conjunction of an assertion and a boolean expression, which doesn't typecheck. To fix this, we need a way of formally "lifting" any bexp b to an assertion. We'll write bassn b for the assertion "b evaluates to true in (a given state)."
Strictly speaking, what we've written -- P /\ b -- is the conjunction of an assertion and a boolean expression, which doesn't typecheck. To fix this, we need a way of formally "lifting" any bexp b to an assertion. We'll write bassn b for the assertion "b evaluates to true in (a given state)."
A couple of facts about bassn will be useful in proving the
validity of the proof rule for conditionals.
Lemma bexp_eval_true : forall b st,
beval st b = true -> (bassn b) st.
Proof.
intros b st Hbe.
unfold bassn. assumption. Qed.
Lemma bexp_eval_false : forall b st,
beval st b = false -> ~ ((bassn b) st).
Proof.
intros b st Hbe contra.
unfold bassn in contra.
rewrite -> contra in Hbe. inversion Hbe. Qed.
Now we can formalize the Hoare proof rule for conditionals
(and prove it correct).
Theorem hoare_if : forall P Q b c1 c2,
{{fun st => P st /\ bassn b st}} c1 {{Q}} ->
{{fun st => P st /\ ~(bassn b st)}} c2 {{Q}} ->
{{P}} (CIf b c1 c2) {{Q}}.
Proof.
unfold hoare_triple.
intros P Q b c1 c2 HTrue HFalse st st' HE HP.
inversion HE; subst.
Case "b is true".
apply (HTrue st st').
assumption.
split. assumption.
apply bexp_eval_true. assumption.
Case "b is false".
apply (HFalse st st').
assumption.
split. assumption.
apply bexp_eval_false. assumption. Qed.
Example if_exercise :
{{fun st => True}}
(CIf (!X === A0) (Y ::= A2) (Y ::= !X +++ A1))
{{fun st => st X <= st Y}}.
Proof.
apply hoare_if.
Case "Then".
eapply hoare_consequence_pre. apply hoare_asgn.
unfold bassn, assn_sub, override. simpl. intros.
inversion H.
assert (st X = 0) as Heqz
by (apply beq_nat_eq; apply sym_eq; assumption).
rewrite -> Heqz.
apply le_S; apply le_S; apply le_n.
Case "Else".
eapply hoare_consequence_pre. apply hoare_asgn.
unfold assn_sub, override; simpl; intros.
rewrite plus_comm; simpl; apply le_S; apply le_n.
Qed.
Finally, we need a rule for reasoning about while loops.
Suppose we have a loop
WHILE b DO c LOOP
and we want to find a pre-condition P and a post-condition Q such that
{{P}} WHILE b DO c LOOP {{Q}}
is a valid triple.
First of all, let's think about the case where b is false at the beginning, so that the loop body never executes at all. In this case, the loop behaves like SKIP, so we might be tempted to write
{{P}} WHILE b DO c LOOP {{P}}.
But, as we remarked above for the conditional, we know a little more at the end -- not just P, but also the fact that b is false in the current state. So we can enrich the postcondition a little:
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
What about the case where the loop body DOES get executed? In order to ensure that P holds when the loop finally exits, we certainly need to make sure that the command c guarantees that P holds whenever c is finished. Moreover, since P holds at the beginning of the first execution of c, and since each execution of c re-establishes P when it finishes, we can always assume that P holds at the beginning of c. This leads us to the following rule:
{{P}} c {{P}}
-----------------------------------
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
This is almost the rule we want, but again it can be improved a little: at the beginning of the loop body, we know not only that P holds, but also that the guard b is true in the current state. This gives us a little more information to use in reasoning about c. Here is the final version of the rule:
{{P /\ b}} c {{P}}
-----------------------------------
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
The proposition P is called an invariant.
Suppose we have a loop
WHILE b DO c LOOP
and we want to find a pre-condition P and a post-condition Q such that
{{P}} WHILE b DO c LOOP {{Q}}
is a valid triple.
First of all, let's think about the case where b is false at the beginning, so that the loop body never executes at all. In this case, the loop behaves like SKIP, so we might be tempted to write
{{P}} WHILE b DO c LOOP {{P}}.
But, as we remarked above for the conditional, we know a little more at the end -- not just P, but also the fact that b is false in the current state. So we can enrich the postcondition a little:
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
What about the case where the loop body DOES get executed? In order to ensure that P holds when the loop finally exits, we certainly need to make sure that the command c guarantees that P holds whenever c is finished. Moreover, since P holds at the beginning of the first execution of c, and since each execution of c re-establishes P when it finishes, we can always assume that P holds at the beginning of c. This leads us to the following rule:
{{P}} c {{P}}
-----------------------------------
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
This is almost the rule we want, but again it can be improved a little: at the beginning of the loop body, we know not only that P holds, but also that the guard b is true in the current state. This gives us a little more information to use in reasoning about c. Here is the final version of the rule:
{{P /\ b}} c {{P}}
-----------------------------------
{{P}} WHILE b DO c LOOP {{P /\ ~b}}
The proposition P is called an invariant.
Formalizing the rule in Coq...
Lemma hoare_while : forall P b c,
{{fun st => P st /\ bassn b st}} c {{P}} ->
{{P}} (WHILE b DO c LOOP) {{fun st => P st /\ ~ (bassn b st)}}.
Proof.
unfold hoare_triple.
intros P b c Hhoare st st' He HP.
(* Like we've seen before, we need to reason by induction on He,
because in the loop case its hypotheses talk about the
whole loop instead of just c *)
remember (WHILE b DO c LOOP) as wcom.
(ceval_cases (induction He) Case); try (inversion Heqwcom); subst.
Case "CEWhileEnd".
split. assumption. apply bexp_eval_false in H. apply H.
Case "CEWhileLoop".
(* Putting both assumption and reflexivity here seems to make
the proof work in both Coq 8.1 and Coq 8.2 *)
destruct IHHe2 as [HP' Hbfalse']; try assumption; try reflexivity.
apply (Hhoare st st'); try assumption.
split. assumption. apply bexp_eval_true. assumption.
split; assumption. Qed.
There are only a few places where we actually need to think
when performing a correctness proof in Hoare Logic. (In
fact, the last two are the only places where real creativity
is required.)
Concretely, a decorated program has one of the following forms:
- the intermediate assertion for a sequential composition
c1;c2
- the "inner" assertions P' and Q' in the consequence rules
- the invariant of a while loop
- the outermost pre- and post-conditions
Concretely, a decorated program has one of the following forms:
- null program
{{ P }}
SKIP
{{ P }}
- sequential composition:
{{ P }}
c1;
{{ Q }}
c2
{{ R }}
(where c1 and c2 are decorated programs)
- assignment:
{{ P where a is substituted for V }}
V ::= a
{{ P }}
- conditional:
{{ P }}
IFB b THEN
{{ P /\ b }}
c1
{{ Q }}
ELSE
{{ P /\ ~b }}
c2
{{ Q }}
- loop:
{{ P }}
WHILE b DO
{{ P /\ b }}
c1
{{ P }}
LOOP
{{ P /\ ~b }}
- consequence:
{{ P }} =>
{{ P' }}
c
{{ Q }}
(where P st implies P' st for all states st)
{{ P }}
c
{{ Q' }} =>
{{ Q }}
(where Q' st implies Q st for all states st)
A simple example:
{{ True }} =>
{{ x = x }}
X ::= ANum x;
{{ X = x }} =>
{{ X = x /\ z = z }}
Z ::= ANum z;
{{ X = x /\ Z = z }} =>
{{ Z - X = z - x }}
WHILE BNot (!X === A0) DO
{{ Z - X = z - x /\ X <> 0 }} =>
{{ (Z - 1) - (X - 1) = z - x }}
Z ::= !Z --- A1;
{{ Z - (X - 1) = z - x }}
X ::= !X --- A1
{{ Z - X = z - x }}
LOOP;
{{ Z - X = z - x /\ ~ (X <> 0) }} =>
{{ Z = z - x }} =>
{{ Z + 1 = z - x + 1 }}
Z ::= !Z + A1
{{ Z = z - x + 1 }}
{{ True }} =>
{{ x = x }}
X ::= ANum x;
{{ X = x }} =>
{{ X = x /\ z = z }}
Z ::= ANum z;
{{ X = x /\ Z = z }} =>
{{ Z - X = z - x }}
WHILE BNot (!X === A0) DO
{{ Z - X = z - x /\ X <> 0 }} =>
{{ (Z - 1) - (X - 1) = z - x }}
Z ::= !Z --- A1;
{{ Z - (X - 1) = z - x }}
X ::= !X --- A1
{{ Z - X = z - x }}
LOOP;
{{ Z - X = z - x /\ ~ (X <> 0) }} =>
{{ Z = z - x }} =>
{{ Z + 1 = z - x + 1 }}
Z ::= !Z + A1
{{ Z = z - x + 1 }}
In this exercise, we've added a new construct to our language
of commands: CRepeat. You will write the evaluation rule
for repeat and add a new hoare logic theorem to the
language for programs involving it.
We recommend that you do this exercise before the ones that follow - it should help solidify your understanding of the material.
We recommend that you do this exercise before the ones that follow - it should help solidify your understanding of the material.
Inductive com : Type :=
| CSkip : com
| CAss : id -> aexp -> com
| CSeq : com -> com -> com
| CIf : bexp -> com -> com -> com
| CWhile : bexp -> com -> com
| CRepeat : com -> bexp -> com
.
repeat behaves like while, but the loop guard is checked
after each execution of the body, with the loop repeating as
long as the guard stays false. Because of this, the body
will always execute at least once. The concrete syntax for
CRepeat c b often looks like:
repeat c until b
repeat c until b
Tactic Notation "com_cases" tactic(first) tactic(c) :=
first;
[ c "CSkip" | c "CAss" | c "CSeq" | c "CIf" | c "CWhile" | c "CRepeat"].
Notation "'SKIP'" :=
CSkip.
Notation "c1 ; c2" :=
(CSeq c1 c2) (at level 80, right associativity).
Notation "V '::=' a" :=
(CAss V a) (at level 60).
Notation "'WHILE' b 'DO' c 'LOOP'" :=
(CWhile b c) (at level 80, right associativity).
Notation "'IF' e1 'THEN' e2 'ELSE' e3" :=
(CIf e1 e2 e3) (at level 80, right associativity).
Notation "'REPEAT' e1 'UNTIL' b2 'END'" :=
(CRepeat e1 b2) (at level 80, right associativity).
first;
[ c "CSkip" | c "CAss" | c "CSeq" | c "CIf" | c "CWhile" | c "CRepeat"].
Notation "'SKIP'" :=
CSkip.
Notation "c1 ; c2" :=
(CSeq c1 c2) (at level 80, right associativity).
Notation "V '::=' a" :=
(CAss V a) (at level 60).
Notation "'WHILE' b 'DO' c 'LOOP'" :=
(CWhile b c) (at level 80, right associativity).
Notation "'IF' e1 'THEN' e2 'ELSE' e3" :=
(CIf e1 e2 e3) (at level 80, right associativity).
Notation "'REPEAT' e1 'UNTIL' b2 'END'" :=
(CRepeat e1 b2) (at level 80, right associativity).
Add new rules for repeat to ceval below. You can use the
rules for while as a guide, but remember that the body of a
repeat should always execute at least once, and that the
loop ends when the guard becomes true. Then update the
ceval_cases tactic to handle these added cases.
Inductive ceval : state -> com -> state -> Prop :=
| CESkip : forall st,
ceval st CSkip st
| CEAss : forall st a1 n V,
aeval st a1 = n ->
ceval st (CAss V a1) (override st V n)
| CESeq : forall c1 c2 st st' st'',
ceval st c1 st' ->
ceval st' c2 st'' ->
ceval st (CSeq c1 c2) st''
| CEIfTrue : forall st st' b1 c1 c2,
beval st b1 = true ->
ceval st c1 st' ->
ceval st (CIf b1 c1 c2) st'
| CEIfFalse : forall st st' b1 c1 c2,
beval st b1 = false ->
ceval st c2 st' ->
ceval st (CIf b1 c1 c2) st'
| CEWhileEnd : forall b1 st c1,
beval st b1 = false ->
ceval st (CWhile b1 c1) st
| CEWhileLoop : forall st st' st'' b1 c1,
beval st b1 = true ->
ceval st c1 st' ->
ceval st' (CWhile b1 c1) st'' ->
ceval st (CWhile b1 c1) st''
(* FILL IN HERE *)
.
Tactic Notation "ceval_cases" tactic(first) tactic(c) := first; [
c "CESkip" | c "CEAss" | c "CESeq" | c "CEIfTrue" | c "CEIfFalse"
| c "CEWhileEnd" | c "CEWhileLoop"
(* FILL IN HERE *)
].
A couple of definitions from above, copied so they use the
new ceval.
Notation "c1 '/' st '~~>' st'" := (ceval st c1 st') (at level 40).
Definition hoare_triple (P:Assertion) (c:com) (Q:Assertion) : Prop :=
forall st st', (c / st ~~> st') -> P st -> Q st'.
Notation "{{ P }} c {{ Q }}" := (hoare_triple P c Q) (at level 90).
Now state and prove a theorem, hoare_repeat, that expresses
an appropriate proof rule for repeat commands. Use
hoare_while as a model.
☐
Formal Hoare-logic proofs often require a bunch of small
algebraic facts, most of which we wouldn't even bother
writing down in an informal proof. This is a bit clunky in
Coq (though it can be made better using some of the
automation features we haven't discussed yet). The best way
to deal with it is to break them out as separate lemmas so
that they don't clutter the main proofs. Here we collect a
bunch of properties that we'll use in the examples below.
Lemma minus_0_n : forall n,
minus 0 n = 0.
Proof.
intros n. induction n as [| n'].
Case "n = 0".
reflexivity.
Case "n = S n'".
simpl. apply IHn'.
Qed.
Lemma neq0_not_true__eq_0 : forall x,
negb (beq_nat x 0) <> true -> x = 0.
Proof.
intros x H. destruct x as [| x'].
Case "x = 0".
reflexivity.
Case "x = S x'".
simpl in H. destruct H. reflexivity.
Qed.
Lemma neq0_true__eq_S : forall n,
negb (beq_nat n 0) = true
-> exists n', n = S n'.
Proof.
intros n H. destruct n as [| n'].
Case "n = 0 (contradiction)".
inversion H.
Case "n = S n'".
exists n'. reflexivity.
Qed.
Lemma S_le__S: forall n m,
S n <= m
-> exists m', n <= m /\ m = S m'.
Proof.
intros n m H. destruct m as [| m'].
Case "m = 0 (contradiction)".
inversion H.
Case "m = S m'".
exists m'. split; try reflexivity.
apply le_S. apply le_S_n. apply H. Qed.
Lemma not_ev_ev_S_gen: forall n,
(~ ev n -> ev (S n)) /\
(~ ev (S n) -> ev (S (S n))).
Proof.
induction n as [| n'].
Case "n = 0".
split; intros H.
SCase "->".
destruct H. apply ev_0.
SCase "<-".
apply ev_SS. apply ev_0.
Case "n = S n'".
destruct IHn' as [Hn HSn]. split; intros H.
SCase "->".
apply HSn. apply H.
SCase "<-".
apply ev_SS. apply Hn. intros contra.
destruct H. apply ev_SS. apply contra. Qed.
Lemma not_ev_ev_S : forall n,
~ ev n -> ev (S n).
Proof.
intros n.
destruct (not_ev_ev_S_gen n) as [H _].
apply H.
Qed.
A very simple while loop that needs no invariant:
A decorated version of this program, establishing a simple
specification:
{{ True }}
WHILE BNot (!X === A0) DO
{{ True /\ X <> 0 }} =>
{{ True }}
X ::= !X --- A1.
{{ True }}
LOOP
{{ True /\ X = 0 }} =>
{{ X = 0 }}
Your job: translate this proof to Coq. (It may help to look at the slow_subtraction proof below for ideas.)
{{ True }}
WHILE BNot (!X === A0) DO
{{ True /\ X <> 0 }} =>
{{ True }}
X ::= !X --- A1.
{{ True }}
LOOP
{{ True /\ X = 0 }} =>
{{ X = 0 }}
Your job: translate this proof to Coq. (It may help to look at the slow_subtraction proof below for ideas.)
Theorem reduce_to_zero_correct :
{{fun st => True}}
reduce_to_zero
{{fun st => st X = 0}}.
Proof.
(* FILL IN HERE *) Admitted.
{{fun st => True}}
reduce_to_zero
{{fun st => st X = 0}}.
Proof.
(* FILL IN HERE *) Admitted.
☐
Definition subtract_slowly : com :=
WHILE BNot (!X === A0) DO
Z ::= !Z --- A1;
X ::= !X --- A1
LOOP.
The decorated program:
{{ X = x /\ Z = z }} =>
{{ Z - X = z - x }}
WHILE BNot (!X === A0) DO
{{ Z - X = z - x /\ X <> 0 }} =>
{{ (Z - 1) - (X - 1) = z - x }}
Z ::= !Z --- A1;
{{ Z - (X - 1) = z - x }}
X ::= !X --- A1
{{ Z - X = z - x }}
LOOP
{{ Z - X = z - x /\ ~ (X <> 0) }} =>
{{ Z = z - x }}
{{ X = x /\ Z = z }} =>
{{ Z - X = z - x }}
WHILE BNot (!X === A0) DO
{{ Z - X = z - x /\ X <> 0 }} =>
{{ (Z - 1) - (X - 1) = z - x }}
Z ::= !Z --- A1;
{{ Z - (X - 1) = z - x }}
X ::= !X --- A1
{{ Z - X = z - x }}
LOOP
{{ Z - X = z - x /\ ~ (X <> 0) }} =>
{{ Z = z - x }}
Definition subtract_slowly_invariant x z :=
fun st => minus (st Z) (st X) = minus z x.
Theorem subtract_slowly_correct : forall x z,
{{fun st => st X = x /\ st Z = z}}
subtract_slowly
{{fun st => st Z = minus z x}}.
Proof.
(* Note that we do NOT unfold the definition of hoare_triple
anywhere in this proof! The goal is to use only the hoare
rules. Your proofs should do the same. *)
intros x z. unfold subtract_slowly.
(* First we need to transform the pre and postconditions so
that hoare_while will apply. In particular, the
precondition should be the loop invariant. *)
eapply hoare_consequence with (P' := subtract_slowly_invariant x z).
apply hoare_while.
Case "Loop body preserves invariant".
(* Split up the two assignments with hoare_seq - using eapply
lets us solve the second one immediately with hoare_asgn *)
eapply hoare_seq. apply hoare_asgn.
(* Now for the first assignment, transform the precondition
so we can use hoare_asgn *)
eapply hoare_consequence_pre. apply hoare_asgn.
(* Finally, we need to justify the implication generated by
hoare_consequence_pre (this bit of reasoning is elided in the
informal proof) *)
unfold subtract_slowly_invariant, assn_sub, override, bassn. simpl.
intros st [Inv GuardTrue]. apply neq0_true__eq_S in GuardTrue.
destruct GuardTrue as [n' Hn'].
rewrite Hn'. simpl. rewrite Hn' in Inv.
omega.
Case "Initial state satisfies invariant".
(* This is the subgoal generated by the precondition part of our
first use of hoare_consequence. It's the first implication
written in the decorated program (though we elided the actual
proof there). *)
unfold subtract_slowly_invariant.
intros st [HX HZ].
rewrite <- HX. rewrite HZ. reflexivity.
Case "Invariant and negated guard imply postcondition".
(* This is the subgoal generated by the postcondition part of
out first use of hoare_consequence. This implication is
the one written after the while loop in the informal proof. *)
intros st [Inv GuardFalse].
unfold subtract_slowly_invariant in Inv.
unfold bassn in GuardFalse. simpl in GuardFalse.
apply neq0_not_true__eq_0 in GuardFalse.
Here's another good use of the omega tactic.
omega. Qed.
Your turn...
The following program adds the variable X into the variable Z
by repeatedly decrementing X and incrementing Z.
Following the pattern of the subtract_slowly example above,
pick a precondition and postcondition that give an
appropriate specification of add_slowly, then decorate the
program accordingly.
(* FILL IN HERE *)
☐
Now write down your specification of add_slowly as a Coq
Hoare_triple and prove it valid.
(* FILL IN HERE *)
☐
Here's another, slightly trickier example. Make sure you
understand the decorated program completely. (Understanding
all the details of the Coq proof is not required, though it
is not actually very hard -- all the required ideas are
already in the informal version, and all the required
miscellaneous facts about arithmetic are recorded above.)
Lemma minus_pred : forall n m,
(exists m', m = S m') ->
m <= n ->
minus n (pred m) = S (minus n m).
Proof.
intros n m Hnot0.
destruct Hnot0 as [m' Heq]. rewrite Heq in *.
clear Heq m. revert n.
induction m' as [| m'']; intros n Hle.
Case "m' = 0".
simpl. destruct n. inversion Hle. omega.
Case "m' = S m''".
simpl.
assert (S m'' <= pred n).
SCase "proof of assertion".
destruct n. inversion Hle. apply le_S_n in Hle.
simpl. assumption.
apply IHm'' in H. clear IHm''.
This looks pretty awful until you just wave the omega wand.
omega.
Qed.
Definition find_parity : com :=
Y ::= A0;
WHILE (BNot (!X === A0)) DO
Y ::= A1 --- !Y;
X ::= !X --- A1
LOOP.
Qed.
Definition find_parity : com :=
Y ::= A0;
WHILE (BNot (!X === A0)) DO
Y ::= A1 --- !Y;
X ::= !X --- A1
LOOP.
Proof of correctness:
{{ X = x }} =>
{{ X = x /\ 0 = 0 }}
Y ::= A0;
{{ X = x /\ Y = 0 }} =>
{{ (Y=0 <-> ev (x-X)) /\ X<=x }}
WHILE BNot (!X === A0) DO
{{ (Y=0 <-> ev (x-X)) /\ X<=x /\ X<>0 }} =>
{{ (1-Y)=0 <-> ev (x-(X-1)) /\ X-1<=x }}
Y ::= A1 --- !Y;
{{ Y=0 <-> ev (x-(X-1)) /\ X-1<=x }}
X ::= !X --- A1
{{ Y=0 <-> ev (x-X) /\ X<=x }}
LOOP
{{ (Y=0 <-> ev (x-X)) /\ X<=x /\ ~(X<>0) }} =>
{{ Y=0 <-> ev x }}
{{ X = x }} =>
{{ X = x /\ 0 = 0 }}
Y ::= A0;
{{ X = x /\ Y = 0 }} =>
{{ (Y=0 <-> ev (x-X)) /\ X<=x }}
WHILE BNot (!X === A0) DO
{{ (Y=0 <-> ev (x-X)) /\ X<=x /\ X<>0 }} =>
{{ (1-Y)=0 <-> ev (x-(X-1)) /\ X-1<=x }}
Y ::= A1 --- !Y;
{{ Y=0 <-> ev (x-(X-1)) /\ X-1<=x }}
X ::= !X --- A1
{{ Y=0 <-> ev (x-X) /\ X<=x }}
LOOP
{{ (Y=0 <-> ev (x-X)) /\ X<=x /\ ~(X<>0) }} =>
{{ Y=0 <-> ev x }}
Definition find_parity_invariant x :=
fun st =>
st X <= x
/\ (st Y = 0 /\ev (x - st X) \/ st Y = 1 /\ ~ev (x - st X)).
Theorem find_parity_correct : forall x,
{{fun st => st X = x}}
find_parity
{{fun st => st Y = 0 <-> ev x}}.
Proof.
intros x. unfold find_parity.
apply hoare_seq with (Q := find_parity_invariant x).
eapply hoare_consequence.
apply hoare_while with (P := find_parity_invariant x).
Case "Loop body preserves invariant".
eapply hoare_seq.
apply hoare_asgn.
eapply hoare_consequence_pre.
apply hoare_asgn.
intros st [[Inv1 Inv2] GuardTrue].
apply neq0_true__eq_S in GuardTrue.
destruct GuardTrue as [X' HX'].
unfold find_parity_invariant, bassn, assn_sub, aeval in *.
rewrite override_eq.
rewrite (override_neq Y X); auto.
rewrite (override_neq X Y); auto.
rewrite override_eq.
rewrite HX' in *; clear HX'.
split. omega.
destruct Inv2 as [[H1 H2] | [H1 H2]]; [right|left]; (split; [omega |]).
replace (S X' - 1) with X' by omega.
apply ev_not_ev_S in H2.
replace (S (x - S X')) with (x-X') in H2 by omega.
auto.
replace (S X' - 1) with X' by omega.
apply not_ev_ev_S in H2.
replace (S (x - S X')) with (x - X') in H2 by omega.
assumption.
Case "Precondition implies invariant".
intros st H. assumption.
Case "Invariant implies postcondition".
unfold bassn, find_parity_invariant. simpl.
intros st [[Inv1 Inv2] GuardFalse].
apply neq0_not_true__eq_0 in GuardFalse.
rewrite GuardFalse in Inv2.
rewrite <- minus_n_O in *.
intuition. (* what a nice tactic! *)
rewrite H0 in H. inversion H.
Case "invariant established before loop".
eapply hoare_consequence_pre.
apply hoare_asgn.
intros st H.
unfold assn_sub, find_parity_invariant, override. simpl.
subst.
split.
apply le_n.
rewrite (minus_diag (st X)).
left. split. reflexivity.
apply ev_0. Qed.
My first attempt at stating an invariant for find_parity
was the following.
Why didn't this work? (Hint: Don't waste time trying to
answer this exercise by attempting a formal proof and seeing
where it goes wrong. Just think about whether the loop body
actually preserves the property.)
(* FILL IN HERE *)
☐
Definition sqrt_loop : com :=
WHILE ((A1 +++ !Z) *** (A1 +++ !Z) <<= !X) DO
Z ::= A1 +++ !Z
LOOP.
Definition sqrt_com : com :=
Z ::= ANum 0;
sqrt_loop.
Definition sqrt_spec (x : nat) : Assertion :=
fun st =>
(mult (st Z) (st Z)) <= x
/\ ~ ((mult (S (st Z)) (S (st Z))) <= x).
Definition sqrt_inv (x : nat) : Assertion :=
fun st =>
st X = x
/\ (mult (st Z) (st Z)) <= x.
Theorem random_fact_1 : forall st,
mult (S (st Z)) (S (st Z)) <= st X
-> bassn (A1 +++ !Z *** A1 +++ !Z <<= !X) st.
Proof.
intros st Hle. unfold bassn. simpl.
destruct (st X) as [|x'].
Case "st X = 0".
inversion Hle.
Case "st X = S x'".
simpl in Hle. apply le_S_n in Hle.
remember (ble_nat (plus (st Z) (mult (st Z) (S (st Z)))) x')
as ble.
destruct ble. reflexivity.
apply sym_eq in Heqble. apply ble_nat_false in Heqble.
unfold not in Heqble. apply Heqble in Hle. destruct Hle.
Qed.
Theorem random_fact_2 : forall st,
bassn (A1 +++ !Z *** A1 +++ !Z <<= !X) st
-> mult (aeval st (A1 +++ !Z)) (aeval st (A1 +++ !Z)) <= st X.
Proof.
intros st Hble. unfold bassn in Hble. simpl in *.
destruct (st X) as [| x'].
Case "st X = 0".
inversion Hble.
Case "st X = S x'".
apply ble_nat_true in Hble. omega. Qed.
Theorem sqrt_com_correct : forall x,
{{fun st => True}} (X ::= ANum x; sqrt_com) {{sqrt_spec x}}.
Proof.
intros x.
apply hoare_seq with (Q := fun st => st X = x).
Case "sqrt_com".
unfold sqrt_com.
apply hoare_seq with (Q := fun st => st X = x /\ st Z = 0).
SCase "sqrt_loop".
unfold sqrt_loop.
eapply hoare_consequence.
apply hoare_while with (P := sqrt_inv x).
SSCase "loop preserves invariant".
eapply hoare_consequence_pre.
apply hoare_asgn. intros st H.
unfold assn_sub. unfold sqrt_inv in *.
destruct H as [[HX HZ] HP]. split.
SSSCase "X is preserved".
rewrite override_neq; auto.
SSSCase "Z is updated correctly".
rewrite (override_eq (aeval st (A1 +++ !Z)) Z st).
subst. apply random_fact_2. assumption.
SSCase "invariant is true initially".
intros st H. destruct H as [HX HZ].
unfold sqrt_inv. split. assumption.
rewrite HZ. simpl. omega.
SSCase "after loop, spec is satisfied".
intros st H. unfold sqrt_spec.
destruct H as [HX HP].
unfold sqrt_inv in HX. destruct HX as [HX Harith].
split. assumption.
intros contra. apply HP. clear HP.
simpl. simpl in contra.
apply random_fact_1. subst x. assumption.
SCase "Z set to 0".
eapply hoare_consequence_pre. apply hoare_asgn.
intros st HX.
unfold assn_sub. split.
rewrite override_neq; auto.
apply override_eq.
Case "assignment of X".
eapply hoare_consequence_pre. apply hoare_asgn.
intros st H.
unfold assn_sub. apply override_eq. Qed.
Write a decorated program corresponding to the above
correctness proof.
(* FILL IN HERE *)
☐
Module Factorial.
Fixpoint real_fact (n:nat) {struct n} : nat :=
match n with
| O => 1
| S n' => mult n (real_fact n')
end.
Recall the factorial Imp program:
Definition fact_body : com :=
Y ::= !Y *** !Z;
Z ::= !Z --- A1.
Definition fact_loop : com :=
WHILE BNot (!Z === A0) DO
fact_body
LOOP.
Definition fact_com : com :=
Z ::= !X;
Y ::= ANum 1;
fact_loop.
Decorate the fact_com program to show that it satisfies the
specification given by the pre and postconditions below.
Just as we have done above, you may omit the algebraic
reasoning about arithmetic, <=, etc., that is (formally)
required by the rule of consequence:
(* FILL IN HERE *)
{{ X = x }}
Z ::= !X;
Y ::= ANum 1;
WHILE BNot (!Z === A0) DO
Y ::= !Y *** !Z;
Z ::= !Z --- A1
LOOP
{{ Y = real_fact x }}
(* FILL IN HERE *)
{{ X = x }}
Z ::= !X;
Y ::= ANum 1;
WHILE BNot (!Z === A0) DO
Y ::= !Y *** !Z;
Z ::= !Z --- A1
LOOP
{{ Y = real_fact x }}
☐
Prove formally that fact_com satisfies this specification,
using your informal proof as a guide. You may want to state
the loop invariant separately (as we did in the examples).
Theorem fact_com_correct : forall x,
{{fun st => st X = x}} fact_com
{{fun st => st Y = real_fact x}}.
Proof.
(* FILL IN HERE *) Admitted.
☐