`std::pow()` is very slow: extract `2.0 * beta * RcpPi` into a variable and use two multiplications instead?