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