13447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein/*
23447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein [The "BSD license"]
33447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein Copyright (c) 2005-2009 Terence Parr
43447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein All rights reserved.
53447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
63447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein Redistribution and use in source and binary forms, with or without
73447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein modification, are permitted provided that the following conditions
83447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein are met:
93447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein 1. Redistributions of source code must retain the above copyright
103447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein     notice, this list of conditions and the following disclaimer.
113447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein 2. Redistributions in binary form must reproduce the above copyright
123447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein     notice, this list of conditions and the following disclaimer in the
133447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein     documentation and/or other materials provided with the distribution.
143447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein 3. The name of the author may not be used to endorse or promote products
153447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein     derived from this software without specific prior written permission.
163447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
173447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
183447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
193447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
203447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
213447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
223447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
233447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
243447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
253447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
263447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
273447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein */
283447a5916aa62f44de24cc441fc9987116ddff52Andrew Sappersteinpackage org.antlr.runtime.misc;
293447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
303447a5916aa62f44de24cc441fc9987116ddff52Andrew Sappersteinimport java.io.*;
313447a5916aa62f44de24cc441fc9987116ddff52Andrew Sappersteinimport java.util.List;
323447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
333447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein/** Stats routines needed by profiler etc...
343447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
353447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein // note that these routines return 0.0 if no values exist in the X[]
363447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein // which is not "correct", but it is useful so I don't generate NaN
373447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein // in my output
383447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
393447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein */
403447a5916aa62f44de24cc441fc9987116ddff52Andrew Sappersteinpublic class Stats {
413447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static final String ANTLRWORKS_DIR = "antlrworks";
423447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
433447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	/** Compute the sample (unbiased estimator) standard deviation following:
443447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *
453447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  Computing Deviations: Standard Accuracy
463447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  Tony F. Chan and John Gregg Lewis
473447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  Stanford University
483447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  Communications of ACM September 1979 of Volume 22 the ACM Number 9
493447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *
503447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  The "two-pass" method from the paper; supposed to have better
513447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  numerical properties than the textbook summation/sqrt.  To me
523447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  this looks like the textbook method, but I ain't no numerical
533447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 *  methods guy.
543447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	 */
553447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static double stddev(int[] X) {
563447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.length;
573447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m<=1 ) {
583447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
593447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
603447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		double xbar = avg(X);
613447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		double s2 = 0.0;
623447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
633447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			s2 += (X[i] - xbar)*(X[i] - xbar);
643447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
653447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		s2 = s2/(m-1);
663447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return Math.sqrt(s2);
673447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
683447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
693447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	/** Compute the sample mean */
703447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static double avg(int[] X) {
713447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		double xbar = 0.0;
723447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.length;
733447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
743447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
753447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
763447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
773447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			xbar += X[i];
783447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
793447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( xbar>=0.0 ) {
803447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return xbar / m;
813447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
823447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return 0.0;
833447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
843447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
853447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static int min(int[] X) {
863447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int min = Integer.MAX_VALUE;
873447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.length;
883447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
893447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
903447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
913447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
923447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			if ( X[i] < min ) {
933447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein				min = X[i];
943447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			}
953447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
963447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return min;
973447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
983447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
993447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static int max(int[] X) {
1003447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int max = Integer.MIN_VALUE;
1013447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.length;
1023447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
1033447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
1043447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1053447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
1063447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			if ( X[i] > max ) {
1073447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein				max = X[i];
1083447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			}
1093447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1103447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return max;
1113447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1123447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1133447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	/** Compute the sample mean */
1143447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static double avg(List<Integer> X) {
1153447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		double xbar = 0.0;
1163447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.size();
1173447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
1183447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
1193447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1203447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
1213447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			xbar += X.get(i);
1223447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1233447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( xbar>=0.0 ) {
1243447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return xbar / m;
1253447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1263447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return 0.0;
1273447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1283447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1293447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static int min(List<Integer> X) {
1303447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int min = Integer.MAX_VALUE;
1313447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.size();
1323447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
1333447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
1343447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1353447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
1363447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			if ( X.get(i) < min ) {
1373447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein				min = X.get(i);
1383447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			}
1393447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1403447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return min;
1413447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1423447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1433447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static int max(List<Integer> X) {
1443447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int max = Integer.MIN_VALUE;
1453447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.size();
1463447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
1473447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
1483447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1493447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
1503447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			if ( X.get(i) > max ) {
1513447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein				max = X.get(i);
1523447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			}
1533447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1543447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return max;
1553447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1563447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1573447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static int sum(int[] X) {
1583447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int s = 0;
1593447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		int m = X.length;
1603447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		if ( m==0 ) {
1613447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			return 0;
1623447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1633447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		for (int i=0; i<m; i++){
1643447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein			s += X[i];
1653447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		}
1663447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return s;
1673447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1683447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1693447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static void writeReport(String filename, String data) throws IOException {
1703447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		String absoluteFilename = getAbsoluteFileName(filename);
1713447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		File f = new File(absoluteFilename);
1723447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		File parent = f.getParentFile();
1733447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		parent.mkdirs(); // ensure parent dir exists
1743447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		// write file
1753447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		FileOutputStream fos = new FileOutputStream(f, true); // append
1763447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		BufferedOutputStream bos = new BufferedOutputStream(fos);
1773447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		PrintStream ps = new PrintStream(bos);
1783447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		ps.println(data);
1793447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		ps.close();
1803447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		bos.close();
1813447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		fos.close();
1823447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1833447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein
1843447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	public static String getAbsoluteFileName(String filename) {
1853447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein		return System.getProperty("user.home")+File.separator+
1863447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein					ANTLRWORKS_DIR +File.separator+
1873447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein					filename;
1883447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein	}
1893447a5916aa62f44de24cc441fc9987116ddff52Andrew Sapperstein}
190