Eu estou tentando usar oscipy.optimize pacote para otimizar um problema de otimização discreta (otimização global). Acc para o doc, simulated annealing implementado emscipy.optimize.anneal deve ser uma boa escolha para o mesmo. Mas não tenho certeza de como forçar o otimizador a pesquisar somente valores inteiros do espaço de pesquisa. Alguém pode ajudar?

Um exemplo ilustrativo:

f(x1,x2) = (1-0.4*x1)^2 + 100*(0.6*x2 -0.4*x1^2)^2

where, $x1, x2 \in I$