sig
  type params = {
    iters_fixed_T : int;
    step_size : float;
    k : float;
    t_initial : float;
    mu_t : float;
    t_min : float;
  }
  val solve :
    Gsl_rng.t ->
    'a ->
    energ_func:('a -> float) ->
    step_func:(Gsl_rng.t -> 'a -> float -> 'a) ->
    ?print_func:('a -> unit) -> Gsl_siman.params -> 'a
end