Bird–Meertens formalism

The Bird–Meertens formalism (BMF) is a calculus for deriving programs from program specifications (in a functional programming setting) by a process of equational reasoning. It was devised by Richard Bird and Lambert Meertens as part of their work within IFIP Working Group 2.1.

It is sometimes referred to in publications as BMF, as a nod to Backus–Naur form. Facetiously it is also referred to as Squiggol, as a nod to ALGOL, which was also in the remit of WG 2.1, and because of the "squiggly" symbols it uses. A less-used variant name, but actually the first one suggested, is SQUIGOL. Martin and Nipkow provided automated support for Squiggol development proofs, using the Larch Prover.^[1]

Basic examples and notations

Map is a well-known second-order function that applies a given function to every element of a list; in BMF, it is written ${\displaystyle *}$:

${\displaystyle f*[e_{1},\dots ,e_{n}]=[f\ e_{1},\dots ,f\ e_{n}].}$

Likewise, reduce is a function that collapses a list into a single value by repeated application of a binary operator. It is written / in BMF. Taking ${\displaystyle \oplus }$ as a suitable binary operator with neutral element e, we have

${\displaystyle \oplus /[e_{1},\dots ,e_{n}]=e\oplus e_{1}\oplus \dots \oplus e_{n}.}$

Using those two operators and the primitives ${\displaystyle +}$ (as the usual addition), and ${\displaystyle +\!\!\!+}$ (for list concatenation), we can easily express the sum of all elements of a list, and the flatten function, as ${\displaystyle {\rm {sum}}=+/}$ and ${\displaystyle {\rm {flatten}}=+\!\!\!+/}$, in point-free style. We have:

${\displaystyle {\rm {sum}}\ [e_{1},\dots ,e_{n}]=+/[e_{1},\dots ,e_{n}]=0+e_{1}+\dots +e_{n}=\sum _{k}e_{k}.}$

${\displaystyle {\rm {flatten}}\ [l_{1},\dots ,l_{n}]=+\!\!\!+/[l_{1},\dots ,l_{n}]=[\,]+\!\!\!+\;l_{1}+\!\!\!+\dots +\!\!\!+\;l_{n}={\text{ the concatenation of all lists }}l_{k}.}$

Example instances of used laws
Map promotion law[3]
Fold promotion law[3]
Generalized Horner's rule[4]
Scan lemma[5]
Fold-scan fusion law[5]

Derivation of Kadane's algorithm^[2]

Similarly, writing ${\displaystyle \cdot }$ for functional composition and ${\displaystyle \land }$ for conjunction, it is easy to write a function testing that all elements of a list satisfy a predicate p, simply as ${\displaystyle {\rm {all}}\ p=(\land /)\cdot (p*)}$:

${\displaystyle {\begin{aligned}{\rm {all}}\ p\ [e_{1},\dots ,e_{n}]&=(\land /)\cdot (p*)\ [e_{1},\dots ,e_{n}]\\&=\land /(p*[e_{1},\dots ,e_{n}])\\&=\land /[p\ e_{1},\dots ,p\ e_{n}]\\&=p\ e_{1}\land \dots \land p\ e_{n}\\&=\forall k\ .\ p\ e_{k}.\end{aligned}}}$

Bird (1989) transforms inefficient easy-to-understand expressions ("specifications") into efficient involved expressions ("programs") by algebraic manipulation. For example, the specification "${\displaystyle \mathrm {max} \cdot \mathrm {map} \;\mathrm {sum} \cdot \mathrm {segs} }$" is an almost literal translation of the maximum segment sum problem,^[6] but running that functional program on a list of size ${\displaystyle n}$ will take time ${\displaystyle {\mathcal {O}}(n^{3})}$ in general. From this, Bird computes an equivalent functional program that runs in time ${\displaystyle {\mathcal {O}}(n)}$, and is in fact a functional version of Kadane's algorithm.

The derivation is shown in the picture, with computational complexities^[7] given in blue, and law applications indicated in red. Example instances of the laws can be opened by clicking on [show]; they use lists of integer numbers, addition, minus, and multiplication. The notation in Bird's paper differs from that used above: ${\displaystyle \mathrm {map} }$, ${\displaystyle \mathrm {concat} }$, and ${\displaystyle \mathrm {foldl} }$ correspond to ${\displaystyle *}$, ${\displaystyle \mathrm {flatten} }$, and a generalized version of ${\displaystyle /}$ above, respectively, while ${\displaystyle \mathrm {inits} }$ and ${\displaystyle \mathrm {tails} }$ compute a list of all prefixes and suffixes of its arguments, respectively. As above, function composition is denoted by "${\displaystyle \cdot }$", which has lowest binding precedence. In the example instances, lists are colored by nesting depth; in some cases, new operations are defined ad hoc (grey boxes).

