AdaCore Blog

SPARK 2014 Rationale: Loop Invariants

by Yannick Moy

Formal verification tools like GNATprove rely on two main inputs from programmers:

  • subprogram contracts (preconditions and postconditions)
  • loop invariants

While the first ones are easy to understand (based on the "contract" analogy, in which a subprogram and its caller have mutual obligations), the second ones are not so simple to grasp. The need is the same though: like calls are "opaque" to the formal verification tool, hence the need for contracts on subprograms, loops are also "opaque" to the formal verification tool, hence the need for loop invariants.

Note that static analysis tools on the contrary do not require either contracts or loop invariants. It is due to the difference in technology between static analysis tools and formal verification tools: the first ones do not require annotations, but they are less powerful, leading to more a posteriori manual work (the review of false positives) if one wants to use them as verification tools instead of simply bug-finding tools.

A loop invariant is a special assertion, expressed with a pragma, that is true at each iteration of the loop. It is executed as a regular assertion, but used differently from assertions by the formal verification tool. For example, here is a function Get_Prime searching for the smallest prime number between Low and High, and the loop invariant giving the range of values of J, and expressing that no integer between Low and the current value J is prime:

function Get_Prime (Low, High : Positive) return Natural is
    J : Positive := Low;
  begin
    while J <= High loop
      if Is_Prime (J) then
        return J;
      end if;
      pragma Loop_Invariant
        (J in Low .. High
          and
        (for all K in Low .. J => not Is_Prime (K)));
      J := J + 1;
    end loop;
    return 0;
  end Get_Prime;

The loop invariant states here properties related to the loop index J: because the value of J changes during the loop, the formal verification tool knows about J only the properties that are stated in the loop invariant. On the code above, GNATprove proves the loop invariant in two stages: 

  • It proves first that the loop invariant is true at the first iteration.
  • It proves then that, assuming the loop invariant held at the previous iteration, it still holds at the next iteration.

This strategy looks a lot like the proof by induction that all students learn in school. Here are the two corresponding checks that GNATprove proves:

loopinv.adb:33:7: info: loop invariant initialization proved
loopinv.adb:33:7: info: loop invariant preservation proved

Now, loop invariants in SPARK 2014 are a bit different from the classical "Hoare" loop invariants (invented by C.A.R. Hoare in 1969), as implemented in EiffelJMLACSL or Spec#.

In all these languages, the loop invariant must be true when reaching a loop, each time the loop resumes, and at loop end. If we had adopted this style of loop invariants in SPARK 2014, the code above would have to be written:

function Get_Prime (Low, High : Positive) return Natural is
    J : Positive := Low;
  begin
    while J <= High loop
      pragma Loop_Invariant
        ((if Low <= High then J in Low .. High + 1)
           and
        (for all K in Low .. J - 1 => not Is_Prime (K)));
      if Is_Prime (J) then
        return J;
      end if;
      J := J + 1;
    end loop;
    return 0;
  end Get_Prime;

You can see immediately that the loop invariant gets more complex, because:

  • Low might be greater than High, hence the guard "if Low <= High" before stating the range of J, for the loop invariant to hold when reaching the loop
  • J may end up being greater than High by 1, hence the range for J "Low .. High + 1", to account for the possible highest value at loop end
  • J has been increased before resuming the loop, hence the range for J from Low to J - 1, due to the fact the loop invariant is checked at the beginning of an iteration.

Hence the decision in SPARK 2014 to allow loop invariants anywhere in the loop, so that the user can put it where it is most natural to express. The loop invariant thus needs not hold when reaching the loop, if the loop is never entered, nor does it need to hold when exiting the loop.

In case you wonder if the loop invariant given previously is useful, it allows you to prove the following contract automatically with GNATprove (expressed with contract cases):

function Get_Prime (Low, High : Positive) return Natural with
     Contract_Cases =>
       --  case 1: there is a prime between Low and High
       ((for some J in Low .. High => Is_Prime (J)) =>
          --  the smallest prime greater or equal to Low is returned
          Get_Prime'Result in Low .. High
            and
          Is_Prime (Get_Prime'Result)
            and
          (for all J in Low .. Get_Prime'Result -1 => not Is_Prime (J)),
       --  case 2: there is no prime between Low and High
        (for all J in Low .. High => not Is_Prime (J)) =>
          --  zero is returned
          Get_Prime'Result = 0);

Posted in #Language    #Formal Verification    #SPARK   

About Yannick Moy

Yannick Moy

Yannick Moy is Head of the Static Analysis Unit at AdaCore. Yannick contributes to the development of SPARK, a software source code analyzer aiming at verifying safety/security properties of programs. He frequently talks about SPARK in articles, conferences, classes and blogs (in particular blog.adacore.com). Yannick previously worked on source code analyzers for PolySpace (now The MathWorks) and at Université Paris-Sud.