Lines Matching +full:x +full:- +full:z
8 | Input: Double-extended number X in location pointed to
11 | Output: The value sinh(X) returned in floating-point register Fp0.
23 | 1. If |X| > 16380 log2, go to 3.
25 | 2. (|X| <= 16380 log2) Sinh(X) is obtained by the formulae
26 | y = |X|, sgn = sign(X), and z = expm1(Y),
27 | sinh(X) = sgn*(1/2)*( z + z/(1+z) ).
30 | 3. If |X| > 16480 log2, go to 5.
32 | 4. (16380 log2 < |X| <= 16480 log2)
33 | sinh(X) = sign(X) * exp(|X|)/2.
34 | However, invoking exp(|X|) may cause premature overflow.
35 | Thus, we calculate sinh(X) as follows:
36 | Y := |X|
37 | sgn := sign(X)
39 | Y' := Y - 16381 log2
40 | sinh(X) := sgnFact * exp(Y').
43 | 5. (|X| > 16480 log2) sinh(X) must overflow. Return
44 | sign(X)*Huge*Huge to generate overflow and an infinity with
70 |--SINH(X) = X FOR DENORMALIZED X
85 |--THIS IS THE USUAL CASE, |X| < 16380 LOG2
86 |--Y = |X|, Z = EXPM1(Y), SINH(X) = SIGN(X)*(1/2)*( Z + Z/(1+Z) )
88 fabsx %fp0 | ...Y = |X|
90 moveml %a1/%d1,-(%sp)
91 fmovemx %fp0-%fp0,(%a0)
93 bsr setoxm1 | ...FP0 IS Z = EXPM1(Y)
98 fadds #0x3F800000,%fp1 | ...1+Z
99 fmovex %fp0,-(%sp)
100 fdivx %fp1,%fp0 | ...Z/(1+Z)
105 movel %d0,-(%sp)
108 fmuls (%sp)+,%fp0 |last fp inst - possible exceptions set
116 fsubd T1(%pc),%fp0 | ...(|X|-16381LOG2_LEAD)
117 movel #0,-(%sp)
118 movel #0x80000000,-(%sp)
122 movel %d0,-(%sp) | ...EXTENDED FMT
123 fsubd T2(%pc),%fp0 | ...|X| - 16381 LOG2, ACCURATE
125 movel %d1,-(%sp)
127 fmovemx %fp0-%fp0,(%a0)