1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
import unittest
from polymatrix.expressionstate.init.initexpressionstate import init_expression_state
import polymatrix.expression.initexpressionbase
class TestQuadraticIn(unittest.TestCase):
def test_1(self):
underlying_terms = {
(0, 0): {
((0, 1),): 1.0, # x1
((0, 1), (2, 1)): 2.0, # x1
((0, 2), (3, 1)): 3.0, # x1 x1
((0, 2), (1, 2), (4, 1)): 4.0, # x1 x1 x2 x2
((0, 2), (1, 1), (5, 1)): 5.0, # x1 x1 x2
}
}
monomial_terms = {
(0, 0): {
tuple(): 1.0,
},
(1, 0): {
((0, 1),): 1.0,
},
(2, 0): {
((0, 1), (1, 1)): 1.0,
},
}
variable_terms = {
(0, 0): {((0, 1),): 1},
(1, 0): {((1, 1),): 1},
}
expr = polymatrix.expression.initexpressionbase.init_quadratic_in_expr(
underlying=polymatrix.expression.initexpressionbase.init_from_terms_expr(terms=underlying_terms, shape=(1, 1)),
monomials=polymatrix.expression.initexpressionbase.init_from_terms_expr(terms=monomial_terms, shape=(3, 1)),
variables=polymatrix.expression.initexpressionbase.init_from_terms_expr(terms=variable_terms, shape=(2, 1),),
)
state = init_expression_state(n_param=2)
state, val = expr.apply(state)
data = val.get_poly(0, 1)
self.assertDictContainsSubset({
tuple(): 1.0,
((2, 1),): 2.0,
}, data)
data = val.get_poly(1, 1)
self.assertDictContainsSubset({
((3, 1),): 3.0,
}, data)
data = val.get_poly(2, 2)
self.assertDictContainsSubset({
((4, 1),): 4.0,
}, data)
data = val.get_poly(1, 2)
self.assertDictContainsSubset({
((5, 1),): 5.0,
}, data)
|
