docs: def _get_fba_problem(model, tfba, solver):
      """Convenience function for returning the right FBA problem instance"""
      p = FluxBalanceProblem(model, solver)
      if tfba:
          p.add_thermodynamic()
      return p