BigDiscrepancy.java

/*
 * Class:        BigDiscrepancy
 * Description:  Base class of all discrepancy classes programmed with
                 multi-precision floating-point numbers
 * Environment:  Java
 * Software:     SSJ
 * Copyright (C) 2001  Pierre L'Ecuyer and Universite de Montreal
 * Organization: DIRO, Universite de Montreal
 * @author       Richard Simard
 * @since        January 2009
 
 * SSJ is free software: you can redistribute it and/or modify it under
 * the terms of the GNU General Public License (GPL) as published by the
 * Free Software Foundation, either version 3 of the License, or
 * any later version.
 
 * SSJ 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.
 
 * A copy of the GNU General Public License is available at
   <a href="http://www.gnu.org/licenses">GPL licence site</a>.
 */
package umontreal.ssj.discrepancy;
 
import umontreal.ssj.hups.PointSet;
import java.math.*;

public abstract class BigDiscrepancy extends Discrepancy {
 
   protected BigDecimal[] C1Big; // functions of gamma[i]
   protected BigDecimal[] C2Big; // functions of gamma[i]
   protected BigDecimal[] C3Big; // functions of gamma[i]
   protected BigDecimal[] UBig; // i/n
   protected BigDecimal[][] FactorBig;
 
   protected void reserveCBig(int s) {
      // Reserve memory for gamma factors
      C1Big = new BigDecimal[s];
      C2Big = new BigDecimal[s];
      C3Big = new BigDecimal[s];
 
      for (int i = 0; i < s; i++) {
         C1Big[i] = new BigDecimal(0);
         C2Big[i] = new BigDecimal(0);
         C3Big[i] = new BigDecimal(0);
      }
   }
 
   protected void reserveFactorBig(int n, int s) {
      // Reserve memory for factors
      FactorBig = new BigDecimal[n][s];
      for (int i = 0; i < n; i++) {
         FactorBig[i] = new BigDecimal[s];
         for (int j = 0; j < s; j++)
            FactorBig[i][j] = new BigDecimal(0);
      }
   }
 
   protected void setUBig(int n) {
      // Precompute all U[i] = i/n, with i = 0, 1, 2, ..., n-1
      UBig = new BigDecimal[n];
      BigDecimal Ninv = new BigDecimal(1); // = 1/n
      Ninv = Ninv.divide(new BigDecimal(n), MathContext.DECIMAL128);
 
      for (int i = 0; i < n; i++) {
         UBig[i] = new BigDecimal(i);
         UBig[i] = UBig[i].multiply(Ninv);
      }
   }

   public BigDiscrepancy(double[][] points, int n, int s) {
      super(points, n, s);
   }

   public BigDiscrepancy(double[][] points, int n, int s, double[] gamma) {
      super(points, n, s, gamma);
   }

   public BigDiscrepancy(int n, int s, double[] gamma) {
      super(n, s, gamma);
   }

   public BigDiscrepancy(PointSet set) {
      super(set);
   }

   public BigDiscrepancy() {
   }
 
}