DataMuseum.dk Presents historical artifacts from the history of: Rational R1000/400

This is an automatic "excavation" of a thematic subset of artifacts from Datamuseum.dk's BitArchive. See our Wiki for more about Rational R1000/400 Excavated with: AutoArchaeologist - Free & Open Source Software.

top - download

⟦6f60d0749⟧ Ada Source

    Length: 3072 (0xc00)
    Types: Ada Source
    Notes: 03_class, FILE, R1k_Segment, e3_tag, function Arcsin, seg_0130bf, separate Generic_Elementary_Functions

Derivation

└─⟦8527c1e9b⟧ Bits:30000544 8mm tape, Rational 1000, Arrival backup of disks in PAM's R1000
    └─ ⟦cfc2e13cd⟧ »Space Info Vol 2« 
        └─⟦this⟧ 

separate (Generic_Elementary_Functions)

function Arcsin (X : Float_Type) return Float_Type is

   Y, Sign_X, U, V, One_M_Y, Z1, Z2, Result : Common_Float;

   Zero : constant := 0.0;
   One  : constant := 1.0;
   Two  : constant := 2.0;
   Half : constant := 0.5;

   Pi_By_2 : constant := 1.57079_63267_94896_61923_13216_91639_75144_2098;

   Piby2_Lead  : constant Common_Float := 16#1.921#;
   Piby2_Trail : constant Common_Float :=
      16#0.FB544_42D18_46989_8CC51_701B8_39A25_2049C#E-3;

begin

-- Filter out exceptional cases.
   Y      := Common_Float (X);
   Sign_X := Copy_Sign (One, Y);
   Y      := abs (Y);

   if Y > One then
      raise Argument_Error;
   end if;

   V := Y;
   if Y < Half then
      U := One - Y * Y;
   else
      One_M_Y := Half + (Half - Y);
      U       := Two * One_M_Y - One_M_Y * One_M_Y;
   end if;

   U := Kf_Sqrt (U);

   if U >= V then
      Kp_Atn (V, U, Z1, Z2);
      Result := Z1 + Z2;
   else
      Kp_Atn (U, V, Z1, Z2);
      Result := (Piby2_Lead - Z1) + (Piby2_Trail - Z2);
   end if;

   return (Float_Type (Copy_Sign (Result, Sign_X)));

end Arcsin;

    nblk1=2
    nid=0
    hdr6=4
        [0x00] rec0=2a rec1=00 rec2=01 rec3=028
        [0x01] rec0=09 rec1=00 rec2=02 rec3=000
    tail 0x2170e741a82b1516c4d09 0x42a00066462061e03