Как сделать код gekko более эффективным ⇐ Python

[Anonymous]

Мне нужно решить следующую проблему


Я переписал задачу в Gekko на непрерывную оптимизацию вот так:

Код: Выделить всё

pt_values = dataframe.select("pt").to_numpy().ravel()
rt_values= dataframe.select("rt").to_numpy().ravel()
number_periods = df.shape[0]
# initialize gekko
m = GEKKO(remote=True)
# set time periods
m.time = np.linspace(0, number_periods-1, number_periods )
# Variables
kt =m.SV(value=ks,lb=0,ub=k_overline,integer=True)
rev_t =m.SV(value=revs)
zt = m.Var(value=zs, lb=-k_overline, ub=r_overline, integer=False)
pt = m.Param(value=pt_value)
rt = m.Param(value=rt_value)

#equations
m.Equation(kt.dt()==m.if3(rt-zt,m.max2(rt-zt,-hr),m.min2(m.min2(k_overline-kt+kt.dt(), rt-zt), hr)))
m.Equation(rev_t.dt()==zt*pt)
#restrictions on zt
m.Equation(kt-kt.dt()+rt>=zt)
p = np.zeros(number_periods)
p[-1] = 1.0
final = m.Param(value=p)

m.Maximize(rev_t*final)
m.options.IMODE = 6
m.solve(disp=False, debug=0)

Проблема в том, что мой df.shape[0]=8000. Это очень много. По этой причине я перешел на сплошные выражения. Дискретная настройка очень быстро выходит из строя. Однако это все еще проблема. Жду предложений, как улучшить. Я готов перейти на другие пакеты оптимизации.

Подробнее здесь: https://stackoverflow.com/questions/788 ... -efficient