  Let $(\mathbb{X},d,\mu)$ be a space of homogeneous type in the sense of R. R.
Coifman and G. Weiss, and $X(\mathbb{X})$ a ball quasi-Banach function space on
$\mathbb{X}$. In this article, the authors introduce the weak Hardy space
$WH_X(\mathbb{X})$ associated with $X(\mathbb{X})$ via the grand maximal
function, and characterize $WH_X(\mathbb{X})$ by other maximal functions and
atoms. The authors then apply these characterizations to obtain the real
interpolation and the boundedness of Calder\'{o}n--Zygmund operators in the
critical case. The main novelties of this article exist in that the authors use
the Aoki--Rolewicz theorem and both the dyadic system and the exponential decay
of approximations of the identity on $\mathbb{X}$, which closely connect with
the geometrical properties of $\mathbb{X}$, to overcome the difficulties caused
by the absence of both the triangle inequality of $\|\cdot\|_{X(\mathbb{X})}$
and the reverse doubling assumption of the measure $\mu$ under consideration,
and also use the relation between the convexification of $X(\mathbb{X})$ and
the weak space $WX(\mathbb{X})$ associated with $X(\mathbb{X})$ to prove that
the infinite summation of atoms converges in the space of distributions on
$\mathbb{X}$. Moreover, all these results have a wide range of generality and,
particularly, even when they are applied to the weighted Lebesgue space, the
Orlicz space, and the variable Lebesgue space, the obtained results are also
new and, actually, some of them are new even on RD-spaces (namely, spaces of
homogeneous type satisfying the additional reverse doubling condition).
