Function optimize_da 

pub fn optimize_da(
    obj: &impl Objective,
    guess: &[f64],
    lower: Vec<f64>,
    upper: Vec<f64>,
    p: &DaParams,
) -> DaResult

Run Dual Annealing. lower/upper empty ⇒ unbounded.