Using types to keep yourself honest

Matthijs Hollemans
by Matthijs Hollemans
25 March 2016

Table of contents

This post explains how you can take advantage of Swift’s type system to make your programs more expressive and robust.

For the last week or so I’ve been playing with machine learning algorithms in Swift. You can often implement these algorithms really succinctly and efficiently by using matrices.

In case you forgot your linear algebra, just think of a matrix as a table of numbers.

When we say, “M is a 4-by-3 matrix”, we mean that M is a table of 4 rows and 3 columns. For the purposes of this article, that’s all you really need to know about matrices.

This is an example of a 4×3 matrix:

A 4x3 matrix

So I wrote some code and created a struct Matrix:

struct Matrix {
  let rows: Int
  let columns: Int
  ...
}

One thing you often need to do with matrices is multiply them, and I created a function for that:

func multiply(m1: Matrix, _ m2: Matrix) -> Matrix {
  // bunch of math...
}

This may all seem very straightforward but there’s something here that bothers me.

Even though m1 and m2 are both Matrix objects, they may actually have different numbers of rows and columns. And that could be a problem.

For example, with matrix multiplication the sizes of the two matrices have to match up in a particular way:

Matrix multiplication

The number of columns in the first matrix has to be the same as the number of rows in the second matrix. If the size of the first matrix is U × V, then the size of the second matrix has to be V × W. That’s how the math works.

The result is a new matrix of size U × W. If the sizes of the matrices don’t match up in this particular way, we can’t multiply them.

For example, the following will work fine:

let A = Matrix(rows: 4, columns: 3)
let B = Matrix(rows: 3, columns: 2)
let C = multiply(A, B)                // gives a 4×2 matrix

Note: In math, matrices are often denoted with capitals, and I’m following that convention here for the variable names.

Multiplying A with B is allowed because A.columns == B.rows. On the other hand, the following is not a valid operation:

let D = multiply(B, A)

The number of columns in matrix B doesn’t match up with the number of rows in matrix A. Here, we have B.columns != A.rows. Mathematically, doing B times A makes no sense.

Currently, the only way to catch these sorts of mistakes is by tripping an assertion during runtime:

func multiply(m1: Matrix, _ m2: Matrix) -> Matrix {
  // do the matrices have the correct sizes?
  precondition(m1.columns == m2.rows)
  
  // bunch of math...
}

Certainly doable, but I don’t like it. The whole point of Swift’s static typing is that the compiler can catch as many programming errors during compile time as possible. It would be nice if we could make the compiler catch this kind of error too.

It turns out we can! In this article I’ll explore how to use Swift’s type system to make such mistakes impossible.

Not a great solution

The naive way to approach this problem is to create different structs for matrices of different sizes:

let A = Matrix_4x3()
let B = Matrix_3x2()

But then you also need a multiply() method that takes these specific types as parameters:

func multiply(m1: Matrix_4x3, _ m2: Matrix_3x2) -> Matrix_4x2

That seems a bit silly and leads to a lot of duplicate code.

What’s worse, you may not actually know the sizes of your matrices at compile time. In a machine learning problem you often need to load a dataset from a file but you won’t know in advance how many rows it has.

So this isn’t a workable solution. However, the idea of declaring different types for matrices of different sizes is promising…

Generics to the rescue

We want to incorporate the dimensions of the matrix into the type of Matrix somehow, without requiring that any of the matrix code, such as multiply(), knows anything about specific sizes.

Let’s define Matrix as follows:

struct Matrix<R,C> {
  let rows: Int
  let columns: Int
  ...
}

It now has two generic parameters, R and C, where R stands for the number of rows and C for the number of columns.

We can then define multiply() like this:

func multiply<U,V,W>(m1: Matrix<U,V>, _ m2: Matrix<V,W>) -> Matrix<U,W> {
  // bunch of math...
  return Matrix(rows: m1.rows, columns: m2.columns)
}

Notice how this captures the rule for matrix multiplication: a matrix of size U × V times a matrix of size V × W gives a new matrix of size U × W.

Here’s an example of how to use this new Matrix:

struct NumExamples {}
struct NumFeatures {}
struct OneDimensional {}

let A = Matrix<NumExamples, NumFeatures>(rows: 20, columns: 10)
let B = Matrix<NumFeatures, OneDimensional>(rows: 10, columns: 1)

We’ve created three new types — named NumExamples, NumFeatures, and OneDimensional — to represent the possible dimensions of our matrices. Note how I’ve given these types descriptive names, so it’s easier to tell their purpose.

The names NumExamples and NumFeatures come from machine learning, since that’s what I’m going to use these matrices for. NumExamples is the number of objects in your dataset and NumFeatures is the number of attributes that each example has. (Of course, if you were going to use the matrices for something else, you’d use different names.)

