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