September 23 2004 11:02:10.636 AM HAMMERSLEY_DATASET (FORTRAN90 version) Generate a Hammersley 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 sequence element corresponding to STEP = 0 * LEAP(1:NDIM), the successive jumps in the sequence. * BASE(1:NDIM), the bases, usually distinct primes or -N (to generate values like J/N). The program generates the data, writes it to the file hammersley_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 have no preference.) (0 or any negative value terminates execution). User input NDIM = 2 Enter N, the number of points to generate: (Try "25" if you have no preference.) (0 or any negative value terminates execution). User input N = 20 Enter STEP, the index of the first subsequence element: (Try "0" or "1" if you have no preference.) (Any negative value terminates execution). User input STEP = 0 Enter SEED(1:NDIM), the sequence index corresponding to STEP = 0 (Try "0 0 ... 0" if you have no preference.) (a negative value terminates execution.) User input 1 0 Enter LEAP(1:NDIM), the jumps in the sequence. (Try "1 1 ... 1" if you have no 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 bases, usually distinct primes, but any NEGATIVE base generates values J/|BASE|. (Try "-N 2 3 5 7 11 ..." if you have no preference.) (any value of 0 or 1 terminates execution.) User input -10 2 The data was written to the file "hammersley_02_00020.txt". * * * Ready to generate a new dataset: * * Enter NDIM, the spatial dimension: (Try "2" if you have no preference.) (0 or any negative value terminates execution). User input NDIM = 3 Enter N, the number of points to generate: (Try "25" if you have no 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 have no preference.) (Any negative value terminates execution). User input STEP = 0 Enter SEED(1:NDIM), the sequence index corresponding to STEP = 0 (Try "0 0 ... 0" if you have no preference.) (a negative value terminates execution.) User input 0 0 0 Enter LEAP(1:NDIM), the jumps in the sequence. (Try "1 1 ... 1" if you have no preference.) (another choice is any prime larger than all bases.) (any value less than 1 terminates execution.) User input 1 5 1 Enter BASE(1:NDIM), the bases, usually distinct primes, but any NEGATIVE base generates values J/|BASE|. (Try "-N 2 3 5 7 11 ..." if you have no preference.) (any value of 0 or 1 terminates execution.) User input 2 -100 3 The data was written to the file "hammersley_03_00010.txt". * * * Ready to generate a new dataset: * * Enter NDIM, the spatial dimension: (Try "2" if you have no preference.) (0 or any negative value terminates execution). User input NDIM = 7 Enter N, the number of points to generate: (Try "25" if you have no 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 have no preference.) (Any negative value terminates execution). User input STEP = 0 Enter SEED(1:NDIM), the sequence index corresponding to STEP = 0 (Try "0 0 ... 0" if you have no preference.) (a negative value terminates execution.) User input 0 0 0 0 3 0 0 Enter LEAP(1:NDIM), the jumps in the sequence. (Try "1 1 ... 1" if you have no 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 2 Enter BASE(1:NDIM), the bases, usually distinct primes, but any NEGATIVE base generates values J/|BASE|. (Try "-N 2 3 5 7 11 ..." if you have no preference.) (any value of 0 or 1 terminates execution.) User input 2 -1 4 -2 2 -20 2 The data was written to the file "hammersley_07_00010.txt". * * * Ready to generate a new dataset: * * Enter NDIM, the spatial dimension: (Try "2" if you have no preference.) (0 or any negative value terminates execution). User input NDIM = -1 HALHAM_NDIM_CHECK - Fatal error! NDIM < 1. NDIM = -1 HAMMERSLEY_DATASET The input value of NDIM = -1 is interpreted as a request for termination. Normal end of execution. September 23 2004 11:02:10.680 AM