|
DataMuseum.dkPresents historical artifacts from the history of: Rational R1000/400 Tapes |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about Rational R1000/400 Tapes Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - downloadIndex: ┃ B T ┃
Length: 1109 (0x455)
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⟧
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;