def cost(rv): """compute the diameter as a calculation of cost Input: - rv -- 1-d array of model parameters Output: - diameter -- scale * | F(x) - F(x')|**2 """ from surrogate import marc_surr as model # prepare x and xprime rv = list(rv) params = rv[:-1] #XXX: assumes Xi' is at rv[-1] params_prime = rv[:i]+rv[-1:]+rv[i+1:-1] #XXX: assumes Xi' is at rv[-1] # get the F(x) response Fx = model(params) # get the F(x') response Fxp = model(params_prime) # compute diameter scale = 1.0 return -scale * (Fx - Fxp)**2
def cost(rv): """compute the diameter as a calculation of cost Input: - rv -- 1-d array of model parameters Output: - diameter -- scale * | F(x) - F(x')|**2 """ from surrogate import marc_surr as model # prepare x and xprime rv = list(rv) params = rv[:-1] #XXX: assumes Xi' is at rv[-1] params_prime = rv[:i] + rv[-1:] + rv[ i + 1:-1] #XXX: assumes Xi' is at rv[-1] # get the F(x) response Fx = model(params) # get the F(x') response Fxp = model(params_prime) # compute diameter scale = 1.0 return -scale * (Fx - Fxp)**2