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;
  }
}