#include <msp/linal/matrix.h>
#include <msp/linal/squarematrix.h>
#include <msp/test/test.h>
#include <msp/linal/matrix.h>
#include <msp/linal/squarematrix.h>
#include <msp/test/test.h>
template<typename T, unsigned N>
bool is_identity(const LinAl::Matrix<T, N, N> &m)
{
template<typename T, unsigned N>
bool is_identity(const LinAl::Matrix<T, N, N> &m)
{