def test_logaddexp_values(self):
x = [1, 2, 3, 4, 5]
y = [5, 4, 3, 2, 1]
z = [6, 6, 6, 6, 6]
for ds, dec in zip(['float32', 'float64'], [6, 15]):
xf = blaze.log(blaze.array(x, dshape=ds))
yf = blaze.log(blaze.array(y, dshape=ds))
zf = blaze.log(blaze.array(z, dshape=ds))
result = blaze.logaddexp(xf, yf)
assert_almost_equal(result, zf, decimal=dec)
def test_log_values(self):
x = [1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024]
y = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
for ds in ['float32', 'float64']:
log2_ = 0.69314718055994530943
xf = blaze.array(x, dshape=ds)
yf = blaze.array(y, dshape=ds)*log2_
result = blaze.log(xf)
assert_almost_equal(result, yf)
def test_log1p(self):
assert_almost_equal(blaze.log1p(0.2), blaze.log(1.2))
assert_almost_equal(blaze.log1p(1e-6), blaze.log(1+1e-6))
