Run-time algorithm specialization

In computer science, run-time algorithm specialization is a methodology for creating efficient algorithms for costly computation tasks of certain kinds. The methodology originates in the field of automated theorem proving and, more specifically, in the Vampire theorem prover project.

The idea is inspired by the use of partial evaluation in optimising program translation. Many core operations in theorem provers exhibit the following pattern. Suppose that we need to execute some algorithm [math]\displaystyle{ \mathit{alg}(A,B) }[/math] in a situation where a value of [math]\displaystyle{ A }[/math] is fixed for potentially many different values of [math]\displaystyle{ B }[/math]. In order to do this efficiently, we can try to find a specialization of [math]\displaystyle{ \mathit{alg} }[/math] for every fixed [math]\displaystyle{ A }[/math], i.e., such an algorithm [math]\displaystyle{ \mathit{alg}_A }[/math], that executing [math]\displaystyle{ \mathit{alg}_A(B) }[/math] is equivalent to executing [math]\displaystyle{ \mathit{alg}(A,B) }[/math].

The specialized algorithm may be more efficient than the generic one, since it can exploit some particular properties of the fixed value [math]\displaystyle{ A }[/math]. Typically, [math]\displaystyle{ \mathit{alg}_A(B) }[/math] can avoid some operations that [math]\displaystyle{ \mathit{alg}(A,B) }[/math] would have to perform, if they are known to be redundant for this particular parameter [math]\displaystyle{ A }[/math]. In particular, we can often identify some tests that are true or false for [math]\displaystyle{ A }[/math], unroll loops and recursion, etc.

Difference from partial evaluation

The key difference between run-time specialization and partial evaluation is that the values of [math]\displaystyle{ A }[/math] on which [math]\displaystyle{ \mathit{alg} }[/math] is specialised are not known statically, so the specialization takes place at run-time.

There is also an important technical difference. Partial evaluation is applied to algorithms explicitly represented as codes in some programming language. At run-time, we do not need any concrete representation of [math]\displaystyle{ \mathit{alg} }[/math]. We only have to imagine [math]\displaystyle{ \mathit{alg} }[/math] when we program the specialization procedure. All we need is a concrete representation of the specialized version [math]\displaystyle{ \mathit{alg}_A }[/math]. This also means that we cannot use any universal methods for specializing algorithms, which is usually the case with partial evaluation. Instead, we have to program a specialization procedure for every particular algorithm [math]\displaystyle{ \mathit{alg} }[/math]. An important advantage of doing so is that we can use some powerful ad hoc tricks exploiting peculiarities of [math]\displaystyle{ \mathit{alg} }[/math] and the representation of [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math], which are beyond the reach of any universal specialization methods.

Specialization with compilation

The specialized algorithm has to be represented in a form that can be interpreted. In many situations, usually when [math]\displaystyle{ \mathit{alg}_A(B) }[/math] is to be computed on many values [math]\displaystyle{ B }[/math] in a row, we can write [math]\displaystyle{ \mathit{alg}_A }[/math] as a code of a special abstract machine, and we often say that [math]\displaystyle{ A }[/math] is compiled. Then the code itself can be additionally optimized by answer-preserving transformations that rely only on the semantics of instructions of the abstract machine.

Instructions of the abstract machine can usually be represented as records. One field of such a record stores an integer tag that identifies the instruction type, other fields may be used for storing additional parameters of the instruction, for example a pointer to another instruction representing a label, if the semantics of the instruction requires a jump. All instructions of a code can be stored in an array, or list, or tree.

Interpretation is done by fetching instructions in some order, identifying their type and executing the actions associated with this type. In C or C++ we can use a switch statement to associate some actions with different instruction tags. Modern compilers usually compile a switch statement with integer labels from a narrow range rather efficiently by storing the address of the statement corresponding to a value [math]\displaystyle{ i }[/math] in the [math]\displaystyle{ i }[/math]-th cell of a special array. One can exploit this by taking values for instruction tags from a small interval of integers.

Data-and-algorithm specialization

There are situations when many instances of [math]\displaystyle{ A }[/math] are intended for long-term storage and the calls of [math]\displaystyle{ \mathit{alg}(A,B) }[/math] occur with different [math]\displaystyle{ B }[/math] in an unpredictable order. For example, we may have to check [math]\displaystyle{ \mathit{alg}(A_1,B_1) }[/math] first, then [math]\displaystyle{ \mathit{alg}(A_2,B_2) }[/math], then [math]\displaystyle{ \mathit{alg}(A_1,B_3) }[/math], and so on. In such circumstances, full-scale specialization with compilation may not be suitable due to excessive memory usage. However, we can sometimes find a compact specialized representation [math]\displaystyle{ A^{\prime} }[/math] for every [math]\displaystyle{ A }[/math], that can be stored with, or instead of, [math]\displaystyle{ A }[/math]. We also define a variant [math]\displaystyle{ \mathit{alg}^{\prime} }[/math] that works on this representation and any call to [math]\displaystyle{ \mathit{alg}(A,B) }[/math] is replaced by [math]\displaystyle{ \mathit{alg}^{\prime}(A^{\prime},B) }[/math], intended to do the same job faster.

See also

• Psyco, a specializing run-time compiler for Python
• multi-stage programming

References

• A. Voronkov, "The Anatomy of Vampire: Implementing Bottom-Up Procedures with Code Trees", Journal of Automated Reasoning, 15(2), 1995 (original idea)

Further reading

• A. Riazanov and A. Voronkov, "Efficient Checking of Term Ordering Constraints", Proc. IJCAR 2004, Lecture Notes in Artificial Intelligence 3097, 2004 (compact but self-contained illustration of the method)
• A. Riazanov and A. Voronkov, Efficient Instance Retrieval with Standard and Relational Path Indexing, Information and Computation, 199(1-2), 2005 (contains another illustration of the method)
• A. Riazanov, "Implementing an Efficient Theorem Prover", PhD thesis, The University of Manchester, 2003 (contains the most comprehensive description of the method and many examples)