/// /// This file is part of ILNumerics Community Edition. /// /// ILNumerics Community Edition - high performance computing for applications. /// Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net /// /// ILNumerics Community Edition is free software: you can redistribute it and/or modify /// it under the terms of the GNU General Public License version 3 as published by /// the Free Software Foundation. /// /// ILNumerics Community Edition is distributed in the hope that it will be useful, /// but WITHOUT ANY WARRANTY; without even the implied warranty of /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the /// GNU General Public License for more details. /// /// You should have received a copy of the GNU General Public License /// along with ILNumerics Community Edition. See the file License.txt in the root /// of your distribution package. If not, see . /// /// In addition this software uses the following components and/or licenses: /// /// ================================================================================= /// The Open Toolkit Library License /// /// Copyright (c) 2006 - 2009 the Open Toolkit library. /// /// Permission is hereby granted, free of charge, to any person obtaining a copy /// of this software and associated documentation files (the "Software"), to deal /// in the Software without restriction, including without limitation the rights to /// use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of /// the Software, and to permit persons to whom the Software is furnished to do /// so, subject to the following conditions: /// /// The above copyright notice and this permission notice shall be included in all /// copies or substantial portions of the Software. /// /// ================================================================================= /// using System; using System.Collections.Generic; using System.Text; using System.Runtime.InteropServices; using ILNumerics.Storage; using ILNumerics.Misc; using ILNumerics.Native; using ILNumerics.Exceptions; namespace ILNumerics { public partial class ILMath { [ThreadStatic] private static ILNRandom m_nrandomGenerator; internal static ILNRandom NRandomGenerator { get { if (m_nrandomGenerator == null) m_nrandomGenerator = new ILNRandom(Environment.TickCount + System.Threading.Thread.CurrentThread.ManagedThreadId); return m_nrandomGenerator; } } /// /// Normal randomly-distributed N-dimensional array elements. /// /// int array or single int paramters specifying /// dimensions for new array to be created. /// N-dimensional array filled with random numbers. /// The elements lie within the range 0.0 ... 1.0 and are choosen to be normally /// distributed. public static ILRetArray randn(params int[] dimensions) { ILSize dim; if (dimensions.Length == 1) dim = new ILSize(dimensions[0], dimensions[0]); else dim = new ILSize(dimensions); double[] data = ILMemoryPool.Pool.New(dim.NumberOfElements); if (dim.NumberOfElements < ILNumerics.Settings.s_minParallelElement1Count) { unsafe { fixed (double* pRetArray = data) { double* pCurIdx = pRetArray; double* pLastElement = pCurIdx + dim.NumberOfElements; while (pCurIdx < pLastElement) *pCurIdx++ = NRandomGenerator.NextDouble(); } } } else { System.Threading.Tasks.Parallel.For(0, dim.NumberOfElements, i => { data[i] = NRandomGenerator.NextDouble(); }); } return new ILRetArray(data, dim); } } }