constraint 4 idea

constraint4.py

#this is the constraint that makes Pyhy[i,i] > Pyhy[i,j]
def con4(inputArray):
    Pyhy_i = inputArray.reshape(n,n)
    # get diagonal vector (Pyhy[i,i])
    Pyhy_diag = np.diag(Pyhy_i)
    # get max of each row
    Pyhy_max = Pyhy.max(axis=1)
    # we want the difference to be positive or zero
    diff = Pyhy_diag - Pyhy_max
    # we want all of the differences to be positive, so we can just return the most negative one
    worst_case = diff.min()
    return worst_case