September 23 2004 10:39:00.472 AM
HALTON_DATASET (FORTRAN90 version)
Generate a Halton dataset.
This program is meant to be used interactively.
It is also possible to prepare a simple input
file beforehand and use it in batch mode.
The program requests input values from the user:
* NDIM, the spatial dimension,
* N, the number of points to generate,
* STEP, the index of the first subsequence element.
* SEED(1:NDIM), the Halton sequence element
corresponding to STEP = 0
* LEAP(1:NDIM), the successive jumps in the
Halton sequence.
* BASE(1:NDIM), the Halton bases,
usually distinct primes.
The program generates the data, writes it to the file
halton_NDIM_N.txt
where "NDIM" and "N" are the numeric values specified
by the user, and then asks the user for more input.
To indicate that no further computations are
desired, it is enough to input a nonsensical value,
such as -1.
*
*
* Ready to generate a new dataset:
*
*
Enter NDIM, the spatial dimension:
(Try "2" if you do not have a preference.)
(0 or any negative value terminates execution).
User input NDIM = 2
Enter N, the number of points to generate:
(Try "25" if you do not have a preference.)
(0 or any negative value terminates execution).
User input N = 10
Enter STEP, the index of the first subsequence element:
(Try "0" or "1" if you do not have a preference.)
(Any negative value terminates execution).
User input STEP = 0
Enter SEED(1:NDIM), the Halton sequence index
corresponding to STEP = 0
(Try "0 0 ... 0" if you do not have a preference).
(a negative value terminates execution.)
User input 1 1
Enter LEAP(1:NDIM), the successive jumps in
the Halton sequence.
(Try "1 1 ... 1" if you do not have a preference).
(another choice is any prime larger than all bases.)
(any value less than 1 terminates execution.)
User input 1 1
Enter BASE(1:NDIM), the Halton bases,
usually distinct primes.
(Try "2 3 5 7 11 13 ..." if you do not have a preference).
(any value 1 or less terminates execution.)
User input 2 3
The data was written to the file "halton_02_00010.txt".
*
*
* Ready to generate a new dataset:
*
*
Enter NDIM, the spatial dimension:
(Try "2" if you do not have a preference.)
(0 or any negative value terminates execution).
User input NDIM = 3
Enter N, the number of points to generate:
(Try "25" if you do not have a preference.)
(0 or any negative value terminates execution).
User input N = 100
Enter STEP, the index of the first subsequence element:
(Try "0" or "1" if you do not have a preference.)
(Any negative value terminates execution).
User input STEP = 0
Enter SEED(1:NDIM), the Halton sequence index
corresponding to STEP = 0
(Try "0 0 ... 0" if you do not have a preference).
(a negative value terminates execution.)
User input 0 1 2
Enter LEAP(1:NDIM), the successive jumps in
the Halton sequence.
(Try "1 1 ... 1" if you do not have a preference).
(another choice is any prime larger than all bases.)
(any value less than 1 terminates execution.)
User input 1 2 3
Enter BASE(1:NDIM), the Halton bases,
usually distinct primes.
(Try "2 3 5 7 11 13 ..." if you do not have a preference).
(any value 1 or less terminates execution.)
User input 5 3 7
The data was written to the file "halton_03_00100.txt".
*
*
* Ready to generate a new dataset:
*
*
Enter NDIM, the spatial dimension:
(Try "2" if you do not have a preference.)
(0 or any negative value terminates execution).
User input NDIM = 7
Enter N, the number of points to generate:
(Try "25" if you do not have a preference.)
(0 or any negative value terminates execution).
User input N = 1000
Enter STEP, the index of the first subsequence element:
(Try "0" or "1" if you do not have a preference.)
(Any negative value terminates execution).
User input STEP = 0
Enter SEED(1:NDIM), the Halton sequence index
corresponding to STEP = 0
(Try "0 0 ... 0" if you do not have a preference).
(a negative value terminates execution.)
User input 0 0 0 0 0
0 0
Enter LEAP(1:NDIM), the successive jumps in
the Halton sequence.
(Try "1 1 ... 1" if you do not have a preference).
(another choice is any prime larger than all bases.)
(any value less than 1 terminates execution.)
User input 1 1 1 1 1
1 1
Enter BASE(1:NDIM), the Halton bases,
usually distinct primes.
(Try "2 3 5 7 11 13 ..." if you do not have a preference).
(any value 1 or less terminates execution.)
User input 2 3 4 5 6
7 8
The data was written to the file "halton_07_01000.txt".
*
*
* Ready to generate a new dataset:
*
*
Enter NDIM, the spatial dimension:
(Try "2" if you do not have a preference.)
(0 or any negative value terminates execution).
User input NDIM = -1
HALHAM_NDIM_CHECK - Fatal error!
NDIM < 1.
NDIM = -1
HALTON_DATASET
The input value of NDIM = -1
is interpreted as a request for termination.
Normal end of execution.
September 23 2004 10:39:00.537 AM