By downloading this code you agree to the
Source Code Use License (PDF). |
{ version = 1.00; (* of twogau.p 2011 Nov 10} (* begin module describe.twogau *) (* name twogau: 2 dimensional Rayleigh function synopsis twogau(data: in, twogaup: in, color: output, output: out) files data: set of Gaussianly distributed variables from the program gentst. twogaup: parameters: first line: total dimensionality D. MUST BE 2 second line: number of points to do. If the end of the data is reached, the actual number of points generated is reported to output. third line: number of steps to generate the fD(r) graph. (The range for this is always -2.5 to +2.5 on both x and y axes.) If the number of steps is less than 1, then no smooth graph is done. fourth line: a real number, "0 <= partial <= 1" by which to multiply the actual fD(r) density by to obtain the density reported to the color file. This allows one to tone down the gray scale, or to avoid having the highest density of color equal the lowest (as when the hue is used and a hue of 1 is the same as a hue of 0). fifth line: printing of data on plot (one character): d=dimension,p=dimension+point,a=all,n=none color: a xyin file for input to the xyplo or riden program. The columns are: 1 symbols: f=from fD(r), s = simulated point'); 2 x: x coordinate 3 y: y coordinate 4 xwidth: width of symbol on x axis 5 ywidth: width of symbol on y axis 6 density: density 7 inverse: 1 - density (for inverse plotting) 8 maximum: MAXimum density 9 minimum: MINimum density 10 maximum: MAXimum density 11 minimum: MINimum density 12 partial: partial density for grey tones Partial is the largest density allowed. When plotted in color, hues come from a color wheel in which the highest color is almost identical to the lowest color. That is the color of hue=1 is almost identical to the color of hue = 0. To avoid this effect, make partial less than 1.0. A partial less than 1.0 also avoids completely black gray scale plots. output: messages to the user, number of points generated. description Plot a Gaussian on x and an independent one on y. old from ring: Simulate mapping from many-dimensional to 2-dimensional Z space. Sets of D Gaussian values are read from the data file, squared, summed and square rooted. The x and y value in Z space is determined from an angle and a radius. The angle is found from the last two Gaussian values, while the radius is determined by the noise (rms) for all dimensions. The statistical function fD(r) is to be graphed in color or gray scale using xyplo, while the simulated points are graphed as points on top of the smooth fD(r) function. The program output is ready to read into the xyplo plotting program. examples twogaup used for generating figures: 16 total dimensionality 100 number of points to do 128 steps for plotting smooth fD(r) graph 0.50 partial d d=dimension,p=dimension+point,a=all,n=none xyplop used for generating figures: 2 2 zerox zeroy graph coordinate center x -2.5 2.5 zx min max (character, real, real) if zx='x' then set xaxis y -2.5 2.5 zy min max (character, real, real) if zy='y' then set yaxis 10 10 xinterval yinterval number of intervals on axes to plot 4 4 xwidth ywidth width of numbers in characters 1 1 xdecimal ydecimal number of decimal places 5 5 xsize ysize size of axes in inches x y c zc if zc='c' then a crosshairs put on zero of x and y n 2 zxl base if zxl='l' then make x axis log to the given base n 2 zyl base if zyl='l' then make y axis log to the given base ********************************************************************* 2 3 xcolumn ycolumn columns of xyin that determine plot location 1 symbol column the xyin column to read symbols from 4 5 xscolumn yscolumn columns of xyin that determine the symbol size 10 8 7 hue saturation brightness columns for color manipulation ********************************************************************* r symbol to plot c(circle)bd(dotted box)x+Ifgpr(rectangle) b symbol flag character in xyin that indicates that this symbol -1.0 symbol sizex side in inches on the x axis of the symbol. -1.0 symbol sizey as for the x axis, get size from yscolumn n no connection (example for connection is c- 0.05 for dashed 0.05 inch) n 0.05 linetype size linetype l.-in and size of dashes or dots ********************************************************************* r symbol to plot c(circle)bd(dotted box)x+Ifgpr(rectangle) f symbol flag character in xyin that indicates that this symbol -1.0 symbol sizex side in inches on the x axis of the symbol. -1.0 symbol sizey as for the x axis, get size from yscolumn n no connection (example for connection is c- 0.05 for dashed 0.05 inch) n 0.05 linetype size linetype l.-in and size of dashes or dots ********************************************************************* c symbol to plot c(circle)bd(dotted box)x+Ifgpr(rectangle) s symbol flag character in xyin that indicates that this symbol 0.0858 symbol sizex side in inches on the x axis of the symbol. 0.0858 symbol sizey as for the x axis, get size from yscolumn n no connection (example for connection is c- 0.05 for dashed 0.05 inch) n 0.05 linetype size linetype l.-in and size of dashes or dots ********************************************************************* g symbol to plot c(circle)bd(dotted box)x+Ifgpr(rectangle) g symbol flag character in xyin that indicates that this symbol -1.0 symbol sizex side in inches on the x axis of the symbol. -1.0 symbol sizey as for the x axis, get size from yscolumn n no connection (example for connection is c- 0.05 for dashed 0.05 inch) n 0.05 linetype size linetype l.-in and size of dashes or dots ********************************************************************* . ********************************************************************* Useful color parameters are: 8 6 10 Light density plot, printable on a black and white device (best). 8 7 10 Dark density plot, printable on a black and white device. 6 8 10 Color plot, red background. 7 8 10 Color plot, purple background (neat). 6 7 10 Color and density varying to make the simulated points easy to see. (red background) 7 6 10 Color and density varying to make the simulated points easy to see. (white background - lovely!) Warning: since the program has changed, these may no longer be correct. documentation ccmm see also gentst.p xyplo.p riden.p author Thomas Dana Schneider bugs none known. Confirm that the density distribution is correct by using program riden. technical notes *) (* end module describe.twogau *) {This manual page was created by makman 1.45}{created by htmlink 1.62}