Глобальная невыпуклая оптимизация, посоветуйте софт

[Xaositect]

Имеется NP-трудная задача, которую хочется в маленьких размерностях попытаться все-таки поподбирать численно.
Численно она получается очень неприятная, нужно оптимизировать невыпуклые полиномы, причем хочется глобальный минимум. Полиномы 6 степени, на интересующей нас задаче получаются от 594 переменных. При использовании обычных методов очень часто попадаем в локальные минимумы.
Существуют ли готовые программные пакеты для нахождения глобальных минимумов неотрицательных полиномов 6 степени от ~500 переменных или я слишком многого хочу?
Полиномы связаны с разложениями тензоров и имеют вид $\sum_{i,j,k} (\sum_r x_{0ir}x_{1jr}x_{2ir} - a_{ijk})^2$.