|
DataMuseum.dkPresents historical artifacts from the history of: Rational R1000/400 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about Rational R1000/400 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - download
Length: 3072 (0xc00)
Types: Ada Source
Notes: 03_class, FILE, R1k_Segment, e3_tag, function Kf_Asncy, seg_0130c9, separate Generic_Elementary_Functions
└─⟦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 Kf_Asncy (X, Cycle : Float_Type) return Float_Type is
Y, Sign_X, U, V, One_M_Y, Z1, Z2, Result, Cy : Common_Float;
Zero : constant := 0.0;
One : constant := 1.0;
Two : constant := 2.0;
Four : constant := 4.0;
Half : constant := 0.5;
Two_Pi_Inv : constant := 0.15915_49430_91895_33576_88837_63372_51436_2034;
begin
-- Filter out exceptional cases.
Cy := Common_Float (Cycle);
if Cy <= 0.0 then
raise Argument_Error;
end if;
Y := Common_Float (X);
Sign_X := Copy_Sign (One, Y);
Y := abs (Y);
if Y > One then
raise Argument_Error;
end if;
if Y = One then
Result := Cy / Four;
return (Float_Type (Copy_Sign (Result, Sign_X)));
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) * Two_Pi_Inv) * Cy;
else
Kp_Atn (U, V, Z1, Z2);
Result := Cy * 0.25 - ((Z1 + Z2) * Two_Pi_Inv) * Cy;
end if;
return (Float_Type (Copy_Sign (Result, Sign_X)));
end Kf_Asncy;
nblk1=2
nid=0
hdr6=4
[0x00] rec0=2b rec1=00 rec2=01 rec3=008
[0x01] rec0=10 rec1=00 rec2=02 rec3=000
tail 0x2170e742e82b1517c8e4d 0x42a00066462061e03