gaussRand

[gottfriedleibniz]

Per #1115 (comment), the implementation looks to be multiplying by variance instead of stddev (source).

Since this method computes two random values at a time, I am wondering whether this library should also handle this inefficiency in some thread-safe or user-defined way. For example, microsoft/STL includes a defaulted "keep" argument.