casa6
Public
Commits
| 252 252 | ref = os.path.join(refdir, ref) |
| 253 253 | self.assertTrue(mytb.open(ref), "Table open failed for " + ref) |
| 254 254 | [ |
| 255 255 | etimes, ewt, ewtsp, eflag, efrow, edata, esigma, esisp |
| 256 256 | ] = _get_table_cols(mytb) |
| 257 257 | mytb.done() |
| 258 258 | self.assertTrue(np.allclose(gwt, ewt), 'WEIGHT comparison failed') |
| 259 259 | if type(gwtsp) != type(None) and type(ewtsp) != type(None): |
| 260 260 | self.assertTrue( |
| 261 261 | np.allclose(gwtsp, ewtsp), 'WEIGHT_SPECTRUM comparison failed' |
| 262 + | f"max diff {np.max(np.abs(gwtsp - ewtsp))}" |
| 262 263 | ) |
| 263 264 | self.assertTrue((gflag == eflag).all(), 'FLAG comparison failed') |
| 264 265 | self.assertTrue((gfrow == efrow).all(), 'FLAG_ROW comparison failed') |
| 265 266 | # all flags must be True where wtsp = 0 |
| 266 267 | self.assertTrue(np.extract(gwtsp == 0, gflag).all()) |
| 267 268 | self.assertTrue(np.allclose(gsigma, esigma), 'SIGMA comparison failed') |
| 268 269 | if type(gsisp) != type(None) and type(esisp) != type(None): |
| 269 270 | self.assertTrue(np.allclose( |
| 270 271 | gsisp, esisp), 'SIGMA_SPECTRUM comparison failed' |
| 271 272 | ) |