Implement evalCPU for p_rs, write function stub for dbeta on which it depends.
This method is included for completeness - it's not actually needed for the
sampling step due to the conjugate beta prior on these terms.