Next: VAXStation 3100 (cvsa3) Up: Linux 1.1.68 (niespx Previous: FORTRAN Version (using

PASCAL Version (using gpc)


Diagnosis resumes after milestone no          7               Page:          3

Program is now RUNNING tests on small integers:
 -1, 0, 1/2, 1, 2, 3, 4, 5, 9, 27, 32 & 240 are O.K.

Searching for Radix and Precision.
Radix =  2.0000000e+00
Closest relative separation found is U1 =  1.1102230e-16

Recalculating radix and precision
 confirms closest relative separation U1 .
Radix confirmed.
The number of significant digits of the Radix is  5.3000000e+01
Some subexpressions appear to be calculated
extra precisely with about 
 1.1000000e+01extra B-digits i.e. 
roughly  3.3113300e+00 extra significant decimals.
That feature is not tested further by this program.
To continue, press any key and newline:


Diagnosis resumes after milestone no         30               Page:          4

Subtraction appears to be normalized, as it should be.

Checking for guard digit on *, /, and -.
 *, /, and - have guard digits, as they should.
To continue, press any key and newline:


Diagnosis resumes after milestone no         40               Page:          5

Checking rounding on multiply, divide and add/subtract.
* is neither chopped nor correctly rounded.
/ is neither chopped nor correctly rounded.
Addition/Subtraction neither rounds nor chops.
Sticky bit used incorrectly or not at all.
FLAW:  lack(s) of guard digits or failure(s) to
correctly round or chop (noted above) count as one
flaw in the final tally below.

Does Multiplication commute? Testing on         20 random pairs.
No failures found in         20 integer pairs.

Running test of square root(x).
Testing if sqrt(X * X) = X for         20 integers X.
Test for Sqrt Monotonicity.
Sqrt has passed a test for Monotonicity.
Testing whether sqrt is rounded or chopped: 
Square root is neither chopped nor correctly rounded.
Observed errors run from -5.0000000e-01
to  5.0000000e-01 ulps.
To continue, press any key and newline:


Diagnosis resumes after milestone no         90               Page:          6

Testing powers Z^i for small integers Z and i.
DEFECT:  computed ( 2.0000000e+00) ^ ( 3.0000000e+00) = 
 8.0000000e+00, which compares unequal to correct  8.0000000e+00;
 they differ by -1.7763568e-15
Error like this may invalidate financial 
calculations involving interest rates.
Similar discrepancies have occurred  2.4700000e+02 times.

To continue, press any key and newline:


Diagnosis resumes after milestone no        100               Page:          7

Seeking Underflow thresholds UnderflowThreshold and E0
Smallest strictly positive number found is E0 = 4.9406565e-324
Since comparison denies Z = 0, evaluating
(Z + Z) / Z should be safe.
What the machine gets for (Z + Z) / Z is:  2.0000000e+00
This is O.K., provided Over/Underflow has NOT just been signaled.
Underflow is gradual; it incurs Absolute Error =
(roundoff in UnderflowThreshold) < E0.
The Underflow threshold is  2.2250739e-308 below which
calculation may suffer larger Relative error then merely roundoff.

Testing X^((X + 1) / (X - 1)) vs. exp(2) =  7.3890561e+00 as X -> 1.
Accuracy seems adequate.
Testing powers Z^Q at four nearly extreme values.
DEFECT:  computed ( 2.0000000e+00) ^ ( 1.0070000e+03) = 
 1.3715310e+303, which compares unequal to correct  1.3715310e+303;
 they differ by  1.8272464e+288
Similar discrepancies have occurred  4.0000000e+00 times.

To continue, press any key and newline:


Diagnosis resumes after milestone no        150               Page:          8

To continue, press any key and newline:


Diagnosis resumes after milestone no        160               Page:          9

Searching for Overflow threshold:  This may generate an error.
Try a few values for N, starting with a large one,
and take the one that just does not stop the machine.
Did you find the correct value for N yet?
N = 
O.K.
Can "Z = -Y" overflow?  Trying it on Y = -8.9884657e+307
Seems O.K.
Overflow threshold is V =  1.7976931e+308
There is no saturation value because
the system traps on overflow.
No Overflow should get signaled for V * 1 =  1.7976931e+308
                            nor for V / 1 =  1.7976931e+308
Any overflow signal separating this * from the one
above is a DEFECT.

To continue, press any key and newline:


Diagnosis resumes after milestone no        190               Page:         10


What message and/or values does Division by Zero produce?
This can interupt your program. You can skip this part if you wish.
Do you wish to compute 1 / 0? 
O.K.
Do you wish to compute 0 / 0?
O.K.
To continue, press any key and newline:


Diagnosis resumes after milestone no        220               Page:         11


The number of  DEFECTs  encountered =                  2
The number of  FLAWs  encountered =                    1

The arithmetic diagnosed may be Acceptable
despite inconvenient Defects.
END OF TEST.


peter@physik3.gwdg.de
Tue Aug 8 13:05:00 GMT+0200 1995