randnum/Normal.cpp File Reference

#include "Normal.h"
#include "randnum.h"
#include "../utility/numutil.h"
#include <cmath>
#include <iostream>

Include dependency graph for Normal.cpp:

Functions
unsigned long	genrand_int32 (void)
double	dummy ()
double	algorithmF ()
Variables
const unsigned long	y []
const long	a []
const unsigned long	q []
const double	c = 0.004996971959878404
const double	d = 1.861970434352886050

---

Function Documentation

double algorithmF

(

)

See Knuth, TAOCP Vol 2 Sec 3.4.1 Algorithm M: The Rectangle-wedge-tail method discovered by G Marsaglia TODO: implement it by filling the gaps TODO: Not yet implemented.

double dummy

(

)

Referenced by Cinfo::getDestFinfo(), Cinfo::getFieldElementFinfo(), Cinfo::getLookupFinfo(), Ksolve::getNvec(), Gsolve::getNvec(), Stoich::getProxyPools(), Cinfo::getSharedFinfo(), Cinfo::getValueFinfo(), Function::getVar(), NeuroMesh::innerBuildDefaultMesh(), NeuroMesh::insertSingleDummy(), and Cinfo::makeCinfoElements().

unsigned long genrand_int32

(

void

)

---

Variable Documentation

const long a[]

Initial value:

 
{
    31,     30,     29,     28,     27,     26,     25,     24,     23,     22,     21,     20,     27,     26,     25,     24,     23,     22,     21,     -1,     19,     19,     21,     20,     19,     18,     17,     16,     15,     14,     13,     12,     19,     18,     17,     16,     15,     14,     13,     12,     11,     10,     9,     8,     7,     6,     5,     4,     3,     2,     1,     0,     23,     22,     21,     20,     19,     18,     17,     16,     15,     14,     13,     12,     11,     10,     9,     8,     7,     6,     5,     4,     3,     2,     1,     0,     51,     50,     49,     48,     47,     46,     45,     44,     43,     42,     41,     40,     39,     38,     37,     36,     35,     34,     33,     32,     31,     30,     29,     28,     27,     26,     25,     24,     23,     22,     21,     20,     19,     18,     17,     16,     15,     14,     13,     12,     11,     10,     9,     8,     7,     6,     5,     4,     3,     2,     1,     0       
}

Referenced by NeuroMesh::indexToSpace(), IzhikevichNrn::initCinfo(), LookupTable::lookup(), SpineMesh::nearest(), PsdMesh::nearest(), CylBase::nearest(), testIsA(), and testSetGetExtField().

const double c = 0.004996971959878404

Referenced by estT(), Gsolve::fillMathDep(), fillPointsOnCircle(), fix(), IzhikevichNrn::initCinfo(), CylBase::nearest(), printSparseMatrix(), Clock::reportClock(), testHopFunc(), testSetGetExtField(), testSparseMatrix(), and testSyncSynapseSize().

const double d = 1.861970434352886050

Referenced by SrcFinfo::addMsg(), SrcFinfo::checkTarget(), Shell::initCinfo(), IzhikevichNrn::initCinfo(), SpineMesh::nearest(), PsdMesh::nearest(), and testSetGetSynapse().

const unsigned long q[]

Initial value:

 
{
    28,     32,   33,   36,   40,   43,   45,   47,   51,   53,   57,   59,   27,   37,   43,   54,   28,   45,   60,   63,   49,   61,   52,   34,   63,   46,   34,   18,   47,   40,   36,   29,   61,   57,   53,   51,   47,   43,   40,   37,   33,   31,   27,   25,   21,   19,   17,   14,   11,   9 ,   8,   5,   54,   51,   49,   46,   44,   41,   39,   37,   35,   33,   31,   29,   28,   26,   24,   23,   21,   20,   19,   18,   16,   15,   14,   13,   12,   12,   11,   10,   9,   9,   8,   7,   7,   6,   6,   5,   5,   5,   4,   4,   4,   3,   3,   3,   3,   3,   2,   2,   2,   2,   2,   2,   2,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1   
}

Referenced by CubeMesh::buildStencil(), CubeMesh::fillSpaceToMeshLookup(), GetHopFunc< A >::getLocalFieldVec(), VoxelPools::gslFunc(), HopFunc1< A >::localFieldOpVec(), HopFunc1< A >::localOpVec(), HopFunc2< A1, A2 >::opVec(), and SparseMatrix< unsigned int >::reorderColumns().

const unsigned long y[]

Initial value:

 
{
    200,     199,     199,     198,     197,     196,     195,     193,     192,     190,     188,     186,     184,     181,     179,     176,     173,     170,     167,     164,     161,     157,     154,     151,     147,     143,     140,     136,     132,     128,     125,     121,     117,     113,     110,     106,     102,     98,     95,     91,     88,     84,     81,     77,     74,     71,     68,     64,     61,     59,     56,     53,     50,     48,     45,     43,     40,     38,     36,     34,     32,     30,     28,     27,     25,     23,     22,     20,     19,     18,     17,     15,     14,     13,     12,     11,     11,     10,     9,     8,     8,     7,     6,     6,     5,     5,     4,     4,     4,     3,     3,     3,     2,     2,     2,     2,     2,     1,     1,     1,     1,     1,     1,     1,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0,     0
}

Refer to: Eine Alias-Methode zur Stichprohenentnahme aus Normalverteilungen. JH Ahrens and U Dieter, 19 89, Computing

We are assuming size of long to be 32 bit

Referenced by CubeMesh::fillSpaceToMeshLookup(), genrand_int32(), Annotator::initCinfo(), TableBase::initCinfo(), Interpol::initCinfo(), Function::initCinfo(), Func::initCinfo(), moose::CompartmentBase::initCinfo(), NeuroMesh::insertDummyNodes(), TableBase::interpolate(), NeuroNode::NeuroNode(), StochSecondOrderSingleSubstrate::operator()(), StochNOrder::operator()(), HopFunc2< A1, A2 >::opVec(), StimulusTable::process(), SteadyState::randomizeInitialCondition(), StimulusTable::reinit(), testWildcard(), writeGroup(), and writeKkit().