Delila Program: twogau

# twogau program

## By downloading this code you agree to the Source Code Use License (PDF). Pascal source code: twogau.p (wget instructions) Instructions on compiling MacOS binary: twogau Alphabetic List of Delila Programs Delila Programs by Most Recent Update Please report broken links delilabundle.zip = All Programs and MacOS Binaries Copyright Statement for Delila Programs

### Documentation for the twogau program is below, with links to related programs in the "see also" section.

```{   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

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}

```