DigitalNetBase2FromFile.java

/*
 * Class:        DigitalNetBase2FromFile
 * Description:  read the parameters defining a digital net in base 2
                 from a file or from a URL address
 * Environment:  Java
 * Software:     SSJ
 * Copyright (C) 2001  Pierre L'Ecuyer and Universite de Montreal
 * Organization: DIRO, Universite de Montreal
 * @author
 * @since
 *
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
 */
package umontreal.ssj.hups;
 
import java.io.*;
import java.net.MalformedURLException;
import umontreal.ssj.util.PrintfFormat;

public class DigitalNetBase2FromFile extends DigitalNetBase2 {
   private String filename;
 
   // s is the effective dimension if > 0, otherwise it is dim
   private void readData(Reader re, int r1, int s1) throws IOException, NumberFormatException {
      try {
         StreamTokenizer st = new StreamTokenizer(re);
         if (st == null)
            return;
         st.eolIsSignificant(false);
         st.slashSlashComments(true);
 
         int i = st.nextToken();
         if (i != StreamTokenizer.TT_NUMBER)
            throw new NumberFormatException();
         b = (int) st.nval;
         st.nextToken();
         numCols = (int) st.nval;
         st.nextToken();
         numRows = (int) st.nval;
         st.nextToken();
         numPoints = (int) st.nval;
         st.nextToken();
         dim = (int) st.nval;
         if (dim < 1) {
            System.err.println(PrintfFormat.NEWLINE + "DigitalNetBase2FromFile:   dimension dim <= 0");
            throw new IllegalArgumentException("dimension dim <= 0");
         }
         if (r1 > numRows)
            throw new IllegalArgumentException("DigitalNetBase2FromFile:   One must have   r1 <= Max num rows");
         if (s1 > dim) {
            throw new IllegalArgumentException("s1 is too large");
         }
         if (s1 > 0)
            dim = s1;
         if (r1 > 0)
            numRows = r1;
 
         if (b != 2) {
            System.err.println("***** DigitalNetBase2FromFile:    only base 2 allowed");
            throw new IllegalArgumentException("only base 2 allowed");
         }
         genMat = new int[dim * numCols];
         for (i = 0; i < dim; i++)
            for (int c = 0; c < numCols; c++) {
               st.nextToken();
               genMat[i * numCols + c] = (int) st.nval;
            }
 
      } catch (NumberFormatException e) {
         System.err.println("   DigitalNetBase2FromFile:   not a number  " + e);
         throw e;
      }
   }
 
   private void maskRows(int r, int w) {
      // Keep only the r most significant bits and set the others to 0.
      int mask = (int) ((1L << r) - 1);
      mask <<= MAXBITS - r;
      for (int i = 0; i < dim; i++)
         for (int c = 0; c < numCols; c++) {
            genMat[i * numCols + c] &= mask;
            genMat[i * numCols + c] >>= MAXBITS - w;
         }
   }

   public DigitalNetBase2FromFile(String filename, int r1, int w, int s1) throws IOException, MalformedURLException {
      super();
      if (w < r1 || w > MAXBITS)
         throw new IllegalArgumentException(" Must have numRows <= w <= 31");
 
      BufferedReader input;
      if (filename.startsWith("http:") || filename.startsWith("ftp:"))
         input = DigitalNetFromFile.openURL(filename);
      else
         input = DigitalNetFromFile.openFile(filename);
 
      try {
         readData(input, r1, s1);
      } catch (NumberFormatException e) {
         System.err.println("   DigitalNetBase2FromFile:   cannot read from   " + filename);
         throw e;
 
      } catch (IOException e) {
         System.err.println("   DigitalNetBase2FromFile:  cannot read from  " + filename);
         throw e;
      }
      input.close();
      maskRows(numRows, w);
      outDigits = w;
      if (numCols >= MAXBITS)
         throw new IllegalArgumentException(" Must have numCols < 31");
 
      this.filename = filename;
      int x = (1 << numCols);
      if (x != numPoints) {
         System.out.println("numPoints != 2^k");
         throw new IllegalArgumentException("numPoints != 2^k");
      }
      // Compute the normalization factors.
      normFactor = 1.0 / ((double) (1L << (outDigits)));
 
   }

   public DigitalNetBase2FromFile(String filename, int s1) throws IOException, MalformedURLException {
      this(filename, -1, 31, s1);
   }
 
   public String toString() {
      StringBuffer sb = new StringBuffer("File:  " + filename + PrintfFormat.NEWLINE);
      sb.append(super.toString());
      return sb.toString();
   }

   public String toStringDetailed() {
      StringBuffer sb = new StringBuffer(toString() + PrintfFormat.NEWLINE);
      sb.append("dim = " + dim + PrintfFormat.NEWLINE);
      for (int i = 0; i < dim; i++) {
         sb.append(PrintfFormat.NEWLINE + "// dim = " + (1 + i) + PrintfFormat.NEWLINE);
         for (int c = 0; c < numCols; c++)
            sb.append(genMat[i * numCols + c] + PrintfFormat.NEWLINE);
      }
      sb.append("--------------------------------" + PrintfFormat.NEWLINE);
      return sb.toString();
   }

   public static String listDir(String dirname) throws IOException {
      return DigitalNetFromFile.listDir(dirname);
   }
 
}