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
