We are using numpy.random_intel with the Intel Anaconda distribution. We need to avoid all uses of Mersenne Twister (which frighteningly is included in random_intel) and rely exclusively on the RDRAND instruction as the randomness source. We then want a random float from the double-tailed geometric distribution.

Unfortunately, the documentation, while including lots of information about distributions, is opaque for how to configure numpy to only use RDRAND as the randomness source. Can you provide me with this information and update the documentation?