| -- CXA5A05.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 the functions Arcsin and Arcsinh provide correct |
| -- results. |
| -- |
| -- TEST DESCRIPTION: |
| -- This test examines both the version of Arcsin and Arcsinh |
| -- the instantiation of the Ada.Numerics.Generic_Elementary_Functions |
| -- with a type derived from type Float, as well as the preinstantiated |
| -- version of this package for type Float. |
| -- Prescribed results, including instances prescribed to raise |
| -- exceptions, are examined in the test cases. In addition, |
| -- certain evaluations are performed where the actual function result |
| -- is compared with the expected result (within an epsilon range of |
| -- accuracy). |
| -- |
| -- TEST FILES: |
| -- The following files comprise this test: |
| -- |
| -- FXA5A00.A (foundation code) |
| -- CXA5A05.A |
| -- |
| -- |
| -- CHANGE HISTORY: |
| -- 20 Mar 95 SAIC Initial prerelease version. |
| -- 06 Apr 95 SAIC Corrected errors in context clause reference and |
| -- use of Cycle parameter. |
| -- 13 Jun 95 SAIC Incorporated use of Dont_Optimize procedure, and |
| -- use of Result_Within_Range function overloaded for |
| -- FXA5A00.New_Float_Type. |
| -- 28 Feb 97 PWB.CTA Removed checks with explict Cycle => 2.0*Pi |
| -- |
| -- CHANGE NOTE: |
| -- According to Ken Dritz, author of the Numerics Annex of the RM, |
| -- one should never specify the cycle 2.0*Pi for the trigonometric |
| -- functions. In particular, if the machine number for the first |
| -- argument is not an exact multiple of the machine number for the |
| -- explicit cycle, then the specified exact results cannot be |
| -- reasonably expected. The affected checks in this test have been |
| -- marked as comments, with the additional notation "pwb-math". |
| -- Phil Brashear |
| --! |
| |
| with Ada.Numerics.Elementary_Functions; |
| with Ada.Numerics.Generic_Elementary_Functions; |
| with FXA5A00; |
| with Report; |
| |
| procedure CXA5A05 is |
| begin |
| |
| Report.Test ("CXA5A05", "Check that the functions Arcsin and Arcsinh " & |
| "provide correct results"); |
| |
| Test_Block: |
| declare |
| |
| use Ada.Numerics; |
| use FXA5A00; |
| |
| package GEF is new Ada.Numerics.Generic_Elementary_Functions(New_Float); |
| package EF renames Ada.Numerics.Elementary_Functions; |
| |
| The_Result : Float; |
| New_Float_Result : New_Float; |
| |
| procedure Dont_Optimize_Float is new Dont_Optimize(Float); |
| procedure Dont_Optimize_New_Float is new Dont_Optimize(New_Float); |
| |
| begin |
| |
| -- Testing of Function Arcsin, both instantiated and pre-instantiated |
| -- versions. |
| |
| -- Check that Argument_Error is raised by the Arcsin function when |
| -- the absolute value of the parameter X is greater than 1.0. |
| |
| begin |
| New_Float_Result := GEF.Arcsin(New_Float(FXA5A00.One_Plus_Delta)); |
| Report.Failed("Argument_Error not raised by Arcsin function " & |
| "when provided a parameter value larger than 1.0"); |
| Dont_Optimize_New_Float(New_Float_Result, 1); |
| exception |
| when Argument_Error => null; -- OK, expected exception. |
| when others => |
| Report.Failed("Unexpected exception raised by Arcsin function " & |
| "when provided a parameter value larger than 1.0"); |
| end; |
| |
| begin |
| The_Result := EF.Arcsin(FXA5A00.Minus_Large); |
| Report.Failed("Argument_Error not raised by Arcsin function " & |
| "when provided a large negative parameter value"); |
| Dont_Optimize_Float(The_Result, 2); |
| exception |
| when Argument_Error => null; -- OK, expected exception. |
| when others => |
| Report.Failed("Unexpected exception raised by Arcsin function " & |
| "when provided a large negative parameter value"); |
| end; |
| |
| |
| -- Check the prescribed result of function Arcsin with parameter 0.0. |
| |
| if GEF.Arcsin(X => 0.0) /= 0.0 or |
| EF.Arcsin(0.0) /= 0.0 |
| then |
| Report.Failed("Incorrect result from Function Arcsin when the " & |
| "value of the parameter X is 0.0"); |
| end if; |
| |
| |
| -- Check the results of the Arcsin function with various input |
| -- parameters. |
| |
| if not Result_Within_Range(GEF.Arcsin(1.0), 1.571, 0.001) or |
| not Result_Within_Range( EF.Arcsin(0.62), 0.669, 0.001) or |
| not Result_Within_Range(GEF.Arcsin(0.01), 0.010, 0.001) or |
| not Result_Within_Range( EF.Arcsin(-0.29), -0.294, 0.001) or |
| not Result_Within_Range(GEF.Arcsin(-0.50), -0.524, 0.001) or |
| not Result_Within_Range( EF.Arcsin(-1.0), -1.571, 0.001) |
| then |
| Report.Failed("Incorrect result from Function Arcsin with " & |
| "various input parameters"); |
| end if; |
| |
| |
| -- Testing of Function Arcsin with specified Cycle parameter. |
| |
| --pwb-math -- Check that Argument_Error is raised by the Arcsin function with |
| --pwb-math -- specified cycle, whenever the absolute value of the parameter X |
| --pwb-math -- is greater than 1.0. |
| --pwb-math |
| --pwb-math begin |
| --pwb-math New_Float_Result := GEF.Arcsin(New_Float(FXA5A00.Large), 2.0*Pi); |
| --pwb-math Report.Failed("Argument_Error not raised by Function Arcsin " & |
| --pwb-math "with specified cycle, when provided a large " & |
| --pwb-math "positive input parameter"); |
| --pwb-math Dont_Optimize_New_Float(New_Float_Result, 3); |
| --pwb-math exception |
| --pwb-math when Argument_Error => null; -- OK, expected exception. |
| --pwb-math when others => |
| --pwb-math Report.Failed("Unexpected exception raised by Function Arcsin " & |
| --pwb-math "with specified cycle, when provided a large " & |
| --pwb-math "positive input parameter"); |
| --pwb-math end; |
| --pwb-math |
| --pwb-math begin |
| --pwb-math The_Result := EF.Arcsin(FXA5A00.Minus_One_Minus_Delta, 2.0*Pi); |
| --pwb-math Report.Failed("Argument_Error not raised by Function Arcsin " & |
| --pwb-math "with specified cycle, when provided an input " & |
| --pwb-math "parameter less than -1.0"); |
| --pwb-math Dont_Optimize_Float(The_Result, 4); |
| --pwb-math exception |
| --pwb-math when Argument_Error => null; -- OK, expected exception. |
| --pwb-math when others => |
| --pwb-math Report.Failed("Unexpected exception raised by Function Arcsin " & |
| --pwb-math "with specified cycle, when provided an input " & |
| --pwb-math "parameter less than -1.0"); |
| --pwb-math end; |
| --pwb-math |
| -- Check that Argument_Error is raised by the Arcsin function with |
| -- specified cycle, whenever the Cycle parameter is zero or negative. |
| |
| begin |
| New_Float_Result := GEF.Arcsin(2.0, 0.0); |
| Report.Failed("Argument_Error not raised by Function Arcsin " & |
| "with specified cycle of 0.0"); |
| Dont_Optimize_New_Float(New_Float_Result, 5); |
| exception |
| when Argument_Error => null; -- OK, expected exception. |
| when others => |
| Report.Failed("Unexpected exception raised by Function Arcsin " & |
| "with specified cycle of 0.0"); |
| end; |
| |
| begin |
| The_Result := EF.Arcsin(2.0, -2.0*Pi); |
| Report.Failed("Argument_Error not raised by Function Arcsin " & |
| "with specified negative cycle parameter"); |
| Dont_Optimize_Float(The_Result, 6); |
| exception |
| when Argument_Error => null; -- OK, expected exception. |
| when others => |
| Report.Failed("Unexpected exception raised by Function Arcsin " & |
| "with specified negative cycle parameter"); |
| end; |
| |
| |
| --pwb-math -- Check the prescribed result of function Arcsin with specified Cycle |
| --pwb-math -- parameter, when the value of parameter X is 0.0. |
| --pwb-math |
| --pwb-math if GEF.Arcsin(X => 0.0, Cycle => 2.0*Pi) /= 0.0 or |
| --pwb-math EF.Arcsin(0.0, 2.0*Pi) /= 0.0 |
| --pwb-math then |
| --pwb-math Report.Failed("Incorrect result from Function Arcsin with " & |
| --pwb-math "specified Cycle parameter, when the value " & |
| --pwb-math "of parameter X is 0.0"); |
| --pwb-math end if; |
| --pwb-math |
| --pwb-math |
| --pwb-math -- Test of the Arcsin function with specified Cycle parameter with |
| --pwb-math -- various input parameters. |
| --pwb-math |
| --pwb-math if not FXA5A00.Result_Within_Range(GEF.Arcsin( 0.01, 2.0*Pi), |
| --pwb-math 0.010, |
| --pwb-math 0.001) or |
| --pwb-math not FXA5A00.Result_Within_Range( EF.Arcsin( 0.14, 2.0*Pi), |
| --pwb-math 0.141, |
| --pwb-math 0.001) or |
| --pwb-math not FXA5A00.Result_Within_Range(GEF.Arcsin( 0.37, 2.0*Pi), |
| --pwb-math 0.379, |
| --pwb-math 0.001) or |
| --pwb-math not FXA5A00.Result_Within_Range( EF.Arcsin( 0.55, 2.0*Pi), |
| --pwb-math 0.582, |
| --pwb-math 0.001) or |
| --pwb-math not FXA5A00.Result_Within_Range(GEF.Arcsin(-0.22, 2.0*Pi), |
| --pwb-math -0.222, |
| --pwb-math 0.001) or |
| --pwb-math not FXA5A00.Result_Within_Range( EF.Arcsin(-0.99, 2.0*Pi), |
| --pwb-math -1.43, |
| --pwb-math 0.01) or |
| --pwb-math not FXA5A00.Result_Within_Range( EF.Arcsin(1.0, 360.0), |
| --pwb-math 90.0, |
| --pwb-math 0.1) or |
| --pwb-math not FXA5A00.Result_Within_Range( EF.Arcsin(1.0, 100.0), |
| --pwb-math 25.0, |
| --pwb-math 0.1) |
| --pwb-math then |
| --pwb-math Report.Failed("Incorrect result from Arcsin with specified " & |
| --pwb-math "cycle parameter with various input parameters"); |
| --pwb-math end if; |
| |
| -- Testing of Arcsinh Function, both instantiated and pre-instantiated |
| -- version. |
| |
| -- Check that no exception occurs on computing the Arcsinh with very |
| -- large (positive and negative) input values. |
| |
| begin |
| New_Float_Result := GEF.Arcsinh(New_Float(FXA5A00.Large)); |
| Dont_Optimize_New_Float(New_Float_Result, 7); |
| exception |
| when others => |
| Report.Failed("Unexpected exception on Arcsinh with large " & |
| "positive value"); |
| end; |
| |
| begin |
| The_Result := EF.Arcsinh(FXA5A00.Minus_Large); |
| Dont_Optimize_Float(The_Result, 8); |
| exception |
| when others => |
| Report.Failed("Unexpected exception on Arcsinh with large " & |
| "negative value"); |
| end; |
| |
| |
| -- Check that no exception occurs on computing the Arcsinh with very |
| -- small (positive and negative) input values. |
| |
| begin |
| New_Float_Result := GEF.Arcsinh(New_Float(FXA5A00.Small)); |
| Dont_Optimize_New_Float(New_Float_Result, 9); |
| exception |
| when others => |
| Report.Failed("Unexpected exception on Arcsinh with small " & |
| "positive value"); |
| end; |
| |
| begin |
| The_Result := EF.Arcsinh(-FXA5A00.Small); |
| Dont_Optimize_Float(The_Result, 10); |
| exception |
| when others => |
| Report.Failed("Unexpected exception on Arcsinh with small " & |
| "negative value"); |
| end; |
| |
| |
| -- Check function Arcsinh for prescribed result with parameter 0.0. |
| |
| if GEF.Arcsinh(X => 0.0) /= 0.0 or |
| EF.Arcsinh(X => 0.0) /= 0.0 |
| then |
| Report.Failed("Incorrect result from Function Arcsinh when " & |
| "provided a 0.0 input parameter"); |
| end if; |
| |
| |
| -- Check the results of the Arcsinh function with various input |
| -- parameters. |
| |
| if not Result_Within_Range(GEF.Arcsinh(0.15), 0.149, 0.001) or |
| not Result_Within_Range( EF.Arcsinh(0.82), 0.748, 0.001) or |
| not Result_Within_Range(GEF.Arcsinh(1.44), 1.161, 0.001) or |
| not Result_Within_Range(GEF.Arcsinh(6.70), 2.601, 0.001) or |
| not Result_Within_Range( EF.Arcsinh(Pi), 1.862, 0.001) or |
| not Result_Within_Range( EF.Arcsinh(-Pi), -1.862, 0.001) or |
| not Result_Within_Range(GEF.Arcsinh(-1.0), -0.881, 0.001) or |
| not Result_Within_Range( EF.Arcsinh(-5.5), -2.406, 0.001) |
| then |
| Report.Failed("Incorrect result from Function Arcsin with " & |
| "various input parameters"); |
| end if; |
| |
| |
| exception |
| when others => Report.Failed ("Exception raised in Test_Block"); |
| end Test_Block; |
| |
| Report.Result; |
| |
| end CXA5A05; |