ARMA implementation of diffraction operators with inverse-root singularities