DataMuseum.dk

top - download

    Length: 857 (0x359)
    Types: TextFile
    Names: »B«

└─⟦5f3412b64⟧ Bits:30000745 8mm tape, Rational 1000, ENVIRONMENT 12_6_5 TOOLS 
    └─ ⟦91c658230⟧ »DATA« 
        └─⟦458657fb6⟧ 
            └─⟦1472c4407⟧ 
                └─⟦this⟧ 
└─⟦d10a02448⟧ Bits:30000409 8mm tape, Rational 1000, ENVIRONMENT, D_12_7_3
    └─ ⟦fc9b38f02⟧ »DATA« 
        └─⟦9b46a407a⟧ 
            └─⟦2e03b931c⟧ 
                └─⟦this⟧ 

TextFile

separate (Generic_Elementary_Functions)

function Arccos (X : Float_Type) return Float_Type is

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

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

   Pi_Lead  : constant Common_Float := 16#3.243#;
   Pi_Trail : constant Common_Float :=
      16#0.F6A88_85A30_8D313_198A2_E0370_7344A_40938#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 := (One - Y) / (One + Y);
   V := Kf_Sqrt (V);
   U := One;
   Kp_Atn (V, U, Z1, Z2);
   Z1 := Z1 + Z1;
   Z2 := Z2 + Z2;


   if Sign_X = One then
      Result := Z1 + Z2;
   else
      Result := (Pi_Lead - Z1) + (Pi_Trail - Z2);
   end if;

   return (Float_Type (Result));

end Arccos;