I guess it boils down to the question of whether the default (atol=0, rtol=1e-7 according to numpy docs) is enough, or do we need better accuracy here. Or, more generally, is there a goal of hunting for corner cases/precision loss/etc in tests.