def dG_to_flux_bounds(dG_min, dG_max, infinity=1000, abstol=1e-06):
    """ Convert standard Gibbs energy range to reaction flux bounds.

    Args:
        dG_min (float): minimum standard Gibbs energy
        dG_max (float): maximum standard Gibbs energy
        infinity (float): value to represent infinity (default: 1000)
        abstol (float): absolute tolerance to consider a value larger than zero (default: 1e-6)

    Returns:

    """
    if dG_min > abstol:
        lb, ub = -infinity, 0
    elif dG_max < -abstol:
        lb, ub = 0, infinity
    else:
        lb, ub = -infinity, infinity
    return lb, ub
