FIX: Player.is_dominated: Fix warnings

Fix #501

Author: oyamad

quantecon/game_theory/normal_form_game.py

@@ -472,17 +472,18 @@ def is_dominated(self, action, tol=None, method=None):
         elif method in ['simplex', 'interior-point']:
             from scipy.optimize import linprog
             m, n = D.shape
-            A = np.empty((n+2, m+1))
-            A[:n, :m] = -D.T
-            A[:n, -1] = 1  # Slack variable
-            A[n, :m], A[n+1, :m] = 1, -1  # Equality constraint
-            A[n:, -1] = 0
-            b = np.empty(n+2)
-            b[:n] = 0
-            b[n], b[n+1] = 1, -1
+            A_ub = np.empty((n, m+1))
+            A_ub[:, :m] = -D.T
+            A_ub[:, -1] = 1  # Slack variable
+            b_ub = np.zeros(n)
+            A_eq = np.empty((1, m+1))
+            A_eq[:, :m] = 1  # Equality constraint
+            A_eq[:, -1] = 0
+            b_eq = np.ones(1)
             c = np.zeros(m+1)
             c[-1] = -1
-            res = linprog(c, A_ub=A, b_ub=b, method=method)
+            res = linprog(c, A_ub=A_ub, b_ub=b_ub, A_eq=A_eq, b_eq=b_eq,
+                          method=method)
             if res.success:
                 return res.x[-1] > tol
             elif res.status == 2:  # infeasible

quantecon/game_theory/tests/test_normal_form_game.py

@@ -162,7 +162,7 @@ def test_player_corner_cases():
         for method in LP_METHODS:
             eq_(player.is_dominated(action, method=method), False)
 
-    e = 1e-8
+    e = 1e-8 * 2
     player = Player([[-e, -e], [1, -1], [-1, 1]])
     action = 0
     eq_(player.is_best_response(action, [1/2, 1/2], tol=e), True)