We report the bivariate $HI$ - and $H2$ -stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for $HI$ , $H2$ , cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the $HI$ and $H2$ conditional distributions and the galaxy stellar mass function (GSMF). For the conditional distributions we use the results from the compilation presented in Calette et al. [(2018) RMxAA, 54, 443.]. For determining the GSMF from $M_{*}∼3×107$ to $3×10^{12} M_{\odot}$ , we combine two spectroscopic samples from the Sloan Digital Sky Survey at the redshift range $0.0033<z<0.2$ . We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $α≈-1.4$ , consistent with previous determinations. The obtained $HI $ MFs agree with radio blind surveys. Similarly, the $H2 $ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $∼\! 96%$ of the $HI$ and $H2$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of $H2$ and total neutral H in late-type galaxies of $∼\!1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting $H2$ into stars and in transforming $HI$ gas into $H2$ . Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas-stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.