1324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver/* 2324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * [The "BSD licence"] 3324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Copyright (c) 2005-2008 Terence Parr 4324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * All rights reserved. 5324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 6324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Conversion to C#: 7324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Copyright (c) 2008-2009 Sam Harwell, Pixel Mine, Inc. 8324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * All rights reserved. 9324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 10324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Redistribution and use in source and binary forms, with or without 11324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * modification, are permitted provided that the following conditions 12324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * are met: 13324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 1. Redistributions of source code must retain the above copyright 14324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * notice, this list of conditions and the following disclaimer. 15324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 2. Redistributions in binary form must reproduce the above copyright 16324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * notice, this list of conditions and the following disclaimer in the 17324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * documentation and/or other materials provided with the distribution. 18324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 3. The name of the author may not be used to endorse or promote products 19324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * derived from this software without specific prior written permission. 20324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 21324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR 22324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 23324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. 24324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, 25324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT 26324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 27324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY 28324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 29324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 30324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 31324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver */ 32324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 33324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruvernamespace Antlr.Runtime.Misc { 34324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver using System.Collections.Generic; 35324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver using Antlr.Runtime.JavaExtensions; 36324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 37324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver using Math = System.Math; 38324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 39324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver /** <summary>Stats routines needed by profiler etc...</summary> 40324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 41324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * <remarks> 42324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * note that these routines return 0.0 if no values exist in the X[] 43324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * which is not "correct", but it is useful so I don't generate NaN 44324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * in my output 45324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * </remarks> 46324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver */ 47324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public class Stats { 48324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public const string ANTLRWORKS_DIR = "antlrworks"; 49324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 50324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver /** <summary>Compute the sample (unbiased estimator) standard deviation following:</summary> 51324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 52324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * <remarks> 53324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Computing Deviations: Standard Accuracy 54324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Tony F. Chan and John Gregg Lewis 55324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Stanford University 56324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * Communications of ACM September 1979 of Volume 22 the ACM Number 9 57324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * 58324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * The "two-pass" method from the paper; supposed to have better 59324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * numerical properties than the textbook summation/sqrt. To me 60324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * this looks like the textbook method, but I ain't no numerical 61324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * methods guy. 62324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver * </remarks> 63324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver */ 64324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static double Stddev(int[] X) { 65324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver int m = X.Length; 66324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver if (m <= 1) { 67324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return 0; 68324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 69324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver double xbar = EnumerableExtensions.Average(X); 70324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver double s2 = 0.0; 71324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver for (int i = 0; i < m; i++) { 72324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver s2 += (X[i] - xbar) * (X[i] - xbar); 73324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 74324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver s2 = s2 / (m - 1); 75324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return Math.Sqrt(s2); 76324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 77324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static double Stddev(List<int> X) { 78324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver int m = X.Count; 79324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver if (m <= 1) { 80324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return 0; 81324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 82324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver double xbar = EnumerableExtensions.Average(X); 83324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver double s2 = 0.0; 84324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver for (int i = 0; i < m; i++) { 85324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver s2 += (X[i] - xbar) * (X[i] - xbar); 86324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 87324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver s2 = s2 / (m - 1); 88324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return Math.Sqrt(s2); 89324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 90324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 91324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver#if DEBUG 92324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver /** <summary>Compute the sample mean</summary> */ 93324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver [System.Obsolete] 94324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static double avg(int[] X) { 95324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return EnumerableExtensions.Average(EnumerableExtensions.DefaultIfEmpty(X, 0)); 96324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 97324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 98324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver [System.Obsolete] 99324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static int min(int[] X) { 100324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return EnumerableExtensions.Min(EnumerableExtensions.DefaultIfEmpty(X, int.MaxValue)); 101324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 102324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 103324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver [System.Obsolete] 104324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static int max(int[] X) { 105324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return EnumerableExtensions.Max(EnumerableExtensions.DefaultIfEmpty(X, int.MinValue)); 106324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 107324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 108324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver [System.Obsolete] 109324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static int sum(int[] X) { 110324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return EnumerableExtensions.Sum(X); 111324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 112324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver#endif 113324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 114324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static void WriteReport(string filename, string data) { 115324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver string absoluteFilename = GetAbsoluteFileName(filename); 116324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 117324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver System.IO.Directory.CreateDirectory(System.IO.Path.GetDirectoryName(absoluteFilename)); 118324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver System.IO.File.AppendAllText(absoluteFilename, data); 119324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 120324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver 121324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver public static string GetAbsoluteFileName(string filename) { 122324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver string personalFolder = System.Environment.GetFolderPath(System.Environment.SpecialFolder.Personal); 123324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver return personalFolder + System.IO.Path.DirectorySeparatorChar + 124324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver ANTLRWORKS_DIR + System.IO.Path.DirectorySeparatorChar + 125324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver filename; 126324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 127324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver } 128324c4644fee44b9898524c09511bd33c3f12e2dfBen Gruver} 129