У меня есть две ошибки: «RuntimeWarning: деление на ноль, встречающееся в double_scalars»; «RuntimeWarning: недопустимое значение, встречающееся в вычитании» - PullRequest
Новая
       54

У меня есть две ошибки: «RuntimeWarning: деление на ноль, встречающееся в double_scalars»; «RuntimeWarning: недопустимое значение, встречающееся в вычитании»

1 голос
/ 28 апреля 2020

Я не понимаю, как python делится на ноль. Это весь мой код: 

from scipy import optimize
import numpy as np
from sympy import *

from tkinter import *
#import ipdb;ipdb.set_trace()
y=[100,200]
alp=[0.5,0.5]
end=[1,2]
e, t, a, w = symbols('e t a w')
u= ((a ** a) * (24 * (1 - t) * w + e) * ((1 - a) ** (1 - a))) * (((1 - t) * w) ** (a-1))

print(latex(u))
def main(x, n):

    r = 0
    p = 0
    #j=-1 * ((alp[i] ** alp[i]) * (24 * (1 - x[i]) * y[i] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) / ((1 - x[1] * y[i]) ** (1 - alp[i]))
    for i in range(len(y)):
        if i >= n:

            r += -1 * ((alp[i] ** alp[i]) * (24 * (1 - x[1]) * y[i] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) * (((1 - x[1]) * y[i]) ** (alp[i]-1))


        elif i < n:

            p += -1 * ((alp[i] ** alp[i]) * (24 * y[i] + x[0] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) * ((y[i]) ** (alp[i]-1))
    return r + p

def tax_rev(n):
    r = 0
    for i in range(n, len(y)):
        r += y[i]

    return r

bounds = [(0, np.inf), (0, 1)]

mat = np.zeros((len(y), 4))


for i in range(len(y)):
    # print(i)
    def constraint1(x):
        return 1 * (x[1] * tax_rev(i) - x[0] * (len(y) -i))

    cons1 = {'type': 'ineq', 'fun': constraint1}

    def co(x):
        return main(x, i)

    max = optimize.shgo(co, bounds=bounds, constraints=cons1)

    mat[i, 0] = round(max.x[0] ,5)
    mat[i, 1] = round(max.x[1] ,5)
    mat[i, 2] = -1 * round(max.fun ,5)
    mat[i, 3] = max.success

tax_from= np.unravel_index(np.argmax(mat, axis=None), mat.shape)

fin_tex = f'{100 * round(float(mat[tax_from[0] ,1]) ,15)}  % from those endowed with {y[tax_from[0]]}  or more. Transfer {round(float(mat[tax_from[0] ,0]) ,15)} units of money to those endowed with less than  {y[tax_from[0]]} .This would maximize utility in the society.'
print(fin_tex)
print(mat)

Я использую этот код для оптимизации уравнения, которое можно увидеть в латексе (уравнение в переменной u). Я не понимаю, почему существует деление на ноль, когда я не делю что-либо на 0. Я думаю, что это вызывает проблему, когда я также оптимизирую.

Деление на ноль, согласно сообщению об ошибке I я получаю здесь: 

def main(x, n):

    r = 0
    p = 0
    #j=-1 * ((alp[i] ** alp[i]) * (24 * (1 - x[i]) * y[i] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) / ((1 - x[1] * y[i]) ** (1 - alp[i]))
    for i in range(len(y)):
        if i >= n:

            r += -1 * ((alp[i] ** alp[i]) * (24 * (1 - x[1]) * y[i] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) * (((1 - x[1]) * y[i]) ** (alp[i]-1))


        elif i < n:

            p += -1 * ((alp[i] ** alp[i]) * (24 * y[i] + x[0] + end[i]) * ((1 - alp[i]) ** (1 - alp[i]))) * ((y[i]) ** (alp[i]-1))
    return r + p

Это происходит в переменной, которую я назвал r.

Я новичок в Python, и любая помощь / совет по моему коду приветствуется :)

1 Ответ

0 голосов
/ 28 апреля 2020

Полное предупреждение: 

RuntimeWarning: divide by zero encountered in double_scalars
  if __name__ == '__main__':
/usr/local/lib/python3.6/dist-packages/scipy/optimize/slsqp.py:63: RuntimeWarning: invalid value encountered in subtract
  jac[i] = (func(*((x0+dx,)+args)) - f0)/epsilon
/usr/local/bin/ipython3:9: RuntimeWarning: invalid value encountered in double_scalars

Таким образом, оптимизатор (slsqp) оценивает якобиан и делится на epsilon. И делает вычитание.

В этом вопросе много ненужного багажа. Например, что sympy/latex бизнес. Вы не используете это в оптимизации.

Основываясь на контексте предупреждения, вам нужно написать код, который тестирует оптимизацию, без этого i l oop.

...