Abstract
We study two key issues militating against the use of the anisotropic three-point correlation function (3PCF) for cosmological parameter inference: difficulties with its computational estimation and high-dimensionality. We show how high-dimensionality may be reduced significantly by multipole decompositions of all angular dependence. This allows deriving the full expression for the multipole moments of the anisotropic 3PCF and its covariance matrix in a basis where the dimensionality reduces from nine to two at each multipole in the plane-parallel limit. We use 2D FFTLog formalism to show how the multipole moments with double momentum integrals over the product of bispectrum and two highly oscillating spherical Bessel functions and its covariance with double momentum integrals over the product of three galaxy power spectra and a combination of four highly oscillating spherical Bessel functions may be computed optimally.