We study the flow forced by precession in rigid non-axisymmetric ellipsoidal containers. To do so, we revisit the inviscid and viscous analytical models which have been previously developed for the spheroidal geometry by respectively Poincaré [Bull. Astron. 27, 321 (1910)] and Busse [J. Fluid Mech. 33, 739 (1968)], and, we report the first numerical simulations of flows in such a geometry. In strong contrast with axisymmetric spheroids where the forced flow is systematically stationary in the precessing frame, we show that the forced flow is unsteady and periodic. Comparisons of the numerical simulations with the proposed theoretical model show excellent agreement for both axisymmetric and nonaxisymmetric containers. Finally, since the studied configuration corresponds to a tidally locked celestial body such as the Earth's Moon, we use our model to investigate the challenging but planetary relevant limit of very small Ekman numbers and the particular case of our Moon.