invalid_matrix in product spaces now allows very small values

This commit is contained in:
Zeno Rogue 2019-12-06 13:06:32 +01:00
parent e9a0d7f817
commit c8db8a328d
1 changed files with 1 additions and 1 deletions

View File

@ -1011,7 +1011,7 @@ EX bool invalid_matrix(const transmatrix T) {
if(std::isnan(T[i][j]) || T[i][j] > 1e8 || T[i][j] < -1e8 || std::isinf(T[i][j]))
return true;
if(prod) {
for(int i=0; i<GDIM; i++) for(int j=0; j<GDIM; j++) if(abs(T[i][j]) > 1e-6) return false;
for(int i=0; i<GDIM; i++) for(int j=0; j<GDIM; j++) if(abs(T[i][j]) > 1e-60) return false;
}
else
for(int i=0; i<GDIM; i++) for(int j=0; j<GDIM; j++) if(T[i][j] > .5 || T[i][j] < -.5) return false;