#include <random.h>

Collaboration diagram for Athena_test::normal_distribution< T >:

Public Member Functions
	normal_distribution (T mu, T s)

T	operator() (URNG &rng) const

Private Attributes
T	m_mu

T	m_s

Detailed Description

template<class T>
class Athena_test::normal_distribution< T >

Definition at line 101 of file random.h.

Constructor & Destructor Documentation

◆ normal_distribution()

template<class T >

Athena_test::normal_distribution< T >::normal_distribution	(	T	mu,
		T	s
	)

inline

Definition at line 104 of file random.h.

104 : m_mu(mu), m_s(s) {}

Member Function Documentation

◆ operator()()

template<class T >

T Athena_test::normal_distribution< T >::operator() ( URNG & rng ) const

Definition at line 114 of file random.h.

 {
   // Stamdard Box-Muller transform a la Numerical Recipes, except that for simplicity
   // i don't try to remember the second generated number between calls.  I don't
   // really care about performance here.
   T v1, v2, r2;
   do {
     v1 = 2*(static_cast<T>(rng()) / URNG::max())-1;
     v2 = 2*(static_cast<T>(rng()) / URNG::max())-1;
     r2 = v1*v1 + v2*v2;
   } while (r2 >= 1);
  
   T fac = sqrt(-2*log(r2) / r2);
   return v1*fac*m_s + m_mu;
 }