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