def __init__(self, worker_id, data, response_surface, maxeval, nsamples,
exp_design=None, search_procedure=None, extra=None):
"""Initialize the optimization strategy.
:param worker_id: Start ID in a multistart setting
:param data: Problem parameter data structure
:param response_surface: Surrogate model object
:param maxeval: Function evaluation budget
:param nsamples: Number of simultaneous fevals allowed
:param exp_design: Experimental design
:param search_procedure: Search procedure for finding
points to evaluate
:param extra: Points to be added to the experimental design
"""
if self.__class__.__name__ == "SyncStrategyNoConstraints":
assert not hasattr(data, "eval_ineq_constraints"), "Objective function has constraints,\n" \
"SyncStrategyNoConstraints can't handle constraints"
self.worker_id = worker_id
self.data = data
self.fhat = response_surface
if self.fhat is None:
self.fhat = RBFInterpolant(surftype=CubicRBFSurface, eta=1e-8, maxp=maxeval)
self.maxeval = maxeval
self.nsamples = nsamples
self.extra = extra
# Default to generate sampling points using Symmetric Latin Hypercube
self.design = exp_design
if self.design is None:
if self.data.dim > 50:
self.design = LatinHypercube(data.dim, data.dim+1)
else:
self.design = SymmetricLatinHypercube(data.dim, 2*(data.dim+1))
self.xrange = np.asarray(data.xup - data.xlow)
# algorithm parameters
self.sigma_min = 0.005
self.sigma_max = 0.2
self.sigma_init = 0.2
self.failtol = max(5, data.dim)
self.succtol = 3
self.numeval = 0
self.status = 0
self.sigma = 0
self.resubmitter = RetryStrategy()
self.xbest = None
self.fbest = np.inf
self.fbest_old = None
# Set up search procedures and initialize
self.search = search_procedure
if self.search is None:
self.search = CandidateDYCORS(data)
# Start with first experimental design
self.sample_initial()
