| -- C490001.A |
| -- |
| -- Grant of Unlimited Rights |
| -- |
| -- Under contracts F33600-87-D-0337, F33600-84-D-0280, MDA903-79-C-0687, |
| -- F08630-91-C-0015, and DCA100-97-D-0025, the U.S. Government obtained |
| -- unlimited rights in the software and documentation contained herein. |
| -- Unlimited rights are defined in DFAR 252.227-7013(a)(19). By making |
| -- this public release, the Government intends to confer upon all |
| -- recipients unlimited rights equal to those held by the Government. |
| -- These rights include rights to use, duplicate, release or disclose the |
| -- released technical data and computer software in whole or in part, in |
| -- any manner and for any purpose whatsoever, and to have or permit others |
| -- to do so. |
| -- |
| -- DISCLAIMER |
| -- |
| -- ALL MATERIALS OR INFORMATION HEREIN RELEASED, MADE AVAILABLE OR |
| -- DISCLOSED ARE AS IS. THE GOVERNMENT MAKES NO EXPRESS OR IMPLIED |
| -- WARRANTY AS TO ANY MATTER WHATSOEVER, INCLUDING THE CONDITIONS OF THE |
| -- SOFTWARE, DOCUMENTATION OR OTHER INFORMATION RELEASED, MADE AVAILABLE |
| -- OR DISCLOSED, OR THE OWNERSHIP, MERCHANTABILITY, OR FITNESS FOR A |
| -- PARTICULAR PURPOSE OF SAID MATERIAL. |
| --* |
| -- |
| -- OBJECTIVE: |
| -- Check that, for a real static expression that is not part of a larger |
| -- static expression, and whose expected type T is a floating point type |
| -- that is not a descendant of a formal scalar type, the value is rounded |
| -- to the nearest machine number of T if T'Machine_Rounds is true, and is |
| -- truncated otherwise. Check that if rounding is performed, and the value |
| -- is exactly halfway between two machine numbers, one of the two machine |
| -- numbers is used. |
| -- |
| -- TEST DESCRIPTION: |
| -- The test obtains a machine number M1 for a floating point subtype S by |
| -- passing a real literal to S'Machine. It then obtains an adjacent |
| -- machine number M2 by using S'Succ (or S'Pred). It then constructs |
| -- values which lie between these two machine numbers: one (A) which is |
| -- closer to M1, one (B) which is exactly halfway between M1 and M2, and |
| -- one (C) which is closer to M2. This is done for both positive and |
| -- negative machine numbers. |
| -- |
| -- Let M1 be closer to zero than M2. Then if S'Machine_Rounds is true, |
| -- C must be rounded to M2, A must be rounded to M1, and B must be rounded |
| -- to either M1 or M2. If S'Machine_Rounds is false, all the values must |
| -- be truncated to M1. |
| -- |
| -- A, B, and C are constructed using the following static expressions: |
| -- |
| -- A: constant S := M1 + (M2 - M1)*Z; -- Z slightly less than 0.5. |
| -- B: constant S := M1 + (M2 - M1)*Z; -- Z equals 0.5. |
| -- C: constant S := M1 + (M2 - M1)*Z; -- Z slightly more than 0.5. |
| -- |
| -- Since these are static expressions, they must be evaluated exactly, |
| -- and no rounding may occur until the final result is calculated. |
| -- |
| -- The checks for equality between the members of (A, B, C) and (M1, M2) |
| -- are performed at run-time within the body of a subprogram. |
| -- |
| -- The test performs additional checks that the rounding performed on |
| -- real literals is consistent for a floating point subtype. A literal is |
| -- assigned to a constant of a floating point subtype S. The same literal |
| -- is then passed to a subprogram, along with the constant, and an |
| -- equality check is performed within the body of the subprogram. |
| -- |
| -- |
| -- CHANGE HISTORY: |
| -- 25 Sep 95 SAIC Initial prerelease version. |
| -- 25 May 01 RLB Repaired to work with the repeal of the round away |
| -- rule by AI-268. |
| -- |
| --! |
| |
| with System; |
| package C490001_0 is |
| |
| type My_Flt is digits System.Max_Digits; |
| |
| procedure Float_Subtest (A, B: in My_Flt; Msg: in String); |
| |
| procedure Float_Subtest (A, B, C: in My_Flt; Msg: in String); |
| |
| |
| -- |
| -- Positive cases: |
| -- |
| |
| -- |----|-------------|-----------------|-------------------|-----------| |
| -- | | | | | | |
| -- 0 P_M1 Less_Pos_Than_Half Pos_Exactly_Half More_Pos_Than_Half P_M2 |
| |
| |
| Positive_Float : constant My_Flt := 12.440193950021943; |
| |
| -- The literal value 12.440193950021943 is rounded up or down to the |
| -- nearest machine number of My_Flt when Positive_Float is initialized. |
| -- The value of Positive_Float should therefore be a machine number, and |
| -- the use of 'Machine in the initialization of P_M1 will be redundant for |
| -- a correct implementation. It's done anyway to make certain that P_M1 is |
| -- a machine number, independent of whether an implementation correctly |
| -- performs rounding. |
| |
| P_M1 : constant My_Flt := My_Flt'Machine(Positive_Float); |
| P_M2 : constant My_Flt := My_Flt'Succ(P_M1); |
| |
| -- P_M1 and P_M2 are adjacent machine numbers. Note that because it is not |
| -- certain whether 12.440193950021943 is a machine number, nor whether |
| -- 'Machine rounds it up or down, 12.440193950021943 may not lie between |
| -- P_M1 and P_M2. The test does not depend on this information, however; |
| -- the literal is only used as a "seed" to obtain the machine numbers. |
| |
| |
| -- The following entities are used to verify that rounding is performed |
| -- according to the value of 'Machine_Rounds. If language rules are |
| -- obeyed, the intermediate expressions in the following static |
| -- initialization expressions will not be rounded; all calculations will |
| -- be performed exactly. The final result, however, will be rounded to |
| -- a machine number (either P_M1 or P_M2, depending on the value of |
| -- My_Flt'Machine_Rounds). Thus, the value of each constant below will |
| -- equal that of P_M1 or P_M2. |
| |
| Less_Pos_Than_Half : constant My_Flt := P_M1 + ((P_M2 - P_M1)*2.9/6.0); |
| Pos_Exactly_Half : constant My_Flt := P_M1 + ((P_M2 - P_M1)/2.0); |
| More_Pos_Than_Half : constant My_Flt := P_M1 + ((P_M2 - P_M1)*4.6/9.0); |
| |
| |
| -- |
| -- Negative cases: |
| -- |
| |
| -- -|-------------|-----------------|-------------------|-----------|----| |
| -- | | | | | | |
| -- N_M2 More_Neg_Than_Half Neg_Exactly_Half Less_Neg_Than_Half N_M1 0 |
| |
| |
| -- The descriptions for the positive cases above apply to the negative |
| -- cases below as well. Note that, for N_M2, 'Pred is used rather than |
| -- 'Succ. Thus, N_M2 is further from 0.0 (i.e. more negative) than N_M1. |
| |
| Negative_Float : constant My_Flt := -0.692074550952117; |
| |
| |
| N_M1 : constant My_Flt := My_Flt'Machine(Negative_Float); |
| N_M2 : constant My_Flt := My_Flt'Pred(N_M1); |
| |
| More_Neg_Than_Half : constant My_Flt := N_M1 + ((N_M2 - N_M1)*4.1/8.0); |
| Neg_Exactly_Half : constant My_Flt := N_M1 + ((N_M2 - N_M1)/2.0); |
| Less_Neg_Than_Half : constant My_Flt := N_M1 + ((N_M2 - N_M1)*2.4/5.0); |
| |
| end C490001_0; |
| |
| |
| --==================================================================-- |
| |
| |
| with TCTouch; |
| package body C490001_0 is |
| |
| procedure Float_Subtest (A, B: in My_Flt; Msg: in String) is |
| begin |
| TCTouch.Assert (A = B, Msg); |
| end Float_Subtest; |
| |
| procedure Float_Subtest (A, B, C: in My_Flt; Msg: in String) is |
| begin |
| TCTouch.Assert (A = B or A = C, Msg); |
| end Float_Subtest; |
| |
| end C490001_0; |
| |
| |
| --==================================================================-- |
| |
| |
| with C490001_0; -- Floating point support. |
| use C490001_0; |
| |
| with Report; |
| procedure C490001 is |
| begin |
| Report.Test ("C490001", "Rounding of real static expressions: " & |
| "floating point subtypes"); |
| |
| |
| -- Check that rounding direction is consistent for literals: |
| |
| Float_Subtest (12.440193950021943, P_M1, "Positive Float: literal"); |
| Float_Subtest (-0.692074550952117, N_M1, "Negative Float: literal"); |
| |
| |
| -- Now check that rounding is performed correctly for values between |
| -- machine numbers, according to the value of 'Machine_Rounds: |
| |
| if My_Flt'Machine_Rounds then |
| Float_Subtest (Pos_Exactly_Half, P_M1, P_M2, "Positive Float: = half"); |
| Float_Subtest (More_Pos_Than_Half, P_M2, "Positive Float: > half"); |
| Float_Subtest (Less_Pos_Than_Half, P_M1, "Positive Float: < half"); |
| |
| Float_Subtest (Neg_Exactly_Half, N_M1, N_M2, "Negative Float: = half"); |
| Float_Subtest (More_Neg_Than_Half, N_M2, "Negative Float: > half"); |
| Float_Subtest (Less_Neg_Than_Half, N_M1, "Negative Float: < half"); |
| else |
| Float_Subtest (Pos_Exactly_Half, P_M1, "Positive Float: = half"); |
| Float_Subtest (More_Pos_Than_Half, P_M1, "Positive Float: > half"); |
| Float_Subtest (Less_Pos_Than_Half, P_M1, "Positive Float: < half"); |
| |
| Float_Subtest (Neg_Exactly_Half, N_M1, "Negative Float: = half"); |
| Float_Subtest (More_Neg_Than_Half, N_M1, "Negative Float: > half"); |
| Float_Subtest (Less_Neg_Than_Half, N_M1, "Negative Float: < half"); |
| end if; |
| |
| |
| Report.Result; |
| end C490001; |