Adapt c-function to directly output pairwise terms

[FabianRathke]

Currently the c-implementation of sum-product only returns singleton terms of $q_c$. When calculating the terms of the objective function, a slow implementation of Matlab is required that explicitly relies on the $\Omega$-Terms. Adapting the c-function could speed up this calculation significantly.