We extend the public CPPTRANSPORT code to calculate the statistical properties of fluctuations in multiple-field inflationary models with curved field space. Our implementation accounts for all physical effects at tree-level in the `in-in' diagrammatic expansion. This includes particle production due to time-varying masses, but excludes scenarios where the curvature perturbation is generated by averaging over the decay of more than one particle. We test our implementation by comparing results in Cartesian and polar field-space coordinates, showing excellent numerical agreement and only minor degradation in compute time. We compare our results with the PYTRANSPORT 2.0 code, which uses the same computational approach but a different numerical implementation, finding good agreement. Finally, we use our tools to study a class of gelaton-like models which could produce an enhanced non-Gaussian signal on equilateral configurations of the Fourier bispectrum. We show this is difficult to achieve using hyperbolic field-space manifolds and simple inflationary potentials.