MAINT: optimize : migration to sparray pass 1 changes (scipy#21785)

* pass 1 changes for migration to sparray

* some more changes to allow sparrays to work

Author: dschult

scipy/optimize/_linprog_util.py

@@ -1215,7 +1215,7 @@ def vstack(blocks):
     if sparse:
         b = b.reshape(-1, 1)
         A = A.tocsc()
-        b -= (A[:, i_shift] * sps.diags(lb_shift)).sum(axis=1)
+        b -= (A[:, i_shift] @ sps.diags(lb_shift)).sum(axis=1)
         b = b.ravel()
     else:
         b -= (A[:, i_shift] * lb_shift).sum(axis=1)
@@ -1249,15 +1249,15 @@ def _autoscale(A, b, c, x0):
             R = R.toarray().flatten()
         R[R == 0] = 1
         R = 1/_round_to_power_of_two(R)
-        A = sps.diags(R)*A if sps.issparse(A) else A*R.reshape(m, 1)
+        A = sps.diags(R)@A if sps.issparse(A) else A*R.reshape(m, 1)
         b = b*R
 
         C = np.max(np.abs(A), axis=0)
         if sps.issparse(A):
             C = C.toarray().flatten()
         C[C == 0] = 1
         C = 1/_round_to_power_of_two(C)
-        A = A*sps.diags(C) if sps.issparse(A) else A*C
+        A = A@sps.diags(C) if sps.issparse(A) else A*C
         c = c*C
 
     b_scale = np.max(np.abs(b)) if b.size > 0 else 1

scipy/optimize/tests/test_linprog.py

@@ -521,7 +521,7 @@ def f(c, A_ub=None, b_ub=None, A_eq=None, b_eq=None, bounds=None):
         ub = x_valid + np.random.rand(n)
         lb = x_valid - np.random.rand(n)
         bounds = np.column_stack((lb, ub))
-        b_eq = A_eq * x_valid
+        b_eq = A_eq @ x_valid
 
         if self.method in {'simplex', 'revised simplex'}:
             # simplex and revised simplex should raise error