void compute_bounds(NumericVector<Number> & lower, NumericVector<Number> & upper, NonlinearImplicitSystem & sys) { FEProblemBase * p = sys.get_equation_systems().parameters.get<FEProblemBase *>("_fe_problem_base"); p->computeBounds(sys, lower, upper); }