def f1(x):
score = interior_score(radius, density, x[0], x[1])
objective = -score + equality_penalty(
x[0] + x[1] + x[-1] - 1, tolerance) + inequality_penalty(
-x[0], error) + inequality_penalty(
x[0] - 1, error) + inequality_penalty(
-x[1], error) + inequality_penalty(x[1] - 1, error)
return objective
def f2(x):
score = surface_score(escape_velocity, surface_temperature, x[0], x[1])
objective = -score + inequality_penalty(
(x[0] + x[1] - 1),
error) + inequality_penalty(-x[0], error) + inequality_penalty(
x[0] - 1, error) + inequality_penalty(
-x[1], error) + inequality_penalty(x[1] - 1, error)
return objective
def f1(x):
score = interior_score(radius, density, x[0], x[1])
objective = -score + modified_inequality_penalty(
(x[0] + x[1] - 1)) + inequality_penalty(
-x[0], error) + inequality_penalty(
x[0] - 1, error) + inequality_penalty(
-x[1], error) + inequality_penalty(x[1] - 1, error)
return objective
def f1(x):
score = ceesa_score(escape_velocity, surface_temperature, radius,
density, orb_eccentricity, x[0], x[1], x[2], x[3],
x[4], x[5], x[6])
objective = -score + inequality_penalty(x[0]-1,error) + inequality_penalty(-x[1],error) + inequality_penalty(-x[2],error) + inequality_penalty(-x[3],error) \
+ inequality_penalty(-x[4],error) + inequality_penalty(-x[5],error) + inequality_penalty(-x[6],error) \
+ inequality_penalty(x[1]-1,error) + inequality_penalty(x[2]-1,error) + inequality_penalty(x[3]-1,error) + inequality_penalty(x[4]-1,error) + inequality_penalty(x[5]-1,error) + inequality_penalty(x[6]-1,error) \
+ equality_penalty(x[2] + x[3] + x[4] + x[5] + x[6] -1,tolerance) + equality_penalty(1-(x[6]-x[5]-x[4]-x[3]-x[2]),tolerance)
return objective
