Static Program Analysis
Static Program Analysis
Static Program Analysis
December 8, 2022
Copyright © 2008–2022 Anders Møller and Michael I. Schwartzbach
Preface v
1 Introduction 1
1.1 Applications of Static Program Analysis . . . . . . . . . . . . . . 1
1.2 Approximative Answers . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Undecidability of Program Correctness . . . . . . . . . . . . . . 6
3 Type Analysis 19
3.1 Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.2 Type Constraints . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.3 Solving Constraints with Unification . . . . . . . . . . . . . . . . 25
3.4 Record Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.5 Limitations of the Type Analysis . . . . . . . . . . . . . . . . . . 33
4 Lattice Theory 37
4.1 Motivating Example: Sign Analysis . . . . . . . . . . . . . . . . . 37
4.2 Lattices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4.3 Constructing Lattices . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.4 Equations, Monotonicity, and Fixed Points . . . . . . . . . . . . . 43
i
ii CONTENTS
6 Widening 79
6.1 Interval Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
6.2 Widening and Narrowing . . . . . . . . . . . . . . . . . . . . . . 82
8 Interprocedural Analysis 99
8.1 Interprocedural Control Flow Graphs . . . . . . . . . . . . . . . 99
8.2 Context Sensitivity . . . . . . . . . . . . . . . . . . . . . . . . . . 103
8.3 Context Sensitivity with Call Strings . . . . . . . . . . . . . . . . 104
8.4 Context Sensitivity with the Functional Approach . . . . . . . . . 107
Bibliography 197
Preface
Static program analysis is the art of reasoning about the behavior of computer
programs without actually running them. This is useful not only in optimizing
compilers for producing efficient code but also for automatic error detection
and other tools that can help programmers. A static program analyzer is a pro-
gram that reasons about the behavior of other programs. For anyone interested
in programming, what can be more fun than writing programs that analyze
programs?
As known from Turing and Rice, all nontrivial properties of the behavior
of programs written in common programming languages are mathematically
undecidable. This means that automated reasoning of software generally must
involve approximation. It is also well known that testing, i.e. concretely running
programs and inspecting the output, may reveal errors but generally cannot
show their absence. In contrast, static program analysis can – with the right kind
of approximations – check all possible executions of the programs and provide
guarantees about their properties. One of the key challenges when developing
such analyses is how to ensure high precision and efficiency to be practically
useful. For example, nobody will use an analysis designed for bug finding if
it reports many false positives or if it is too slow to fit into real-world software
development processes.
These notes present principles and applications of static analysis of pro-
grams. We cover basic type analysis, lattice theory, control flow graphs, dataflow
analysis, fixed-point algorithms, widening and narrowing, path sensitivity, rela-
tional analysis, interprocedural analysis, context sensitivity, control flow analysis,
several flavors of pointer analysis, and key concepts of semantics-based abstract
interpretation. A tiny imperative programming language with pointers and
first-class functions is subjected to numerous different static analyses illustrating
the techniques that are presented.
We take a constraint-based approach to static analysis where suitable constraint
systems conceptually divide the analysis task into a front-end that generates
constraints from program code and a back-end that solves the constraints to
produce the analysis results. This approach enables separating the analysis
v
vi Preface
specification, which determines its precision, from the algorithmic aspects that
are important for its performance. In practice when implementing analyses, we
often solve the constraints on-the-fly, as they are generated, without representing
them explicitly.
We focus on analyses that are fully automatic (i.e., not involving program-
mer guidance, for example in the form of loop invariants or type annotations)
and conservative (sound but incomplete), and we only consider Turing com-
plete languages (like most programming languages used in ordinary software
development).
The analyses that we cover are expressed using different kinds of constraint
systems, each with their own constraint solvers:
• term unification constraints, with an almost-linear union-find algorithm,
• conditional subset constraints, with a cubic-time algorithm, and
Introduction
Static program analysis aims to automatically answer questions about the possi-
ble behaviors of programs. In this chapter, we explain why this can be useful
and interesting, and we discuss the basic characteristics of analysis tools.
Analysis for program correctness The most successful analysis tools that have
been designed to detect errors (or verify absence of errors) target generic cor-
rectness properties that apply to most or all programs written in specific pro-
gramming languages. In unsafe languages like C, such errors sometimes lead to
critical security vulnerabilities. In more safe languages like Java, such errors are
typically less severe, but they can still cause program crashes. Examples of such
properties are:
• Does there exist an input that leads to a null pointer dereference, division-
by-zero, or arithmetic overflow?
• Are all variables initialized before they are read?
• Are arrays always accessed within their bounds?
• Can there be dangling references, i.e., use of pointers to memory that has
been freed?
• Does the program terminate on every input? Even in reactive systems such
as operating systems, the individual software components, for example
device driver routines, are expected to always terminate.
Other correctness properties depend on specifications provided by the program-
mer for the individual programs (or libraries), for example:
• Are all assertions guaranteed to succeed? Assertions express program
specific correctness properties that are supposed to hold in all executions.
• Is function hasNext always called before function next, and is open always
called before read? Many libraries have such so-called typestate correctness
properties.
• Does the program throw an ActivityNotFoundException or a
SQLiteException for some input?
With web and mobile software, information flow correctness properties have
become extremely important:
• Can input values from untrusted users flow unchecked to file system
operations? This would be a violation of integrity.
• Can secret information become publicly observable? Such situations are
violations of confidentiality.
The increased use of concurrency (parallel or distributed computing) and event-
driven execution models gives rise to more questions about program behavior:
while (n > 1) {
if (n % 2 == 0) // if n is even, divide it by two
n = n / 2;
else // if n is odd, multiply by three and add one
n = 3 * n + 1;
}
In 1937, Collatz conjectured that the answer is “yes”. As of 2020, the conjecture
has been checked for all inputs up to 268 , but nobody has been able to prove it
for all inputs [Roo19].
Even straight-line programs can be difficult to reason about. Does the fol-
lowing program output true for some integer inputs?
This was an open problem since 1954 until 2019 when the answer was found
after over a million hours of computing [BS19].
Rice’s theorem [Ric53] is a general result from 1953 which informally states
that all interesting questions about the input/output behavior of programs
(written in Turing-complete programming languages1 ) are undecidable. This is
easily seen for any special case. Assume for example the existence of an analyzer
that decides if a variable in a program has a constant value in any execution. In
other words, the analyzer is a program A that takes as input a program T , one
of T ’s variables x, and some value k, and decides whether or not x’s value is
always equal to k whenever T is executed.
A yes
(T, x, k)
Is the value of variable x
always equal to k when
T is executed?
no
We could then exploit this analyzer to also decide the halting problem by
using as input the following program where TM(j) simulates the j’th Turing
machine on empty input:
x = 17; if (TM(j)) x = 18;
Here x has a constant value 17 if and only if the j’th Turing machine does not
halt on empty input. If the hypothetical constant-value analyzer A exists, then
we have a decision procedure for the halting problem, which is known to be
impossible [Tur37].
At first, this seems like a discouraging result, however, this theoretical result
does not prevent approximative answers. While it is impossible to build an
analysis that would correctly decide a property for any analyzed program, it is
often possible to build analysis tools that give useful answers for most realistic
programs. As the ideal analyzer does not exist, there is always room for building
more precise approximations (which is colloquially called the full employment
theorem for static program analysis designers).
Approximative answers may be useful for finding bugs in programs, which
may be viewed as a weak form of program verification. As a case in point,
consider programming with pointers in the C language. This is fraught with
dangers such as null dereferences, dangling pointers, leaking memory, and
unintended aliases. Ordinary compilers offer little protection from pointer errors.
Consider the following small program which may perform every kind of error:
int main(int argc, char *argv[]) {
if (argc == 42) {
char *p,*q;
p = NULL;
printf("%s",p);
1 From this point on, we only consider Turing complete languages.
1.2 APPROXIMATIVE ANSWERS 5
q = (char *)malloc(100);
p = q;
free(q);
*p = ’x’;
free(p);
p = (char *)malloc(100);
p = (char *)malloc(100);
q = p;
strcat(p,q);
assert(argc > 87);
}
}
Standard compiler tools such as gcc -Wall detect no errors in this program.
Finding the errors by testing might miss the errors (for this program, no errors
are encountered unless we happen to have a test case that runs the program
with exactly 42 arguments). However, if we had even approximative answers
to questions about null values, pointer targets, and branch conditions then
many of the above errors could be caught statically, without actually running
the program.
Exercise 1.1: Describe all the pointer-related errors in the above program.
Ideally, the approximations we use are conservative (or safe), meaning that all
errors lean to the same side, which is determined by our intended application.
As an example, approximating the memory usage of programs is conservative if
the estimates are never lower than what is actually possible when the programs
are executed. Conservative approximations are closely related to the concept
of soundness of program analyzers. We say that a program analyzer is sound if
it never gives incorrect results (but it may answer maybe). Thus, the notion of
soundness depends on the intended application of the analysis output, which
may cause some confusion. For example, a verification tool is typically called
sound if it never misses any errors of the kinds it has been designed to detect, but
it is allowed to produce spurious warnings (also called false positives), whereas
an automated testing tool is called sound if all reported errors are genuine, but
it may miss errors.
Program analyses that are used for optimizations typically require soundness.
If given false information, the optimization may change the semantics of the
program. Conversely, if given trivial information, then the optimization fails to
do anything.
Consider again the problem of determining if a variable has a constant value.
If our intended application is to perform constant propagation optimization,
then the analysis may only answer yes if the variable really is a constant and
must answer maybe if the variable may or may not be a constant. The trivial
solution is of course to answer maybe all the time, so we are facing the engineering
challenge of answering yes as often as possible while obtaining a reasonable
6 1 INTRODUCTION
analysis performance.
A yes, definitely!
(T, x, k)
Is the value of variable x
always equal to k when
T is executed?
maybe, don’t know
given Turing machine is correct, and it halts in the reject state otherwise. Our
goal is to show that P cannot exist.
If P exists, then we can also build another Turing machine, let us call it M ,
that takes as input the encoding e(T ) of a Turing machine T and then builds the
encoding e(ST ) of yet another Turing machine ST , which behaves as follows:
ST is essentially a universal Turing machine that is specialized to simulate T on
input e(T ). Let w denote the input to ST . Now ST is constructed such that it
simulates T on input e(T ) for at most |w| moves. If the simulation ends in T ’s
accept state, then ST goes to its fail state. It is obviously possible to create ST in
such a way that this is the only way it can reach its fail state. If the simulation does
not end in T ’s accept state (that is, |w| moves have been made, or the simulation
reaches T ’s reject or fail state), then ST goes to its accept state or its reject state
(which one we choose does not matter). This completes the explanation of how
ST works relative to T and w. Note that ST never diverges, and it reaches its fail
state if and only if T accepts input e(T ) after at most |w| moves. After building
e(ST ), M passes it to our hypothetical program analyzer P . Assuming that P
works as promised, it ends in accept if ST is correct, in which case we also let M
halt in its accept state, and in reject otherwise, in which case M similarly halts
in its reject state.
M
accept
accept
e(T ) construct e(ST ) e(ST ) P
from e(T )
reject
reject
We now ask: Does M accept input e(M )? That is, what happens if we run
M with T = M ? If M does accept input e(M ), it must be the case that P
accepts input e(ST ), which in turn means that ST is correct, so its fail state is
unreachable. In other words, for any input w, no matter its length, ST does
not reach its fail state. This in turn means that T does not accept input e(T ).
However, we have T = M , so this contradicts our assumption that M accepts
input e(M ). Conversely, if M rejects input e(M ), then P rejects input e(ST ), so
the fail state of ST is reachable for some input v. This means that there must
exist some w such that the fail state of ST is reached in |w| steps on input v, so
T must accept input e(T ), and again we have a contradiction. By construction
M halts in either accept or reject on any input, but neither is possible for input
e(M ). In conclusion, the ideal program correctness analyzer P cannot exist.
Exercise 1.2: In the above proof, the hypothetical program analyzer P is only
required to correctly analyze programs that always halt. Show how the proof
can be simplified if we want to prove the following weaker property: There
exists no Turing machine P that can decide whether or not the fail state is
reachable in a given Turing machine. (Note that the given Turing machine is
now not assumed to be total.)
Chapter 2
A Tiny Imperative
Programming Language
Basic Expressions
The basic expressions all denote integer values:
Int → 0 | 1 | -1 | 2 | -2 | . . .
Id → x | y | z | . . .
Exp → Int
| Id
| Exp + Exp | Exp - Exp | Exp * Exp | Exp / Exp | Exp > Exp | Exp == Exp
| ( Exp )
10 2 A TINY IMPERATIVE PROGRAMMING LANGUAGE
| input
Expressions Exp include integer literals Int and variables (identifiers) Id. The
input expression reads an integer from the input stream. The comparison
operators yield 0 for false and 1 for true. Function calls, pointer operations, and
record operations will be added later.
Statements
The simple statements Stm are familiar:
Stm → Id = Exp ;
| output Exp ;
| Stm Stm
|
?
| if ( Exp ) { Stm } else { Stm }
| while ( Exp ) { Stm }
?
We use the notation . . . to indicate optional parts. In the conditions we
interpret 0 as false and all other values as true. The output statement writes an
integer value to the output stream.
Functions
A function declaration F contains a function name, a list of parameters, local
variable declarations, a body statement, and a return expression:
?
Fun → Id ( Id , . . . , Id ) { var Id , . . . , Id ; Stm return Exp; }
Function names and parameters are identifiers, like variables. The var block
declares a collection of uninitialized local variables. Function calls are an extra
kind of expression:
Functions as Values
We also allow functions as first-class values. The name of a function can be used
as a kind of variable that refers to the function, and such function values can be
assigned to ordinary variables, passed as arguments to functions, and returned
from functions.
We add a generalized form of function calls (sometimes called computed or
indirect function calls, in contrast to the simple direct calls described earlier):
Unlike simple function calls, the function being called is now an expression
that evaluates to a function value. Function values allow us to illustrate the
main challenges that arise with methods in object-oriented languages and with
higher-order functions in functional languages.
The following example program contains three functions:
twice(f, x) {
return f(f(x));
}
inc(y) {
return y+1;
}
main(z) {
return twice(inc, z);
}
In the main function, the inc function is passed as argument to the twice func-
tion, which calls the given function twice.
Pointers
To be able to build data structures and dynamically allocate memory, we intro-
duce pointers:
Exp → . . .
| alloc Exp
| & Id
| * Exp
| null
The first expression allocates a new cell in the heap initialized with the value of
the given expression and results in a pointer to the cell. The second expression
creates a pointer to a program variable, and the third expression dereferences
a pointer value (this is also called a load operation). In order to assign values
through pointers we allow another form of assignment (called a store operation):
Stm → . . . | * Exp = Exp;
In such an assignment, if the expression on the left-hand-side evaluates to a
pointer to a cell, then the value of the right-hand-side expression is stored in
that cell. Pointers and integers are distinct values, and pointer arithmetic is not
possible.
The following example illustrates the various pointer operations:
x = alloc null;
y = &x;
*x = 42;
z = **y;
12 2 A TINY IMPERATIVE PROGRAMMING LANGUAGE
The first line allocates a cell initially holding the value null, after the second line
y points to the variable x, the third line assigns the value 42 to the cell allocated
in the first line (thereby overwriting the null value), and the fourth line reads
the new value of that cell via two pointer dereferences.
Records
A record is a collection of fields, each having a name and a value. The syntax for
creating records and for reading field values looks as follows:
Exp → . . .
| { Id : Exp , . . . , Id : Exp }
| Exp . Id
Here is an example:
x = {f: 1, g: 2};
y = x.f;
The first line creates a record with two fields: one with name f and value 1, and
one with name g and value 2. The second line reads the value of the f field.
To update the values of record fields we allow two other forms of assign-
ment, for writing directly to a record held by variable or indirectly via a pointer,
respectively:
Stm → . . .
| Id . Id = Exp ;
| ( * Exp ) . Id = Exp;
x = {f: 1, g: 2};
y = &x;
x.f = 3;
(*y).g = 4;
Here, x holds a record, y holds a pointer to the same record, and the two last
assignments update the values of the fields f and g of the record.
Records are passed by value, so, for example, if x holds a record then a simple
assignment z = x; copies the record to z. For simplicity, the values of record
fields cannot themselves be records (but they can be, for example, pointers to
records).
Programs
A complete program is just a collection of functions:
Exercise 2.1: Identify some of the under-specified parts of the TIP language,
and propose meaningful choices to make it more well-defined.
iterate(n) {
var f;
f = 1;
while (n>0) {
f = f*n;
n = n-1;
}
return f;
}
recurse(n) {
var f;
if (n==0) { f=1; }
else { f=n*recurse(n-1); }
return f;
}
foo(p,x) {
var f,q;
if (*p==0) { f=1; }
else {
q = alloc 0;
*q = (*p)-1;
f=(*p)*(x(q,x));
14 2 A TINY IMPERATIVE PROGRAMMING LANGUAGE
}
return f;
}
main() {
var n;
n = input;
return foo(&n,foo);
}
2.3 Normalization
A rich and flexible syntax is useful when writing programs, but when describing
and implementing static analyses, it is often convenient to work with a syntacti-
cally simpler language. For this reason we sometimes normalize programs by
transforming them into equivalent but syntactically simpler ones. A particu-
larly useful normalization is to flatten nested pointer expressions, such that
pointer dereferences are always of the form *Id rather than the more general
*Exp, and similarly, function calls are always of the form Id(Id,. . . ,Id) rather
than Exp(Exp,. . . ,Exp). It may also be useful to flatten arithmetic expressions,
arguments to direct calls, branch conditions, and return expressions.
As an example,
x = f(y+3)*5;
can be normalized to
t1 = y+3;
t2 = f(t1);
x = t2*5;
where t1 and t2 are fresh variables, whereby each statement performs only one
operation.
Exercise 2.2: Argue that any TIP program can be normalized so that all
expressions, with the exception of right-hand side expressions of assignments,
are variables. (This is sometimes called A-normal form [FSDF93].)
TIP uses lexical scoping, however, we make the notationally simplifying as-
sumption that all declared variable and function names are unique in a program,
i.e. that no identifiers is declared more than once.
Exercise 2.5: Argue that any TIP program can be normalized so that all
declared identifiers are unique.
ite
n return
f = ... f
while
var
1 > f = ... n = ...
f
n 0 * −
f n n 1
With this representation, it is easy to extract the set of statements and their
structure for each function in the program.
16 2 A TINY IMPERATIVE PROGRAMMING LANGUAGE
For the sequence S1 S2 , we eliminate the exit node of S1 and the entry node of
S2 and glue the statements together:
S1
S2
Similarly, the other control structures are modeled by inductive graph construc-
tions (sometimes with branch edges labeled with true and false):
2.5 CONTROL FLOW GRAPHS 17
E E E
false true true false false true
S S1 S2 S
Using this systematic approach, the iterative factorial function results in the
following CFG:
var f
f=1
n>0
false true
f=f*n
n=n−1
return f
Exercise 2.6: Draw the AST and the CFG for the rec program from Section 2.2.
Type Analysis
The TIP programming language does not have explicit type declarations, but of
course the various operations are intended to be applied only to certain kinds of
values. Specifically, the following restrictions seem reasonable:
We assume that their violation results in runtime errors. Thus, for a given
program we would like to know that these requirements hold during execution.
Since this is a nontrivial question, we immediately know (Section 1.3) that it is
undecidable.
We resort to a conservative approximation: typability. A program is typable if
it satisfies a collection of type constraints that is systematically derived, typically
from the program AST. The type constraints are constructed in such a way
that the above requirements are guaranteed to hold during execution, but the
converse is not true. Thus, our type analysis will be conservative and reject some
programs that in fact will not violate any requirements during execution.
In most mainstream programming languages with static type checking, the
programmer must provide type annotations for all declared variables and func-
tions. Type annotations serve as useful documentation, and they also make it
20 3 TYPE ANALYSIS
easier to design and implement type systems. TIP does not have type annota-
tions, so our type analysis must infer all the types, based on how the variables
and functions are being used in the program.
Exercise 3.1: Type checking also in mainstream languages like Java may reject
programs that cannot encounter runtime type errors. Give an example of
such a program. To make the exercise nontrivial, every instruction in your
program should be reachable by some input.
Exercise 3.2: Even popular programming languages may have static type
systems that are unsound. Inform yourself about Java’s covariant typing of
arrays. Construct an example Java program that passes all of javac’s type
checks but generates a runtime error due to this covariant typing. (Note that,
because you do receive runtime errors, Java’s dynamic type system is sound,
which is important to avert malicious attacks.)
3.1 Types
We first define a language of types that will describe possible values:
Type → int
| Type
| ( Type , . . . , Type ) → Type
trees defined over the above constructors. (A possibly infinite tree is regular if it
contains only finitely many different subtrees.)
For example, we need infinite types to describe the type of the foo function
from Section 2.2, since the second parameter x may refer to the foo function
itself:
( int,( int,( int,( int,...)→int)→int)→int)→int
To express such recursive types concisely, we add the µ operator and type
variables to the language of types:
Type → . . .
| µ TypeVar . Type
| TypeVar
TypeVar → t | u | . . .
We use meta-variables τ ∈ Type and α ∈ TypeVar, often with subscripts, ranging
over types and type variables. A type of the form µα.τ is considered identical to
the type τ [µα.τ /α].1 With this extra notation, the type of the foo function can
be expressed like this:
µt.( int,t)→int
Exercise 3.3: Explain how regular types can be represented by finite automata
so that two types are equal if their automata accept the same language. Show
an automaton that represents the type µt.( int,t)→int.
We allow free type variables (i.e., type variables that are not bound by an
enclosing µ). Such type variables are implicitly universally quantified, meaning
that they represent any type. Consider for example the following function:
store(a,b) {
*b = a;
return 0;
}
It has type (t, t)→int where t is a free type variable meaning that it can be any
type, which corresponds to the polymorphic behavior of the function. Note
that such type variables are not necessarily entirely unconstrained: the type of a
may be anything, but it must match the type of whatever b points to. The more
restricted type (int, int)→int is also a valid type for the store function, but we
are usually interested in the most general solutions.
Exercise 3.4: What are the types of rec, f, and n in the recursive factorial
program from Section 2.2?
1 Think of a term µα.τ as a quantifier that binds the type variable α in the sub-term τ . An
occurrence of α in a term τ is free if it is not bound by an enclosing µα. The notation τ1 [τ2 /α]
denotes a copy of τ1 where all free occurrences of α have been substituted by τ2 .
22 3 TYPE ANALYSIS
Exercise 3.5: Write a TIP program that contains a function with type
((int)→int)→(int,int)→int.
Type variables are not only useful for expressing recursive types; we also use
them in the following section to express systems of type constraints.
not yield any constraints and that parenthesized expression are not present in
the abstract syntax.
For the program
short() {
var x, y, z;
x = input;
y = alloc x;
*y = x;
z = *y;
return z;
}
[[short]] = ()→[[z]]
[[input]] = int
[[x]] = [[input]]
[[alloc x]] = [[x]]
[[y]] = [[alloc x]]
[[y]] = [[x]]
[[z]] = [[*y]]
[[y]] = [[*y]]
Most of the constraint rules are straightforward. For example, for any syntac-
tic occurrence of E1 ==E2 in the program being analyzed, the two sub-expressions
E1 and E2 must have the same type, and the result is always of type integer.
Exercise 3.6: Explain each of the above type constraint rules, most importantly
those involving functions and pointers.
For a complete program, we add constraints to ensure that the types of the
parameters and the return value of the main function are int:
where c and c0 are term constructors and each τi and τi0 is a sub-term. In the
previous example two of the constraints are [[y]] = [[x]] and [[y]] = [[*y]], so by the
term equality axiom we also have [[x]] = [[*y]].
24 3 TYPE ANALYSIS
τ1 = τ1
τ1 = τ2 =⇒ τ2 = τ1
τ1 = τ2 ∧ τ2 = τ3 =⇒ τ1 = τ3
for all terms τ1 , τ2 , and τ3 .
A solution assigns a type to each type variable, such that all equality con-
straints are satisfied.2 The correctness claim for the type analysis is that the
existence of a solution implies that the specified runtime errors cannot occur
during execution. A solution for the identifiers in the short program is the
following:
[[short]] = ()→int
[[x]] = int
[[y]] = int
[[z]] = int
Exercise 3.8: Give a reasonable definition of what it means for one solution
to be “more general” than another. (See page 21 for an example of two types
where one is more general than the other.)
Exercise 3.9: This exercise demonstrates the importance of the term equality
axiom. First explain what the following TIP code does when it is executed:
var x,y;
x = alloc 1;
y = alloc (alloc 2);
x = y;
Then generate the type constraints for the code, and apply the unification
algorithm (by hand).
Exercise 3.10: Extend TIP with procedures, which, unlike functions, do not
return anything. Show how to extend the language of types and the type
constraint rules accordingly.
2 We can define precisely what we mean by “solution” as follows. A substitution is a mapping σ
from type variables to types. Applying a substitution σ to a type τ , denoted τ σ, means replacing
each free type variable α in τ by σ(α). A solution to a set of type constraints is a substitution σ where
τ1 σ is identical to τ2 σ for each of the type constraints τ1 = τ2 .
3.3 SOLVING CONSTRAINTS WITH UNIFICATION 25
In pseudo-code:
procedure MakeSet(x)
x.parent := x
end procedure
procedure Find(x)
if x.parent 6= x then
x.parent := Find(x.parent)
end if
return x.parent
end procedure
procedure Union(x, y)
xr := Find(x)
y r := Find(y)
if xr 6= y r then
xr .parent := y r
end if
end procedure
3 We here consider a simple version of union-find without union-by-rank; for a description of the
full version with almost-linear worst case time complexity see a textbook on data structures.
26 3 TYPE ANALYSIS
Exercise 3.11: Argue that the unification algorithm works correctly, in the
sense that it finds a solution to the given constraints if one exists. Additionally,
argue that if multiple solutions exist, the algorithm finds the uniquely most
general one (cf. Exercise 3.8).
(The most general solution, when one exists, for a program expression is also
called the principal type of the expression.)
The unification solver only needs to process each constraint once. This means
that although we conceptually first generate the constraints and then solve them,
in an implementation we might as well interleave the two phases and solve the
constraints on-the-fly, as they are being generated.
The complicated factorial program from Section 2.2 generates the following
constraints (duplicates omitted):
[[foo]] = ([[p]],[[x]])→[[f]] [[*p==0]] = int
[[*p]] = int [[f]] = [[1]]
[[1]] = int [[0]] = int
[[p]] = [[*p]] [[q]] = [[alloc 0]]
[[alloc 0]] = [[0]] [[q]] = [[(*p)-1]]
[[q]] = [[*q]] [[*p]] = int
[[f]] = [[(*p)*(x(q,x))]] [[(*p)*(x(q,x))]] = int
[[x(q,x)]] = int [[x]] = ([[q]],[[x]])→[[x(q,x)]]
[[input]] = int [[main]] = ()→[[foo(&n,foo)]]
[[n]] = [[input]] [[&n]] = [[n]]
[[foo]] = ([[&n]],[[foo]])→[[foo(&n,foo)]] [[*p]] = [[0]]
[[(*p)-1]] = int [[foo(&n,foo)]] = int
These constraints have a solution, where most variables are assigned int, except
these:
[[p]] = int
[[q]] = int
[[alloc 0]] = int
[[x]] = µt.( int,t)→int
[[foo]] = µt.( int,t)→int
[[&n]] = int
[[main]] = ()→int
As mentioned in Section 3.1, recursive types are needed for the foo function and
the x parameter. Since a solution exists, we conclude that our program is type
correct.
Exercise 3.12: Check (by hand or using the Scala implementation) that the
constraints and the solution shown above are correct for the complicated
factorial program.
28 3 TYPE ANALYSIS
Exercise 3.13: Consider this fragment of the example program shown earlier:
x = input;
*y = x;
z = *y;
Explain step-by-step how the unification algorithm finds the solution, includ-
ing how the union-find data structure looks in each step.
Recursive types are also required when analyzing TIP programs that manip-
ulate data structures. The example program
var p;
p = alloc null;
*p = p;
creates these constraints:
[[null]] = t
[[alloc null]] = [[null]]
[[p]] = [[alloc null]]
[[p]] = [[p]]
which for [[p]] has the solution [[p]] = µt. t that can be unfolded to [[p]] = ....
Exercise 3.14: Show what the union-find data structure looks like for the
above example program.
Exercise 3.15: Generate and solve the constraints for the iterate example
program from Section 2.2.
3.4 RECORD TYPES 29
Exercise 3.16: Generate and solve the type constraints for this program:
map(l,f,z) {
var r;
if (l==null) r=z;
else r=f(map(*l,f,z));
return r;
}
foo(i) {
return i+1;
}
main() {
var h,t,n;
t = null;
n = 42;
while (n>0) {
n = n-1;
h = alloc null;
*h = t;
t = h;
}
return map(h,foo,0);
}
What is the output from running the program?
(Try to find the solutions manually; you can then use the Scala implementation
to check that they are correct.)
For example, the record type {a:int,b:int} describes records that have two fields,
a and b, both with integer values. Record types with different sets of field names
are considered as different term constructors.
Our goal is to be able to check conservatively whether the different kinds of
errors listed in the beginning of the chapter may appear in the program being
analyzed. This means that we need to distinguish records from other kinds of
values. Specifically, we want the analysis to check that field accesses are only
30 3 TYPE ANALYSIS
performed on records, not on other types of values, and that the fields being
accessed are present.
As a first attempt, the type constraints for record construction and field
lookup can be expressed as follows, inspired by our treatment of pointers and
dereferences. (Field write operations are handled in Exercise 3.19.)
{ X 1 :E1 ,. . . , X n :En }: [[{ X 1 :E1 ,. . . , X n :En }]] = {X1 :[[E1 ]], . . . ,Xn :[[En ]]}
E.X: [[E]] = { . . . ,X:[[E.X]], . . . }
Intuitively, the constraint rule for field lookup says that the type of E must be
a record type that contains a field named X with the same type as E.X. The
right-hand-side of this constraint rule is, however, not directly expressible in our
language of types. One way to remedy this, without requiring any modifications
of our unification algorithm, is to require that every record type contains all
record fields that exist in the program, and add a special type, , representing
absent fields:
Type → . . . |
Let F = {f1 , f2 , . . . , fm } be the set of all field names. We then use the following
two constraint rules instead of the ones above.
Each αi denotes a fresh type variable. The constraints of the form [[E.X]] 6= can
easily be checked separately after the unification process has been completed.
Exercise 3.17: Explain why it is easier to check the constraints of the form
[[E.X]] 6= after instead of during the unification process (i.e., as soon as the
other constraints related to field access operations are processed).
generate the following constraints, assuming that a, b, and c are the only fields
in the program.
3.4 RECORD TYPES 31
The unconstrained type variables x1 and x2 in the solution for [[x]] reflect the
fact that the fields b and c are irrelevant for the field lookup x.a. Similarly,
appears in the solution for [[{b: 42, c: 87}]], because the a field is absent in
that record.
Exercise 3.18: Generate and solve the type constraints for the following pro-
gram:
var a,b,c,d;
a = {f:3, g:17};
b = a.f;
c = {f:alloc 5, h:15};
d = c.f;
What happens if you change c.f to c.g in the last line?
Exercise 3.19: Specify suitable type constraint rules for both forms of field
write statements, X 1 .X 2 =E; and (*E1 ).X=E2 ;. Then generate and solve the
type constraints for the following program:
var a,b,c,d;
a = {f:null, g:17};
b = alloc {g:42, h:87};
a.f = b;
(*b).g = 117;
c = (*(a.f)).g;
Exercise 3.20: As mentioned in Section 2.1, the values of record fields in TIP
programs cannot themselves be records. How can we extend the type analysis
to check whether a given program satisfies this property?
32 3 TYPE ANALYSIS
Exercise 3.21: Assume we extend the TIP language with array operations.
Array values are constructed using a new form of expressions (not to be
confused with the syntax for records):
Exp → . . . | { Exp , . . . ,Exp }
and individual elements are read and written as follows:
a = { 17, 87 };
a[0] = 42;
x = a[0] + a[1]; // x is now 129
Unlike records, arrays are constructed in the heap and passed by reference, so
in the first line, the contents of the array are not copied, and a is like a pointer
to the array containing the two integers.
The type system is extended accordingly with an array type constructor:
Type → . . . | Type []
As an example, the type int[][] denotes arrays of arrays of integers.
Give appropriate type constraints for array operations. Then use the type
analysis to check that the following program is typable and infer the type of
each program variable:
var x,y,z,t;
x = {2,4,8,16,32,64};
y = x[x[3]];
z = {{},x};
t = z[1];
t[2] = y;
3.5 LIMITATIONS OF THE TYPE ANALYSIS 33
Exercise 3.23: Discuss how TIP could be extended with strings and operations
on strings, and how the type analysis could be extended accordingly to check,
for example, that the string operations are only applied to strings and not to
other types of values.
f() {
var x;
x = alloc 17;
x = 42;
return x + 87;
}
This program has no type errors at runtime, but it is rejected by our type checker
because the analysis is flow-insensitive: the order of execution of the program
instructions is abstracted away by the analysis, so intuitively it does not know
that x must be an integer at the return expression. In the following chapters we
shall see how to perform flow-sensitive analysis that does distinguish between
the different program points.
Another limitation, which is even more significant from a practical point of
view, is the current treatment of polymorphic types. In fact, polymorphic types
are not very useful in their current form. Consider this example program:
f(x) {
34 3 TYPE ANALYSIS
return *x;
}
main() {
return f(alloc 1) + *(f(alloc(alloc 2));
}
It never causes an error at runtime but is not typable since it among others
generates constraints equivalent to
which are clearly unsolvable. For this program, we could analyze the f function
first, resulting in this polymorphic type:
[[f]] = ( x)→ x
When analyzing the main function, at each call to f we could then instantiate the
polymorphic type according to the type of the argument: At the first call, the
argument has type int so in this case we treat f as having the type ( int)→int, and
at the second call, the argument has type int so here we treat f as having the
type ( int)→ int. The key property of the program that enables this technique is
the observation that the polymorphic function is not recursive. This idea is called
let-polymorphism (and this is essentially how Damas-Hindley-Milner-style type
analysis actually works in ML and related languages). In Section 8.2 we shall see
a closely related program analysis mechanism called context sensitivity. The price
of the increased precision of let-polymorphism in type analysis is that the worst-
case complexity increases from almost-linear to exponential [KTU90, Mai90].
Even with let-polymorphism, infinitely many other examples will inevitably
remain rejected. An example:
polyrec(g,x) {
var r;
if (x==0) { r=g; } else { r=polyrec(2,0); }
return r+1;
}
main() {
return polyrec(null,1);
}
With functions that are both polymorphic and recursive, type analysis becomes
undecidable in the general case [Hen93, KTU93].
Exercise 3.24: Explain the runtime behavior of the polyrec program, and
why it is unfairly rejected by our type analysis, and why let-polymorphism
does not help.
3.5 LIMITATIONS OF THE TYPE ANALYSIS 35
Yet another limitation of the type system presented in this chapter is that it
ignores many other kinds of runtime errors, such as dereference of null pointers,
reading of uninitialized variables, division by zero, and the more subtle escaping
stack cell demonstrated by this program:
baz() {
var x;
return &x;
}
main() {
var p;
p = baz();
*p = 1;
return *p;
}
The problem in this program is that *p denotes a stack cell that has “escaped”
from the baz function. As we shall see in Section 11.7, such problems can instead
be handled by other kinds of static analysis.
Chapter 4
Lattice Theory
The technique for static analysis that we will study next is based on the mathe-
matical theory of lattices, which we briefly review in this chapter. The connection
between lattices and program analysis was established in the seminal work by
Kildall, Kam and Ullman [Kil73, KU77].
var a,b,c;
38 4 LATTICE THEORY
a = 42;
b = 87;
if (input) {
c = a + b;
} else {
c = a - b;
}
Here, the analysis could conclude that a and b are positive numbers in all possible
executions at the end of the program. The sign of c is either positive or negative
depending on the concrete execution, so the analysis must report > for that
variable.
For this analysis we have an abstract domain consisting of the five abstract
values {+, -, 0, >, ⊥}, which we can organize as follows with the least precise
information at the top and the most precise information at the bottom:
>
- 0 +
The ordering reflects the fact that ⊥ represents the empty set of integer values
and > represents the set of all integer values. Note that > may arise for different
reasons: (1) In the example above, there exist executions where c is positive and
executions where c is negative, so, for this choice of abstract domain, > is the
only sound option. (2) Due to undecidability, imperfect precision is inevitable,
so no matter how we design the analysis there will be programs where, for
example, some variable can only have a positive value in any execution but the
analysis is not able to show that it could not also have a negative value (recall
the TM(j) example from Chapter 1).
The five-element abstract domain shown above is an example of a so-called
complete lattice. We continue the development of the sign analysis in Section 5.1,
but we first need the mathematical foundation in place.
4.2 Lattices
A partial order is a set S equipped with a binary relation v where the following
conditions are satisfied:
• reflexivity: ∀x ∈ S : x v x
• transitivity: ∀x, y, z ∈ S : x v y ∧ y v z =⇒ x v z
• anti-symmetry: ∀x, y ∈ S : x v y ∧ y v x =⇒ x = y
4.2 LATTICES 39
The abstract domain shown in Section 4.1 is an example of a partial order, where
S = {-, 0, +, >, ⊥} and v specifies the ordering of the elements, for example,
⊥ v + and + v >.
When x v y we say that y is a safe approximation of x, or that x is at least as
precise as y. Formally, a partial order is a pair (S, v), but we often use the same
name for the partial order and its underlying set S. Also, we sometimes write
y w x (“y is greater than or equal x”) instead of x v y (“x is less than or equal
y”).
Let X ⊆ S. We say that y ∈ S is an upper bound for X, written X v y, if we
have ∀x ∈ X : x v y. Similarly, y ∈ S is a lower F bound for X, written y v X, if
∀x ∈ X : y v x. A least upper bound, written X, is defined by:
G G
Xv X ∧ ∀y ∈ S : X v y =⇒ Xvy
G X v X ∧ ∀y ∈ S : y v X =⇒ y v G X
Exercise 4.3: Argue that the abstract domain presented in Section 4.1 is a
complete lattice.
Exercise 4.4: Prove that if (S, v) is a nonempty finite lattice (i.e., a lattice
where S is nonempty and finite), then (S, v) is also a complete lattice.
40 4 LATTICE THEORY
It follows from the following exercise that to show that a partial order is
a complete lattice, it suffices to argue that a least upper bound exists for each
subset of its elements; the presence of greatest lower bounds then comes for free
(and the dual property also holds).
Exercise 4.7: Prove that if S is a partially ordered set, then every subset of S
has a least upper boundF if and only if every subset of S has a greatest lower
bound. (Hint: X = {y ∈ S | y v X}.)
F
Every complete lattice has a unique largest element denoted > (pronounced
top) and a unique smallest element denoted ⊥ (pronounced bottom).
F
Exercise
F 4.8: Assume S is a partial order where S and S exist. Prove that
F
S and S are the unique largest element and F
F the unique smallest element,
respectively, in S. In other words, we have > = S and ⊥ = S. F
F F
Exercise 4.9: F where S, S, ∅ and ∅ exist.
FAssume S is a partial order F F
Prove that S = ∅ and that
F F S = ∅. (Together with Exercise 4.8 we
F
then have > = ∅ and ⊥ = ∅.)
F
4.3 CONSTRUCTING LATTICES 41
The height of a lattice is defined to be the length of the longest path from ⊥
to >. As an example, the height of the sign analysis lattice from Section 4.1 is 2.
For some lattices the height is infinite (see Section 6.1).
{0, 1, 2, 3}
The above powerset lattice has height 4. In general, the lattice (P(A), ⊆) has
height |A|. We use powerset lattices in Chapter 5 to represent sets of variables
or expressions.
The reverse powerset lattice for a set A is the complete lattice (P(A), ⊇).
Exercise 4.10: Draw the Hasse diagram of the reverse powerset lattice for the
set {foo, bar, baz}.
>
a1 a2 ···
is a complete lattice with height 2. As an example, the set Sign = {+, -, 0, >, ⊥}
with the ordering described in Section 4.1 forms a complete lattice that can also
be expressed as flat({+, 0, -}).
42 4 LATTICE THEORY
L1 × L2 × . . . × Ln = {(x1 , x2 , . . . , xn ) | xi ∈ Li }
Exercise 4.11: Show that the t and u operators for a product lattice L1 ×
L2 × . . . × Ln can be computed componentwise (i.e. in terms of the t and u
operators from L1 , L2 , . . . , Lk ).
Exercise 4.13: Show that the t and u operators for a map lattice A → L can
be computed pointwise (i.e. in terms of the t and u operators from L).
Exercise 4.14: Show that if A is finite and L has finite height then the height
of the map lattice A → L is height(A → L) = |A| · height(L).
1 We often abuse notation by using the same symbol v for many different order relations, in this
case from the n + 1 different lattices, but it should always be clear from the context which lattice it
belongs to. The same applies to the other operators w, t, u and the top/bottom symbols >, ⊥.
2 The notation [a 7→ x , a 7→ x , . . .] means the function that maps a to x , a to x , etc.
1 1 2 2 1 1 2 2
4.4 EQUATIONS, MONOTONICITY, AND FIXED POINTS 43
var a,b; // 1
a = 42; // 2
b = a + input; // 3
a = a - b; // 4
a1 = >
b1 = >
a2 = +
3 We use the term constraint variable to denote variables that appear in mathematical constraint
systems, to avoid confusion with program variables that appear in TIP programs.
44 4 LATTICE THEORY
b2 = b1
a3 = a2
b3 = a2 + >
a4 = a3 - b3
b4 = b3
For example, a2 denotes the abstract value of a at the program point immediately
after line 2. The operators + and - here work on abstract values, which we return
to in Section 5.1. In this constraint system, the constraint variables have values
from the abstract value lattice Sign defined in Section 4.3. We can alternatively
derive the following equivalent constraint system where each constraint variable
instead has a value from the abstract state lattice StateSigns from Section 4.3:4
x1 = [a 7→ >, b 7→ >]
x2 = x1 [a 7→ +]
x3 = x2 [b 7→ x2 (a) + >]
x4 = x3 [a 7→ x3 (a) - x3 (b)]
Here, each constraint variable models the abstract state at a program point;
for example, x1 models the abstract state at the program point immediately
after line 1. Notice that each equation only depends on preceding ones for this
example program, so in this case the solution can be found by simple substitution.
However, mutually recursive equations may appear, for example for programs
that contain loops (see Section 5.1).
Also notice that it is important for the analysis of this simple program that the
order of statements is taken into account, which is called flow-sensitive analysis.
Specifically, when a is read in line 3, the value comes from the assignment to a
in line 2, not from the one in line 4.
Exercise 4.16: Give a solution to the constraint system above (that is, values
for x1 , . . . , x4 that satisfy the four equations).
Exercise 4.17: Why is the unification solver from Chapter 3 not suitable for
this kind of constraints?
4 The notation f [a 7→ x , . . . , a
1 n n 7→ xn ] means the function that maps ai to xi , for each
i = 1, . . . , n and for all other inputs gives the same output as the function f .
4.4 EQUATIONS, MONOTONICITY, AND FIXED POINTS 45
Exercise 4.23: The operators t and u can be viewed as functions. For example,
x1 t x2 takes as input x1 , x2 ∈ L and returns an element from L as output.
Show that t and u are monotone.
Exercise 4.25: Show that set difference, X\Y , as a function with two argu-
ments over a powerset lattice is monotone in the first argument X but not in
the second argument Y .
We say that x ∈ L is a fixed point for f if f (x) = x. A least fixed point x for f
is a fixed point for f where x v y for every fixed point y for f .
Let L be a complete lattice. An equation system5 over L is of the form
x1 = f1 (x1 , . . . , xn )
x2 = f2 (x1 , . . . , xn )
..
.
xn = fn (x1 , . . . , xn )
where x1 , . . . , xn are variables and f1 , . . . , fn : Ln → L are functions which we
call constraint functions. A solution to an equation system provides a value from
L for each variable such that all equations are satisfied.
We can combine the n functions into one, f : Ln → Ln ,
f (x1 , . . . , xn ) = f1 (x1 , . . . , xn ), . . . , fn (x1 , . . . , xn )
where all constraint are equalities. On page 49 we discuss other forms of constraints.
4.4 EQUATIONS, MONOTONICITY, AND FIXED POINTS 47
x3 = x2 [b 7→ x2 (a) + >]
x4 = x3 [a 7→ x3 (a) - x3 (b)]
Exercise 4.28: Show that the four constraint functions f1 , . . . , f4 are monotone.
(Hint: see Exercise 4.26.)
Exercise 4.29: Argue that your solution from Exercise 4.16 is the least fixed
point of the function f defined by
f (x1 , . . . , x4 ) = f1 (x1 , . . . , x4 ), . . . , f4 (x1 , . . . , x4 ) .
(Note that when applying this theorem to the specific equation system shown
above, f is a function over the product lattice Ln .)
The proof of this theorem is quite simple. Observe that ⊥ v f (⊥) since ⊥
is the least element. Since f is monotone, it follows that f (⊥) v f 2 (⊥) and by
induction that f i (⊥) v f i+1 (⊥) for any i. Thus, we have an increasing chain:
⊥ v f (⊥) v f 2 (⊥) v . . .
Since L is assumed to have finite height, we must for some k have that f k (⊥) =
f k+1 (⊥), i.e. f k (⊥) is a fixed point for f . By Exercise 4.2, f k (⊥) must be the least
upper bound of all elements in the chain, so lfp(f ) = f k (⊥). Assume now that
x is another fixed point. Since ⊥ v x it follows that f (⊥) v f (x) = x, since f is
6 Theordinary Kleene fixed-point theorem requires f to be continuous (see page 168), on the
other hand it does not require L to be a complete lattice but it can be any complete partial order
(possibly with infinite height); see also Exercise 12.8.
48 4 LATTICE THEORY
procedure NaiveFixedPointAlgorithm(f )
x := ⊥
while x 6= f (x) do
x := f (x)
end while
return x
end procedure
(Instead of computing f (x) both in the loop condition and in the loop body, a
trivial improvement is to just compute it once in each iteration and see if the
result changes.) The computation of a fixed point can be illustrated as a walk
up the lattice starting at ⊥:
11111111111
00000000000
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
This algorithm is called “naive” because it does not exploit the special structures
that are common in analysis lattices. We shall see various less naive fixed-point
algorithms in Section 5.3.
The least fixed point is the most precise possible solution to the equation
system, but the equation system is (for a sound analysis) merely a conservative
approximation of the actual program behavior (again, recall the TM(j) example
from Chapter 1). This means that the semantically most precise possible (while
still correct) answer is generally below the least fixed point in the lattice. We shall
see examples of this in Chapter 5 and study the connection to semantics in more
detail in Chapter 12.
4.4 EQUATIONS, MONOTONICITY, AND FIXED POINTS 49
Exercise 4.30: Explain step-by-step how the naive fixed-point algorithm com-
putes the solution to the equation system from Exercise 4.16.
The time complexity of computing a fixed point with this algorithm depends
on
• the height of the lattice, since this provides a bound for the number of
iterations of the algorithm, and
• the cost of computing f (x) and testing equality, which are performed in
each iteration.
Exercise 4.32: Does the fixed-point theorem also hold without the assumption
that the lattice has finite height? If yes, give a proof; if no, give a counterex-
ample. (Hint: see [CC79a].)
x1 w f1 (x1 , . . . , xn )
x2 w f2 (x1 , . . . , xn )
..
.
xn w fn (x1 , . . . , xn )
x1 = x1 t f1 (x1 , . . . , xn )
x2 = x2 t f2 (x1 , . . . , xn )
..
.
xn = xn t fn (x1 , . . . , xn )
x1 v f1 (x1 , . . . , xn )
x2 v f2 (x1 , . . . , xn )
..
.
xn v fn (x1 , . . . , xn )
x1 = x1 u f1 (x1 , . . . , xn )
x2 = x2 u f2 (x1 , . . . , xn )
..
.
xn = xn u fn (x1 , . . . , xn )
Classical dataflow analysis starts with a CFG and a complete lattice with finite
height. The lattice describes abstract information we wish to infer for the different
CFG nodes. It may be fixed for all programs, or it may be parameterized based
on the given program. To every node v in the CFG, we assign a constraint
variable1 [[v]] ranging over the elements of the lattice. For each node we then
define a dataflow constraint that relates the value of the variable of the node to
those of other nodes (typically the neighbors), depending on what construction
in the programming language the node represents. If all the constraints for the
given program happen to be equations or inequations with monotone right-hand
sides, then we can use the fixed-point algorithm from Section 4.4 to compute
the analysis result as the unique least solution.
The combination of a complete lattice and a space of monotone functions
is called a monotone framework [KU77]. For a given program to be analyzed, a
monotone framework can be instantiated by specifying the CFG and the rules
for assigning dataflow constraints to its nodes.
An analysis is sound if all solutions to the constraints correspond to correct
information about the program. The solutions may be more or less imprecise, but
computing the least solution will give the highest degree of precision possible.
We return to the topic of analysis correctness and precision in Chapter 12.
Throughout this chapter we use the subset of TIP without function calls,
pointers, and records; those language features are studied in Chapters 10 and 11
and in Exercise 5.10.
1 As for type analysis, we will ambiguously use the notation [[S]] for [[v]] if S is the syntax associated
with node v. The meaning will always be clear from the context.
52 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
- 0 +
We want an abstract value for each program variable, so we define the map
lattice
States = Vars → Sign
where Vars is the set of variables occurring in the given program. Each element
of this lattice can be thought of as an abstract state, hence its name. For each CFG
node v we assign a constraint variable [[v]] denoting an abstract state that gives
the sign values for all variables at the program point immediately after v. The
lattice States n , where n is the number of CFG nodes, then models information
for all the CFG nodes.
The dataflow constraints model the effects of program execution on the
abstract states. For simplicity, we here focus on a subset of TIP that does not
contain pointers or records, so integers are the only type of values we need to
consider.
First, we define an auxiliary function JOIN (v) that combines the abstract
states from the predecessors of a node v:
G
JOIN (v) = [[w]]
w∈pred(v)
Note that JOIN (v) is a function of all the constraint variables [[v1 ]], . . . , [[vn ]] for
the program. For example, with the following CFG, we have JOIN ([[a=c+2]]) =
[[c=b]] t [[c=-5]].
b > 5
true false
c=b c=−5
a=c+2
5.1 SIGN ANALYSIS, REVISITED 53
The most interesting constraint rule for this analysis is the one for assignment
statements, that is, nodes v of the form X = E:
This constraint rule models the fact that the abstract state after an assignment
X = E is equal to the abstract state immediately before the assignment, except
that the abstract value of X is the result of abstractly evaluating the expression
E. The eval function performs an abstract evaluation of expression E relative to
an abstract state σ:
The function sign gives the sign of an integer constant, and opb is an abstract
evaluation of the given operator,2 defined by the following tables:
b+ ⊥ 0 - + > b- ⊥ 0 - + >
⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥
0 ⊥ 0 - + > 0 ⊥ 0 + - >
- ⊥ - - > > - ⊥ - > - >
+ ⊥ + > + > + ⊥ + + > >
> ⊥ > > > > > ⊥ > > > >
b* ⊥ 0 - + > b/ ⊥ 0 - + >
⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥
0 ⊥ 0 0 0 0 0 ⊥ ⊥ 0 0 >
- ⊥ 0 + - > - ⊥ ⊥ > > >
+ ⊥ 0 - + > + ⊥ ⊥ > > >
> ⊥ 0 > > > > ⊥ ⊥ > > >
b> ⊥ 0 - + > ==
b ⊥ 0 - + >
⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥ ⊥
0 ⊥ 0 + 0 > 0 ⊥ + 0 0 >
- ⊥ 0 > 0 > - ⊥ 0 > 0 >
+ ⊥ + + > > + ⊥ 0 0 > >
> ⊥ > > > > > ⊥ > > > >
For the subset of TIP we have chosen to focus on, no other kinds of CFG
nodes affect the values of variables, so for the remaining nodes we have this
trivial constraint rule:
[[v]] = JOIN (v)
Exercise 5.1: In the CFGs we consider in this chapter (for TIP without function
calls), entry nodes have no predecessors.
(a) Argue that the constraint rule [[v]] = JOIN (v) for such nodes is equiva-
lent to defining [[v]] = ⊥.
(b) Argue that removing all equations of the form [[v]] = ⊥ from an equation
system does not change its least solution.
Exercise 5.3: Check that the above tables indeed define monotone operators
on the Sign lattice.
Exercise 5.4: Argue that these tables are the most precise possible for the
Sign lattice, given that soundness must be preserved. (An informal argument
suffices for now; we shall see a more formal approach to stating and proving
this property in Section 12.4.)
5.1 SIGN ANALYSIS, REVISITED 55
Using the fixed-point algorithm from Section 4.4, we can now obtain the ana-
lysis result for the given program by computing lfp(af ) where af (x1 , . . . , xn ) =
af v1 (x1 , . . . , xn ), . . . , af vn (x1 , . . . , xn ) .
Recall the example program from Section 4.1:
var a,b,c;
a = 42;
b = 87;
if (input) {
c = a + b;
} else {
c = a - b;
}
v1
var a,b,c v2
a = 42 v3
b = 87 v4
input v5
true false
c=a+b v6 c=a−b v7
v8
Exercise 5.6: Generate the equation system for this example program. Then
solve the equations using the fixed-point algorithm from Section 4.4.
(Notice that the least upper bound operation is exactly what we need to model
the merging of information at v8 !)
56 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
Exercise 5.7: Write a small TIP program where the sign analysis leads to an
equation system with mutually recursive constraints. Then explain step-by-
step how the fixed-point algorithm from Section 4.4 computes the solution.
We lose some information in the above analysis, since for example the expres-
sions (2>0)==1 and x-x are analyzed as >, which seems unnecessarily coarse.
(These are examples where the least fixed point of the analysis equation system
is not identical to the semantically best possible answer.) Also, + divided by
+ results in > rather than + since e.g. 1/2 is rounded down to zero. To handle
some of these situations more precisely, we could enrich the sign lattice with
element 1 (the constant 1), 0+ (positive or zero), and 0- (negative or zero) to
keep track of more precise abstract values:
>
0- 0+
- 0 +
Exercise 5.8: Define the six operators on the extended Sign lattice (shown
above) by means of 8 × 8 tables. Check that they are monotone. Does this
new lattice improve precision for the expressions (2>0)==1, x-x, and 1/2?
Exercise 5.9: Show how the eval function could be improved to make the sign
analysis able to show that the final value of z cannot be a negative number in
the following program:
var x,y,z;
x = input;
y = x*x;
z = (x-x+1)*y;
5.2 CONSTANT PROPAGATION ANALYSIS 57
Exercise 5.10: Explain how to extend the sign analysis to handle TIP programs
that use records (see Chapter 2).
One approach, called field insensitive analysis, simply mixes together the
different fields of each record. Another approach, field sensitive analysis,
instead uses a more elaborate lattice that keeps different abstract values for
the different field names.
The results of a sign analysis could in theory be used to eliminate division-
by-zero errors by rejecting programs in which denominator expressions have
sign 0 or >. However, the resulting analysis will probably unfairly reject too
many programs to be practical. Other more powerful analysis techniques, such
as interval analysis (Section 6.1) and path sensitivity (Chapter 7) would be more
useful for detecting such errors. F
Notice that in this analysis we use the operation (in the definition of JOIN ),
but we never use the operation. In fact, when implementing analyses with
F
monotone frameworks, it is common that is ignored entirely even though it
F
mathematically exists.
>
· · · −2 −1 0 1 2 ···
var x,y,z;
x = 27;
y = input;
z = 2*x+y;
if (x < 0) { y = z-3; } else { y = 12; }
output y;
into
var x,y,z;
x = 27;
y = input;
z = 54+y;
if (0) { y = z-3; } else { y = 12; }
output y;
which, following a reaching definitions analysis and dead code elimination (see
Section 5.7), can be reduced to this shorter and more efficient program:
var y;
y = input;
output 12;
This kind of optimization was among the first uses of static program ana-
lysis [Kil73].
procedure RoundRobin(f1 , . . . , fn )
(x1 , . . . , xn ) := (⊥, . . . , ⊥)
while (x1 , . . . , xn ) 6= f (x1 , . . . , xn ) do
for i := 1 . . . n do
xi := fi (x1 , . . . , xn )
end for
end while
return (x1 , . . . , xn )
end procedure
Exercise 5.13: Prove that the round-robin algorithm computes the least fixed
point of f . (Hint: see the proof of the fixed-point theorem, and consider
the ascending chain that arises from the sequence of xi := fi (x1 , . . . , xn )
operations.)
Exercise 5.14: Continuing Exercise 4.30, how many iterations are required by
the naive fixed-point algorithm and the round-robin algorithm, respectively,
to reach the fixed point?
procedure ChaoticIteration(f1 , . . . , fn )
60 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
(x1 , . . . , xn ) := (⊥, . . . , ⊥)
while (x1 , . . . , xn ) 6= f (x1 , . . . , xn ) do
choose i nondeterministically from {1, . . . , n}
xi := fi (x1 , . . . , xn )
end while
return (x1 , . . . , xn )
end procedure
This is not a practical algorithm, because its efficiency and termination depend
on how i is chosen in each iteration. Additionally, naively computing the loop
condition may now be more expensive than executing the loop body. However,
if it terminates, the algorithm produces the right result.
Exercise 5.15: Prove that the chaotic-iteration algorithm computes the least
fixed point of f , if it terminates. (Hint: see your solution to Exercise 5.13.)
In the general case, every constraint variable [[vi ]] may depend on all other
variables. Most often, however, an actual instance of fi will only read the values
of a few other variables, as in the examples from Exercise 4.28 and Exercise 5.6.
We represent this information as a map
which for each node v tells us the subset of other nodes for which [[v]] occurs in
a nontrivial manner on the right-hand side of their dataflow equations. That is,
dep(v) is the set of nodes whose information may depend on the information of
v. We also define its inverse: dep −1 (v) = {w | v ∈ dep(w)}.
For the example from Exercise 5.6, we have, in particular, dep(v5 ) = {v6 , v7 }.
This means that whenever [[v5 ]] changes its value during the fixed-point compu-
tation, only f6 and f7 may need to be recomputed.
Armed with this information, we can present a simple work-list algorithm:
procedure SimpleWorkListAlgorithm(f1 , . . . , fn )
(x1 , . . . , xn ) := (⊥, . . . , ⊥)
W := {v1 , . . . , vn }
while W 6= ∅ do
vi := W.removeNext()
y := fi (x1 , . . . , xn )
if y 6= xi then
xi := y
for each vj ∈ dep(vi ) do
W.add(vj )
end for
end if
5.3 FIXED-POINT ALGORITHMS 61
end while
return (x1 , . . . , xn )
end procedure
The set W is here called the work-list with operations ‘add’ and ‘removeNext’
for adding and (nondeterministically) removing an item. The work-list initially
contains all nodes, so each fi is applied at least once. It is easy to see that the
work-list algorithm terminates on any input: In each iteration, we either move
up in the Ln lattice, or the size of the work-list decreases. As usual, we can
only move up in the lattice finitely many times as it has finite height, and the
while-loop terminates when the work-list is empty. Correctness follows from
observing that each iteration of the algorithm has the same effect on (x1 , . . . , xn )
as one iteration of the chaotic-iteration algorithm for some nondeterministic
choice of i.
Exercise 5.16: Argue that a sound, but probably not very useful choice for the
dep map is one that always returns the set of all CFG nodes.
Exercise 5.17: As stated above, we can choose dep(v5 ) = {v6 , v7 } for the
example equation system from Exercise 5.6. Argue that a good strategy for the
sign analysis is to define dep = succ. (We return to this topic in Section 5.8.)
Exercise 5.18: Explain step-by-step how the work-list algorithm computes the
solution to the equation system from Exercise 5.6. (Since the ‘removeNext’
operation is nondeterministic, there are many correct answers!)
(b) Now let us modify the CFG construction slightly, such that a dummy
“no-op” node is inserted at the merge point after the two branches of
each if block. This will increase the number of CFG nodes by at most a
constant factor. Argue that we now have |pred (v)| ≤ 2 and |succ(v)| ≤ 2
for all nodes v.
Assuming that |dep(v)| and |dep −1 (v)| are bounded by a constant for all
nodes v, the worst-case time complexity of the simple work-list algorithm can
be expressed as
O(n · h · k)
62 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
where n is the number of CFG nodes in the program being analyzed, h is the
height of the lattice L for abstract states, and k is the worst-case time required to
compute a constraint function fi (x1 , . . . , xn ).
Exercise 5.20: Prove the above statement about the worst-case time complexity
of the simple work-list algorithm. (It is reasonable to assume that the work-list
operations ‘add’ and ‘removeNext’ take constant time.)
Exercise 5.22: Estimate the worst-case time complexity of the sign analysis
with the simple work-list algorithm, using the formula above. (As this for-
mula applies to any dataflow analysis implemented with the simple work-list
algorithm, the actual worst-case complexity of this specific analysis may be
asymptotically better!)
var x,y,z;
x = input;
5.4 LIVE VARIABLES ANALYSIS 63
while (x>1) {
y = x/2;
if (y>3) x = x-y;
z = x-4;
if (z>0) x = x/2;
z = z-1;
}
output x;
the lattice modeling abstract states is thus:4
States = (P({x,y,z}), ⊆)
z>0 x = x/2
z = z−1
output x
For every CFG node v we introduce a constraint variable [[v]] denoting the subset
of program variables that are live at the program point before that node. The
analysis will be conservative, since the computed set may be too large. We use
the auxiliary definition [
JOIN (v) = [[w]]
w∈succ(v)
Unlike the JOIN function from sign analysis, this one combines abstract states
from the successors instead of the predecessors. We have defined the order
relation as v=⊆, so t = ∪.
As in sign analysis, the most interesting constraint rule is the one for assign-
ments:
4 A word of caution: For historical reasons, some textbooks and research papers, e.g. [Hec77,
ALSU06], describe dataflow analyses using the lattices “upside down”. This makes no difference
whatsoever to the analysis results (because of the lattice dualities discussed in Chapter 4), but it can
be confusing.
64 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
Exercise 5.23: Explain why the constraint rule for assignments, as defined
above, is sound.
where vars(E) denotes the set of variables occurring in E. For variable declara-
tions and exit nodes:
var X1 , . . . ,Xn : [[v]] = JOIN (v) \ {X1 , . . . , Xn }
[[exit]] = ∅
For all other nodes:
[[v]] = JOIN (v)
Exercise 5.24: Argue that the right-hand sides of the constraints define mono-
tone functions.
The example program yields these constraints:
[[entry]] = [[var x,y,z]]
[[var x,y,z]] = [[x=input]] \ {x, y, z}
[[x=input]] = [[x>1]] \ {x}
[[x>1]] = ([[y=x/2]] ∪ [[output x]]) ∪ {x}
[[y=x/2]] = ([[y>3]] \ {y}) ∪ {x}
[[y>3]] = [[x=x-y]] ∪ [[z=x-4]] ∪ {y}
[[x=x-y]] = ([[z=x-4]] \ {x}) ∪ {x,y}
[[z=x-4]] = ([[z>0]] \ {z}) ∪ {x}
[[z>0]] = [[x=x/2]] ∪ [[z=z-1]] ∪ {z}
[[x=x/2]] = ([[z=z-1]] \ {x}) ∪ {x}
[[z=z-1]] = ([[x>1]] \ {z}) ∪ {z}
[[output x]] = [[exit]] ∪ {x}
[[exit]] = ∅
whose least solution is:
[[entry]] = ∅
[[var x,y,z]] = ∅
[[x=input]] = ∅
5.4 LIVE VARIABLES ANALYSIS 65
[[x>1]] = {x}
[[y=x/2]] = {x}
[[y>3]] = {x, y}
[[x=x-y]] = {x, y}
[[z=x-4]] = {x}
[[z>0]] = {x, z}
[[x=x/2]] = {x, z}
[[z=z-1]] = {x, z}
[[output x]] = {x}
[[exit]] = ∅
From this information a clever compiler could deduce that y and z are never live
at the same time, and that the value written in the assignment z=z-1 is never
read. Thus, the program may safely be optimized into the following one, which
saves the cost of one assignment and could result in better register allocation:
var x,yz;
x = input;
while (x>1) {
yz = x/2;
if (yz>3) x = x-yz;
yz = x-4;
if (yz>0) x = x/2;
}
output x;
(a) Show how the live variables analysis can be modified to compute
strongly live variables.
(b) Show for each program point in the program from Exercise 5.25 the set
of strongly live variables, as computed by your new analysis.
We can estimate the worst-case time complexity of the live variables analysis,
with for example the naive fixed-point algorithm from Section 4.4. We first
observe that if the program has n CFG nodes and b variables, then the lattice
P(Vars)n has height b·n, which bounds the number of iterations we can perform.
Each lattice element can be represented as a bitvector of length b · n. Using the
observation from Exercise 5.19 we can ensure that |succ(v)| ≤ 2 for any node v.
For each iteration we therefore have to perform O(n) intersection, difference, or
equality operations on sets of size b, which can be done in time O(b · n). Thus,
we reach a time complexity of O(b2 · n2 ).
Exercise 5.28: Can you obtain an asymptotically better bound on the worst-
case time complexity of live variables analysis with the naive fixed-point
algorithm, if exploiting properties of the structures of TIP CFGs and the
analysis constraints?
Exercise 5.29: Recall from Section 5.3 that the work-list algorithm relies on
a function dep(v) for avoiding recomputation of constraint functions that
are guaranteed not to change outputs. What would be a good strategy for
defining dep(v) in general for live variables analysis of any given program?
Exercise 5.30: Estimate the worst-case time complexity of the live variables
analysis with the simple work-list algorithm, by using the formula from
page 61.
var x,y,z,a,b;
z = a+b;
y = a*b;
while (y > a+b) {
a = a+1;
x = a+b;
}
we have four different nontrivial expressions, so our lattice for abstract states is
{a+b, a*b} {a+b, y>a+b} {a+b, a+1} {a*b, y>a+b} {a*b, a+1} {y>a+b, a+1}
{a+b, a*b, y>a+b} {a+b, a*b, a+1} {a+b, y>a+b, a+1} {a*b, y>a+b, a+1}
The top element of our lattice is ∅, which corresponds to the trivial information
that no expressions are known to be available. The CFG for the above program
looks as follows:
68 5 DATAFLOW ANALYSIS WITH MONOTONE FRAMEWORKS
var x,y,z,a,b
z = a+b
y = a*b
false
y > a+b
true
a = a+1
x = a+b
Here, the function ↓X removes all expressions that contain the variable X, and
exps collects all nontrivial expressions:
exps(X) = ∅
exps(I) = ∅
exps(input) = ∅
exps(E1 op E2 ) = {E1 op E2 } ∪ exps(E1 ) ∪ exps(E2 )
5.5 AVAILABLE EXPRESSIONS ANALYSIS 69
[[entry]] = ∅
For all other kinds of nodes, the collected sets of expressions are simply propa-
gated from the predecessors:
Exercise 5.31: Explain informally why the constraints are monotone and the
analysis is sound.
z = aplusb;
y = a*b;
while (y > aplusb) {
a = a+1;
aplusb = a+b;
x = aplusb;
}
[[exit]] = ∅
The rules for the remaining nodes, include branch conditions and output state-
ments, are the same as for available expressions analysis.
On the example program:
var x,a,b;
x = input;
a = x-1;
b = x-2;
while (x>0) {
output a*b-x;
x = x-1;
}
output a*b;
5.7 REACHING DEFINITIONS ANALYSIS 71
the analysis reveals that a*b is very busy inside the loop. The compiler can
perform code hoisting and move the computation to the earliest program point
where it is very busy. This would transform the program into this more efficient
version:
var x,a,b,atimesb;
x = input;
a = x-1;
b = x-2;
atimesb = a*b;
while (x>0) {
output atimesb-x;
x = x-1;
}
output atimesb;
var x,y,z;
x = input;
while (x>1) {
y = x/2;
if (y>3) x = x-y;
z = x-4;
if (z>0) x = x/2;
z = z-1;
}
output x;
For every CFG node v the variable [[v]] denotes the set of assignments that may
define values of variables at the program point after the node. We define
[
JOIN (v) = [[w]]
w∈pred(v)
where this time the ↓X function removes all assignments to the variable X. For
all other nodes we define:
[[v]] = JOIN (v)
This analysis can be used to construct a def-use graph, which is like a CFG except
that edges go from definitions (i.e. assignments) to possible uses. Here is the
def-use graph for the example program:
x = input
x>1
y = x/2
y>3 x = x−y
z = x−4
z>0 x = x/2
z = z−1
output x
The def-use graph is a further abstraction of the program and is the basis of
widely used optimizations such as dead code elimination and code motion.
Exercise 5.33: Show that the def-use graph is always a subgraph of the transi-
tive closure of the CFG.
A backward analysis is one that for each program point computes information
about the future behavior. Examples of this are live variables analysis and very
busy expressions analysis. They can be characterized by the right-hand sides of
constraints only depending on successors of the CFG node. Thus, the analysis
starts at the exit node and moves backward in the CFG. For such analyses, the
JOIN function is defined using succ, and dep can be defined by pred .
The distinction between forward and backward applies to any flow-sensitive
analysis. For analyses that are based on a powerset lattice, we can also distin-
guish between may and must analysis.
A may analysis is one that describes information that may possibly be true
and, thus, computes an over-approximation. Examples of this are live variables
analysis and reaching definitions analysis. They can be characterized by the
lattice order being ⊆ and constraint functions that use the ∪ operator to combine
information.
Conversely, a must analysis is one that describes information that must defi-
nitely be true and, thus, computes an under-approximation. Examples of this
are available expressions analysis and very busy expressions analysis. They can
be characterized by the use of ⊇ as lattice order and constraint functions that
use ∩ to combine information.
Thus, our four examples that are based on powerset lattices show every
possible combination, as illustrated by this diagram:
Forward Backward
May Reaching Definitions Live Variables
Must Available Expressions Very Busy Expressions
These classifications are mostly botanical, but awareness of them may provide
inspiration for constructing new analyses.
Exercise 5.34: Which among the following analyses are distributive, if any?
(a) Available expressions analysis.
(b) Very busy expressions analysis.
Exercise 5.35: Let us design a flow-sensitive type analysis for TIP. (This exercise
can be seen as an alternative to the approach in Exercise 3.22.) In the simple
version of TIP we focus on in this chapter, we only have integer values at run-
time, but for the analysis we can treat the results of the comparison operators
> and == as a separate type: boolean. The results of the arithmetic operators +,
-, *, / can similarly be treated as type integer. As lattice for abstract states we
choose
States = Vars → P({integer, boolean})
such that the analysis can keep track of the possible types for every variable.
main(a,b) {
var x,y;
x = a+b;
if (x) { // warning: using integer as branch condition
output 17;
}
y = a>b;
return y+3; // warning: using boolean in addition
}
5.9 INITIALIZED VARIABLES ANALYSIS 75
Exercise 5.37: What is the JOIN function for initialized variables analysis?
detection tool could now check for every use of a variable that it is contained
in the computed set of initialized variables, and emit a warning otherwise. A
warning would be emitted for this trivial example program:
main() {
var x;
return x;
}
Exercise 5.39: Write a TIP program where such an error detection tool would
emit a spurious warning. That is, in your program there are no reads from
uninitialized variables in any execution but the initialized variables analysis
is too imprecise to show it.
(a) How should we order the two elements? That is, which one is > and
which one is ⊥?
(b) How should the constraint rule for assignments be modified to fit with
this alternative lattice?
F
In the simple work-list algorithm, JOIN (v) = w∈dep −1 (v) [[w]] is computed
in each iteration of the while-loop. However, often [[w]] has not changed since last
time v was processed, so much of that computation may be redundant. (When
we introduce inter-procedural analysis in Chapter 8, we shall see that dep −1 (v)
may become large.) We now present another work-list algorithm based on trans-
fer functions that avoids some of that redundancy. With this algorithm, for a
forward analysis each variable xi denotes the abstract state for the program point
before the corresponding CFG node vi , in contrast to the other fixed-point solvers
we have seen previously where xi denotes the abstract state for the program
point after vi (and conversely for a backward analysis).
procedure PropagationWorkListAlgorithm(t1 , . . . , tn )
(x1 , . . . , xn ) := (⊥, . . . , ⊥)
W := {v1 , . . . , vn }
while W 6= ∅ do
vi := W.removeNext()
y := tvi (xi )
for each vj ∈ dep(vi ) do
z := xj t y
if xj 6= z then
xj := z
W.add(vj )
end if
end for
end while
return (x1 , . . . , xn )
end procedure
Compared to the simple work-list algorithm, this variant typically avoids many
redundant least-upper-bound computations. In each iteration of the while-loop,
the transfer function of the current node vi is applied, and the resulting abstract
state is propagated (hence the name of the algorithm) to all dependencies. Those
that change are added to the work-list. We thereby have tv ([[v]]) v [[w]] for all
nodes v, w where w ∈ succ(v).
Widening
where
N = {−∞, . . . , −2, −1, 0, 1, 2, . . . , ∞}
is the set of integers extended with infinite endpoints and the order on intervals
is defined by inclusion:
[l1 , h1 ] v [l2 , h2 ] ⇐⇒ l2 ≤ l1 ∧ h1 ≤ h2
[− , ]
8
8
[−2,2]
[− ,0] [0, ]
8
8
[−2,1] [−1,2]
[− ,−1] [1, ]
8
8
[−2,0] [−1,1] [0,2]
[− ,−2] [2, ]
8
8
[−2,−1] [−1,0] [0,1] [1,2]
This lattice does not have finite height, since it contains for example the following
infinite chain:
Before we specify the constraint rules, we define a function eval that performs
an abstract evaluation of expressions:
eval (σ, X) = σ(X)
eval (σ, I) = [I, I]
eval (σ, input) = [−∞, ∞]
eval (σ, E1 op E2 ) = op(eval
b (σ, E1 ), eval (σ, E2 ))
The abstract operators are all defined by:
op([l
b 1 , h1 ], [l2 , h2 ]) = [ min x op y, max x op y]
x∈[l1 ,h1 ],y∈[l2 ,h2 ] x∈[l1 ,h1 ],y∈[l2 ,h2 ]
Exercise 6.1: Explain informally why the definition of eval given above is
a conservative approximation compared to evaluating TIP expressions con-
cretely. Give an example of how the definition of eval could be modified to
make the analysis more precise (and still sound).
6.1 INTERVAL ANALYSIS 81
Exercise 6.2: This general definition of op b looks simple in math, but it is non-
trivial to implement it efficiently. Write pseudo-code for an implementation of
the abstract greater-than operator b >. (To be usable in practice, the execution
time of your implementation should be less than linear in the input numbers!)
It is acceptable to sacrifice optimal precision, but see how precise you can
make it.
In Chapter 12 we provide a more formal treatment of the topics of soundness
and precision.
The JOIN function is the usual one for forward analyses:
G
JOIN (v) = [[w]]
w∈pred(v)
In the preceding chapter, we defined the analysis result for each analysis as
the least solution to the analysis constraints for the given program. Kleene’s
fixed-point theorem (Section 4.4, page 47) tells us that this is well-defined: a
unique least solution always exists for such analyses. However, as the interval
analysis defined in this section uses an infinite-height lattice, that fixed-point
theorem does not apply, so the careful reader may wonder, do the interval
analysis constraints actually have a least solution for any given program? The
answer is affirmative.
One way to see that the least fixed point exists is to use a stronger variant of
the fixed-point theorem that relies on transfinite iteration sequences and holds
without the finite-height assumption [CC79a]; see also Exercise 4.32. Another
approach is to use a different fixed-point theorem, known as Tarski’s fixed-point
theorem [Tar55]:1
fixed-point theorem (page 47) has the advantage (for finite-height lattices) that it is constructive
in the sense that it directly leads to an algorithm for computing the least fixed point. The theorem
presented here is a corollary of Tarski’s more general theorem that the set of fixed points of f forms
a complete lattice, but we do not need that more general property here.
82 6 WIDENING
f i (⊥) for i = 0, 1, . . .
integer constants occurring in the given program, but other heuristics could also
be used. For single intervals we define the function ω 0 : Intervals → Intervals by
Exercise 6.5: Show that interval analysis with widening, with this definition
of ω, always terminates and yields a solution that is a safe approximation of
lfp(f ).
Widening is not only useful for infinite-height lattices; it can also be used as
an acceleration technique for analysis that have finite-height lattices but converge
too slowly and where a loss of precision is tolerable.
Widening generally shoots above the target, but a subsequent technique
called narrowing may improve the result. Let fω denote the result of analysis
with widening. Narrowing simply consists of computing f (fω ). By Exercise 6.5
we have lfp(f ) v fω . Then f (fω ) v ω(f (fω )) = (ω ◦ f )(fω ) = fω since ω is
extensive and fω is a fixed point of ω ◦ f . In other words, f (fω ) is at least as
precise as fω . Furthermore, f (f (fω )) v f (fω ) since f is monotone, so f (fω ) ∈
{x ∈ L | f (x) v x} and thereby lfp(f ) v f (fω ) by Tarski’s fixed-point theorem.
In other words, f (fω ) is a safe approximation of lfp(f ). This technique may in
fact be iterated arbitrarily many times, as seen in the following exercise.
[x 7→ ⊥, y 7→ ⊥]
[x 7→ [8, 8], y 7→ [0, 1]]
[x 7→ [8, 8], y 7→ [0, 2]]
[x 7→ [8, 8], y 7→ [0, 3]]
..
.
[x 7→ ⊥, y 7→ ⊥]
[x 7→ [7, ∞], y 7→ [0, 1]]
[x 7→ [7, ∞], y 7→ [0, 7]]
[x 7→ [7, ∞], y 7→ [0, ∞]]
Exercise 6.7: Exactly how many narrowing steps are necessary to reach this
solution?
This result is really the best we could hope for, for this program. For that
reason, further narrowing has no effect. However, in general, the decreasing
sequence
fω w f (fω ) w f 2 (fω ) w . . .
Exercise 6.8: Give an example of a TIP program where the narrowing se-
quence diverges for the interval analysis, when using widening followed by
narrowing.
∇: L × L → L
The widening operator ∇ (usually written with infix notation) must satisfy
∀x, y ∈ L : x v x∇y ∧y v x∇y (meaning that it is an upper bound operator) and
that for any increasing sequence z0 v z1 v z2 v . . . , the sequence y0 , y1 , y2 . . .
defined by y0 = z0 and yi+1 = yi ∇ zi+1 for i = 0, 1, . . . converges after a finite
6.2 WIDENING AND NARROWING 85
number of steps. With such an operator, we can approximate the least fixed
point of f by computing the following sequence:
x0 = ⊥
xi+1 = xi ∇ f (xi )
This sequence eventually converges, that is, for some k we have xk+1 = xk .
Additionally, the result is a safe approximation of the ordinary fixed point:
lfp(f ) v xk .
This leads us to the following variant of the naive fixed-point algorithm with
(traditional) widening:
procedure NaiveFixedPointAlgorithmWithWidening(f )
x := ⊥
while x 6= f (x) do
x := x ∇ f (x)
end while
return x
end procedure
The other fixed-point algorithms (Section 5.3) can be extended with this form
of widening in a similar manner.
Note that if we choose as a special case ∇ = t, the computation of x0 , x1 , . . .
proceeds exactly as with the ordinary naive fixed-point algorithm.
Exercise 6.10: Show that t is a widening operator (albeit perhaps not a very
useful one) if L has finite height.
and (
h1 if h2 ≤ h1
h3 =
min{i ∈ B | h2 ≤ i} otherwise
Compared to the definition of ω 0 for simple widening (see page 83), we now
coarsen the interval end points only if they are unstable compared to the last
iteration. Intuitively, an interval that does not become larger during an iteration
of the fixed-point computation cannot be responsible for divergence.
Now we can define ∇ based on ∇0 , similarly to how we previously defined
ω pointwise in terms of ω 0 :
(σ1 , . . . , σn ) ∇ (σ10 , . . . , σn0 ) = (σ100 , . . . , σn00 ) where σi00 (X) = σi (X) ∇0 σi0 (X)
for i = 1, . . . , n and X ∈ Vars
Exercise 6.11: Show that this definition of ∇ for the interval analysis satisfies
the requirements for being a widening operator.
With this more advanced form of widening but without using narrowing,
for the small example program from page 83 we obtain the same analysis result
as with the combination of simple widening and narrowing we looked at earlier.
Exercise 6.12: Explain why the “simple” form of widening (using the unary
ω operator) is just a special case of the “traditional” widening mechanism
(using the binary ∇ operator).
With the simple form of widening, the analysis effectively just uses a finite
subset of L. In contrast, the traditional form of widening is fundamentally
more powerful: Although each program being analyzed uses only finitely many
elements of L, no finite-height subset suffices for all programs [CC92].
We can be even more clever by observing that divergence can only appear in
presence of recursive dataflow constraints (see Section 5.1) and apply widening
only at, for example, CFG nodes that are loop heads.2 In the above definition of
∇, this means changing the definition of σi00 to
(
00 σi (X) ∇ σi0 (X) if node i is a loop head
σi (X) =
σi0 (X) otherwise
Exercise 6.13: Argue why applying widening only at CFG loop heads suffices
for guaranteeing convergence of the fixed-point computation.
Then give an example of a program where this improves precision for the
interval analysis, compared to widening at all CFG nodes.
2 Aslong as we ignore function calls and only analyze individual functions, the loop heads are
the while nodes in CFGs for TIP programs. If we also consider interprocedural analysis (Chapter 8)
then recursive function calls must also be taken into account.
6.2 WIDENING AND NARROWING 87
Exercise 6.14: We can define another widening operator for interval analysis
that does not require a set B of integer constants. In the definition of ∇0 and
∇ from page 85, we simply chance l3 and h3 as follows:
(
l1 if l1 ≤ l2
l3 =
−∞ otherwise
and (
h1 if h2 ≤ h1
h3 =
∞ otherwise
Intuitively, this widening coarsens unstable intervals to +/ − ∞.
(a) Argue that after this change, ∇ still satisfies the requirements for being
a widening operator.
(b) Give an example of a program that is analyzed less precisely after this
change.
Chapter 7
Until now, we have ignored the values of branch and loop conditions by simply
treating if- and while-statements as a nondeterministic choice between the two
branches, which is called control insensitive analysis. Such analyses are also path
insensitive, because they do not distinguish different paths that lead to a given
program point. The information about branches and paths can be important for
precision. Consider for example the following program:
x = input;
y = 0;
z = 0;
while (x > 0) {
z = z+x;
if (17 > y) { y = y+1; }
x = x-1;
}
The previous interval analysis (with widening) will conclude that after the
while-loop, the variable x is in the interval [−∞, ∞], y is in the interval [0, ∞],
and z is in the interval [−∞, ∞]. However, in view of the conditionals being
used, this result is too pessimistic.
Exercise 7.1: What would be the ideal (i.e., most precise, yet sound) analysis
result for x, y, and z at the exit program point in the example above, when
using the Intervals lattice to describe abstract values? (Later in this chapter
we shall see an improved interval analysis that obtains that result.)
90 7 PATH SENSITIVITY AND RELATIONAL ANALYSIS
x = input;
y = 0;
z = 0;
while (x > 0) {
assert(x > 0);
z = z+x;
if (17 > y) { assert(17 > y); y = y+1; }
x = x-1;
}
assert(!(x > 0));
(We here also extend TIP with a unary negation operator !.) It is always safe to
ignore the assert statements, which amounts to this trivial constraint rule:
With that constraint rule, no extra precision is gained. It requires insight into the
specific static analysis to define nontrivial and sound constraints for assertions.
For the interval analysis, extracting the information carried by general condi-
tions, or predicates, such as E1 > E2 or E1 == E2 relative to the lattice elements
is complicated and in itself an area of considerable study. For simplicity, let us
consider conditions only of the two kinds X > E and E > X. The former kind
of assertion can be handled by the constraint rule
assert(X > E): [[v]] = JOIN (v)[X 7→ gt(JOIN (v)(X), eval (JOIN (v), E))]
Exercise 7.2: Argue that this constraint for assert is sound and monotone.
Negated conditions are handled in similar fashions, and all other conditions
are given the trivial constraint by default.
With this refinement, the interval analysis of the above example will conclude
that after the while-loop the variable x is in the interval [−∞, 0], y is in the
interval [0, 17], and z is in the interval [0, ∞].
Exercise 7.4: Discuss how more conditions may be given nontrivial constraints
for assert to improve analysis precision further.
As the analysis now takes the information in the branch conditions into
account, this kind of analysis is called control sensitive (or branch sensitive). An al-
ternative approach to control sensitivity that does not involve assert statements
is to model each branch node in the CFG using two constraint variables instead
of just one, corresponding to the two different outcomes of the evaluation of the
branch condition. Another approach is to associate dataflow constraints with
CFG edges instead of nodes. The technical details of such approaches will be
different compared to the approach taken here, but the overall idea is the same.
if (condition) {
open();
flag = 1;
} else {
flag = 0;
}
...
if (flag) {
close();
}
We here assume that open and close are built-in functions for opening and
closing a specific file. (A more realistic setting with multiple files can be handled
using techniques presented in Chapter 11.) The file is initially closed, condition
is some complex expression, and the “. . . ” consists of statements that do not
call open or close or modify flag. We wish to design an analysis that can check
that close is only called if the file is currently open, that open is only called if
the file is currently closed, and that the file is definitely closed at the program
exit. In this example program, the critical property is that the branch containing
the call to close is taken only if the branch containing the call to open was taken
earlier in the execution.
As a starting point, we use this powerset lattice for modeling the open/closed
92 7 PATH SENSITIVITY AND RELATIONAL ANALYSIS
[[open()]] = {open}
[[close()]] = {closed}
For the entry node, we define:
[[entry]] = {closed}
and for every other node, which does not modify the file status, the constraint is
simply
[[v]] = JOIN (v)
where JOIN is defined as usual for a forward, may analysis:
[
JOIN (v) = [[w]]
w∈pred(v)
In the example program, the close function is clearly called if and only if
open is called, but the current analysis fails to discover this.
Exercise 7.5: Write the constraints being produced for the example program
and show that the solution for [[flag]] (the node for the last if condition) is
{open, closed}.
Arguing that the program has the desired property obviously involves the
flag variable, which the lattice above ignores. So, we can try with a slightly
more sophisticated lattice – a product lattice of two powerset lattices that keeps
track of both the status of the file and the value of the flag:
} else {
assert(!condition);
flag = 0;
}
...
if (flag) {
assert(flag);
close();
} else {
assert(!flag);
}
This is still insufficient, though. At the program point after the first if-else
statement, the analysis only knows that open may have been called and flag
may be 0.
Exercise 7.6: Specify the constraints that fit with the L0 lattice. Then show that
the analysis produces the lattice element ({open, closed}, {flag = 0, flag 6= 0})
at the program point after the first if-else statement.
L00 = Paths → L
where Paths is a finite set of path contexts. A path context is typically a predicate
over the program state.1 (For instance, a condition expression in TIP defines
such a predicate.) In general, each statement is then analyzed in |Paths| different
path contexts, each describing a set of paths that lead to the statement, which is
why this kind of analysis is called path sensitive. For the example above, we can
use Paths = {flag = 0, flag 6= 0}.
The constraints for open, close, and entry can now be defined as follows.2
[[open()]] = λp.{open}
[[close()]] = λp.{closed}
[[entry]] = λp.{closed}
The constraints for assignments make sure that flag gets special treatment:
1 Anotherway to select Paths is to use sequences of branch nodes.
2 Wehere use the lambda abstration notation to denote a function: if f = λx.e then f (x) = e.
Thus, λp.{open} is the function that returns {open} for any input p.
94 7 PATH SENSITIVITY AND RELATIONAL ANALYSIS
S
flag = 0: [[v]] = [flag = 0 7→ p∈Paths JOIN (v)(p),
flag 6= 0 7→ ∅]
S
flag = I: [[v]] = [flag 6= 0 7→ p∈Paths JOIN (v)(p),
flag = 0 7→ ∅]
S
flag = E: [[v]] = λq. p∈Paths JOIN (v)(p)
Here, I is an integer constant other than 0 and E is a non-integer-constant
expression. The definition of JOIN follows from the lattice structure and from
the analysis being forward:
[
JOIN (v)(p) = [[w]](p)
w∈pred(v)
The constraint for the case flag = 0 models the fact that flag is definitely 0 after
the statement, so the open/closed information is obtained from the predecessors,
independent of whether flag was 0 or not before the statement. Also, the
open/closed information is set to the bottom element ∅ for flag 6= 0 because that
path context is infeasible at the program point after flag = 0. The constraint for
flag = I is dual, and the last constraint covers the cases where flag is assigned
an unknown value.
For assert statements, we also give special treatment to flag:
assert(flag): [[v]] = [flag 6= 0 7→ JOIN (v)(flag 6= 0),
flag = 0 7→ ∅]
Notice the small but important difference compared to the constraint for
flag = 1 statements. As before, the case for negated expressions is similar.
Finally, for any other node v, including other assert statements, the constraint
keeps the dataflow information for different path contexts apart but otherwise
simply propagates the information from the predecessors in the CFG:
Although this is sound, we could make more precise constraints for assert
nodes by recognizing other patterns that fit into the abstraction given by the
lattice.
For our example program, the following constraints are generated:
[[entry]] = λp.{closed}
[[condition]] = [[entry]]
[[assert(condition)]] = [[condition]]
[[open()]] = λp.{open}
S
[[flag = 1]] = flag 6= 0 7→ p∈Paths [[open()]](p), flag = 0 7→ ∅
[[assert(!condition)]] = [[condition]]
7.2 PATHS AND RELATIONS 95
S
[[flag = 0]] = flag = 0 7→ p∈Paths [[assert(!condition)]](p),
flag 6
= 0 →
7 ∅
[[...]] = λp. [[flag = 1]](p) ∪ [[flag = 0]](p)
[[flag]] = [[...]]
[[assert(flag)]] = [flag 6= 0 7→ [[flag]](flag 6= 0), flag = 0 7→ ∅]
[[close()]] = λp.{closed}
[[assert(!flag)]] = [flag = 0 7→ [[flag]](flag =0), flag 6= 0 7→ ∅]
[[exit]] = λp. [[close()]](p) ∪ [[assert(!flag)]](p)
flag = 0 flag 6= 0
[[entry]] {closed} {closed}
[[condition]] {closed} {closed}
[[assert(condition)]] {closed} {closed}
[[open()]] {open} {open}
[[flag = 1]] ∅ {open}
[[assert(!condition)]] {closed} {closed}
[[flag = 0]] {closed} ∅
[[...]] {closed} {open}
[[flag]] {closed} {open}
[[assert(flag)]] ∅ {open}
[[close()]] {closed} {closed}
[[assert(!flag)]] {closed} ∅
[[exit]] {closed} {closed}
The analysis produces the lattice element [flag = 0 7→ {closed}, flag 6= 0 7→
{open}] for the program point after the first if-else statement. The constraint
for the assert(flag) statement will eliminate the possibility that the file is
closed at that point. This ensures that close is only called if the file is open, as
desired.
Exercise 7.8: For the present example, the basic lattice L is a defined as a
powerset of a finite set A = {open, closed}.
(a) Show that Paths → P(A) is isomorphic to P(Paths × A) for any finite
set A. (This explains why such analyses are called relational: each
element of P(Paths × A) is a (binary) relation between Paths and A.)
(b) Reformulate the analysis using P({flag = 0, flag 6= 0} × {open, closed})
instead of L00 (without affecting the analysis precision).
Exercise 7.9: Describe a variant of the example program above where the
present analysis would be improved if combining it with constant propaga-
tion.
In general, the program analysis designer is left with the choice of Paths.
Often, Paths consists of combinations of predicates that appear in conditionals
in the program. This quickly results in an exponential blow-up: for k predicates,
96 7 PATH SENSITIVITY AND RELATIONAL ANALYSIS
Exercise 7.10: Assume that we change the rule for open from
[[open()]] = λp.{open}
to
[[open()]] = λp. if JOIN (v)(p) = ∅ then ∅ else {open}
Argue that this is sound and for some programs more precise than the original
rule.
Exercise 7.12: Construct yet another variant of the open/close example pro-
gram where the desired property can only be established with a choice of
Paths that includes a predicate that does not occur as a conditional expression
in the program source. (Such a program may be challenging to handle with
iterative refinement techniques.)
7.2 PATHS AND RELATIONS 97
Exercise 7.13: The following TIP code computes the absolute value of x:
if (x < 0) {
sgn = -1;
} else {
sgn = 1;
}
y = x * sgn;
Design an analysis (i.e., define a lattice and describe the relevant constraint
rules) that is able to show that y is always positive or zero after the last
assignment in this program.
Chapter 8
Interprocedural Analysis
So far, we have only analyzed the body of individual functions, which is called
intraprocedural analysis. We now consider interprocedural analysis of whole pro-
grams containing multiple functions and function calls.
We use the subset of the TIP language containing functions, but still ignore
pointers and functions as values. As we shall see, the CFG for an entire program
is then quite simple to obtain. It becomes more complicated when adding
function values, which we discuss in Chapter 10.
First we construct the CFGs for all individual function bodies as usual. All
that remains is then to glue them together to reflect function calls properly.
We need to take care of parameter passing, return values, and values of local
variables across calls. For simplicity we assume that all function calls are
performed in connection with assignments:
X = f (E1 , . . . , ,En );
Exercise 8.1: Show how any program can be normalized (cf. Section 2.3) to
have this form.
In the CFG, we represent each function call statement using two nodes: a
call node representing the connection from the caller to the entry of f, and an
after-call node where execution resumes after returning from the exit of f:
100 8 INTERPROCEDURAL ANALYSIS
= f (E 1,...,En )
X=
result = E
As discussed in Section 2.5, CFGs can be constructed such that there is always a
unique entry node and a unique exit node for each function.
We can now glue together the caller and the callee as follows:
f(b1 , ..., b n )
= f (E 1,...,E n )
X=
result = E
The connection between the call node and its after-call node is represented by
a special edge (not in succ and pred ), which we need for propagating abstract
values for local variables of the caller.
With this interprocedural CFG in place, we can apply the monotone frame-
work. Examples are given in the following sections.
8.1 INTERPROCEDURAL CONTROL FLOW GRAPHS 101
Exercise 8.2: How many edges may the interprocedural CFG contain in a
program with n CFG nodes?
Recall the intraprocedural sign analysis from Sections 4.1 and 5.1. That
analysis models values with the lattice Sign:
+ 0 −
and abstract states are represented by the map lattice States = Vars → Sign. For
any program point, the abstract state only provides information about variables
that are in scope; all other variables can be set to ⊥.
To make the sign analysis interprocedural, we define constraints to model the
information flow at function calls and returns. First, we treat the call nodes as
no-ops that simply collect information from their predecessors. Thus, if v is a call
node, we define [[v]] = JOIN (v). For an entry node v of a function f (b1 , . . . ,bn )
we consider the abstract states for all callers pred (v ) and model the passing of
parameters:
G
[[v]] = sw
w∈pred(v)
where1
sw = ⊥[b1 7→ eval ([[w]], Ew w
1 ), . . . , bn 7→ eval ([[w]], En )]
where Ew i is the i’th argument at the call node w. As discussed in Section 4.4,
constraints can be expressed using inequations instead of equations. The con-
straint rule above can be reformulated as follows, where v is a function entry
node and w ∈ pred (v) is a caller:
sw v [[v]]
Intuitively, this shows how information flows from the call node (the left-hand-
side of v) to the function entry node (the right-hand-side of v). This for-
mulation fits well with algorithms like PropagationWorkListAlgorithm from
Section 5.10.
Exercise 8.3: Explain why these two formulations of the constraint rule for
function entry nodes are equivalent.
1 In this expression, ⊥ denotes the bottom element of the Vars → Sign, that is, it maps every
For the entry node v of the main function with parameters b1 , . . . , bn we have
this special rule that models the fact that main is implicitly called with unknown
arguments:
[[v]] = ⊥[b1 7→ >, . . . , bn 7→ >]
Function exit nodes are modeled as no-ops, just like call nodes, so if v is an
exit node we can define [[v]] = JOIN (v). For an after-call node v that stores
the return value in the variable X and where v 0 is the accompanying call node
and w ∈ pred (v) is the function exit node, the dataflow can be modeled by the
following constraint:
[[v]] = [[v 0 ]][X 7→ [[w]](result)]
The constraint obtains the abstract values of the local variables from the call
node v 0 and the abstract value of result from w.
Notice that in this definition of the analysis constraints we exploit the fact
that the variant of the TIP language we use in this chapter does not have global
variables, a heap, nested functions, or higher-order functions.
Exercise 8.4: Write and solve the constraints that are generated by the inter-
procedural sign analysis for the following program:
inc(a) {
return a+1;
}
main() {
var x,y;
x = inc(17);
y = inc(87);
return x+y;
}
Exercise 8.5: Assume we extend TIP with global variables. Such variables are
declared before all functions and their scope covers all functions. Write a TIP
program with global variables that is analyzed incorrectly (that is, unsoundly)
with the current analysis. Then show how the constraint rules above should
be modified to accommodate this language feature.
Function entry nodes may have many predecessors, and similarly, function
exit nodes may have many successors. For this reason, algorithms like Propaga-
tionWorkListAlgorithm (Section 5.10) are often preferred for interprocedural
dataflow analysis.
Exercise 8.6: For the interprocedural sign analysis, how can we choose dep(v)
when v is a call node, an after-call node, a function entry node, or a function
exit node?
8.2 CONTEXT SENSITIVITY 103
main() {
var x,y;
x = f(0); // call 1
y = f(87); // call 2
return x + y;
}
Due to the first call to f the parameter z may be 0, and due to the second call it
may be a positive number, so in the abstract state at the entry of f, the abstract
value of z is >. That value propagates through the body of f and back to the
callers, so both x and y also become >. This is an example of dataflow along
interprocedurally invalid paths: according to the analysis constraints, dataflow
from one call node propagates through the function body and returns not only
at the matching after-call node but at all after-call nodes. Although the analysis
is still sound, the resulting loss of precision may be unacceptable.
A naive solution to this problem is to use function cloning. In this specific
example we could clone f and let the two calls invoke different but identical func-
tions. A similar effect would be obtained by inlining the function body at each
call. More generally this may, however, increase the program size significantly,
and in case of (mutually) recursive functions it would result in infinitely large
programs. As we shall see next, we can instead encode the relevant information
to distinguish the different calls by the use of more expressive lattices, much
like the path-sensitivity approach in Chapter 7.
As discussed in the previous section, a basic context-insensitive dataflow ana-
lysis can be expressed using a lattice States n where States is the lattice describing
abstract states and n = |Nodes| (or equivalently, using a lattice Nodes → States).
Context-sensitive analysis instead uses a lattice of the form
n
Contexts → lift(States)
where G
JOIN (v, c) = [[w]](c)
w∈pred(v)
to match the new lattice with context sensitivity. Note that information for
different call contexts is kept apart, and that the reachability information is
propagated along. How to model the dataflow at call nodes, after-call nodes,
function entry nodes, and function exit nodes depends on the context sensitivity
strategy, as described in the following sections.
which has different abstract values for z depending on the caller. Notice that the
information for the context is unreachable, since f is not the main function but
is always executed from c1 or c2 .
Parameter passing at function calls is modeled in the same way as in context-
insensitive analysis, but now taking the call contexts into account. Assume w is
a call node and v is the entry node of the function f (b1 , . . . ,bn ) being called.
0
The abstract state scw defined by
(
c0 unreachable if [[w]](c0 ) = unreachable
sw =
⊥[b1 7→ eval ([[w]](c ), E1 ), . . . , bn 7→ eval ([[w]](c ), Ew
0 w 0
n )] otherwise
describes the dataflow from w to v, like the context-insensitive variant but now
parameterized with a context c0 for the call node. This definition also models
the fact that no new dataflow can appear from call node w in context c0 if that
combination of node and context is unreachable (perhaps because the analysis
has not yet encountered any dataflow to that node and context).
The constraint rule for a function entry node v where w ∈ pred (v) is a caller
and c0 ∈ Contexts is a call context can then be written as follows.
0
c
sw v [[v]](c) where c = w
0
Informally, for any call context c0 at the call node w, an abstract state scw is built
by evaluating the function arguments and propagated to call context c at the
106 8 INTERPROCEDURAL ANALYSIS
function entry node v. In this simple case where k = 1, the call node w is directly
used as context c for the function entry node, but for larger values of k it is
necessary to express how the call site is pushed onto the stack (represented by
the call string from the call context).
Alternatively, the constraint rule can be expressed as an equation by collecting
the abstract states for all relevant call nodes and contexts:
G 0
[[v]](c) = scw
w ∈ pred(v) ∧
c=w∧
c0 ∈ Contexts
Exercise 8.8: Verify that this constraint rule for function entry nodes indeed
leads to the lattice element shown above for the example program.
Exercise 8.9: Give a constraint rule for the entry node of the special function
main. (Remember that main is always reachable in context and that the
values of its parameters can be any integers.)
Assume v is an after-call node that stores the return value in the variable X,
and that v 0 is the associated call node and w ∈ pred (v) is the function exit node.
The constraint rule for v merges the abstract state from the v 0 and the return
value from w, now taking the call contexts and reachability into account:
(
unreachable if [[v 0 ]](c) = unreachable ∨ [[w]](v 0 ) = unreachable
[[v]](c) =
[[v 0 ]](c)[X 7→ [[w]](v 0 )(result)] otherwise
Notice that with this kind of context sensitivity, v 0 is both a call node and a call
context, and the abstract value of result is obtained from the exit node w in call
context v 0 .
Exercise 8.10: Write and solve the constraints that are generated by the inter-
procedural sign analysis for the program from Exercise 8.4, this time with
context sensitivity using the call string approach with k = 1. (Even though
this program does not need context sensitivity to be analyzed precisely, it
illustrates the mechanism behind the call string approach.)
8.4 CONTEXT SENSITIVITY WITH THE FUNCTIONAL APPROACH 107
Explain informally what this tells us about the program P . Give an example
of what program P may be.
Exercise 8.12: Write a TIP program that needs the call string bound k = 2 or
higher to be analyzed with optimal precision using the sign analysis. That is,
some variable in the program is assigned the abstract value > by the analysis
if and only if k < 2.
Exercise 8.13: Generalize the constraint rules shown above to work with any
k ≥ 1, not just k = 1.
In summary, the call string approach distinguishes calls to the same function
based on the call sites that appear in the call stack. In practice, k = 1 sometimes
gives inadequate precision, and k ≥ 2 is generally too expensive. For this reason,
it is common to select k individually for each call site, based on heuristics.
main() {
var x,y;
x = f(42); // call 1
y = f(87); // call 2
return x + y;
}
The call string approach with k ≥ 1 will analyze the f function twice, which
is unnecessary because the abstract value of the argument is + at both calls.
Rather than distingushing calls based on information about control flow from
the call stack, the functional approach to context sensitivity distinguishes calls
108 8 INTERPROCEDURAL ANALYSIS
based on the data from the abstract states at the calls. In the most general form,
the functional approach uses
Contexts = States
although a subset often suffices. With this set of call contexts, the analysis lattice
becomes
n
States → lift(States)
which clearly leads to a significant increase of the theoretical worst-case com-
plexity compared to context insensitive analysis.
The idea is that a lattice element for a CFG node v is a map mv : States →
lift(States) such that mv (s) approximates the possible states at v given that the
current function containing v was entered in a state that matches s. The situation
mv (s) = unreachable means that there is no execution of the program where the
function is entered in a state that matches s and v is reached. If v is the exit node
of a function f , the map mv is a summary of f , mapping abstract entry states to
abstract exit states, much like a transfer function (see Section 5.10) models the
effect of executing a single instruction but now for an entire function.
Returning to the example program from Section 8.2 (page 103), we will now
define the analysis constraints such that, in particular, we obtain the following
lattice element at the exit of the function f:3
⊥[z 7→ 0] 7→ ⊥[z 7→ 0, result 7→ 0],
⊥[z 7→ +] 7→ ⊥[z 7→ +, result 7→ +],
all other contexts 7→ unreachable]
Exercise 8.15: Verify that this constraint rule for function entry nodes indeed
leads to the lattice element shown above for the example program.
Exercise 8.16: Give a constraint rule for the entry node of the special function
main. (Remember that main is always reachable and that the values of its
parameters can be any integers.)
Assume v is an after-call node that stores the return value in the variable X,
and that v 0 is the associated call node and w ∈ pred (v) is the function exit node.
The constraint rule for v merges the abstract state from the v 0 and the return
value from w, while taking the call contexts and reachability into account:
(
unreachable if [[v 0 ]](c) = unreachable ∨ [[w]](scv0 ) = unreachable
[[v]](c) =
[[v 0 ]](c)[X 7→ [[w]](scv0 )(result)] otherwise
To find the relevant context for the function exit node, this rule builds the
same abstract state as the one built at the call node.
Exercise 8.17: Assume we have analyzed a program P using context sensitive
interprocedural sign analysis with the functional approach, and the analysis
result contains the following lattice element for the exit node of a function
named foo:
[x 7→ -, y 7→ -, result 7→ ⊥] 7→ [x 7→ +, y 7→ +, result 7→ +],
[x 7→ +, y 7→ +, result 7→ ⊥] 7→[x 7→ -, y 7→ -, result 7→ -],
all other contexts 7→ unreachable
Explain informally what this tells us about the program P . What could foo
look like?
Exercise 8.18: Write and solve the constraints that are generated by the inter-
procedural sign analysis for the program from Exercise 8.4, this time with
context sensitivity using the functional approach.
Exercise 8.19: Show that this claim about the precision of the functional
approach is correct.
parameters but not other parts of the program state. In the version of TIP used
in this chapter, there are no pointers or global variables, so the entire program
state at function entries is defined by the values of the parameters, which means
that the analysis presented in this section coincides with parameter sensitivity.
When analyzing object oriented programs, a popular choice is object sensitivity,
which is essentially a variant of the functional approach that distinguishes calls
not on the entire abstract states at function entries but only on the abstract values
of the receiver objects.
Chapter 9
Distributive Analysis
Frameworks
Recall from Exercise 5.26 that an analysis is distributive if all its constraint func-
tions are distributive as defined in Exercise 4.20: A function f : L1 → L2 where
L1 and L2 are lattices is distributive when ∀x, y ∈ L1 : f (x) t f (y) = f (x t y).
In this chapter, we demonstrate how distributivity enables efficient algorithms
for context- and flow-sensitive analysis. These results build on a combination of
(1) the functional approach to context sensitivity (Section 8.4), and (2) a clever
compact representations of distributive functions. The techniques presented in
this chapter originate from Reps, Horwitz and Sagiv [RHS95, SRH96] but are
presented in a constraint-based style to make the connections to the preceding
chapters more clear.
b = a + c; // b,c
if (input) { // b,c
c = 1; // b
} // b,c
At the program point immediately after the variable declarations, all three
variables are possibly-uninitialized. Notice that b remains possibly-uninitialized
after the assignment b = a + c because its value depends on c (and therefore
this analysis is not simply the dual of the initialized variables analysis from
Section 5.9). At the final program point, c is possibly-uninitialized despite the
assignment c = 1 because the branch containing that assignment may not be
taken.
We shall now design such an analysis that is both flow sensitive and context
sensitive using the functional approach to context sensitivity, following the
methodology from Chapters 5 and 8. As discussed in Section 8.4, the functional
approach to context sensitivity is precise enough to completely avoid the problem
with interprocedurally invalid paths. The results of such an analysis can be used
as follows. When analyzing a program, we want to know for a given variable
X and program point v inside a function f whether f may be called in some
context c such that X is possibly-uninitialized at v in context c – and if so, a
warning can be emitted informing the programmer that the program behavior
is likely unintended if the instruction at v uses X.
First, we choose a suitable lattice of abstract states,
States = P(Vars)
foo(x, y) {
var z;
1 unreachable here denotes the bottom element of the lifted lattice as in Section 8.2.
9.1 POSSIBLY-UNINITIALIZED VARIABLES ANALYSIS 113
z = x - y;
z = z * z;
return z;
}
The entry of the main function is always reachable with the empty context:2
The transfer function for a variable declaration, var X, models the fact that
variables in TIP are uninitialized right after they have been declared:
As in Section 8.2, the analysis constraint for [[v]] when v is a variable decla-
ration or an assignment can be expressed using the transfer function and the
JOIN function:
(
tv (JOIN (v, c)) if JOIN (v, c) ∈ P(Vars)
[[v]](c) =
unreachable if JOIN (v, c) = unreachable
where
G
JOIN (v, c) = [[w]](c)
w∈pred(v)
Exercise 9.1: Explain intuitively why the above definition of the transfer func-
tion for assignments is correct for possibly-uninitialized variables analysis.
Exercise 9.2: Specify suitable constraint rules for other kinds of CFG nodes
than variable declarations and assignments (most importantly, non-main
function entry nodes and after-call nodes).
Finally, we can F
collect the set of possibly-uninitialized variables at each
program point v as c∈Contexts [[v]](c).
2 In TIP, the parameters of the main function obtain values from the program input, and other
Exercise 9.3: At each program point in the following program, which variables
are possibly-uninitialized according to the context sensitive analysis described
above? What is the result if instead using context insensitive analysis?
main() {
var x,y,z;
x = input;
z = p(x,y);
return z;
}
p(a,b) {
if (a > 0) {
b = input;
a = a - b;
b = p(a,b);
output(a);
output(b);
}
return b;
}
In Section 8.4 we saw that the jump function mv for a function exit node v has
the interesting property that it “summarizes” the entire function for all possible
contexts the function may be called in, by mapping entry abstract states to exit
abstract states. If, for example v is the exit node of the foo function on page 112,
then mv is this function:
(
s ∪ {z} if x ∈ s ∨ y ∈ s
mv (s) =
s\{z} otherwise
This means that once mv (s) has been computed for some abstract state s, then
the result can be reused for all calls to foo where the corresponding abstract state
for the function entry is s, without revisiting the function body. The work-list
algorithm from Section 5.3 automatically performs such reuse.
Another important observation from Section 8.4 is that the functional ap-
proach to context sensitivity can be computationally expensive when the number
of contexts is large. The following exercise shows that a basic implementation
of the analysis is not scalable, despite the reuse effect for calls with the same
call context and the use of unreachable for avoiding analysis of functions in
unreachable contexts.
9.2 AN ALTERNATIVE FORMULATION 115
Interestingly, both phases are context insensitive (not context sensitive!), and
yet the resulting analysis achieves the same precision as the original formulation
above for computing which variables may be possibly-uninitialized at a given
program point.
The pre-analysis phase We first describe the pre-analysis. It uses this lattice
for a program with n CFG nodes:
n
lift(P(Vars) → P(Vars))
Even though we think of this as a context insensitive analysis, the analysis lattice
looks very similar to the one from page 112. In the old analysis, an element of
the lattice for individual CFG nodes was a set of possibly-uninitialized program
variables for each context. In this context insensitive pre-analysis, an element
of the lattice for individual CFG nodes is instead either unreachable or a jump
function mv : P(Vars) → P(Vars) with almost the same meaning as before,
relating the set of possibly-uninitialized variables at the entry of the function
containing v with the set of possibly-uninitialized variables at v. The lattice lift
operation is used slightly differently; this new analysis does not track which
116 9 DISTRIBUTIVE ANALYSIS FRAMEWORKS
contexts each function may be called in, but only which functions may be called
in some context, thereby providing a coarser notion of reachability.
With this choice of analysis lattice, the constraint rules for the different kinds
of CFG nodes can be defined as follows.
The entry node entry main of the main function is always reachable, and for
any given set s of possibly-uninitialized variables at the entry of the function, s
is trivially the set of possibly-uninitialized variables at this node since it is the
same node:
The subscript at [[ · ]]1 indicates that the constraint variables here are for analysis
phase 1. For the entry node entry f of any other function f , the constraint is the
same, except that we account for reachability of the function:
Notice that the only information that is carried from the caller to the callee in this
pre-analysis is whether the function is reachable – there is no flow of information
between the functions about which variables are possibly-uninitialized, unlike
the first formulation of the possibly-uninitialized variables analysis.
G
[[X = E]]1 = tX=E ◦ [[w]]1 if [[w]]1 6= unreachable
w ∈ pred(X=E) ∧ for some w ∈ pred (X=E)
[[w]]1 6= unreachable
For each reachable predecessor w in the CFG, [[w]]1 is a jump function for w. By
computing the least upper bound of those functions and composing with the
transfer function, we take one step further, thereby obtaining the jump function
for v as illustrated with the thick dashed edge below.
9.2 AN ALTERNATIVE FORMULATION 117
entry
[[w]]1
w ...
[[v]]1 v tv
[[X = ]]1 (s) = [[v 0 ]]1 (s)\{X} ∪ (tX=result ◦[[w]]1 ◦tX1 ,...,Xn =E1 ,...,En ◦[[v 0 ]]1 )(s)∩{X}
3 As explained in Section 8.2, an interprocedurally valid path is a path in the CFG where all calls
f entry f
entry
s
[[v 0 ]]1
v0
tX1 ,...,Xn =E1 ,...,En [[w]]1
[[v]]1
v
tX=result
Exercise 9.6: What are suitable constraint rules for the remaining kinds of
CFG nodes (for example, if E)?
The main analysis phase After the pre-analysis has been executed, jump func-
tions are available for all CFG nodes. This makes it easy to compute the sets
of possibly-uninitialized variables in the main analysis phase. This analysis
uses the same lattice for abstract states as the first formulation of the analysis in
Section 9.1:
States = P(Vars)
The analysis is flow sensitive but context insensitive, and we want to compute
information only for functions that may be reachable, so we use this analysis
lattice for a program with n CFG nodes:
n
lift(P(Vars))
By using a lifted lattice, we can distinguish between nodes that belong to un-
reachable functions and nodes that are in possibly reachable functions but with
the empty set of possibly-uninitialized variables.
The entry of the main function is always reachable:
[[entry main ]]2 6= unreachable
(The subscript at [[ · ]]2 indicates that these constraint variables belong to analysis
phase 2.)
For the entry node entry f of a function f with parameters X1 , . . . , Xn , the
set of possibly-uninitialized variables is obtained from all the calls to f while
taking reachability into account:
G
[[entry f ]]2 = sw
w∈pred(entry f )
9.2 AN ALTERNATIVE FORMULATION 119
where (
unreachable if [[w]]2 = unreachable
sw =
tw w w
X1 ,...,Xn =E1 ,...,En ([[w]] 2 ) otherwise
w w
and tw X1 ,...,Xn =Ew w models parameter passing for call node w : f (E1 , . . . ,En )
1 ,...,En
as before.
For all non-entry nodes we can now simply use the jump functions computed
by the pre-analysis. If v is a non-entry node and v 0 is the entry of the function
containing v, we have the following constraint that applies v’s jump function
[[v 0 ]]1 to the set of possibly-uninitialized variables at v 0 unless v 0 is unreachable:
(
unreachable if [[v 0 ]]2 = unreachable
[[v]]2 =
[[v 0 ]]1 ([[v 0 ]]2 ) otherwise
v0
[[w]]1
Exercise 9.8: As a continuation of Exercise 9.4, prove that the alternative two-
phase formulation of possibly-uninitialized variables analysis has the same
worst-case complexity as the original formulation. (Hint: The bottleneck is
the pre-analysis, not the main analysis.)
Despite the apparent lack of progress suggested by Exercises 9.7 and 9.8, this
alternative formulation of the analysis paves the way for a trick presented in the
next section.
120 9 DISTRIBUTIVE ANALYSIS FRAMEWORKS
1. functions of the form f : P(D) → P(D) that are distributive and where D
is a finite set,
2. functions of the form f : (D → L) → (D → L) that are distributive and
where D is a finite set and L is a complete lattice.
The first one includes the functions used in the alternative formulation of
possibly-uninitialized variables analysis in Section 9.2 and forms the foundation
of the IFDS framework presented in Section 9.4, while the second one is used in
the IDE framework in Section 9.6.
The first family of functions is a special case of the second one in the following
sense. The lattices P(D) and D → L are isomorphic (see page 43) if L is set to
the two-element lattice {>, ⊥}. An element s ∈ P(D) is then represented by its
characteristic function in D → L:
(
> if d ∈ s
ms (d) =
⊥ if d 6∈ s
• d1 d2 d3
• d1 d2 d3
represents the function g where g(•) = {d1 }, g(d1 ) = ∅, and g(d2 ) = g(d3 ) =
{d3 }, and thereby the function f where f (S) = {d1 , d3 } if d2 ∈ S or d3 ∈ S and
f (S) = {d1 } otherwise. In general, the edges are defined by
{• •} ∪ {• y | y ∈ f (∅)} ∪ {x y | y ∈ f ({x}) ∧ y ∈
/ f (∅)}
where x y denotes an edge from x in the top row to y in the bottom row.
(The condition y ∈ / f (∅) is safe because f (∅) is covered by the • y edges.)
Intuitively, the edges describe how the output depends on the input. In possibly-
uninitialized variables analysis where D = Vars, the edge d2 d3 means: if d2
is possibly-uninitialized at input, then d3 is possibly-uninitialized at output. The
symbol • can be thought of as a dataflow fact that always holds. The edge • d1
means that d1 is unconditionally possibly-uninitialized at output.
Exercise 9.11: Assume Vars = {a, b, c, d}. Draw the bipartite graph represen-
tations of the transfer functions in possibly-uninitialized variables analysis
for each of the following two statements:
var b,d;
a = b + d;
Then give a general description of how to produce the bipartite graph repre-
sentations of the transfer functions for variable declarations and assignments
in possibly-uninitialized variables analysis.
Exercise 9.12: Assume fA : P(D) → P(D) and fB : P(D) → P(D) are dis-
tributive and D is a finite set. Prove that fA ◦ fB and fA t fB are distributive.4
This can be illustrated graphically with examples, here is one for composition:
• d1 d2 d3 • d1 d2 d3
fA
• d1 d2 d3 fA ◦ fB
fB
• d1 d2 d3 • d1 d2 d3
The edges d2 d1 and d3 d1 could be omitted in the resulting graph; they are
unnecessary because of the edge • d1 , but it is simpler to include them.
Exercise 9.13: Draw a similar example illustrating least upper bound of two
distributive functions represented as bipartite graphs.
To conclude this part, all the transfer functions and jump functions that
appear in possibly-uninitialized variables analysis and other distributive flow-
sensitive analyses with a similar analysis lattice can be represented compactly
and constructed efficiently using bipartite graphs.
The bipartite graph representation generalizes to distributive functions of the
form f : (D → L) → (D → L) where D is a finite set and L is a complete lattice.
Assume f is such a function and define g : (D ∪ {•}) × (D ∪ {•}) → (L → L) as
follows.
function g:
• d1 d2 d3
• d1 d2 d3
hv 0 , •i hv 0 , xi hv 0 , yi hv 0 , zi
The exploded CFG has the property that a dataflow fact d ∈ D may hold
at CFG node v if and only if hv, di is reachable from hentry main , •i along an
interprocedurally valid path in the exploded CFG (see Exercise 9.23).
Exercise 9.17: Draw the exploded CFG for the following program using the
transfer functions for possibly-uninitialized variables analysis.
main() {
var a, b, c;
b = input;
a = plus(a, b);
return a - c;
}
plus(d, e) {
return d + e;
}
The IFDS algorithm proceeds in two phases similar to the alternative formu-
lation of possibly-uninitialized variables analysis in Section 9.2. The first phase
(corresponding to the pre-analysis phase in Section 9.2) is called the tabulation
5 Exploded CFGs are called exploded supergraphs in [RHS95, SRH96].
9.4 THE IFDS FRAMEWORK 125
Phase 1 The path edges are constructed as the least solution to the following
constraints for the different kinds of CFG nodes. Thus, the lattice used in this
analysis phase is the powerset lattice of path edges.
The first constraint rule expresses that the program entry is always reachable:
hentry main , •i hentry main , •i ∈ P
Second, if v is a function entry node, v1 ∈ pred (v) is a call node that calls
the function containing v, and v0 is the entry node of the function containing
v1 , this constraint rule models the flow of reachable dataflow facts from the call
node to the function entry node:
v0 f
w
d1
d3
v0
d2
v
d5
w0
d4
Notice how this rule uses the path edges of the exit node w0 of the function being
called as a summary of the function, and that the dataflow follows interproce-
durally valid paths only, which is essential for obtaining the desired context
sensitivity.
The fourth constraint rule models the flow of function-local state that is
transferred from a call node v 0 to its after-call node v, where v0 is the entry node
of the function that contains v and v 0 :
Exercise 9.18: Explain how these constraint rules correspond to those in the
pre-analysis phase in Section 9.2.
Phase 2 The second phase (corresponding to the main analysis phase in Sec-
tion 9.2) that computes the final analysis results is remarkably simple and entirely
intra-procedural: The abstract state [[v]] ∈ P(D) for any CFG node v is obtained
directly from v’s path edges, as expressed by this single rule where v0 is the
entry node in the function containing v:
∀d1 , d2 : hv0 , d1 i hv, d2 i ∈ P ∧ d2 ∈ D =⇒ d2 ∈ [[v]]
This works because a path edge hv0 , d1 i hv, d2 i ∈ P has been added in the
first analysis phase only if the dataflow fact d2 may hold at v. Because of the
9.4 THE IFDS FRAMEWORK 127
Exercise 9.20: Show that the worst-case time complexity of IFDS analysis is
O(|E| · |D|3 ) where E is the set of CFG edges in the program being analyzed.
(Compare this result with Exercises 9.4 and 9.8.)
Exercise 9.22: Which of the five analyses in Exercise 5.34 fit into the IFDS
framework? For the analyses that do, describe how their transfer functions
can be expressed as bipartite graphs.
6 The IFDS and IDE papers [RHS95, SRH96] use upside-down lattices compared to the presenta-
Exercise 9.23: As mentioned earlier in this section, the IFDS framework has
an interesting connection to a certain kind of graph reachability: Dataflow
fact d ∈ D may hold at CFG node v (i.e., d ∈ [[v]]2 ) if and only if there exists
an interprocedurally valid path from hentry main , •i to hv, di in the exploded
CFG. Explain why this property holds.
In constant propagation analysis from Section 5.2, the lattice of abstract states
is Vars → flat(Z). With D = Vars and L = flat(Z), the transfer functions are
of the form f : (D → L) → (D → L) that we studied in Section 9.3. Constant
propagation analysis is, however, not distributive, as seen in Exercise 5.34. Copy-
constant propagation analysis is a less precise variant that is distributive, and it is
suitable for motivating the IDE framework presented in Section 9.6.
(
⊥ if a = ⊥ or b = ⊥
a op
b b=
> otherwise
For any composite expression, such as, x*y+3, the analysis simply yields >. This
means that the transfer functions for assignments, X = E, can be defined by
considering only three cases of right-hand-side expressions: (1) E is an integer
literal, Int, (2) E is a single variable, Y (an assignment of the form X = Y is called
a copy assignment), and (3) E is any other expression.
Exercise 9.24: Describe the labeled bipartite graph representation of the trans-
fer functions for the three kinds of expressions.
9.6 THE IDE FRAMEWORK 129
Exercise 9.25: Describe the result of performing context insensitive (and path
insensitive) copy-constant propagation analysis on the following program.
main() {
var x;
x = p(7);
return x;
}
p(a) {
var y,z;
if (a > 0) {
y = a - 1;
y = p(y);
}
z = -2 * y + 5;
return z;
}
Exercise 9.28: Draw the exploded CFG for the program from Exercise 9.25
using the transfer functions from copy-constant propagation analysis.
In IDE, path edges (called jump functions in [SRH96]) are similarly labeled
with functions of the form L → L. Recall from Section 9.4 that IFDS builds a set
of path edges P, starting with the empty set of edges and in each step adding
an edge. IDE instead builds the edge labels, starting with ⊥ everywhere and in
each step moving upwards in the L → L lattice at some edge. For this reason,
we introduce a different notation. For each pair of exploded CFG nodes, hv1 , d1 i
and hv2 , d2 i, the constraint variable [[hv1 , d1 i hv2 , d2 i]]P : L → L denotes the
label of the corresponding edge.
Phase 1 The first phase of IDE builds the path edges, much like phase one in
IFDS but using the more general lattice consisting of L → L functions. Each of
the following constraint rules directly corresponds to one of the constraint rules
in IFDS.
The first constraint rule expresses that the program entry is always reachable
and has the identity jump function:
The composition of the edge functions express how the abstract value of hv0 , d1 i
influences the abstract value of hv, d5 i. This rule can be illustrated like the
corresponding rule in IFDS but now with edge labels (where m = m4 ◦ m3 ◦
m2 ◦ m1 ):
v0 f
w
d1
d3
m1
m
v0 m2
m3
d2
v
d5
w0
m4 d4
The fourth constraint rule applies to after-call nodes v where v 0 is the asso-
ciated call node and v0 is the entry node of the function that contains v and
v0 :
m2
∀d1 , d2 , d3 : m1 = [[hv0 , d1 i hv 0 , d2 i]]P ∧ hv 0 , d2 i hv, d3 i ∈ E
=⇒ m2 ◦ m1 v [[hv0 , d1 i hv, d3 i]]P
Similar to IFDS, this constraint rule also applies to any other kind of node v
with a predecessor v 0 ∈ pred (v) where v0 is the entry node of the function that
contains v and v 0 .
Phase 2 The second phase of IDE analysis uses the path edges (representing
jump functions) that have been produced in the first phase to compute an
abstract value [[hv, di]] ∈ lift(L) for each exploded CFG node hv, di. We here use
a lifted lattice with the bottom element unreachable representing nodes that are
unreachable from the program entry. This analysis phase closely resembles the
main phase of the analysis described in Section 9.2.
As usual, we specify the analysis in a constraint-based style. First, the pro-
gram entry is always reachable:
Second, at any CFG node v, the abstract value of d ∈ D can be obtained from the
abstract values at the entry v0 of the function containing v and the path edges
132 9 DISTRIBUTIVE ANALYSIS FRAMEWORKS
The least solution to the resulting constraints for each analysis phase for a
given program can be computed using a work-list algorithm, much like for IFDS
but with a little more effort for the second phase (pseudo-code can be found in
[SRH96]).
Exercise 9.29: Explain why the second phase of IDE does not need separate
constraints for modeling the flow of return values from function exit nodes to
after-call nodes.
The efficiency of IDE depends on the edge label functions being efficiently
representable. This means that functions that appear on the edges E in the
exploded CFG and on the path edges P can be represented in constant space (i.e.,
independently of the size of the program), and function composition, least upper
bound, equality, and function application can be computed in constant time.
For copy-constant propagation analysis, we only need the identity function and
constant functions, which trivially have this property (assuming that operations
on constants can be performed in constant time).
2. The worst-case time complexity of IDE analysis is the same as for IFDS:
O(|E|·|D|3 ) where E is the set of CFG edges in the program being analyzed,
provided that the edge functions are efficiently representable and the
height of L is constant.
Exercise 9.33: Prove the statement above about the worst-case time complexity
of IDE analysis.
Exp → λId.Exp
| Id
| Exp Exp
(In Section 10.3 we demonstrate this analysis technique on the TIP language.)
The concrete syntax also allows parentheses.1 For simplicity we assume that all
λ-bound variables are distinct. To construct a CFG for a term in this calculus,
we need to approximate for every expression E the set of closures to which it may
evaluate. A closure can be modeled by a symbol of the form λX that identifies
1 As customary, application is left-associative, and application has higher precedence than ab-
straction.
136 10 CONTROL FLOW ANALYSIS
λX ∈ [[λX.E]]
for every abstraction λX.E, which models that the actual argument E2 may flow
into the formal argument X and that the value of the function body E is among
the possible results of the function call E1 E2 .
As an example, for the small program (λs.λz.sz)(λn.λt.λe.e)(λr.λp.r), which
contains 7 λ-abstractions and 3 applications, the following constraints are pro-
duced:
λs ∈ [[λs.λz.sz]]
λz ∈ [[λz.sz]]
λn ∈ [[λn.λt.λe.e]]
λt ∈ [[λt.λe.e]]
λe ∈ [[λe.e]]
λr ∈ [[λr.λp.r]]
λp ∈ [[λp.r]]
λs ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[s]] ∧ [[λz.sz]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λz ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[z]] ∧ [[sz]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λn ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[n]] ∧ [[λt.λe.e]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λt ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[t]] ∧ [[λe.e]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λe ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[e]] ∧ [[e]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λr ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[r]] ∧ [[λp.r]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λp ∈ [[λs.λz.sz]] =⇒
([[λn.λt.λe.e]] ⊆ [[p]] ∧ [[r]] ⊆ [[(λs.λz.sz)(λn.λt.λe.e)]])
λs ∈ [[s]] =⇒
([[z]] ⊆ [[s]] ∧ [[λz.sz]] ⊆ [[sz]])
λz ∈ [[s]] =⇒
([[z]] ⊆ [[z]] ∧ [[sz]] ⊆ [[sz]])
λn ∈ [[s]] =⇒
10.1 CLOSURE ANALYSIS FOR THE λ-CALCULUS 137
[[(λs.λz.sz)(λn.λt.λe.e)]] = {λz}
[[(λs.λz.sz)(λn.λt.λe.e)(λr.λp.r)]] = {λt}
In particular, we see from this solution that s in the application sz can only be the
abstraction λn.λt.λe.e (according to the solution for [[s]]), and at the outermost
application (corresponding to the entire program), the abstraction being applied
can only be λz.sz (according to the solution for [[(λs.λz.sz)(λn.λt.λe.e)]]).
Exercise 10.1: Show that the resulting constraints for any λ-calculus term are
monotone and can be solved by a fixed-point computation, with an appropri-
ate choice of lattice.
Exercise 10.2: Show that a unique least solution exists, since solutions are
closed under intersection.
The following algorithm can find the least solution to the collection of con-
straints in time O(n3 ) for a program of size n (for example, measured by the
number of AST nodes). It maintains a directed graph where nodes correspond to
constraint variables and edges reflect inclusion constraints. For each constraint
variable x,
• x.sol ⊆ T holds the solution for x,
• x.succ ⊆ V is the set of successors of x (i.e., the edges of the graph), and
• x.cond (t) ⊆ V × V represents a set of conditional constraints for x and t.
Additionally we have
• W ⊆ T × V is a worklist.
All the sets are initially empty. The constraints are processed one by one using
the following helper functions for adding tokens and edges and propagating
tokens to satisfy the constraints that have been seen so far.
procedure addToken(t, x)
if t 6∈ x.sol then
x.sol .add(t)
W.add(t, x)
10.2 THE CUBIC ALGORITHM 139
end if
end procedure
procedure addEdge(x, y)
if x 6= y ∧ y 6∈ x.succ then
x.succ.add(y)
for t in x.sol do
addToken(t, y)
end for
end if
end procedure
procedure propagate()
while W 6= ∅ do
(t, x) := W.removeNext()
for (y, z) in x.cond (t) do
addEdge(y, z)
end for
for y in x.succ do
addToken(t, y)
end for
end while
end procedure
t ∈ x:
addToken(t, x)
propagate()
t ∈ x =⇒ y ⊆ z:
if t ∈ x.sol then
addEdge(y, z)
propagate()
else
x.cond (t).add(y, z)
end if
Exercise 10.3: Show, step by step, how this algorithm finds the smallest solu-
tion for the constraints from the example program in Section 10.1.
Exercise 10.4: Explain how the algorithm can be extended to also support
unconditional inclusion constraints of the form x ⊆ y.
140 10 CONTROL FLOW ANALYSIS
To analyze the time complexity of this algorithm, we assume that the numbers
of tokens and variables are both O(n). This is clearly the case in closure analysis
of a program of size n. The number of constraints generated is O(n) for the first
kind and O(n2 ) for the second kind.
Each pair of a token and a variable can be added at most once to the worklist,
so the number of iterations of propagate is O(n2 ). For the same reason, W can be
implemented with constant time operations using an array. The addEdge func-
tion is called O(n2 ) times in total because there are O(n2 ) conditional constraints,
and each of them is processed either immediately or at most once by propagate.
The function addToken takes time O(1) if representing each x.sol as a bit vector.
The addEdge function calls addToken O(n3 ) times in total because there are
O(n2 ) edges and O(n) tokens. If representing x.succ using, for example, arrays
or hashtables, each operation on x.succ takes time O(n). This means that all calls
to addEdge take time O(n3 ) in total. Each iteration of propagate takes time O(n)
if excluding the time for addEdge, assuming that x.cond (t) is also implemented
using arrays or hashtables (and duplicates in x.cond (t) do not affect correctness).
Adding up, the total cost for the algorithm is O(n3 ). The fact that this seems like
a lower bound as well is referred to as the cubic time bottleneck.2
Numerous algorithmic variations and improvements are possible, including:
• cycle elimination (collapsing nodes if there is a cycle of inclusion con-
straints) [FFSA98],
• maintaining the worklist in topological order [PKH07],
• interleaving solution propagation and constraint processing [PKH07, PB09],
• using a shared bit vector representation to save memory [HT01b, SCD+ 13],
• type filtering [LH03, SCD+ 13],
• on-demand processing [HT01a],
• difference propagation [LH03, PKH04], and
• subsumed node compaction [RC00].
These techniques have mostly been studied in connection to points-to analysis
(see Chapter 11). For a control flow analysis perspective, see the survey by
Midtgaard [Mid12].
assuming that any function with the right number of arguments could be called.
However, we can do much better by performing a control flow analysis.
Our lattice is the powerset of the set of tokens containing X for every function
name X, ordered by subset inclusion. For every syntax tree node v we introduce
a constraint variable [[v]] denoting the set of functions v could point to. For a
function named f we have the constraint
f ∈ [[f ]]
[[E]] ⊆ [[X]]
and, finally, for computed function calls E(E1 ,. . . ,En ) we have for every defi-
nition of a function f with arguments a1f , . . . , anf and return expression Ef0 this
constraint:
foo(n,f) {
var r;
if (n==0) { f = ide; }
r = f(n);
return r;
}
main() {
var x,y;
x = input;
142 10 CONTROL FLOW ANALYSIS
The control flow analysis (with the special rule for direct calls) generates the
following constraints:
inc ∈ [[inc]]
dec ∈ [[dec]]
ide ∈ [[ide]]
[[ide]] ⊆ [[f]]
[[f(n)]] ⊆ [[r]]
inc ∈ [[f]] =⇒ [[n]] ⊆ [[i]] ∧ [[i+1]] ⊆ [[f(n)]]
dec ∈ [[f]] =⇒ [[n]] ⊆ [[j]] ∧ [[j-1]] ⊆ [[f(n)]]
ide ∈ [[f]] =⇒ [[n]] ⊆ [[k]] ∧ [[k]] ⊆ [[f(n)]]
[[input]] ⊆ [[x]]
[[foo(x,inc)]] ⊆ [[y]]
[[foo(x,dec)]] ⊆ [[y]]
foo ∈ [[foo]]
foo ∈ [[foo]] =⇒ [[x]] ⊆ [[n]] ∧ [[inc]] ⊆ [[f]] ∧ [[r]] ⊆ [[foo(x,inc)]]
foo ∈ [[foo]] =⇒ [[x]] ⊆ [[n]] ∧ [[dec]] ⊆ [[f]] ∧ [[r]] ⊆ [[foo(x,dec)]]
main ∈ [[main]]
[[inc]] = {inc}
[[dec]] = {dec}
[[ide]] = {ide}
[[f]] = {inc, dec, ide}
[[foo]] = {foo}
On this basis, we can construct the following interprocedural CFG for the pro-
gram, which then can be used as basis for subsequent interprocedural dataflow
analyses.
10.3 TIP WITH FIRST-CLASS FUNCTION 143
main foo
entry entry
var r
var x,y
if (n == 0)
x = input true
f = ide false
... = f(n)
true false
y = ... y = ...
return y
exit
144 10 CONTROL FLOW ANALYSIS
(a) Write the constraints that are produced by the control flow analysis for
this program.
(b) Compute the least solution to the constraints.
mon object-oriented languages, such as Java or C#, it is also useful, but the added
structure provided by the class hierarchy and the type system permits some
simpler alternatives. In the object-oriented setting the question is which method
implementations may be executed at a given method invocation site:
x.m(a,b,c)
The simplest solution is to scan the class library and select any method named m
whose signature accepts the types of the actual arguments. A better choice, called
Class Hierarchy Analysis (CHA), is to consider only the part of the class hierarchy
that is spanned by the declared type of x. A further refinement, called Rapid Type
Analysis (RTA), is to restrict further to the classes of which objects are actually
allocated. Yet another technique, called Variable Type Analysis (VTA), performs
intraprocedural control flow analysis while making conservative assumptions
about the remaining program.
These techniques are of course much faster than full-blown control flow
analysis, and for real-life programs they are often sufficiently precise.
Chapter 11
Pointer Analysis
The last TIP language feature we consider is pointers and dynamic memory
allocation (see the syntax in Section 2.1). For simplicity, we ignore records
(except in the last part of Section 11.2).
To illustrate the problem with pointers, assume we wish to perform a sign
analysis or constant propagation analysis of TIP code like this:
...
*x = 42;
*y = -87;
z = *x;
Here, the value of z depends on whether or not x and y are aliases, meaning that
they point to the same cell. Without knowledge of such aliasing information, it
quickly becomes impossible to produce useful dataflow and control flow analysis
results.
The first analyses that we shall study are flow-insensitive. The end result
of such an analysis is a function pt : Vars → P(Cells) that for each pointer
variable X returns the set pt(X) of abstract cells it may point to. This class of
analyses is called points-to analysis. We wish to perform a conservative analysis
that computes sets that may be too large but never too small.
Given such points-to information, many other facts can be approximated. If
we wish to know whether pointer variables x and y may be aliases, then a safe
answer is obtained by checking whether pt(x) ∩ pt(y) is nonempty.
The initial values of local variables are undefined in TIP programs, however,
for these flow-insensitive points-to analyses we assume that all the variables are
initialized before they are used. (In other words, these analyses are sound only
for programs that never read from uninitialized variables; see also Section 5.9.)
An almost-trivial analysis, called address taken, is to simply return all possible
abstract cells, except that any variable X is only included if the expression &X
occurs in the given program. This only suffices for very simple applications, so
more ambitious approaches are usually preferred. If we restrict ourselves to
typable programs, then any points-to analysis could be improved by removing
those cells whose types do not match the type of the pointer variable.
• X 1 = &X 2 ;
• X1 = X2;
• X 1 = *X 2 ;
• *X 1 = X 2 ;
• X = null;
The last two constraints are generated for every abstract cell c ∈ Cells, but
it is safe to skip the abstract cells that represent program variables X where
&X does not occur in the program. The null assignment is ignored, since it
corresponds to the trivial constraint ∅ ⊆ [[X]]. Notice that these constraints
match the requirements of the cubic algorithm from Section 10.2. The resulting
points-to function is defined simply as pt(p) = [[p]].
Consider the following example program fragment.
p = alloc null;
x = y;
x = z;
*p = z;
p = q;
q = &y;
x = *p;
p = &z;
alloc-1 ∈ [[p]]
[[y]] ⊆ [[x]]
[[z]] ⊆ [[x]]
c ∈ [[p]] =⇒ [[z]] ⊆ [[c]] for each c ∈ Cells
[[q]] ⊆ [[p]]
y ∈ [[q]]
c ∈ [[p]] =⇒ [[c]] ⊆ [[x]] for each c ∈ Cells
z ∈ [[p]]
where Cells = {p, q, x, y, z, alloc-1}. The points-to sets for the least solution is
quite precise in this case:
pt(p) = {alloc-1, y, z}
pt(q) = {y}
pt(x) = pt(y) = pt(z) = ∅
Note that although this algorithm is flow insensitive, the directionality of the
set inclusion constraints implies that the dataflow is still modeled with some
accuracy.
150 11 POINTER ANALYSIS
Exercise 11.2: Use Andersen’s algorithm to compute the points-to sets for the
variables in the following program fragment:
a = &d;
b = &e;
a = b;
*a = alloc null;
Exercise 11.3: Use Andersen’s algorithm to compute the points-to sets for the
variables in the following program fragment:
z = &x;
w = &a;
a = 42;
if (a > b) {
*z = &a;
y = &b;
} else {
x = &b;
y = w;
}
The following two exercises show that Andersen’s analysis is, perhaps sur-
prisingly, not the most precise possible flow-insensitive points-to analysis. This
is mostly of theoretical interest, however [BCSS11]. The lack of flow-sensitivity
(see Section 11.6) and path-sensitivity (and context sensitivity for interprocedu-
ral analysis, see Section 11.4) are much more important factors in practice.
11.2 ANDERSEN’S ALGORITHM 151
Exercise 11.4: Use Andersen’s algorithm to compute the points-to sets for the
variables in the following program fragment:
a = &b;
a = &x;
b = &c;
c = &d;
x = &y;
y = &z;
*a = **a;
Notice that the last statement has not been normalized. As observed by
Horwitz [Hor97], the normalization of the statement causes imprecision in
the following sense: According to Andersen’s analysis, which requires the
statement to be normalized using two temporary variables, b and x may be
aliases, but that is not possible with any TIP program that contains the above
set of assignments, no matter which control flow exists between them.
(a) Show that Andersen’s analysis concludes for this code that d may point
to c.
(b) Argue that for any TIP program that has this set of assignments, no
matter which control flow exists between them, d never points to c in
any execution. (Hint: a may point to b, and b may point to c, but not
simultaneously.)
Exercise 11.6: Recall from Exercise 5.26 that an analysis is distributive if all
its constraint functions are distributive. Show that Andersen’s analysis is not
distributive. (Hint: consider the constraint for the statement x=*y or *x=y.)
Let us now consider how to perform pointer analysis for TIP programs
that also use records (see Section 2.1 and Exercise 5.10). One simple approach
called field insensitive analysis is to treat all fields in a record as equivalent,
which is trivial for Andersen’s analysis. Another simple approach called field-
based analysis is to conversely treat each field name as a single global variable
without distinguishing between the different cells that contain the fields. Those
152 11 POINTER ANALYSIS
X = { X 1 :X 01 ,. . . , X k :X 0k }: [[X01 ]] ⊆ [[X.X1 ]] ∧ . . . ∧
[[X0k ]] ⊆ [[X.Xk ]]
X = alloc { X 1 :X 01 ,. . . , X k :X 0k }: alloc-i ∈ [[X]] ∧
[[X01 ]] ⊆ [[alloc-i.X1 ]] ∧ . . . ∧
[[X0k ]] ⊆ [[alloc-i.Xk ]]
X 1 = X 2 .X 3 : [[X2 .X3 ]] ⊆ [[X 1 ]]
X1 = X2: [[X2 ]] ⊆ [[X1 ]] ∧
[[X2 .f ]] ⊆ [[X1 .f ]] for each f ∈ Fields
X 1 = *X 2 : c ∈ [[X 2 ]] =⇒ ([[c]] ⊆ [[X 1 ]] ∧
[[c.f ]] ⊆ [[X 1 .f ]])
for
each c ∈ Cells and f ∈ Fields
X 1 .X 2 = X 3 :
see Exercise 11.7
(*X 1 ).X 2 = X 3 :
As usual, other syntactic forms can be normalized to these basic constructs, and
each constraint models the flow of pointer values.
As an example, for the (already normalized) program
i = null;
x = alloc {f: i}; // alloc-1
y = alloc {h: x}; // alloc-2
t = *y;
z = t.h;
alloc-1 ∈ [[x]]
[[i]] ⊆ [[alloc-1.f]]
alloc-2 ∈ [[y]]
[[x]] ⊆ [[alloc-2.h]]
c ∈ [[y]] =⇒ ([[c]] ⊆ [[t]] ∧ [[c.f ]] ⊆ [[t.f ]]) for each c ∈ Cells and f ∈ Fields
[[t.h]] ⊆ [[z]]
and after computing the least solution we get these points-to sets for the program
variables:
11.3 STEENSGAARD’S ALGORITHM 153
Exercise 11.7: Specify a suitable constraint rule for both forms of field write
statements, X 1 .X 2 =X 3 ; and (*X 1 ).X 2 = X 3 ;. Then explain how the ana-
lysis works on this program:
x = alloc {f: 1, g: 2};
y = alloc {h: x, g: 3};
z = alloc 4;
(*y).h = z;
In Java-like languages, the fields of an object are always accessed via a refer-
ence (i.e., a pointer) to the object. The Java syntax E.X for accessing a field X of
an object pointed to by an expression E corresponds to the TIP (or C) syntax
(*E).X.
Notice that at most two unification constraints are constructed for each pointer
operation. (This is different from Andersen’s analysis where the last two kinds
of statements require a constraint for each abstract cell.)
As usual, term constructors satisfy the general term equality axiom:
α1 = α2 =⇒ α1 = α2
For the example program from Section 11.2, Steensgaard’s algorithm gener-
ates the following constraints:
[[p]] = [[alloc-1]]
[[x]] = [[y]]
[[x]] = [[z]]
[[p]] = α1
[[z]] = α1
[[p]] = [[q]]
[[q]] = [[y]]
[[x]] = α2
[[p]] = α2
[[p]] = [[z]]
We can use the unification algorithm from Section 3.3 to find the resulting
points-to sets:
This result is less precise than what we obtained using Andersen’s algorithm,
but reached using the faster algorithm.
Notice that unification can never fail in Steensgaard’s analysis, since there is
only one kind of term constructor (unlike the type analysis in Chapter 3).
Exercise 11.9: Use Steensgaard’s algorithm to compute the points-to sets for
the two programs from Exercise 11.2 and Exercise 11.3. Are the results less
precise than when using Andersen’s analysis?
11.4 INTERPROCEDURAL POINTER ANALYSIS 155
[[X1 ]] = α ∧ [[X2 ]] = α
to
[[X1 ]] = [[X2 ]]
without affecting the analysis results?
(*x)(x);
X = X 0 (X 1 ,. . . , X n );
so that the involved expressions are all variables. Similarly, all return expressions
are assumed to be just variables, as in return X;.
f ∈ [[f ]]
Exercise 11.14: Continuing Exercise 11.13, can we still use the cubic algorithm
(Section 10.2) to solve the analysis constraints? If so, is the analysis time still
O(n3 ) where n is the size of the program being analyzed?
NN
where the bottom element NN means definitely not null and the top element >
represents values that may be null. We then form the following map lattice for
abstract states:
States = Cells → Null
For every CFG node v we introduce a constraint variable [[v]] denoting an element
from the map lattice. We shall use each constraint variable to describe an abstract
state for the program point immediately after the node.
For all nodes that do not involve pointer operations we have the constraint:
where G
JOIN (v) = [[w]]
w∈pred(v)
where G
load (σ, X1 , X2 ) = σ[X1 7→ σ(α)]
α∈pt(X2 )
Similar reasoning gives constraints for the other operations that affect pointer
variables:
X = alloc P : [[v]] = JOIN (v)[X 7→ NN, alloc-i 7→ >]
X 1 = &X 2 : [[v]] = JOIN (v)[X 1 7→ NN]
X1 = X2: [[v]] = JOIN (v)[X 1 7→ JOIN (v)(X 2 )]
X = null: [[v]] = JOIN (v)[X 7→ >]
Exercise 11.15: Explain why the above four constraints are monotone and
sound.
where
store(σ, X1 , X2 ) = σ [α 7→ σ(α) t σ(X2 ) ]
α∈pt(X1 )
Exercise 11.16: Show an alternative constraint for load operations using weak
update, together with an example program where the modified analysis then
gives a result that is less precise than the analysis presented above.
where x, y, and z are program variables. We will seek to answer questions about
disjointness of the structures contained in program variables. In the example
above, x and y are not disjoint whereas y and z are. Such information may be
useful, for example, to automatically parallelize execution in an optimizing com-
piler. For such analysis, flow-insensitive reasoning is sometimes too imprecise.
However, we can create a flow-sensitive variant of Andersen’s analysis.
We use a lattice of points-to graphs, which are directed graphs in which the
nodes are the abstract cells for the given program and the edges correspond
to possible pointers. Points-to graphs are ordered by inclusion of their sets of
edges. Thus, ⊥ is the graph without edges and > is the completely connected
graph. Formally, our lattice for abstract states is then
ordered by the usual subset inclusion. For every CFG node v we introduce
a constraint variable [[v]] denoting a points-to graph that describes all possible
stores at that program point. For the nodes corresponding to the various pointer
manipulations we have these constraints:
σ ↓ x = {(s, t) ∈ σ | s 6= x}
160 11 POINTER ANALYSIS
p q
alloc−3 alloc−4
x y
alloc−1 alloc−2
From this result we can safely conclude that x and y will always be disjoint.
Note that this analysis also computes a flow sensitive points-to map that for
each program point v is defined by:
pt(p) = {t | (p, t) ∈ [[v]]}
This analysis is more precise than Andersen’s algorithm, but clearly also more
expensive to perform. As an example, consider the program:
x = &y;
x = &z;
After these statements, Andersen’s algorithm would predict that pt(x) = {y, z}
whereas the flow-sensitive analysis computes pt(x) = {z} for the final program
point.
11.7 ESCAPE ANALYSIS 161
main() {
var p;
p=baz(); *p=1;
return *p;
}
Having performed a points-to analysis, we can easily perform an escape analysis
to catch such errors. We just need to check that the possible cells for return
expressions in the points-to graph cannot reach arguments or variables defined
in the function itself, since all other locations must then necessarily reside in
earlier frames on the invocation stack.
Chapter 12
Abstract Interpretation
In the preceding chapters we have used the term soundness of an analysis only
informally: if an analysis is sound, the properties it infers for a given program
hold in all actual executions of the program. The theory of abstract interpreta-
tion provides a solid mathematical foundation for what it means for an analysis
to be sound, by relating the analysis specification to the formal semantics of
the programming language. Another use of abstract interpretation is for un-
derstanding whether an analysis design, or a part of an analysis design, is as
precise as possible relative to a choice of analysis lattice and where imprecision
may arise. The fundamental ideas of abstract interpretation were introduced by
Cousot and Cousot in the 1970s [CC76, CC77, CC79b].
as the “concrete semantics” and the analysis specification is called the “abstract semantics”.
164 12 ABSTRACT INTERPRETATION
ConcreteStates = Vars ,→ Z
For every CFG node v we have a constraint variable that ranges over sets of
concrete states:
{[v]} ⊆ ConcreteStates
The idea is that {[v]} shall denote the set of concrete states that are possible
at the program point immediately after the instruction represented by v, in
some execution of the program. This is called a collecting semantics, because
it “collects” the possible states. In Exercises 12.59–12.64 we shall study other
kinds of collecting semantics that collect relevant information suitable for other
analysis, such as live variables analysis. We choose to focus on the program
point immediately after the instruction of the CFG node, instead of the program
point before, to align with our sign analysis and the other forward analyses from
Chapter 5.
Consider this simple program as an example:
var x;
x = 0;
while (input) {
x = x + 2;
}
Its CFG looks as follows, where the bullets represent the program points that
have associated constraint variables.
entry
var x
x = 0
input
true
false
x = x + 2
exit
The solution we are interested in maps the constraint variable {[x = 0]} to the
single state where x is zero, {[x = x + 2]} is mapped to the set of all states where
x is a positive even number, and similarly for the other program points.
As a first step, we define some useful auxiliary functions, ceval , csucc, and
CJOIN , that have a close connection to the auxiliary functions used in the
2 We use the notation A ,→ B to denote the set of partial functions from A to B.
12.1 A COLLECTING SEMANTICS FOR TIP 165
specification of the sign analysis, but now for concrete executions instead of
abstract executions.
The function ceval : ConcreteStates × Exp → P(Z) gives the semantics of
evaluating an expression E ∈ Exp relative to a concrete state ρ ∈ ConcreteStates,
which results in a set of possible integer values, defined inductively as follows:3
ceval (ρ, X) = {ρ(X)}
ceval (ρ, I) = {I}
ceval (ρ, input) = Z
ceval (ρ, E1 + E2 ) = {z1 + z2 | z1 ∈ ceval (ρ, E1 ) ∧ z2 ∈ ceval (ρ, E2 )}
ceval (ρ, E1 / E2 ) = {z1 / z2 | z1 ∈ ceval (ρ, E1 ) ∧ z2 ∈ ceval (ρ, E2 )}
Evaluation of the other binary operators is defined similarly. In this simple
subset of TIP we consider here, evaluating an expression cannot affect the values
of the program variables. Also note that division by zero simply results in the
empty set of values.
We overload ceval such that it also works on sets of concrete states,
ceval : P(ConcreteStates) × Exp → P(Z):
[
ceval (R, E) = ceval (ρ, E)
ρ∈R
For a CFG node v, CJOIN (v) denotes the set of states at the program point
immediately before the instruction represented by v, relative to the states at the
program points after the relevant other nodes according to the csucc function:4
and for simplicity we do not restrict the numeric computations to, for example, 64 bit signed integers.
4 Note that CJOIN (v) is a function of all the constraint variables {[v ]}, . . . , {[v ]} for the entire
1 n
program, just like JOIN (v) is a function of all the constraint variables [[v1 ]], . . . , [[vn ]].
166 12 ABSTRACT INTERPRETATION
Exercise 12.1: Convince yourself that this definition of CJOIN makes sense,
especially for the cases where v is a node with multiple incoming edges (like
input in the example on page 164) or it is the first node of a branch (like x =
x + 2 in the example).
This rule formalizes the runtime behavior of assignments: For every state ρ
that may appear immediately before executing the assignment, the state after
the assignment is obtained by overwriting the value of X with the result of
evaluating E.
If v is a variable declaration, var X1 , . . . ,Xn , we use this rule:
The only possible initial state at entry nodes is the partial map that is undefined
for all program variables, denoted []:
{[entry]} = {[]}
For all other kinds of nodes, we have this trivial constraint rule:
Notice the resemblance with the analysis constraints from Section 5.1.
Exercise 12.2: Define a suitable constraint rule that expresses the semantics
of assert statements (see Section 7.1) in our collecting semantics. (For this
exercise, think of assert statements as written explicitly by the programmers
anywhere in the programs, not just for use in control sensitive analysis.)
n
The constraint for each CFG node v is a function cf v : P(ConcreteStates) →
P(ConcreteStates), much like the analysis of a program is expressed as an equa-
tion system in Sections 4.4 and 5.1. In the nsame way, we can combine n the n
functions into one, cf : P(ConcreteStates) → P(ConcreteStates) , defined
by
cf (x1 , . . . , xn ) = cf v1 (x1 , . . . , xn ), . . . , cf vn (x1 , . . . , xn )
in which case the equation system looks like
x = cf (x)
n
where x ∈ P(ConcreteStates) . Thus, cf captures the semantics of the pro-
gram as one single function, in the same way af in Chapter 5 (page 54) captures
the analysis of the program as one single function.
Exercise 12.3: Show that the constraints defined by the rules above are mono-
tone. Then show that the combined constraint function cf is also monotone.
Exercise 12.4:
(a) Check that both of the above solutions are indeed solutions to the con-
straints for this particular program.
(b) Give an example of yet another solution.
(c) Argue that solution 1 is the least of all solutions.
var a,b,c;
a = 42;
b = 87;
if (input) {
c = a + b;
} else {
c = a - b;
}
For this program, at the program points immediately after the assignment
b = 87, immediately after the assignment c = a - b (at the end of the else
branch), and the exit, the following sets of concrete states are possible according
to the collecting semantics:
Exercise 12.5: Check that the least fixed point of the semantic constraints for
the program is indeed these sets, for the three program points.
In comparison, the sign analysis specified in Section 5.1 computes the following
abstract states for the same program points:
In this specific case the analysis result is almost the best we could hope for,
with that choice of analysis lattice. Notice that the abstract value of c at the
program point after c = a - b is >, although the only possible value in concrete
executions is −45. This is an example of a conservative analysis result.
As shown above, Tarski’s fixed-point theorem ensures that the collecting
semantics is well defined – even though it is generally non-computable. The
variant of Kleene’s fixed-point theorem from Chapter 4 (page 47) is not strong
enough to ensure this property, as the lattice has infinite height.
Exercise 12.6: Show that the lattice (P(ConcreteStates))n has infinite height.
where A is totally ordered, i.e., ∀x, y ∈ A : x v y ∨ y v x. For a variant of Kleene’s theorem that
only requires monotonicity, see Exercise 4.32.
12.2 ABSTRACTION AND CONCRETIZATION 169
Exercise 12.8: Prove that if L is a complete lattice (not necessarily with finite
height) and the function f : L → L is a complete join morphism, then lfp(f ) =
i
F
i≥0 f (⊥) is a unique least fixed point for f .
Exercise 12.9: Show that the constraints defined by the rules above for the col-
lecting semantics, including the combined constraint function cf , are indeed
complete join morphisms.
For use later in this chapter, let us introduce the notation {[P ]} = lfp(cf ) and
[[P ]] = lfp(af ) where cf is the semantic constraint function and af is the analysis
constraint function for a given program P . In other words, {[P ]} denotes the
semantics of P , and [[P ]] denotes the analysis result for P .
αa : P(Z) → Sign
αb : P(ConcreteStates) →
n States
αc : P(ConcreteStates) → States n
170 12 ABSTRACT INTERPRETATION
As before, P(Z) is the powerset lattice defined over the set of integers ordered
by subset, Sign is the sign lattice from Section 4.1, we define ConcreteStates =
Vars → Z and State = Vars → Sign as in Sections 12.1 and 4.1, respectively, and
n is the number of CFG nodes. The functions are defined as follows, to precisely
capture the informal descriptions given earlier:
⊥ if D is empty
+ if D is nonempty and contains only positive integers
αa (D) = - if D is nonempty and contains only negative integers
0 if D is nonempty and contains only the integer 0
> otherwise
for any D ∈ P(Z)
Exercise 12.12: Argue that the three functions γa , γb , and γc from the sign
analysis example are monotone.
γc γc
y
x
αc αc
n n
P(ConcreteStates) States n P(ConcreteStates) States n
Exercise 12.13: Show that all three pairs of abstraction and concretization
functions (αa , γa ), (αb , γb ), and (αc , γc ) from the sign analysis example are
Galois connections. (See also Exercise 12.29.)
Exercise 12.14: Show that αa ◦γa , αb ◦γb , and αc ◦γc are all equal to the identity
function, for the three pairs of abstraction and concretization functions from
the sign analysis.
Exercise 12.15: Argue that γ ◦ α is typically not the identity function, when
α : L1 → L2 is an abstraction function and γ : L2 → L1 is the associated
concretization function for some analysis. (Hint: consider αa and γa from the
sign analysis example.)
172 12 ABSTRACT INTERPRETATION
if and only if
∀x ∈ L1 , y ∈ L2 : α(x) v y ⇐⇒ x v γ(y)
This theorem is useful for some of the later exercises in this section.
We shall use this result in the soundness argument in Section 12.3. Not
surprisingly, the dual property also holds: γ satisfies γ( B) = b∈B γ(b) for
F F
every B ⊆ L2 when α and γ form a Galois connection.
Exercise 12.19: Show that if α and γ form a Galois connection, then α(⊥) = ⊥
and γ(>) = >.
We have argued that the Galois connection property is natural for any rea-
sonable pair of an abstraction function and a concretization function, including
those that appear in our sign analysis example. The following exercise tells
us that it always suffices to specify either α or γ, then the other is uniquely
determined if requiring that they together form a Galois connection.
12.2 ABSTRACTION AND CONCRETIZATION 173
The result from Exercise 12.20 means that once the analysis designer has
specified the collecting semantics and the analysis lattice and constraint rules,
then the relation between the semantic domain and the analysis domain may
be specified using an abstraction function α (resp. a concretization function γ),
and then the associated concretization function γ (resp. abstraction function α)
is uniquely determined – provided that one exists such that the two functions
form a Galois connection.
The dual property also holds: if γ satisfies γ( B) = b∈B γ(b) for every
F F
B ⊆ L2 then there exists a function α : L1 → L2 such that α and γ form a
Galois connection.
The following exercise demonstrates that the Galois connection property can
be used as a “sanity check” when designing analysis lattices.
174 12 ABSTRACT INTERPRETATION
Exercise 12.22: Instead of using the usual Sign lattice from Section 5.1, assume
we chose to design our sign analysis based on this lattice:
>
0- 0+
At first, this may seem like a reasonable lattice for an analysis, but there is
something strange about it: How should we define eval (σ, 0)? Or equivalently,
how should we define αa ({0})? We could somewhat arbitrarily choose either
0+ or 0-.
Show that, with this choice of lattice, there does not exist an abstraction
function αa such that αa and γa form a Galois connection.
x = input > 3;
if (!x) {
output x;
}
After the first statement, the abstract value of x is 0+ (since > yields 0 or 1), so
at the output statement, the abstract value of x is then 0+. This means that the
analysis is precise enough to be able to conclude that the output value is definitely
not a negative number. Now consider a situation where the analysis designer
wants to increase analysis precision generally by adding control sensitivity to the
analysis, using the approach presented in Chapter 7. With the representation of
12.2 ABSTRACTION AND CONCRETIZATION 175
This rule soundly models the fact that only executions where x has the value
0 satisfy the branch condition. However, with this change of the analysis, the
abstract value of x at the output statement is now 0-, so the analysis is no longer
able to conclude that the output value is definitely not a negative number. In
other words, even though the analysis has apparently become more precise by
adding a simple form of control sensitivity, it is less precise for some programs,
which may be counterintuitive.
Exercise 12.23: Continuing Exercise 12.22, let us add a lattice element 0, below
0- and 0+ and above ⊥, with γa (0) = {0}. Show that with this modification,
an abstraction function αa exists such that αa and γa form a Galois connection.
Exercise 12.24: In interval analysis (Section 6.1) we use the domain P(Z)
representing sets of concrete values and the domain Intervals representing
abstract values. What are the abstraction function and the concretization
function between these domains? Do they form a Galois connection?
For the sign analysis (with the ordinary Sign lattice) we can also specify the
connection between concrete values in Z and abstract values in Sign using a
representation function, β : Z → Sign:
+ if d > 0
β(d) = - if d < 0
0 if d = 0
The abstraction function αa can then be induced from the representation func-
tion: G
αa (D) = {β(d) | d ∈ D}
Exercise 12.25: Show that this definition of αa coincides with the one from
page 169.
The next two exercises show that we can build Galois connections using
product lattices and map lattices (see Section 4.3).
Exercise 12.29: Use the results of Exercises 12.27 and 12.28 to give a simpler
solution to Exercise 12.13.
12.3 Soundness
In other words, soundness means that the analysis result over-approximates the
abstraction of the semantics of the program. For the sign analysis, the property
can be illustrated like this:
12.3 SOUNDNESS 177
[[P ]]
αc
{[P ]}
n
P(ConcreteStates) States n
For the simple TIP example program considered in Section 12.1, the sound-
ness property is indeed satisfied (here showing the information just for the
program point immediately after the c = a - b statement):
αc (. . . , {[c = a - b]}, . . . ) =
αc (. . . , [a 7→ 42, b 7→ 87, c 7→ −45], . . . ) v
(. . . , [a 7→ +, b 7→ +, c 7→ >], . . . ) =
(. . . , [[c = a - b]], . . . )
If we specify the relation between the two domains using concretization
functions instead of using abstraction functions, we may dually define soundness
as the property that the concretization of the analysis result over-approximates
the semantics of the program:
γc
[[P ]]
{[P ]}
n
P(ConcreteStates) States n
Exercise 12.30: Show that if α and γ form a Galois connection, then the two
definitions of soundness stated above are equivalent.
(Hint: see Exercise 12.17.)
178 12 ABSTRACT INTERPRETATION
Exercise 12.31: Prove that eval is a sound abstraction of ceval , in the sense
defined above.
(Hint: Use induction in the structure of the TIP expression. As part of the
proof, you need to show that each abstract operator is a sound abstraction
the corresponding concrete operator, for example for the addition operator:
αa ({z1 + z2 | z1 ∈ D1 ∧ z2 ∈ D2 }) v αa (D1 ) b
+ αa (D2 ) for all sets D1 , D2 ⊆
Z.)
Exercise 12.32: Prove that succ is a sound abstraction of csucc and that JOIN
is a sound abstraction of CJOIN , in the sense defined above.
n
Let cfv : P(ConcreteStates) → P(ConcreteStates) and af v : States n →
States denote v’s constraint function from the semantics and the analysis, re-
spectively, for every CFG node v. For example, if v represents an assignment
statement X = E, we have:
cfv ({[v1 ]}, . . . , {[vn ]}) = ρ[X 7→ z] ρ ∈ CJOIN (v) ∧ z ∈ ceval (ρ, E)
af v ([[v1 ]], . . . , [[vn ]]) = σ[X 7→ eval (σ, E)] where σ = JOIN (v)
The function af v is a sound abstraction of cfv if the following property holds for
all R1 , . . . , Rn ⊆ ConcreteStates:
αb (cfv (R1 , . . . , Rn )) v af v (αb (R1 ), . . . , αb (Rn ))
12.3 SOUNDNESS 179
cf αc af
αc
n
P(ConcreteStates) States n
Exercise 12.33: Prove that each kind of CFG node v, the sign analysis con-
straint af v is a sound abstraction of the semantic constraint cfv . (The most
interesting case is the one where v is an assignment node.) Then use that
result to prove that af is a sound abstraction of cf .
α0 ◦ cg v ag ◦ α
As an example for our sign analysis, eval being a sound abstraction of ceval
means that αa ◦ceval v eval ◦αb (for a fixed expression), which can be illustrated
as follows:
180 12 ABSTRACT INTERPRETATION
αb
ceval
αa
P(Z) Sign
Intuitively, when starting from a set of concrete states, if we first abstract the
states and then evaluate abstractly with eval we get an abstract value that over-
approximates the one we get if we first evaluate concretely with ceval and then
abstract the values.
With the result of Exercise 12.33, it follows that the sign analysis is sound
with respect to the semantics and the abstraction/concretization functions, as
shown next.
Recall that the analysis result for a given program P is computed as [[P ]] =
lfp(af ) (if not using widening, which we discuss later in this section). By the
fixed-point theorems from Section 12.1, the semantics of P is similarly given
by {[P ]} = lfp(cf ). According to the definition of soundness, we thus need to
show that α(lfp(cf )) v lfp(af ) (if the connection between the lattices is specified
using an abstraction function α) or lfp(cf ) v γ(lfp(af )) (if using a concretization
function γ, cf. Exercise 12.17).
The central result we need is the following soundness theorem:
We know from Tarski’s fixed point theorem that lfp(cf ) and lfp(af ) are well-
defined. As lfp(af ) is a fixed point of af we have af (lfp(af )) = lfp(af ), so
γ(af (lfp(af ))) = γ(lfp(af )). Since af is a sound abstraction of cf we then have
cf (γ(lfp(af ))) v γ(lfp(af )). This means that γ(lfp(af )) ∈ {x ∈ L1 | cf (x) v x},
so by Tarski’s fixed point theorem, lfp(cf ) v γ(lfp(af )).
To summarize, a general recipe for specifying and proving soundness of an
analysis consists of the following steps:
1. Specify the analysis, i.e. the analysis lattice and the constraint generation
rules, and check that all the analysis constraint functions are monotone
(as we did for the sign analysis example in Section 5.1).
2. Specify the collecting semantics, and check that the semantic constraint
functions are monotone (as we did for the sign analysis example in Sec-
tion 12.1). The collecting semantics must capture the desired aspects of
concrete execution, such that it formalizes the ideal properties that the
analysis is intended to approximate. For the sign analysis example, we
designed the collecting semantics such that it collects the reachable states
for every program point; other analyses may need other kinds of collecting
semantics (see Section 12.6).
3. Establish the connection between the semantic lattice and the analysis
lattice (as we did for the sign analysis example in Section 12.2), either
by an abstraction function or by a concretization function. Alternatively,
one may use a representation function (Exercise 12.26) and induce the
abstraction and concretization functions. For the sign analysis example, the
lattices are defined in three layers, leading to the definitions of αa , αb , αc
and γa , γb , γc . By ensuring that the Galois connection property is satisfied,
it is only a matter of taste whether one chooses to specify this connection
using abstraction functions or using concretization functions, thanks to the
Galois connection dualities we have seen in Section 12.2. (However, notice
that the soundness theorem does not require the presence of a Galois
connection, but merely that the abstraction and concretization functions
are monotone.)
5. Soundness of the analysis then follows from the soundness theorem stated
above.
For analyses that use widening, we know from Exercise 6.9 that the ana-
lysis result is always a safe approximation of the least solution to the analysis
constraints. This means that the soundness theorem can also be used for such
analyses.
182 12 ABSTRACT INTERPRETATION
Exercise 12.35: Prove that the interval analysis (Section 6.1) with widening
(using the definition of ∇ from page 86) is sound with respect to the collecting
semantics from Section 12.1.
12.4 Optimality
Assume we are developing a new analysis, and that we have chosen an analysis
lattice and the rules for generating analysis constraints for the various program-
ming language constructs. To enable formal reasoning about the soundness and
precision of the analysis, we have also provided a suitable collecting semantics for
the programming language (as in Section 12.1) and abstraction/concretization
functions that define the meaning of the analysis lattice elements (as in Sec-
tion 12.2). Furthermore, assume we have proven that the analysis is sound using
the approach from Section 12.3. We may now ask: Are our analysis constraints as
precise as possible (yet sound), relative to the chosen analysis lattice?
As in the previous section, let α : L1 → L2 be an abstraction function where
L1 is the lattice for a collecting semantics and L2 is the lattice for an analysis, such
that α and γ form a Galois connection, and consider two functions cf : L1 → L1
and af : L2 → L2 that represent, respectively, the semantic constraints and
the analysis constraints for a given program. We say that af is the optimal6
abstraction of cf if
af = α ◦ cf ◦ γ
(which can also be written: af (b) = α(cf (γ(b))) for all b ∈ L2 ). Using the lattices
and abstraction/concretization functions from the sign analysis example, this
property can be illustrated as follows.
6 In the literature on abstract interpretation, the term “best” is sometimes used instead of “opti-
mal”.
12.4 OPTIMALITY 183
cf αc af
γc
n
P(ConcreteStates) States n
(Compare this with the illustration of soundness from page 179.) To see that α ◦
cf ◦ γ is indeed the most precise monotone function that is a sound abstraction of
cf , assume g : L2 → L2 is some monotone function that is a sound abstraction of
cf , that is, α(cf (a)) v g(α(a)) for all a ∈ L1 . Then for all a0 ∈ L1 , α(cf (γ(a0 ))) v
g(α(γ(a0 ))) v g(a0 ). The last inequality holds because α ◦ γ is reductive and g is
monotone. Thus, α ◦ cf ◦ γ v g, meaning that α ◦ cf ◦ γ is the most precise such
function.
Notice what the optimality condition tells us: We can obtain the best possible
(i.e., most precise yet sound) abstraction of cf for a given analysis lattice element
y by first concretizing y, then applying the semantic function cf , and finally
abstracting. Unfortunately this observation does not automatically give us
practical algorithms for computing optimal abstraction, but it enables us to
reason about the precision of our manually specified analysis constraints.
The above definition of optimality focuses on af and cf , but it can be gen-
eralized to all levels of the analysis as follows. Assume L1 and L01 are lattices
used in a collecting semantics and L2 and L02 are lattices used for an analysis.
Let α : L1 → L2 and α0 : L01 → L02 be abstraction functions, and let γ : L2 → L1
and γ 0 : L02 → L01 be concretization functions such that α, γ, α0 , and γ 0 form two
Galois connections. Consider two functions cg : L1 → L01 and ag : L2 → L02 . We
say that ag is the optimal abstraction of cg if the following property holds:
ag = α0 ◦ cg ◦ γ
Let us look at some examples from the sign analysis. First, it is easy to see that
* (page 53) is the optimal abstraction of
our definition of abstract multiplication b
concrete multiplication, denoted “·”:
*s2 = αa γa (s1 ) · γa (s2 )
s1 b
for any s1 , s2 ∈ Sign, where we overload the · operator to work on sets of integers:
D1 · D2 = {z1 · z2 | z1 ∈ D1 ∧ z2 ∈ D2 } for any D1 , D2 ⊆ Z.
Despite the result of the previous exercise, the eval function from Section 5.1
is not the optimal abstraction of ceval . Here is a simple counterexample: Let
σ ∈ States such that σ(x) = > and consider the TIP expression x - x. We then
have
eval (σ, x - x) = >
while
αa ceval (γb (σ), x - x) = 0
(This is essentially the same observation as the one in Exercise 5.9, but this time
stated more formally.) Interestingly, defining the eval function inductively and
compositionally in terms of optimal abstractions does not make the function
itself optimal.
Exercise 12.37: Assume we only work with normalized TIP programs (as
in Exercise 2.2). Give an alternative computable definition of eval for sign
analysis (i.e., an algorithm for computing eval (σ, E) for any abstract state σ
and normalized TIP expression E), such that eval is the optimal abstraction
of ceval .
Recall the definition of the abstract operators in interval analysis from Chap-
ter 6 (page 80):
op([l
b 1 , h1 ], [l2 , h2 ]) = [ min x op y, max x op y]
x∈[l1 ,h1 ],y∈[l2 ,h2 ] x∈[l1 ,h1 ],y∈[l2 ,h2 ]
Exercise 12.39: Which of the other abstractions used in interval analysis (Sec-
tion 6.1) are optimal?
To be able to reason about optimality of the abstractions used in, for example,
live variables analysis or reaching definitions analysis, we first need a style of
collecting semantics that is suitable for those analyses, which we return to in
Section 12.6.
12.5 Completeness
As usual in logics, the dual of soundness is completeness. In Section 12.3 we
defined soundness of an analysis for a program P as the property α({[P ]}) v [[P ]].
Consequently, it is natural to define that an analysis is complete for P if the
following property holds:
[[P ]] v α({[P ]})
12.5 COMPLETENESS 185
x = input;
y = x;
z = x - y;
Let σ denote the abstract state after the statement y = x such that σ(x) = σ(y) =
>. Any sound abstraction of the semantics of the single statement z = x - y
will result in an abstract state that maps z to >, but the answer 0 would be
more precise and still sound in the analysis result for the final program point.
Intuitively, the analysis does not know about the correlation between x and y.
For this specific example program, we could in principle improve analysis
precision by changing the constraint generation rules to recognize the special
pattern consisting of the statement y = x followed by z = x - y. Instead of
such an ad hoc approach to gain precision, relational analysis (see Chapter 7) is
usually a more viable solution.
In Section 12.3 we observed (Exercise 12.30) that analysis soundness could
equivalently be defined as the property {[P ]} v γ([[P ]]). However, a similar
equivalence does not hold for completeness, as shown by the following two
exercises.
Exercise 12.40: We know from Exercise 12.17 that if α and γ form a Galois
connection then α(x) v y ⇐⇒ x v γ(y) for all x, y. Prove (by showing a
counterexample) that the converse property does not hold, i.e. α(x) w y ⇐⇒
6
x w γ(y).
(Hint: consider αa and γa from the sign analysis.)
Exercise 12.41: Give an example of a program P such that [[P ]] v α({[P ]}) and
γ([[P ]]) 6v {[P ]} for sign analysis.
Thus, {[P ]} = γ([[P ]]) is a (much) stronger property than α({[P ]}) = [[P ]].
If {[P ]} = γ([[P ]]) is satisfied, the analysis captures exactly the semantics of
P without any approximation; we say that the analysis is exact for P . Every
7 The literature on abstract interpretation often uses the term “complete” for what we call “sound
Exercise 12.42:
(a) Explain why the sign analysis is both complete and exact for this pro-
gram:
var x;
x = 0;
(b) Explain why the sign analysis is complete but not exact for this program:
var x;
x = 1;
ag ◦ α v α0 ◦ cg
(Compare this with the definition of soundness of abstractions from page 179.)
Again, let us consider sign analysis as example. In Section 12.4 we saw that
abstract multiplication is optimal. In fact, it is also complete:
for any D1 , D2 ⊆ Z. This tells us that the analysis, perhaps surprisingly, never
loses any precision at multiplications.
For abstract addition, the situation is different, as shown in the next exercise.
8 Theseobservations show that we could in principle have chosen define the concept of complete-
ness using the concretization function γ instead of using the abstraction function α, but that would
have been much less useful.
12.5 COMPLETENESS 187
Exercise 12.44: For which of the operators -, /, >, and == is sign analysis
complete? What about input?
(This result follows directly from Exercise 12.18.) Intuitively this means, perhaps
surprisingly, that joining abstract information, for example when the branches
merge after an if statement, is never to blame for precision losses.
Exercise 12.45: In sign analysis, is the analysis constraint function for assign-
ments af X=E a complete abstraction of the corresponding semantic constraint
function cf X=E , given that E is an expression for which eval is complete?
For abstractions that are sound, completeness implies optimality (but not
vice versa, cf. exercises 12.36 and 12.43):
Exercise 12.46: Prove that if ag is sound and complete with respect to cg, it is
also optimal.
We have seen in Section 12.4 that for any Galois connection, there exists an
optimal (albeit often non-computable) abstraction of every concrete operation.
The same does not hold for soundness and completeness, as shown in the
following exercise.
Exercise 12.47: Prove that there exists a sound and complete abstraction ag
of a given concrete operation cg if and only if the optimal abstraction of cg is
sound and complete.
(We have seen examples of abstractions that are optimal but not sound and
complete, so this result implies that sound and complete abstractions do not
always exist.)
Exercise 12.48: Prove that if af is sound and complete with respect to cf then
α({[P ]}) = [[P ]], where cf is the semantic constraint function and af is the
analysis constraint function for a given program P , and α is the abstraction
function.
This theorem relies on the fact that sound and complete abstractions are
closed under composition, unlike optimal abstractions, as seen in the following
188 12 ABSTRACT INTERPRETATION
exercise.
Exercise 12.49:
(a) Prove that if af is sound and complete with respect to cf then af ◦ af is
sound and complete with respect to cf ◦ cf . (Hint: see your solution to
Exercise 12.48.)
(b) Prove that the following statement does not hold in general: If af is
the optimal abstraction of cf then af ◦ af is the optimal abstraction of
cf ◦ cf . (Hint: consider the example program on page 185.)
Since [[P ]] = i≥0 af i (⊥) by the fixed-point theorem, part (b) of this result
F
shows that even when the analysis constraint af for a given program P is
the most precise possible (relative to the chosen analysis lattice), the analysis
result [[P ]] may not be the most precise possible (that can be expressed in the
analysis lattice). In contrast, part (a) is the central property in the proof of
the soundness and completeness theorem in Exercise 12.48.
Exercise 12.50: Which of the abstractions used in interval analysis (Section 6.1)
are complete? In particular, is abstract addition complete?
There are not many programs for which our simple sign analysis is complete
and gives a nontrivial analysis result, so to be able to demonstrate how these
observations may be useful, let us modify the analysis to use the following
slightly different lattice instead of the usual Sign lattice from Section 5.1.
>
0- 0+
∅ if s = ⊥
{0} if s = 0
γa (s) = {0, 1, 2, 3, . . . } if s = 0+
{0, −1, −2, −3, . . . }
if s = 0-
Z if s = >
From Exercise 12.22 we know that an abstraction function αa exists such that αa
and γa form a Galois connection.
Exercise 12.52: Give a definition of eval that is optimal for expressions of the
form t * t where t is any program variable (recall Exercise 12.37).
The remaining abstract operators can be defined similarly, and the rest of
the eval function and other analysis constraints can be reused from the ordinary
sign analysis.
The modified sign analysis concludes that output is 0+ for the following
small program:
x1 = input;
x2 = input;
y1 = x1 * x1;
y2 = x2 * x2;
output y1 + y2;
Exercise 12.53: Explain why the modified sign analysis is sound and complete
for this program.
Assume the analysis is built to raise an alarm if the output of the analyzed
program is a negative value for some input. In this case, it will not raise an
alarm for this program, and because we know the analysis is sound (over-
approximating all possible behaviors of the program), this must be the correct
result. Now assume instead that the analysis is built to raise an alarm if the
output of the analyzed program is a positive value for some input. In this case,
it does raise an alarm, and because we know the analysis is complete for this
program (though not for all programs), this is again the correct result – there
must exist an execution of the program that outputs a positive value; we can
trust that the alarm is not a false positive.
Exercise 12.54: Design a relational sign analysis that is sound and complete
for the three-line program from page 185.
Exercises 12.50 and 12.54 might suggest that increasing analysis precision
generally makes an analysis complete for more programs, but that is not the
case: The trivial analysis that uses a one-element analysis lattice is sound and
complete for all programs, but it is obviously useless because its abstraction
discards all information about any analyzed program.
190 12 ABSTRACT INTERPRETATION
The semantics of variable declaration nodes can be defined similarly. All other
kinds of nodes do not change the state:
that the traces reach the program exit.) More formally, we define h[P ]i as the
least solution to the following two constraints (similarly to how we defined [[P ]]
and {[P ]} earlier):
(entry, []) ∈ h[P ]i
which contains the information that the value of x is the same at the entry and
the exit, in any execution.
Exercise 12.55: Describe the trace semantics of the program from page 164.
Exercise 12.56: Explain why it makes sense to define the trace semantics as
the least solution to the constraints. (Hint: see the discussion for the reachable
states collecting semantics on page 167.)
defined by
αt (T ) = (R1 , . . . , Rn )
Exercise 12.57: Show that αt (as defined above) is indeed a complete join
morphism.
192 12 ABSTRACT INTERPRETATION
The existence of this Galois connection shows that the domain of the reach-
able states collecting semantics is in some sense an abstraction of the domain of
the trace semantics.
The following exercise shows that composition of Galois connections leads
to new Galois connections.
Exercise 12.58: Let α1 : L1 → L2 , γ1 : L2 → L1 , α2 : L2 → L3 , and γ2 : L3 →
L2 . Assume both (α1 , γ1 ) and (α2 , γ2 ) are Galois connections. Prove that
(α2 ◦ α1 , γ1 ◦ γ2 ) is then also a Galois connection.
In Section 12.2 we
nhave established a Galois connection between the domain
P(ConcreteStates) of the reachable states collecting semantics and the do-
main States n of the sign analysis, and we have now also established a Galois
connection between the domain
n P(Traces) of the trace semantics and the do-
main P(ConcreteStates) . Applying the result from Exercise 12.58 then gives
us a Galois connection between P(Traces) and States n , which can be illustrated
like this:
γt ◦ γc
γt γc
αt αc
αc ◦ αt
n
P(Traces) P(ConcreteStates) States n
Exercise 12.59: Prove that the reachable states collecting semantics is sound
with respect to the trace semantics. (Even though the collecting semantics is
not a computable analysis, we can still apply the notion of soundness and the
proof techniques from Section 12.3.)
where (
0 δ[X 7→ v 0 ] if v 0 is an assignment X = E
δ =
δ otherwise
In each element of the traces, the first two components are exactly as before.
The reaching definitions map in the third component is initially empty and is
updated whenever an assignment is encountered.
We also define an abstraction function αRD : P(Traces RD ) → (P(Nodes))n
that extracts the sets of semantically reaching definitions from a set of instru-
mented traces for a program with n program points:
αRD (T ) = (R1 , . . . , Rn )
where
where [[P ]]RD denotes the reaching definitions analysis result for P .
Exercise 12.60: Describe h[P ]iRD and [[P ]]RD for P being the example program
from Section 5.7.
Exercise 12.61: With this instrumented trace semantics, use the approach from
Section 12.3 to prove that the reaching definitions analysis from Section 5.7 is
sound.
Reasoning about other analyses generally requires other forms of instru-
mented semantics that capture the semantic properties of interest.
194 12 ABSTRACT INTERPRETATION
Exercise 12.62: Use the approach from Section 12.3 and an appropriate in-
strumented trace semantics to prove that the available expressions analysis
from Section 5.5 is sound. (This is more tricky than Exercise 12.61, because
available expressions analysis is a “must” analysis!)
Exercise 12.63: Use the approach from Section 12.3 and an appropriate in-
strumented trace semantics to prove that the live variables analysis from
Section 5.4 is sound. As part of this, you need to specify an appropriate
collecting semantics that formally captures what it means for a variable to
be live (see the informal definition in Section 5.4). (This is more tricky than
Exercise 12.61, because live variables analysis is a “backward” analysis!)
Exercise 12.64: Investigate for some of the abstractions used in analyses pre-
sented in the preceding chapters (for example, live variables analysis or reach-
ing definitions analysis) whether or not they are optimal and/or complete.
Exercise 12.65: Show that there exists a Galois connection between the instru-
mented trace semantics for reaching definitions αRD (h[P ]iRD ) and the ordinary
trace semantics α(h[P ]iRD ).
Index of Notation
[[·]] Constraint variable (for analysis) 22, 51, 136, 148, 153
{[·]} Constraint variable (for collecting semantics) 164
h[·]i Constraint variable (for trace semantics) 190
A≤k Bounded tuples 104
[] Empty function 166
[a1 7→ x1 , . . .] Function definition 42
F
F Greatest lower bound 39
Least upper bound 39
⊥ Bottom element 37, 40
◦ Function composition 45, 121
↓ Projection (various meanings) 68, 72, 159
λx.e Lambda abstraction 93, 135
Z Integers 57
lfp Least fixed point 47, 81, 169
µα.τ Recursive type 21
u Meet operator (see alsoF ) 39
F
t Join operator (see also ) 39
v Partial order relation 38
w Reverse partial order relation 39
τ1 [τ2 /α] Term substitution 21
(τ1 , τ2 ) → τ3 Function type 20
> Top element 37, 40
ocp Abstract operator 53
f [a 7→ x] Function update 44
Pointer type 20
195
Bibliography
[ALSU06] Alfred V. Aho, Monica S. Lam, Ravi Sethi, and Jeffrey D. Ullman.
Compilers: Principles, Techniques, and Tools. Addison Wesley, 2006.
[And94] Lars Ole Andersen. Program analysis and specialization for the C pro-
gramming language. PhD thesis, University of Copenhagen, 1994.
[BR02] Thomas Ball and Sriram K. Rajamani. The SLAM project: debugging
system software via static analysis. In Conference Record of POPL 2002:
The 29th SIGPLAN-SIGACT Symposium on Principles of Programming
Languages, Portland, OR, USA, January 16-18, 2002, pages 1–3. ACM,
2002.
197
198 BIBLIOGRAPHY
[PB09] Fernando Magno Quintão Pereira and Daniel Berlin. Wave propa-
gation and deep propagation for pointer analysis. In Proceedings of
the CGO 2009, The Seventh International Symposium on Code Generation
and Optimization, Seattle, Washington, USA, March 22-25, 2009, pages
126–135. IEEE Computer Society, 2009.
[PKH04] David J. Pearce, Paul H. J. Kelly, and Chris Hankin. Online cycle de-
tection and difference propagation: Applications to pointer analysis.
Softw. Qual. J., 12(4):311–337, 2004.
[PKH07] David J. Pearce, Paul H. J. Kelly, and Chris Hankin. Efficient field-
sensitive pointer analysis of C. ACM Trans. Program. Lang. Syst.,
30(1):4, 2007.
[RHS95] Thomas W. Reps, Susan Horwitz, and Shmuel Sagiv. Precise interpro-
cedural dataflow analysis via graph reachability. In Conference Record
of POPL’95: 22nd ACM SIGPLAN-SIGACT Symposium on Principles of
Programming Languages, San Francisco, California, USA, January 23-25,
1995, pages 49–61. ACM Press, 1995.
[Ric53] Henry Gordon Rice. Classes of recursively enumerable sets and their
decision problems. Transactions of the American Mathematical Society,
74(2):358–366, 1953.
[SCD+ 13] Manu Sridharan, Satish Chandra, Julian Dolby, Stephen J. Fink, and
Eran Yahav. Alias analysis for object-oriented programs. In Aliasing in
Object-Oriented Programming. Types, Analysis and Verification, volume
7850 of Lecture Notes in Computer Science, pages 196–232. Springer,
2013.
[SP81] Micha Sharir and Amir Pnueli. Two approaches to interprocedural data
flow analysis, chapter 7, pages 189–234. Prentice-Hall, 1981.
[SRH96] Shmuel Sagiv, Thomas W. Reps, and Susan Horwitz. Precise interpro-
cedural dataflow analysis with applications to constant propagation.
Theor. Comput. Sci., 167(1&2):131–170, 1996.
[Wan87] Mitchell Wand. A simple algorithm and proof for type inference.
Fundamenta Informaticae, 10(2):115–121, 1987.