sink_cpp function

A variational approximation of an expected permutation matrix using C++