def test_dtype():
assert theano_code(x, dtypes={x: 'float32'}).type.dtype == 'float32'
assert theano_code(x, dtypes={x: 'float64'}).type.dtype == 'float64'
assert theano_code(x + 1, dtypes={x: 'float32'}).type.dtype == 'float32'
assert theano_code(x + y, dtypes={
x: 'float64',
y: 'float32'
}).type.dtype == 'float64'
def test_symbol():
xt = theano_code(x)
assert isinstance(xt, (tt.TensorVariable, ts.ScalarVariable))
assert xt.name == x.name
assert theano_code(x, broadcastables={
x: (False, )
}).broadcastable == (False, )
assert theano_code(x, broadcastables={x: (False, )}).name == x.name
def test_AppliedUndef():
t = sy.Symbol('t')
f = sy.Function('f')
cache = {}
ft = theano_code(f(t), cache=cache)
assert isinstance(ft, tt.TensorVariable)
assert ft.name == 'f_t'
assert theano_code(f(t), cache=cache) is ft
assert theano_code(f(t), cache={}) is not ft
def test_add():
expr = x + y
comp = theano_code(expr)
assert comp.owner.op == theano.tensor.add
comp = theano_code(expr, broadcastables={x: (False, ), y: (False, )})
assert comp.broadcastable == (False, )
comp = theano_code(expr,
broadcastables={
x: (False, True),
y: (False, False)
})
assert comp.broadcastable == (False, False)
def test_MatrixSlice():
n = diofant.Symbol('n', integer=True)
X = diofant.MatrixSymbol('X', n, n)
Y = X[1:2:3, 4:5:6]
Yt = theano_code(Y)
# assert tuple(Yt.owner.op.idx_list) == (slice(1,2,3), slice(4,5,6))
assert Yt.owner.inputs[0] == theano_code(X)
k = diofant.Symbol('k')
kt = theano_code(k, dtypes={k: 'int32'})
start, stop, step = 4, k, 2
Y = X[start:stop:step]
Yt = theano_code(Y, dtypes={n: 'int32', k: 'int32'})
def test_MatrixSlice_2():
n = diofant.Symbol('n', integer=True)
X = diofant.MatrixSymbol('X', n, n)
Y = X[1:2:3, 4:5:6]
Yt = theano_code(Y)
assert tuple(Yt.owner.op.idx_list) == (slice(1, 2, 3), slice(4, 5, 6))
def test_DenseMatrix():
t = sy.Symbol('theta')
for MatrixType in [sy.Matrix, sy.ImmutableMatrix]:
X = MatrixType([[sy.cos(t), -sy.sin(t)], [sy.sin(t), sy.cos(t)]])
tX = theano_code(X)
assert isinstance(tX, tt.TensorVariable)
assert tX.owner.op == tt.join
def test_Piecewise():
# A piecewise linear
xt, yt = theano_code(x), theano_code(y)
expr = sy.Piecewise((0, x < 0), (x, x < 2), (1, True)) # ___/III
result = theano_code(expr)
assert result.owner.op == tt.switch
expected = tt.switch(xt < 0, 0, tt.switch(xt < 2, xt, 1))
assert theq(result, expected)
expr = sy.Piecewise((x, x < 0))
result = theano_code(expr)
expected = tt.switch(xt < 0, xt, np.nan)
assert theq(result, expected)
expr = sy.Piecewise((0, sy.And(x > 0, x < 2)), (x, sy.Or(x > 2, x < 0)))
result = theano_code(expr)
expected = tt.switch(tt.and_(xt > 0, xt < 2), 0,
tt.switch(tt.or_(xt > 2, xt < 0), xt, np.nan))
assert theq(result, expected)
def test_BlockMatrix():
n = diofant.Symbol('n', integer=True)
A = diofant.MatrixSymbol('A', n, n)
B = diofant.MatrixSymbol('B', n, n)
C = diofant.MatrixSymbol('C', n, n)
D = diofant.MatrixSymbol('D', n, n)
At, Bt, Ct, Dt = map(theano_code, (A, B, C, D))
Block = diofant.BlockMatrix([[A, B], [C, D]])
Blockt = theano_code(Block)
solutions = [
tt.join(0, tt.join(1, At, Bt), tt.join(1, Ct, Dt)),
tt.join(1, tt.join(0, At, Ct), tt.join(0, Bt, Dt))
]
assert any(theq(Blockt, solution) for solution in solutions)
def test_trig():
assert theq(theano_code(diofant.sin(x)), tt.sin(xt))
assert theq(theano_code(diofant.tan(x)), tt.tan(xt))
def test_MatMul():
X = diofant.MatrixSymbol('X', 4, 4)
Y = diofant.MatrixSymbol('X', 4, 4)
Z = diofant.MatrixSymbol('X', 4, 4)
expr = X * Y * Z
assert isinstance(theano_code(expr).owner.op, tt.Dot)
def test_Transpose():
X = diofant.MatrixSymbol('X', 4, 4)
assert isinstance(theano_code(X.T).owner.op, tt.DimShuffle)
def test_Relationals():
xt, yt = theano_code(x), theano_code(y)
assert theq(theano_code(x > y), xt > yt)
assert theq(theano_code(x < y), xt < yt)
assert theq(theano_code(x >= y), xt >= yt)
assert theq(theano_code(x <= y), xt <= yt)
def test_cache():
sx = sy.Symbol('x')
cache = {}
tx = theano_code(sx, cache=cache)
assert theano_code(sx, cache=cache) is tx
assert theano_code(sx, cache={}) is not tx
def test_Integers():
assert theano_code(diofant.Integer(3)) == 3
def test_symbols_are_created_once():
expr = x**x
comp = theano_code(expr)
assert theq(comp, xt**xt)
def test_MatrixSymbol():
X = diofant.MatrixSymbol('X', 4, 5)
Xt = theano_code(X)
assert isinstance(Xt, tt.TensorVariable)
assert Xt.broadcastable == (False, False)
def test_slice():
assert theano_code(slice(1, 2, 3)) == slice(1, 2, 3)
assert str(theano_code(slice(1, x, 3), dtypes={x: 'int32'})) ==\
str(slice(1, xt, 3))
def test_Derivative():
def simp(expr):
return theano_simplify(fgraph_of(expr))
assert theq(simp(theano_code(sy.Derivative(sy.sin(x), x, evaluate=False))),
simp(theano.grad(tt.sin(xt), xt)))
def test_factorial():
n = diofant.Symbol('n')
assert theano_code(diofant.factorial(n))
def test_many():
expr = sy.exp(x**2 + sy.cos(y)) * sy.log(2 * z)
comp = theano_code(expr)
expected = tt.exp(xt**2 + tt.cos(yt)) * tt.log(2 * zt)
def test_MatAdd():
X = diofant.MatrixSymbol('X', 4, 4)
Y = diofant.MatrixSymbol('X', 4, 4)
Z = diofant.MatrixSymbol('X', 4, 4)
expr = X + Y + Z
assert isinstance(theano_code(expr).owner.op, tt.Elemwise)
def test_AppliedUndef():
t = sy.Symbol('t')
f = sy.Function('f')
ft = theano_code(f(t))
assert isinstance(ft, tt.TensorVariable)
assert ft.name == 'f_t'
def test_Rationals():
assert theq(theano_code(diofant.Rational(2, 3)), tt.true_div(2, 3))
assert theq(theano_code(S.Half), tt.true_div(1, 2))
