spline - Man Page

Use a local interpolation algorithm (the cubic Bessel algorithm), so that spline can be used as a real-time filter. The slope of the interpolating curve at each point in a dataset will be chosen by fitting a quadratic function through that point and the two adjacent points in the dataset. If -f is specified then the -t option, otherwise optional, must be used as well. Also, if -f is specified then the -k, -p, and -T options may not be used.

If -f is not specified, then the default (global) interpolation algorithm will be used.

Set the boundary condition parameter for each constructed spline to be k. (The default value is 1.0.) In each of its components, the spline will satisfy the two boundary conditions y"[0]=ky"[1] and y"[n]=ky"[n-1]. Here y[0] and y[1] signify the values of a specified component of the vector-valued dependent variable y at the first two points of a dataset, and y[n-1] and y[n] the values at the last two points. Setting k to zero will yield a "natural" spline, i.e., one that has zero curvature at the two ends of the dataset. The -k option may not be used if -f or -p is specified.

Subdivide the interval over which interpolation occurs into n subintervals. The number of data points computed, and written to the output, will be n+1. The default value for n is 100.

Construct a periodic spline. If this option is specified, the y values for the first and last points in each dataset must be equal. The -f and -k options may not be used if -p is specified.

Each interpolating curve will be a spline under tension. This option sets the tension value (the default is 0.0).

If tension equals zero, the curve will be a piecewise cubic spline. Increasing the tension above zero makes the curve "tighter", and  reduces the likelihood of spurious inflection points. That is because between each pair of successive points in a dataset, the curve will satisfy the fourth-order differential equation  y""=sgn(tension)*(tension^2)y" in each of its components. As tension increases to positive infinity, it will converge to a polygonal line. The -T option may not be used if -f is specified.

For each dataset, set the interval over which interpolation occurs to be the interval between  tmin and tmax. If tspacing is not specified, the interval will be divided into the number of subintervals specified by the -n option.

If the -t option is not used, the interval over which interpolation occurs will be the entire range of the independent variable in the dataset. The -t option must always be used if the -f option is used to request filter-like behavior (see above).

Set the dimensionality of the dependent variable y in the input and output files to be dimension. The default dimension is 1.

Set the data format for the input file(s) to be data-format, which may be one of the following.

a

ASCII format (the default). Each file is a sequence of floating point numbers, interpreted as the  t and y coordinates of the successive data points in a dataset. If y is d-dimensional,  there will be d+1 numbers for each point. The t and  y coordinates of a point need not appear on the same line, and points need not appear on different lines. But if a blank line occurs (i.e., two newlines in succession are seen), it is interpreted as the end of a dataset, and the beginning of the next.

f

Single precision binary format. Each file is a sequence of floating point numbers, interpreted as the t and y coordinates of the successive data points in a dataset. If y is d-dimensional,  there will be d+1 numbers for each point. Successive datasets are separated by a single occurrence of the quantity FLT_MAX, which is the largest possible single precision floating point number. On most machines this is approximately 3.4x10^38.

d

Double precision binary format. Each file is a sequence of double precision floating point numbers, interpreted as the t and y coordinates of the successive data points in a dataset. If y is d-dimensional, there will be d+1 numbers for each point. Successive datasets are separated by a single occurrence of the quantity DBL_MAX, which is the largest possible double precision floating point number. On most machines this is approximately 1.8x10^308.

i

Integer binary format. Each file is a sequence of integers, interpreted as the t and y coordinates of the successive data points  in a dataset. If y is d-dimensional,  there will be d+1 numbers for each point. Successive datasets are separated by a single occurrence of the quantity INT_MAX, which is the largest possible integer. On most machines this is 2^31-1.

Automatically generate values for t, the independent variable (the default values of step_size and lower_limit are 1.0 and 0.0, respectively).

Irrespective of data format (`a', `f', `d', or `i'), this option specifies that the values of t are missing from the input file: the dataset(s) to be read contain only values of y, the dependent variable. So if y is d-dimensional, there will be only d numbers for each point. The increment from each t value to the next will be step_size, and the first  t value will be lower_limit. This option is useful, e.g., when interpolating curves rather than functions.

Automatically generate values for t, the independent variable. This is a variant form of the -a option. The increment from each t value to the next will be the distance in d-dimensional space between the corresponding y values, and the first t value will be 0.0. That is, t will be "polygonal arclength". This option is useful when interpolating curves rather than functions.

Set the data format for the output file to be data-format. The interpretation of data-format is the same as for the -I option. The default is `a', i.e., ASCII format.

Set the numerical precision for the t and y values in the output file to be significant-digits. This takes effect only if the output file is written in `a' format, i.e., in ASCII. significant-digits must be a positive integer (the default is 6).

Omit the independent variable t from the output file; for each point, supply only the dependent variable y. If y is d-dimensional, there will be only d numbers for each point, not d+1. This option is useful when interpolating curves rather than functions.

Print a list of command-line options, and exit.

Print the version number of spline and the plotting utilities package, and exit.