#ifndef MSP_LINAL_SQUAREMATRIX_H_
#define MSP_LINAL_SQUAREMATRIX_H_
+#include <cmath>
#include <stdexcept>
#include "matrix.h"
template<typename T, unsigned S>
SquareMatrix<T, S> &SquareMatrix<T, S>::invert()
{
+ using std::abs;
+
SquareMatrix<T, S> r = identity();
for(unsigned i=0; i<S; ++i)
{