def test_hlog_inv(self): expected = _xall result = trans.hlog_inv(trans.hlog(_xall)) d = (result - expected) / expected assert_almost_equal(d, np.zeros(len(d)), decimal=2) result = trans.hlog_inv(trans.hlog(_xall, b=10), b=10) d = (result - expected) / expected 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 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)