def test_tlog(self):
th = 2
r = 10 ** 4
d = np.log10(2 ** 18)
result = trans.tlog(_xall, th, r, d)
assert_almost_equal(result[_xall < th], np.log10(th) * 1.0 * r / d, decimal=3)
self.assertTrue(np.all(result[_xall > th]))
assert_almost_equal(r, result.max())
def test_tlog(self):
th = 2
r = 10**4
d = np.log10(2**18)
result = trans.tlog(_xall, th, r, d)
assert_almost_equal(result[_xall < th],
np.log10(th) * 1. * r / d,
decimal=3)
self.assertTrue(np.all(result[_xall > th]))
assert_almost_equal(r, result.max())
def test_hlog(self):
hlpos = trans.hlog(_xpos)
hlneg = trans.hlog(_xneg)
assert_almost_equal((hlpos[-1] - _ymax) / _ymax, 0, decimal=2)
assert_almost_equal(hlpos, -hlneg[::-1]) # check symmetry
# test that values get larger as b decreases
hlpos10 = trans.hlog(_xpos, b=10)
self.assertTrue(np.all(hlpos10 >= hlpos))
# check that converges to tlog for large values
tlpos = trans.tlog(_xpos)
i = np.where(_xpos > 1e4)[0]
tlpos_large = tlpos[i]
hlpos_large = hlpos10[i]
d = (hlpos_large - tlpos_large) / hlpos_large
assert_almost_equal(d, np.zeros(len(d)), decimal=2)
# test spline option
result1 = trans.hlog(_xall, use_spln=True)
result2 = trans.hlog(_xall, use_spln=False)
d = (result1 - result2) / result1
assert_almost_equal(d, np.zeros(len(d)), decimal=2)
def test_hlog(self):
hlpos = trans.hlog(_xpos)
hlneg = trans.hlog(_xneg)
assert_almost_equal((hlpos[-1] - _ymax) / _ymax, 0, decimal=2)
assert_almost_equal(hlpos, -hlneg[::-1]) # check symmetry
# test that values get larger as b decreases
hlpos10 = trans.hlog(_xpos, b=10)
self.assertTrue(np.all(hlpos10 >= hlpos))
# check that converges to tlog for large values
tlpos = trans.tlog(_xpos)
i = np.where(_xpos > 1e4)[0]
tlpos_large = tlpos[i]
hlpos_large = hlpos10[i]
d = (hlpos_large - tlpos_large) / hlpos_large
assert_almost_equal(d, np.zeros(len(d)), decimal=2)
# test spline option
transformation = Transformation(transform='hlog', direction='forward')
result1 = transformation(_xall, use_spln=True)
result2 = transformation(_xall, use_spln=False)
d = (result1 - result2) / result1
assert_almost_equal(d, np.zeros(len(d)), decimal=2)
def test_tlog_inv(self):
th = 2
expected = _xall.copy()
expected[_xall <= th] = th
result = trans.tlog_inv(trans.tlog(_xall, th=th), th=th)
assert_almost_equal(result, expected)
def test_tlog_inv(self):
th = 2
expected = _xall.copy()
expected[_xall <= th] = th
result = trans.tlog_inv(trans.tlog(_xall, th=th), th=th)
assert_almost_equal(result, expected)