| 49 | 
  | 
#define MATH_SQUAREMATRIX_HPP  | 
| 50 | 
  | 
 | 
| 51 | 
  | 
#include "math/RectMatrix.hpp" | 
| 52 | 
+ | 
#include "utils/NumericConstant.hpp" | 
| 53 | 
  | 
 | 
| 54 | 
  | 
namespace oopse { | 
| 55 | 
  | 
 | 
| 139 | 
  | 
    bool isSymmetric() const { | 
| 140 | 
  | 
      for (unsigned int i = 0; i < Dim - 1; i++) | 
| 141 | 
  | 
        for (unsigned int j = i; j < Dim; j++) | 
| 142 | 
< | 
          if (fabs(this->data_[i][j] - this->data_[j][i]) > oopse::epsilon)  | 
| 142 | 
> | 
          if (fabs(this->data_[i][j] - this->data_[j][i]) > epsilon)  | 
| 143 | 
  | 
            return false; | 
| 144 | 
  | 
                         | 
| 145 | 
  | 
      return true; | 
| 158 | 
  | 
    bool isDiagonal() const { | 
| 159 | 
  | 
      for (unsigned int i = 0; i < Dim ; i++) | 
| 160 | 
  | 
        for (unsigned int j = 0; j < Dim; j++) | 
| 161 | 
< | 
          if (i !=j && fabs(this->data_[i][j]) > oopse::epsilon)  | 
| 161 | 
> | 
          if (i !=j && fabs(this->data_[i][j]) > epsilon)  | 
| 162 | 
  | 
            return false; | 
| 163 | 
  | 
                         | 
| 164 | 
  | 
      return true; | 
| 170 | 
  | 
        return false; | 
| 171 | 
  | 
                 | 
| 172 | 
  | 
      for (unsigned int i = 0; i < Dim ; i++) | 
| 173 | 
< | 
        if (fabs(this->data_[i][i] - 1) > oopse::epsilon) | 
| 173 | 
> | 
        if (fabs(this->data_[i][i] - 1) > epsilon) | 
| 174 | 
  | 
          return false; | 
| 175 | 
  | 
                     | 
| 176 | 
  | 
      return true; |