OneDimensional tells you that the matrix B only has one column. In linear algebra, we’d call that a column vector instead of a matrix. To make that distinction clearer in our code, it would be cool if we could write something like this:

typealias ColumnVector<Rows> = Matrix<Rows, OneDimensional>
typealias RowVector<Columns> = Matrix<OneDimensional, Columns>

Those typealiases would make ColumnVector and RowVector special cases of Matrix. But unfortunately this syntax is not supported in Swift 2.2. It might be coming in Swift 3.0.

Anyway, back to the example. When you now write,

let C = multiply(A, B)

it gives a new 20×1 matrix, as expected. Unlike before, however, an invalid attempt at multiplication gives a compiler error:

let D = multiply(B, A)

// error: cannot convert value of type 'Matrix<NumFeatures, OneDimensional>'
// to expected argument type 'Matrix<_, _>'

The error message is a bit vague, but what’s great is that we’ve used Swift’s type system to catch this kind of error. Instead of crashing the app at runtime, it is now impossible to multiply two matrices that do not have the correct dimensions.

Or is it? Well, you can still lie to the compiler:

let A = Matrix<NumExamples, NumFeatures>(rows: 20, columns: 10)
let B = Matrix<NumFeatures, OneDimensional>(rows: 500, columns: 1)

By changing B’s number of rows to 500, we still end up in the same situation as before. Now multiply(A, B) is no longer valid.

Just having these extra types is not enough… We need to make sure that the type NumFeatures somehow always refers to the same number, regardless of where it is used.

Protocols to the rescue

We could do something like this:

struct NumExamples { let size = 20 }
struct NumFeatures { let size = 10 }

But this fixes the sizes of these dimensions at compile time. Remember, we want to be able to set the matrix sizes at runtime, for example by reading a dataset from a file — and we may not know beforehand how much data is in that file. Hardcoding the matrix sizes is a no-go.

Instead, let’s define a new protocol:

protocol Dimension {
  static var size: Int { get set }
}

Then Matrix becomes:

struct Matrix<R: Dimension, C: Dimension> {
  let rows: Int
  let columns: Int

  init() {
    self.rows = R.size
    self.columns = C.size
  }
}

Notice there no longer is an init(rows:columns:) function. The size of the matrix comes directly from the types R and C.

The final step is to make our dimension types conform to the new protocol:

struct NumExamples: Dimension { static var size = 20 }
struct NumFeatures: Dimension { static var size = 10 }
struct OneDimensional: Dimension { static var size = 1 }

Now we can write multiply() as follows:

func multiply<U: Dimension, V: Dimension, W: Dimension>
             (m1: Matrix<U,V>, _ m2: Matrix<V,W>) -> Matrix<U,W> {
  // bunch of math...
  return Matrix<U,W>()
}

It is now impossible that matrices m1 and m2 don’t match up. The compiler simply won’t accept it.

let A = Matrix<NumExamples, NumFeatures>()
let B = Matrix<NumFeatures, OneDimensional>()

let C = multiply(A, B)   // yay!

let D = multiply(B, A)   // compiler error

There is no longer a way to make inadvertent mistakes. Of course, you can still cheat the system by doing something like this:

let A = Matrix<NumExamples, NumFeatures>()
NumFeatures.size = 500
let B = Matrix<NumFeatures, OneDimensional>()

Even Swift’s type system can’t stop you if you’re intent on being evil! (It’s probably smart to keep that precondition() inside multiply().)

By the way, you actually need that ability to change NumFeatures.size. But you should use it carefully. Just as we don’t know what this particular size will be until we run the program, there’s no reason why it should remain the same the whole time. For example, you might need to process multiple datasets of different sizes using the same routines.

Of course, you can do more with matrices than just multiply them. Here’s another example of where these dimension types come in useful:

func processData<M: Dimension, N: Dimension>
                (X: Matrix<M, N>, _ y: Matrix<M, OneDimensional>) 
                -> Matrix<OneDimensional, N> {
  // do impressive stuff...
}

let X = Matrix<NumExamples, NumFeatures>()
let y = Matrix<NumExamples, OneDimensional>()
processData(X, y)

This function takes a matrix X and a column vector y and does some work on them; for example, it might train a learning system. The constraint here is that X and y must have the same number of rows. Thanks to our dimension types, the compiler can enforce that constraint.

Conclusion

We used types to better express to the compiler what our program is doing. And that helps the compiler catch mistakes.

Written by Matthijs Hollemans.
First published on Friday, 25 March 2016.
If you liked this post, say hi on Twitter @mhollemans or LinkedIn.
Find the source code on my GitHub.

Code Your Own Synth Plug-Ins With C++ and JUCENew e-book: Code Your Own Synth Plug-Ins With C++ and JUCE
Interested in how computers make sound? Learn the fundamentals of audio programming by building a fully-featured software synthesizer plug-in, with every step explained in detail. Not too much math, lots of in-depth information! Get the book at Leanpub.com