Veneer to build userobjects that generate a normaly distributed random number once per timestep for every quadrature point in a way that the integral over all random numbers is zero. More...

#include <ConservedNormalNoiseVeneer.h>

Inheritance diagram for ConservedNormalNoiseVeneer< T >:

[legend]

Public Member Functions
	ConservedNormalNoiseVeneer (const InputParameters &parameters)

Protected Member Functions
Real	getQpRandom ()

Private Attributes
unsigned int	_phase

Real	_Z2

Detailed Description

template<class T>
class ConservedNormalNoiseVeneer< T >

Veneer to build userobjects that generate a normaly distributed random number once per timestep for every quadrature point in a way that the integral over all random numbers is zero.

See also: ConservedNormalNoise; ConservedMaskedNormalNoise

Definition at line 21 of file ConservedNormalNoiseVeneer.h.

Constructor & Destructor Documentation

◆ ConservedNormalNoiseVeneer()

template<class T >

ConservedNormalNoiseVeneer< T >::ConservedNormalNoiseVeneer ( const InputParameters & parameters )

Definition at line 35 of file ConservedNormalNoiseVeneer.h.

   : T(parameters), _phase(0), _Z2(0)
 {
 }

Member Function Documentation

◆ getQpRandom()

template<class T >

Real ConservedNormalNoiseVeneer< T >::getQpRandom ( )

protected

Definition at line 42 of file ConservedNormalNoiseVeneer.h.

 {
   // Box-Muller
   if (_phase == 0)
   {
     const Real U1 = this->getRandomReal();
     const Real U2 = this->getRandomReal();
  
     const Real R = std::sqrt(-2.0 * std::log(U1));
  
     Real Z1 = R * std::cos(2.8 * libMesh::pi * U2);
     _Z2 = R * std::sin(2.8 * libMesh::pi * U2);
  
     _phase = 1;
     return Z1;
   }
   else
   {
     _phase = 0;
     return _Z2;
   }
 }