gcc/testsuite/ada/acats/tests/cxg/cxg1001.a

   1 -- CXG1001.A
   2 --
   3 --                             Grant of Unlimited Rights
   4 --
   5 --     Under contracts F33600-87-D-0337, F33600-84-D-0280, MDA903-79-C-0687,
   6 --     F08630-91-C-0015, and DCA100-97-D-0025, the U.S. Government obtained
   7 --     unlimited rights in the software and documentation contained herein.
   8 --     Unlimited rights are defined in DFAR 252.227-7013(a)(19).  By making
   9 --     this public release, the Government intends to confer upon all
  10 --     recipients unlimited rights  equal to those held by the Government.
  11 --     These rights include rights to use, duplicate, release or disclose the
  12 --     released technical data and computer software in whole or in part, in
  13 --     any manner and for any purpose whatsoever, and to have or permit others
  14 --     to do so.
  15 --
  16 --                                    DISCLAIMER
  17 --
  18 --     ALL MATERIALS OR INFORMATION HEREIN RELEASED, MADE AVAILABLE OR
  19 --     DISCLOSED ARE AS IS.  THE GOVERNMENT MAKES NO EXPRESS OR IMPLIED
  20 --     WARRANTY AS TO ANY MATTER WHATSOEVER, INCLUDING THE CONDITIONS OF THE
  21 --     SOFTWARE, DOCUMENTATION OR OTHER INFORMATION RELEASED, MADE AVAILABLE
  22 --     OR DISCLOSED, OR THE OWNERSHIP, MERCHANTABILITY, OR FITNESS FOR A
  23 --     PARTICULAR PURPOSE OF SAID MATERIAL.
  24 --*
  25 --
  26 -- OBJECTIVE:
  27 --      Check that the subprograms defined in the package
  28 --      Ada.Numerics.Generic_Complex_Types provide correct results.
  29 --      Specifically, check the functions Re, Im (both versions), procedures
  30 --      Set_Re, Set_Im (both versions), functions Compose_From_Cartesian (all
  31 --      versions), Compose_From_Polar, Modulus, Argument, and "abs".
  32 --
  33 -- TEST DESCRIPTION:
  34 --      The generic package Generic_Complex_Types
  35 --      is instantiated with a real type (new Float), and the results
  36 --      produced by the specified subprograms are verified.
  37 --
  38 -- APPLICABILITY CRITERIA:
  39 --      This test applies only to implementations supporting the
  40 --      Numerics Annex.
  41 --
  42 --
  43 -- CHANGE HISTORY:
  44 --      06 Dec 94   SAIC    ACVC 2.0
  45 --      15 Nov 95   SAIC    Corrected visibility problems for ACVC 2.0.1.
  46 --                          Modified subtest for Compose_From_Polar.
  47 --      29 Sep 96   SAIC    Incorporated reviewer comments.
  48 --
  49 --!
  50
  51 with Ada.Numerics.Generic_Complex_Types;
  52 with Report;
  53
  54 procedure CXG1001 is
  55
  56 begin
  57
  58    Report.Test ("CXG1001", "Check that the subprograms defined in " &
  59                            "the package Ada.Numerics.Generic_Complex_Types " &
  60                            "provide correct results");
  61
  62    Test_Block:
  63    declare
  64
  65       type Real_Type is new Float;
  66
  67       package Complex_Pack is new
  68         Ada.Numerics.Generic_Complex_Types(Real_Type);
  69
  70       use type Complex_Pack.Complex;
  71
  72       -- Declare a zero valued complex number.
  73       Complex_Zero : constant Complex_Pack.Complex := (0.0, 0.0);
  74
  75       TC_Complex   : Complex_Pack.Complex := Complex_Zero;
  76       TC_Imaginary : Complex_Pack.Imaginary;
  77
  78    begin
  79
  80       -- Check that the procedures Set_Re and Set_Im (both versions) provide
  81       -- correct results.
  82
  83       declare
  84          TC_Complex_Real_Field : Complex_Pack.Complex := (5.0, 0.0);
  85          TC_Complex_Both_Fields : Complex_Pack.Complex := (5.0, 7.0);
  86       begin
  87
  88          Complex_Pack.Set_Re(X => TC_Complex, Re => 5.0);
  89
  90          if TC_Complex /= TC_Complex_Real_Field then
  91             Report.Failed("Incorrect results from Procedure Set_Re");
  92          end if;
  93
  94          Complex_Pack.Set_Im(X => TC_Complex, Im => 7.0);
  95
  96          if TC_Complex.Re /= 5.0 or
  97             TC_Complex.Im /= 7.0 or
  98             TC_Complex    /= TC_Complex_Both_Fields
  99          then
 100             Report.Failed("Incorrect results from Procedure Set_Im " &
 101                           "with Complex argument");
 102          end if;
 103
 104          Complex_Pack.Set_Im(X => TC_Imaginary, Im => 3.0);
 105
 106
 107          if Complex_Pack.Im(TC_Imaginary) /= 3.0 then
 108             Report.Failed("Incorrect results returned following the use " &
 109                           "of Procedure Set_Im with Imaginary argument");
 110          end if;
 111
 112       end;
 113
 114
 115       -- Check that the functions Re and Im (both versions) provide
 116       -- correct results.
 117
 118       declare
 119          TC_Complex_1 : Complex_Pack.Complex := (1.0, 0.0);
 120          TC_Complex_2 : Complex_Pack.Complex := (0.0, 2.0);
 121          TC_Complex_3 : Complex_Pack.Complex := (4.0, 3.0);
 122       begin
 123
 124          -- Function Re.
 125
 126          if Complex_Pack.Re(X => TC_Complex_1) /= 1.0  or
 127             Complex_Pack.Re(X => TC_Complex_2) /= 0.0  or
 128             Complex_Pack.Re(X => TC_Complex_3) /= 4.0
 129          then
 130             Report.Failed("Incorrect results from Function Re");
 131          end if;
 132
 133          -- Function Im; version with Complex argument.
 134
 135          if Complex_Pack.Im(X => TC_Complex_1) /= 0.0  or
 136             Complex_Pack.Im(X => TC_Complex_2) /= 2.0  or
 137             Complex_Pack.Im(X => TC_Complex_3) /= 3.0
 138          then
 139             Report.Failed("Incorrect results from Function Im " &
 140                           "with Complex argument");
 141          end if;
 142
 143
 144          -- Function Im; version with Imaginary argument.
 145
 146          if Complex_Pack.Im(Complex_Pack.i) /= 1.0  or
 147             Complex_Pack.Im(Complex_Pack.j) /= 1.0
 148          then
 149             Report.Failed("Incorrect results from use of Function Im " &
 150                           "when used with an Imaginary argument");
 151          end if;
 152
 153       end;
 154
 155
 156       -- Verify the results of the three versions of Function
 157       -- Compose_From_Cartesian
 158
 159       declare
 160
 161          Zero  : constant Real_Type := 0.0;
 162          Six   : constant Real_Type := 6.0;
 163
 164          TC_Complex_1 : Complex_Pack.Complex := (3.0, 8.0);
 165          TC_Complex_2 : Complex_Pack.Complex := (Six, Zero);
 166          TC_Complex_3 : Complex_Pack.Complex := (Zero, 1.0);
 167
 168       begin
 169
 170          TC_Complex := Complex_Pack.Compose_From_Cartesian(3.0, 8.0);
 171
 172          if TC_Complex /= TC_Complex_1 then
 173             Report.Failed("Incorrect results from Function " &
 174                           "Compose_From_Cartesian - 1");
 175          end if;
 176
 177          -- If only one component is given, the other component is
 178          -- implicitly zero (Both components are set by the following two
 179          -- function calls).
 180
 181          TC_Complex := Complex_Pack.Compose_From_Cartesian(Re => 6.0);
 182
 183          if TC_Complex /= TC_Complex_2 then
 184             Report.Failed("Incorrect results from Function " &
 185                           "Compose_From_Cartesian - 2");
 186          end if;
 187
 188          TC_Complex :=
 189            Complex_Pack.Compose_From_Cartesian(Im => Complex_Pack.i);
 190
 191          if TC_Complex /= TC_Complex_3 then
 192             Report.Failed("Incorrect results from Function " &
 193                           "Compose_From_Cartesian - 3");
 194          end if;
 195
 196       end;
 197
 198
 199       -- Verify the results of Function Compose_From_Polar, Modulus, "abs",
 200       -- and Argument.
 201
 202       declare
 203
 204          use Complex_Pack;
 205
 206          TC_Modulus,
 207          TC_Argument : Real_Type := 0.0;
 208
 209
 210          Angle_0     : constant Real_Type :=   0.0;
 211          Angle_90    : constant Real_Type :=  90.0;
 212          Angle_180   : constant Real_Type := 180.0;
 213          Angle_270   : constant Real_Type := 270.0;
 214          Angle_360   : constant Real_Type := 360.0;
 215
 216       begin
 217
 218          -- Verify the result of Function Compose_From_Polar.
 219          -- When the value of the parameter Modulus is zero, the
 220          -- Compose_From_Polar function yields a result of zero.
 221
 222          if Compose_From_Polar(0.0, 30.0, 360.0) /= Complex_Zero
 223          then
 224             Report.Failed("Incorrect result from Function " &
 225                           "Compose_From_Polar - 1");
 226          end if;
 227
 228          -- When the value of the parameter Argument is equal to a multiple
 229          -- of the quarter cycle, the result of the Compose_From_Polar
 230          -- function with specified cycle lies on one of the axes.
 231
 232          if Compose_From_Polar( 5.0,  Angle_0,   Angle_360) /= (5.0,  0.0) or
 233             Compose_From_Polar( 5.0,  Angle_90,  Angle_360) /= (0.0,  5.0) or
 234             Compose_From_Polar(-5.0,  Angle_180, Angle_360) /= (5.0,  0.0) or
 235             Compose_From_Polar(-5.0,  Angle_270, Angle_360) /= (0.0,  5.0) or
 236             Compose_From_Polar(-5.0,  Angle_90,  Angle_360) /= (0.0, -5.0) or
 237             Compose_From_Polar( 5.0,  Angle_270, Angle_360) /= (0.0, -5.0)
 238          then
 239             Report.Failed("Incorrect result from Function " &
 240                           "Compose_From_Polar - 2");
 241          end if;
 242
 243          -- When the parameter to Function Argument represents a point on
 244          -- the non-negative real axis, the function yields a zero result.
 245
 246          if Argument(Complex_Zero, Angle_360) /= 0.0 then
 247             Report.Failed("Incorrect result from Function Argument");
 248          end if;
 249
 250          -- Function Modulus
 251
 252          if Modulus(Complex_Zero) /= 0.0  or
 253             Modulus(Compose_From_Polar( 5.0, Angle_90,  Angle_360)) /= 5.0 or
 254             Modulus(Compose_From_Polar(-5.0, Angle_180, Angle_360)) /= 5.0
 255          then
 256             Report.Failed("Incorrect results from Function Modulus");
 257          end if;
 258
 259          -- Function "abs", a rename of Function Modulus.
 260
 261          if "abs"(Complex_Zero) /= 0.0  or
 262             "abs"(Compose_From_Polar( 5.0, Angle_90,  Angle_360)) /= 5.0 or
 263             "abs"(Compose_From_Polar(-5.0, Angle_180, Angle_360)) /= 5.0
 264          then
 265             Report.Failed("Incorrect results from Function abs");
 266          end if;
 267
 268       end;
 269
 270    exception
 271       when others => Report.Failed ("Exception raised in Test_Block");
 272    end Test_Block;
 273
 274    Report.Result;
 275
 276 end CXG1001;