The BI-spectra and Non-Gaussianity Operator (BINGO) code, written in Fortran, computes the scalar bi-spectrum and the non-Gaussianity parameter fNL in single field inflationary models involving the canonical scalar field. BINGO can calculate all the different contributions to the bi-spectrum and the parameter fNL for an arbitrary triangular configuration of the wavevectors.
https://ui.adsabs.harvard.edu/abs/2013JCAP...05..026H and https://ui.adsabs.harvard.edu/abs/2015JCAP...02..029S