CppTransport solves the 2- and 3-point functions of the perturbations produced during an inflationary epoch in the very early universe. It is implemented for models with canonical kinetic terms, although the underlying method is quite general and could be scaled to handle models with a non-trivial field-space metric or an even more general non-canonical Lagrangian.
http://adsabs.harvard.edu/abs/2016JCAP...12..033D http://adsabs.harvard.edu/abs/2016arXiv160900380S