/// <summary>
/// Normal randomly-distributed N-dimensional array elements.
/// </summary>
/// <param name="dimensions">int array or single int paramters specifying
/// dimensions for new array to be created.</param>
/// <returns>N-dimensional array filled with random numbers.</returns>
/// <remarks>The elements lie within the range 0.0 ... 1.0 and are choosen to be normally
/// distributed.</remarks>
public static ILArray<double> randn(params int[] dimensions) {
ILDimension dim;
if (dimensions.Length == 1)
dim = new ILDimension(dimensions[0],dimensions[0]);
else
dim = new ILDimension(dimensions);
if (m_nrandomGenerator == null)
m_nrandomGenerator = new ILNRandom(Environment.TickCount);
double[] data = ILMemoryPool.Pool.New<double>(dim.NumberOfElements);
unsafe {
fixed (double* pRetArray = data) {
double* pCurIdx = pRetArray;
double* pLastElement = pCurIdx + dim.NumberOfElements;
while (pCurIdx < pLastElement)
*pCurIdx++ = m_nrandomGenerator.NextDouble();
}
}
return new ILArray<double>(data, dim);
}
/// <summary>
/// Set seed to both rand and randn functions
/// </summary>
/// <param name="seed">Any ole' number will do</param>
public static void setseed(int seed)
{
m_nrandomGenerator = new ILNRandom(seed);
m_randomGenerator = new Random(seed);
}