from numpy import *
from openopt import NSP

N = 10

def objFun(x):
    r = 0.0
    for j in range(N):
        for i in range(N):
            r += abs(x[i]) * (cos(i) + sin(j) + 15)
    return r

x0 = cos(1+asfarray(range(N)))

p = NSP(objFun, x0, maxFunEvals = 1e7, xtol = 1e-8)

r = p.solve('ralg', iprint = 10, fEnough = 0.01)

# intel Atom 1.66 GHz:
# 20.6 sec in CPython 2.7
# 9.9 sec in PyPy 1.8