if soln.game_value() > epsilon:
# L should be positive stable
positive_stable = all([eig > -options.ABS_TOL for eig in eigs])
- if not positive_stable:
- print(str(eigs))
self.assertTrue(positive_stable)
elif soln.game_value() < -epsilon:
# L should be negative stable
negative_stable = all([eig < options.ABS_TOL for eig in eigs])
- if not negative_stable:
- print(str(eigs))
self.assertTrue(negative_stable)
# The dual game's value should always equal the primal's.
while G._condition() > MAX_COND:
# Try again until the condition number is satisfactory.
- G = random_orthant_game()
+ G = random_icecream_game()
L = random_lyapunov_like_icecream(G._K.dimension())
G = SymmetricLinearGame(L, G._K, G._e1, G._e2)