# 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
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]
# 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)