I am sure this question has been answered before but I could not find any similar question.
In Diffie-Hellman protocol implementations we can try to find a large primes using safe prime formula (i.e. $p = 2k + 1$ such that $k$ is a prime). So computer generates a relatively large prime $k$ and from that we get a prime (or pseudo-prime) $p$. To find a primitive root quickly, we generate a random number $x$ and test if $x, x^2, x^k \not\equiv 1 \bmod p$ if so then it is a primitive root (because the order of $\mathbb{Z}^*_p$ should divide $\phi(p)$).
My question is what is the success probability to generate a primitive root as there are $\phi(\phi(p))$ primitive roots modulo $p$. Thank you so much in advance.