  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 non-linear 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 non-Gaussianity. 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_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.
