def test_symbols():
a, b = Symbol('a'), Symbol('b')
d = PUSOMatrix()
d[(0, )] -= a
d[(0, 1)] += 2
d[(1, )] += b
assert d == {(0, ): -a, (0, 1): 2, (1, ): b}
assert d.subs(a, 2) == {(0, ): -2, (0, 1): 2, (1, ): b}
assert d.subs(b, 1) == {(0, ): -a, (0, 1): 2, (1, ): 1}
assert d.subs({a: -3, b: 4}) == {(0, ): 3, (0, 1): 2, (1, ): 4}
def test_puso_remove_value_when_zero():
d = PUSOMatrix()
d[(0, )] += 1
d[(0, )] -= 1
assert d == {}
d.refresh()
assert d.degree == 0
assert d.num_binary_variables == 0
assert d.variables == set()
def test_create_var():
d = PUSOMatrix.create_var(0)
assert d == {(0, ): 1}
assert d.name == 0
assert type(d) == PUSOMatrix
d = PUSOMatrix.create_var(1)
assert d == {(1, ): 1}
assert d.name == 1
assert type(d) == PUSOMatrix
def test_puso_to_pubo_to_puso():
puso = {
(0, 1): -4,
(0, 2): 3,
(): -2,
(0, ): 1,
(2, ): -2,
(0, 1, 2): 3,
(0, 2, 3): -1
}
assert puso == pubo_to_puso(puso_to_pubo(puso))
# type asserting
assert type(puso_to_pubo(puso)) == PUBO
assert type(puso_to_pubo(PUSOMatrix(puso))) == PUBOMatrix
assert type(puso_to_pubo(PUSO(puso))) == PUBO
puso = {
('0', 1): -4,
('0', '2'): 3,
(): -2,
('0', ): 1,
('2', ): -2,
('0', 1, '2'): 3,
('0', '2', 3): -1,
('2', 0, 0, '1', 0): -2,
(0, 1, 1, 0, 3, 0, 1, 1, 3, 2, 3): -8
}
# need to reformat qubo so it is sorted with the same hash
assert PUSO(puso) == pubo_to_puso(puso_to_pubo(puso))
# type asserting
assert type(puso_to_pubo(puso)) == PUBO
assert type(puso_to_pubo(PUSO(puso))) == PUBO
def test_round():
d = PUSOMatrix({(0, ): 3.456, (1, ): -1.53456})
assert round(d) == {(0, ): 3, (1, ): -2}
assert round(d, 1) == {(0, ): 3.5, (1, ): -1.5}
assert round(d, 2) == {(0, ): 3.46, (1, ): -1.53}
assert round(d, 3) == {(0, ): 3.456, (1, ): -1.535}
def test_normalize():
temp = {(0, ): 4, (1, ): -2}
d = PUSOMatrix(temp)
d.normalize()
assert d == {k: v / 4 for k, v in temp.items()}
temp = {(0, ): -4, (1, ): 2}
d = PUSOMatrix(temp)
d.normalize()
assert d == {k: v / 4 for k, v in temp.items()}
def test_puso_addition():
temp = PUSOMatrix({(0, ): 1, (0, 1): 2, (2, 1, 0): -1})
temp1 = {(0, ): -1, (1, 0): 3}
temp2 = {(0, 1): 5, (0, 1, 2): -1}
temp3 = {(0, ): 2, (0, 1): -1, (0, 1, 2): -1}
# add constant
d = temp.copy()
d += 5
d[()] -= 2
d == {(0, ): 1, (0, 1): 2, (): 3, (0, 1, 2): -1}
# __add__
d = temp.copy()
g = d + temp1
assert g == temp2
# __iadd__
d = temp.copy()
d += temp1
assert d == temp2
# __radd__
d = temp.copy()
g = temp1 + d
assert g == temp2
# __sub__
d = temp.copy()
g = d - temp1
assert g == temp3
# __isub__
d = temp.copy()
d -= temp1
assert d == temp3
# __rsub__
d = temp.copy()
g = temp1 - d
assert g == PUSOMatrix(temp3) * -1
def test_puso_degree():
d = PUSOMatrix()
assert d.degree == 0
d[(0, )] += 2
assert d.degree == 1
d[(1, )] -= 3
assert d.degree == 1
d[(1, 2)] -= 2
assert d.degree == 2
d[(0, 1, 2)] -= 1
assert d.degree == 3
d[(0, 1, 2, 4, 8)] -= 1
assert d.degree == 5
def test_pretty_str():
def equal(expression, string):
assert expression.pretty_str() == string
z = [PUSOMatrix() + {(i, ): 1} for i in range(3)]
a, b = Symbol('a'), Symbol('b')
equal(z[0], "z(0)")
equal(-z[0], "-z(0)")
equal(z[0] * 0, "0")
equal(2 * z[0] * z[1] - 3 * z[2], "2 z(0) z(1) - 3 z(2)")
equal(0 * z[0] + 1, "1")
equal(0 * z[0] - 1, "-1")
equal(0 * z[0] + a, "(a)")
equal(0 * z[0] + a * b, "(a*b)")
equal((a + b) * (z[0] * z[1] - z[2]), "(a + b) z(0) z(1) + (-a - b) z(2)")
equal(2 * z[0] * z[1] - z[2], "2 z(0) z(1) - z(2)")
equal(-z[2] + z[0] * z[1], "-z(2) + z(0) z(1)")
equal(-2 * z[2] + 2 * z[0] * z[1], "-2 z(2) + 2 z(0) z(1)")
def test_pusomatrix_solve_bruteforce():
H = PUSOMatrix({
(0, 1): 1,
(1, 2): 1,
(1, ): -1,
(2, ): -2,
(3, 4, 5): -1,
(3, ): -1,
(4, ): -1,
(5, ): -1
})
sol = H.solve_bruteforce()
assert sol in (
{
0: -1,
1: 1,
2: 1,
3: 1,
4: 1,
5: 1
},
{
0: 1,
1: -1,
2: 1,
3: 1,
4: 1,
5: 1
},
)
assert allclose(H.value(sol), -7)
H = PUSOMatrix({
(0, ): 0.25,
(1, ): -0.25,
(0, 1): -0.25,
(): 1.25,
(3, 4, 5): -1,
(3, ): -1,
(4, ): -1,
(5, ): -1
})
sols = H.solve_bruteforce(True)
assert (sols == [{
0: 1,
1: 1,
3: 1,
4: 1,
5: 1
}, {
0: -1,
1: 1,
3: 1,
4: 1,
5: 1
}, {
0: -1,
1: -1,
3: 1,
4: 1,
5: 1
}])
assert all(allclose(H.value(s), -3) for s in sols)
def test_num_terms():
d = PUSOMatrix({(0, ): 1, (0, 3): 2, (0, 2): -1})
assert d.num_terms == len(d)
def test_puso_max_index():
d = PUSOMatrix({(0, ): 1, (0, 3): 2, (0, 4, 3): 3})
assert d.max_index == 4
def test_puso_num_binary_variables():
d = PUSOMatrix({(0, ): 1, (0, 3): 2, (0, 4, 3): 3})
assert d.num_binary_variables == 3
def test_puso_update():
d = PUSOMatrix({(0, ): 1, (0, 1): 2})
d.update({(0, ): 0, (1, 0): 1, (1, ): -1, (0, 2, 1): -1})
assert d == {(0, 1): 1, (1, ): -1, (0, 1, 2): -1}
def test_puso_reinitialize_dictionary():
d = PUSOMatrix({(0, ): 1, (1, 0): 2, (2, 0): 0, (0, 1): 1, (2, 0, 1): 1})
assert d == {(0, ): 1, (0, 1): 3, (0, 1, 2): 1}
def test_puso_multiplication():
temp = PUSOMatrix({(0, ): 1, (0, 1): 2, (0, 2, 3): 4})
temp += 2
# __mul__
d = temp.copy()
g = d * 3
assert g == {(0, ): 3, (0, 1): 6, (): 6, (0, 2, 3): 12}
d = temp.copy()
g = d * 0
assert g == {}
# __imul__
d = temp.copy()
d *= 3
assert d == {(0, ): 3, (0, 1): 6, (): 6, (0, 2, 3): 12}
d = temp.copy()
d *= 0
assert d == {}
# __rmul__
d = temp.copy()
g = 3 * d
assert g == {(0, ): 3, (0, 1): 6, (): 6, (0, 2, 3): 12}
d = temp.copy()
g = 0 * d
assert g == {}
# __truediv__
d = temp.copy()
g = d / 2
assert g == {(0, ): .5, (0, 1): 1, (): 1, (0, 2, 3): 2}
# __itruediv__
d = temp.copy()
d /= 2
assert d == {(0, ): .5, (0, 1): 1, (): 1, (0, 2, 3): 2}
# __floordiv__
d = temp.copy()
g = d // 2
assert g == {(0, 1): 1, (): 1, (0, 2, 3): 2}
# __ifloordiv__
d = temp.copy()
d //= 2
assert d == {(0, 1): 1, (): 1, (0, 2, 3): 2}
# __mul__ but with dict
d = temp.copy()
d *= {(1, ): 2, (0, ): -1}
assert d == {
(0, 1): 2,
(): -1,
(0, ): 2,
(1, ): 2,
(0, 1, 2, 3): 8,
(2, 3): -4
}
# __pow__
d = temp.copy()
d **= 2
assert d == {
(): 25,
(1, ): 4,
(2, 3): 8,
(0, ): 4,
(1, 2, 3): 16,
(0, 1): 8,
(0, 2, 3): 16
}
temp = d.copy()
assert d**3 == d * d * d
temp = d.copy()
assert d**4 == d * d * d * d
def test_puso_default_valid():
d = PUSOMatrix()
assert d[(0, )] == 0
d[(0, )] += 1
assert d == {(0, ): 1}