Any linear gradient based search algorithm can be used or if there is possibility of several local minima, a simulated annealing or genetic algorithm can be employed.
the adaptive simulated annealing method is proposed for obtaining the pareto-optimal setpoints.