Reconciling Abstraction with High Performance

Assuming no prior knowledge of MetaOCaml and only a basic familiarity with functional programming, this book explains and illustrates how to implement a simple domain-specific language for linear algebra, with layers of optimizations for sparsity and memory layout of matrices and vectors, and their algebraic properties.