Returns information describing the computer’s floating-point arithmetic.

Synopsis

#include <imsls.h>

float imsls_f_machine (int n)

The type double function is imsls_d_machine.

Required Arguments

int n   (Input)
Index indicating which value is to be returned. The index must be between 1 and 8.

Return Value

The requested value is returned. If n is out of range, NaN is returned.

Description

Function imsls_f_machine returns information describing the computer’s floating-point arithmetic. This can be used to make programs machine independent. In addition, some of the functions are also important in setting missing values.

Assume that float numbers are represented in Nf-digit, base B form as

where σ is the sign; 0 ≤ xk < B for k = 1, 2, …, Nf; and

Note that B = imsls_i_machine(6); Nf = imsls_i_machine(7);

and

The ANSI/IEEE 754-1985 standard for binary arithmetic uses NaN as the result of various otherwise illegal operations, such as computing 0/0. On computers that do not support NaN, a value larger than imsls_d_machine(2) is returned for imsls_f_machine(6). On computers that do not have a special representation for infinity, imsls_f_machine(2) returns the same value as imsls_f_machine(7).

Function imsls_f_machine is defined by the following table:

Function imsls_d_machine retrieves machine constants that define the computer’s double arithmetic. Note that for double B = imsls_i_machine(6), Nd = imsls_i_machine(10),

and

Missing values in functions are always indicated by NaN. This is imsls_f_machine(6) in single precision and imsls_d_machine(6) in double precision. There is no missing-value indicator for integers. Users will almost always have to convert from their missing value indicators to NaN.

Example

In this example, all eight values returned by imsls_f_machine and by imsls_d_machine on a machine with IEEE arithmetic are printed.

#include <imsls.h>

main()
{
    int             n;
    float           fans;
    double          dans;

    for (n = 1;  n <= 8;  n++) {
        fans = imsls_f_machine(n);
        printf("imsls_f_machine(%d) = %g\n", n, fans);
    }

    for (n = 1;  n <= 8;  n++) {
        dans = imsls_d_machine(n);
        printf("imsls_d_machine(%d) = %g\n", n, dans);
    }
}

Output

imsls_f_machine(1) = 1.17549e-38
imsls_f_machine(2) = 3.40282e+38
imsls_f_machine(3) = 5.96046e-08
imsls_f_machine(4) = 1.19209e-07
imsls_f_machine(5) = 0.30103
imsls_f_machine(6) = NaN
imsls_f_machine(7) = Inf
imsls_f_machine(8) = -Inf
imsls_d_machine(1) = 2.22507e-308
imsls_d_machine(2) = 1.79769e+308
imsls_d_machine(3) = 1.11022e-16
imsls_d_machine(4) = 2.22045e-16
imsls_d_machine(5) = 0.30103
imsls_d_machine(6) = NaN
imsls_d_machine(7) = Inf
imsls_d_machine(8) = -Inf

Visual Numerics, Inc. Visual Numerics - Developers of IMSL and PV-WAVE http://www.vni.com/ PHONE: 713.784.3131 FAX:713.781.9260