Ejemplo n.º 1
0
# To change this license header, choose License Headers in Project Properties.
# To change this template file, choose Tools | Templates
# and open the template in the editor.

import OPTutils
import scipy as sp
from matplotlib import pyplot as plt
import GPdc
import ESutils

d=2
lb = sp.array([[-1.]*d])
ub = sp.array([[1.]*d])
[ojf,truexmin] = OPTutils.gensquexpdraw(d,sp.array([-1.]*d),sp.array([1.]*d),ignores=1e-3)

O = OPTutils.opt(ojf,lb,ub)
for i in xrange(20):
    O.step()


para = dict()
para['kindex'] = GPdc.SQUEXPCS
para['mprior'] = sp.array([0.]+[-1.]*d+[-2.])
para['sprior'] = sp.array([1.]*(d+1)+[2.])
#para['s'] = 1e-6
para['ninit'] = 10
#para['maxf'] = 2500
para['volper'] = 1e-7
para['DH_SAMPLES'] = 8
para['DM_SAMPLES'] = 8
para['DM_SUPPORT'] = 400
Ejemplo n.º 2
0
import sys

sys.path.append("./..")
import OPTutils
import scipy as sp
from matplotlib import pyplot as plt
import GPdc
import ESutils
from tqdm import tqdm, tqdm_gui

runn = 60
d = 2
lb = sp.array([[-1.0] * d])
ub = sp.array([[1.0] * d])

[ojf, truexmin, trueymin] = OPTutils.gensquexpdraw(d, sp.array([-1.0] * d), sp.array([1.0] * d))


kindex = GPdc.SQUEXP

mprior = sp.array([0.0] + [-1.0] * d)
sprior = sp.array([1.0] * (d + 1))
maxf = 4000
volper = 1e-7
s = 1e-5
ninit = 10

# para = [kindex,hlb,hub,maxf,s,ninit]
para = [kindex, mprior, sprior, volper, s, ninit]

Ejemplo n.º 3
0
# To change this template file, choose Tools | Templates
# and open the template in the editor.
import sys
sys.path.append("./..")
import OPTutils
import scipy as sp
from matplotlib import pyplot as plt
import GPdc
import ESutils
from tqdm import tqdm, tqdm_gui
runn = 60
d = 2
lb = sp.array([[-1.] * d])
ub = sp.array([[1.] * d])

[ojf, truexmin, trueymin] = OPTutils.gensquexpdraw(d, sp.array([-1.] * d),
                                                   sp.array([1.] * d))

kindex = GPdc.SQUEXP

mprior = sp.array([0.] + [-1.] * d)
sprior = sp.array([1.] * (d + 1))
maxf = 4000
volper = 1e-7
s = 1e-5
ninit = 10

#para = [kindex,hlb,hub,maxf,s,ninit]
para = [kindex, mprior, sprior, volper, s, ninit]

OE = OPTutils.EIMLE(ojf, lb, ub, para)