A Simple and Practical Linear Algebra Library Interface with Static Size Checking

Akinori Abe
Eijiro Sumii

Linear algebra is a major field of numerical computation and is widely applied. Most linear algebra libraries (in most programming languages) do not statically guarantee consistency of the dimensions of vectors and matrices, causing runtime errors. While advanced type systems—specifically, dependent types on natural numbers—can ensure consistency among the sizes of collections such as lists and arrays, such type systems generally require non-trivial changes to existing languages and application programs, or tricky type-level programming.

We have developed a linear algebra library interface that verifies the consistency (with respect to dimensions) of matrix operations by means of generative phantom types, implemented via fairly standard ML types and module system. To evaluate its usability, we ported to it a practical machine learning library from a traditional linear algebra library. We found that most of the changes required for the porting could be made mechanically, and changes that needed human thought are minor.

In Oleg Kiselyov and Jacques Garrigue: Proceedings ML Family/OCaml Users and Developers workshops (ML/OCaml 2014), Gothenburg, Sweden, September 4-5, 2014, Electronic Proceedings in Theoretical Computer Science 198, pp. 1–21.
Published: 5th December 2015.

ArXived at: https://dx.doi.org/10.4204/EPTCS.198.1 bibtex PDF
References in reconstructed bibtex, XML and HTML format (approximated).
Comments and questions to: eptcs@eptcs.org
For website issues: webmaster@eptcs.org