add .cpp for calculating posterior Ps-Pu and log(Ps)-log(Pu)