Question

I’m writing a templated mathematical matrix class to get used to some new c++11 features, the basic declaration is as follows:

template <typename Type, int kNumRows, int kNumCols>
class Matrix { ... };

The class has a member function to return one of its minors (which is later used to calculate the determinant of NxN matrices).

Matrix<Type, kNumRows - 1, kNumCols - 1> minor(const int row, const int col) {
  static_assert(kNumRows > 2, "");
  static_assert(kNumCols > 2, "");

  ...
}

I then created a non-member function to calculate the determinant of any square matrix:

template <typename Type, int kSize>
Type determinant(const Matrix<Type, kSize, kSize>& matrix) {
  switch (kSize) {
  case 2:
    return 0; // For now unimportant
  case 3:
    // Recursively call the determinant function on a minor matrix
    return determinant(matrix.minor(0, 0));
  }
  ...
}

In main() I create a 3×3 matrix and call determinant on it. This will not compile. The compiler effectively moves to case 3, creating a minor matrix and calling determinant on it. It then steps into case 3 again, resulting in a static_assert by trying to create a 1×1 minor.

The question is simple: am I missing something here? Is calling a templated function like this recursively simply not allowed? Is this a compiler fault (I doubt it)?

For completeness’ sake: I’m using Clang++.

Answer 1

Templates determine what to do at compile time, but a switch statement determines what to do at run time. The compiler generates code, or at least verifies the validity, for all switch cases, even if the correct case is “obvious” at compile time.

Instead of using switch, try overloading determinant:

template <typename Type>
Type determinant(const Matrix<Type, 1, 1>& matrix) {
    return matrix(0,0);
}

template <typename Type>
Type determinant(const Matrix<Type, 2, 2>& matrix) {
    return 0; // (incorrect math)
}

template <typename Type, int kSize>
Type determinant(const Matrix<Type, kSize, kSize>& matrix) {
    return determinant(matrix.minor(0,0)); // (incorrect math)
}