eqmatrix now actually takes its 'eps' argument into account

This commit is contained in:
Zeno Rogue 2019-05-04 18:29:52 +02:00
parent 156f56feb0
commit 94ccddb873
2 changed files with 2 additions and 2 deletions

View File

@ -3649,7 +3649,7 @@ hyperpoint xspinpush0(ld alpha, ld x);
transmatrix cspin(int a, int b, ld alpha); transmatrix cspin(int a, int b, ld alpha);
transmatrix cpush(int cid, ld alpha); transmatrix cpush(int cid, ld alpha);
bool eqmatrix(transmatrix A, transmatrix B, ld eps = 1e-6); bool eqmatrix(transmatrix A, transmatrix B, ld eps = 1e-2);
void set_column(transmatrix& T, int i, const hyperpoint& H); void set_column(transmatrix& T, int i, const hyperpoint& H);
#define DF_INIT 0 // always display these #define DF_INIT 0 // always display these

View File

@ -327,7 +327,7 @@ inline hyperpoint xspinpush0(ld alpha, ld x) {
bool eqmatrix(transmatrix A, transmatrix B, ld eps) { bool eqmatrix(transmatrix A, transmatrix B, ld eps) {
for(int i=0; i<MDIM; i++) for(int i=0; i<MDIM; i++)
for(int j=0; j<MDIM; j++) for(int j=0; j<MDIM; j++)
if(std::abs(A[i][j] - B[i][j]) > 1e-2) if(std::abs(A[i][j] - B[i][j]) > eps)
return false; return false;
return true; return true;
} }