The homomorphism lemma and its applications to parallel implementations

A function h on lists is called a list homomorphism if there exists an associative binary operator ${\displaystyle \oplus }$ and neutral element ${\displaystyle e}$ such that the following holds:

${\displaystyle {\begin{aligned}&h\ [\,]&&=\ e\\&h\ (l+\!\!\!+\;m)&&=\ h\ l\oplus h\ m.\end{aligned}}}$

The homomorphism lemma states that h is a homomorphism if and only if there exists an operator ${\displaystyle \oplus }$ and a function f such that ${\displaystyle h=(\oplus /)\cdot (f*)}$.

A point of great interest for this lemma is its application to the derivation of highly parallel implementations of computations. Indeed, it is trivial to see that ${\displaystyle f*}$ has a highly parallel implementation, and so does ${\displaystyle \oplus /}$ — most obviously as a binary tree. Thus for any list homomorphism h, there exists a parallel implementation. That implementation cuts the list into chunks, which are assigned to different computers; each computes the result on its own chunk. It is those results that transit on the network and are finally combined into one. In any application where the list is enormous and the result is a very simple type – say an integer – the benefits of parallelisation are considerable. This is the basis of the map-reduce approach.

See also

• Catamorphism
• Anamorphism
• Paramorphism
• Hylomorphism

References

1. Ursula Martin; Tobias Nipkow (Apr 1990). "Automating Squiggol". In Manfred Broy; Cliff B. Jones (eds.). Proc. IFIP WG 2.2/2.3 Working Conference on Programming Concepts and Methods. North-Holland. pp. 233–247.
2. Bird 1989, Sect.8, p.126r.
3. Bird 1989, Sect.2, p.123l.
4. Bird 1989, Sect.7, Lem.1, p.125l.
5. Bird 1989, Sect.5, p.124r.
6. Where ${\displaystyle \mathrm {max} }$, ${\displaystyle \mathrm {sum} }$, and ${\displaystyle \mathrm {segs} }$ returns the largest value, the sum, and the list of all segments (i.e. sublists) of a given list, respectively.
7. Each expression in a line denotes an executable functional program to compute the maximum segment sum.

• Meertens, Lambert (1986). "Algorithmics: Towards programming as a mathematical activity". In de Bakker, J.W.; Hazewinkel, M.; Lenstra, J.K. (eds.). Mathematics and Computer Science. CWI Monographs. Vol. 1. North-Holland. pp. 289–334.
• Meertens, Lambert; Bird, Richard (1987). "Two Exercises Found in a Book on Algorithmics" (PDF). North-Holland.
• Backhouse, Roland (1988). An Exploration of the Bird-Meertens Formalism (PDF) (Technical report).
• Bird, Richard S. (1989). "Algebraic Identities for Program Calculation" (PDF). The Computer Journal. 32 (2): 122–126. doi:10.1093/comjnl/32.2.122.
• Cole, Murray (1993). "Parallel Programming, List Homomorphisms and the Maximum Segment Sum Problem". Parallel Computing: Trends and Applications, PARCO 1993, Grenoble, France. pp. 489–492.
• Backhouse, Roland; Hoogendijk, Paul (1993). Elements of a Relational Theory of Datatypes (PDF). pp. 7–42. doi:10.1007/3-540-57499-9_15. ISBN 978-3-540-57499-6.
• Bunkenburg, Alexander (1994). O'Donnell, John T.; Hammond, Kevin (eds.). The Boom Hierarchy (PDF). Functional Programming, Glasgow 1993: Proceedings of the 1993 Glasgow Workshop on Functional Programming, Ayr, Scotland, 5–7 July 1993. London: Springer. pp. 1–8. doi:10.1007/978-1-4471-3236-3_1. ISBN 978-1-4471-3236-3.
• Bird, Richard; de Moor, Oege (1997). Algebra of Programming. International Series in Computing Science. Vol. 100. Prentice Hall. ISBN 0-13-507245-X.
• Gibbons, Jeremy (2020). Troy Astarte (ed.). The School of Squiggol: A History of the Bird-Meertens Formalism (PDF). Formal Methods (Workshop on History of Formal Methods). LNCS. Vol. 12233. Springer. doi:10.1007/978-3-030-54997-8_2.