Abstract. We develop a new, efficient code for solving the second-order Einstein-Boltzmann equations, and use it to estimate the intrinsic CMB non-Gaussianity arising from the nonlinear evolution of density perturbations. The full calculation involves contributions from recombination and less tractable contributions from terms integrated along the line of sight. We investigate the bias that this intrinsic bispectrum implies for searches of primordial nonGaussianity. We find that the inclusion or omission of certain line of sight terms can make a large impact. When including all physical effects but lensing and time-delay, we find that the local-type f NL would be biased by f intr NL = 0.5, below the expected sensitivity of the Planck satellite. The speed of our code allows us to confirm the robustness of our results with respect to a number of numerical parameters.