"""
This is a gamma-law equation of state: p = rho e (gamma - 1), where
gamma is the constant ratio of specific heats.
"""



[docs]
def pres(gamma, rho, eint):
    """
    Given the density and the specific internal energy, return the
    pressure

    Parameters
    ----------
    gamma : float
        The ratio of specific heats
    rho : float
        The density
    eint : float
        The specific internal energy

    Returns
    -------
    out : float
       The pressure

     """
    p = rho * eint * (gamma - 1.0)
    return p




[docs]
def dens(gamma, p, eint):
    """
    Given the pressure and the specific internal energy, return
    the density

    Parameters
    ----------
    gamma : float
        The ratio of specific heats
    p : float
        The pressure
    eint : float
        The specific internal energy

    Returns
    -------
    out : float
       The density

    """
    rho = p / (eint * (gamma - 1.0))
    return rho




[docs]
def rhoe(gamma, p):
    """
    Given the pressure, return (rho * e)

    Parameters
    ----------
    gamma : float
        The ratio of specific heats
    p : float
        The pressure

    Returns
    -------
    out : float
       The internal energy density, rho e

    """
    return p / (gamma - 1.0)