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 - metrics - download

⟦d919381a6⟧ Ada Source

    Length: 3072 (0xc00)
    Types: Ada Source
    Notes: 03_class, FILE, R1k_Segment, e3_tag, function Kf_Cotcy, seg_0130cf, 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 Kf_Cotcy (X, Cycle : Float_Type) return Float_Type is

-- On input, X is a floating-point value in Float_Type; Cycle is a
-- user specificed period.
-- On output, the value of cot(X*2pi/Cycle) is returned.

   R, Y, Sign_Y, Z, User_Period, P_Lead, P_Trail, K : Common_Float;
   U, V : Common_Float;
   Zero : Common_Float := 0.0;
   I_Flag, Sc_Flag : Quadrant_Flag;
   Two_Pi : constant := 6.28318_53071_79586_47692_52867_66559_00576_8394;

begin

-- Argument reduction by  "rem"

   User_Period := Common_Float (Cycle);

   if User_Period <= 0.0 then
      raise Argument_Error;
   end if;

   if (X = 0.0) then
      raise Constraint_Error;
   end if;

   Y := Common_Float (X);
   if (Y > 0.0) then
      Sign_Y := 1.0;
   else
      Sign_Y := -1.0;
      Y      := abs (Y);
   end if;

   R := Y rem User_Period;

   K      := Round (4.0 * R / User_Period);
   I_Flag := Quadrant_Flag (Integer (K) mod 4);

   Sc_Flag := I_Flag mod 2;

   P_Lead  := Leading_Part (User_Period, 4);
   P_Trail := User_Period - P_Lead;

   K := K * 0.25;

   R := (R - K * P_Lead) - K * P_Trail;
   R := (R / User_Period) * Two_Pi;

   U := Kf_Sin (R, Zero);
   V := Kf_Cos (R, Zero);

   if (Sc_Flag = 1) then
      if (V = 0.0) then
         raise Constraint_Error;
      end if;
      Z := -U / V;
   else
      if (U = 0.0) then
         raise Constraint_Error;
      end if;
      Z := V / U;
   end if;

   return (Float_Type (Sign_Y * Z));

end Kf_Cotcy;

    nblk1=2
    nid=0
    hdr6=4
        [0x00] rec0=29 rec1=00 rec2=01 rec3=030
        [0x01] rec0=20 rec1=00 rec2=02 rec3=000
    tail 0x2170e743a82b15193d57b 0x42a00066462061e03