template <size_t N>
class Rivet::Matrix;

using Rivet::Matrix< N >::EMatrix =  Eigen::Matrix<double,N,N>;

static inline Matrix< N > mkZero()

static inline Matrix< N > mkDiag(
    Vector< N > diag
)

static inline Matrix< N > mkIdentity()

inline Matrix()

inline Matrix(
    const Matrix< N > & other
)

inline Matrix & set(
    const size_t i,
    const size_t j,
    const double value
)

inline double get(
    const size_t i,
    const size_t j
) const

inline Vector< N > getRow(
    const size_t row
) const

inline Matrix< N > & setRow(
    const size_t row,
    const Vector< N > & r
)

inline Vector< N > getColumn(
    const size_t col
) const

inline Matrix< N > & setColumn(
    const size_t col,
    const Vector< N > & c
)

inline Matrix< N > transpose() const

inline Matrix< N > inverse() const

inline double det() const

inline double trace() const

inline Matrix< N > operator-() const

inline constexpr size_t size() const

inline bool isZero(
    double tolerance =1E-5
) const

inline bool isEqual(
    Matrix< N > other
) const

inline bool isSymm() const

inline bool isDiag() const

inline bool operator==(
    const Matrix< N > & a
) const

inline bool operator!=(
    const Matrix< N > & a
) const

inline Matrix< N > & operator*=(
    const Matrix< N > & m
)

inline Matrix< N > & operator*=(
    const double a
)

inline Matrix< N > & operator/=(
    const double a
)

inline Matrix< N > & operator+=(
    const Matrix< N > & m
)

inline Matrix< N > & operator-=(
    const Matrix< N > & m
)

friend Matrix< M > add(
    const Matrix< M > & ,

    const Matrix< M > & 
);

friend Matrix< M > multiply(
    const double ,

    const Matrix< M > & 
);

friend Matrix< M > multiply(
    const Matrix< M > & ,

    const Matrix< M > & 
);

friend Vector< M > multiply(
    const Matrix< M > & ,

    const Vector< M > & 
);

friend Matrix< M > divide(
    const Matrix< M > & ,

    const double 
);