// Since fW_2TO56_RSH is scaled by 2^56, it must be multiplied by 2^-56// before the shift constant 1.10000 * 2^63 is subtracted to yield fNfloat.// Thus, fNfloat contains the floating point version of N{ .mfi      ldfpd           fP3, fP2  = [rAD_P](p14) fcmp.gt.unc.s1  p14,p0 = fAbsX,fMAX_DBL_NORM_ARG      nop.i           0}{ .mfb      nop.m           0      fms.s1          fNfloat = fW_2TO56_RSH, f2TOM56, fRSHF(p15) br.cond.spnt    SINH_CERTAIN_OVERFLOW};;{ .mfi      getf.sig        rN        = fW_2TO56_RSH      nop.f           0      mov             rExp_bias_minus_1 = 0xfffe};;// rIndex_1 has index_1// rIndex_2_16 has index_2 * 16// rBiased_M has M// rM has true M// r = x - Nfloat * ln2_by_128_hi// f = 1 - Nfloat * ln2_by_128_lo{ .mfi      and             rIndex_1 = 0x0f, rN      fnma.s1         fR   = fNfloat, fLn2_by_128_hi, fNormX      shr             rM = rN,  0x7}{ .mfi      and             rIndex_2_16 = 0x70, rN      fnma.s1         fF   = fNfloat, fLn2_by_128_lo, f1      sub             rN_neg = r0, rN};;{ .mmi      and             rIndex_1_neg = 0x0f, rN_neg      add             rBiased_M = rExp_bias_minus_1, rM      shr             rM_neg = rN_neg,  0x7}{ .mmi      and             rIndex_2_16_neg = 0x70, rN_neg      add             rAD_T2 = rAD_TB2, rIndex_2_16      shladd          rAD_T1 = rIndex_1, 4, rAD_TB1};;// rAD_T1 has address of T1// rAD_T2 has address if T2{ .mmi      setf.exp        f2M = rBiased_M      ldfe            fT2  = [rAD_T2]      nop.i           0}{ .mmi      add             rBiased_M_neg = rExp_bias_minus_1, rM_neg      add             rAD_T2_neg = rAD_TB2, rIndex_2_16_neg      shladd          rAD_T1_neg = rIndex_1_neg, 4, rAD_TB1};;// Create Scale = 2^M// Load T1 and T2{ .mmi      ldfe            fT1  = [rAD_T1]      nop.m           0      nop.i           0}{ .mmf      setf.exp        f2M_neg = rBiased_M_neg      ldfe            fT2_neg  = [rAD_T2_neg]      fma.s1          fF_neg   = fNfloat, fLn2_by_128_lo, f1};;{ .mfi      nop.m           0      fma.s1          fRsq = fR, fR, f0      nop.i           0}{ .mfi      ldfe            fT1_neg  = [rAD_T1_neg]      fma.s1          fP54 = fR, fP5, fP4      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fP32 = fR, fP3, fP2      nop.i           0}{ .mfi      nop.m           0      fnma.s1         fP54_neg = fR, fP5, fP4      nop.i           0};;{ .mfi      nop.m           0      fnma.s1         fP32_neg = fR, fP3, fP2      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fP5432  = fRsq, fP54, fP32      nop.i           0}{ .mfi      nop.m           0      fma.s1          fS2  = fF,fT2,f0      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fS1  = f2M,fT1,f0      nop.i           0}{ .mfi      nop.m           0      fma.s1          fP5432_neg  = fRsq, fP54_neg, fP32_neg      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fS1_neg  = f2M_neg,fT1_neg,f0      nop.i           0}{ .mfi      nop.m           0      fma.s1          fS2_neg  = fF_neg,fT2_neg,f0      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fP     = fRsq, fP5432, fR      nop.i           0}{ .mfi      nop.m           0      fma.s1          fS   = fS1,fS2,f0      nop.i           0};;{ .mfi      nop.m           0      fms.s1          fP_neg     = fRsq, fP5432_neg, fR      nop.i           0}{ .mfi      nop.m           0      fma.s1          fS_neg   = fS1_neg,fS2_neg,f0      nop.i           0};;{ .mfb      nop.m           0      fmpy.s0         fTmp = fLn2_by_128_lo, fLn2_by_128_lo // Force inexact(p14) br.cond.spnt    SINH_POSSIBLE_OVERFLOW};;{ .mfi      nop.m           0      fma.s1          fExp = fS, fP, fS      nop.i           0}{ .mfi      nop.m           0      fma.s1          fExp_neg = fS_neg, fP_neg, fS_neg      nop.i           0};;{ .mfb      nop.m           0      fms.d.s0        f8 = fExp, f1, fExp_neg      br.ret.sptk     b0                  // Normal path exit};;// Here if 0 < |x| < 0.25SINH_SMALL:{ .mfi      add             rAD_T1 = 0x1a0, rAD_TB1      fcmp.lt.s1      p7, p8 = fNormX, f0       // Test sign of x      cmp.gt          p6, p0 = -60, rExp_x      // Test |x| < 2^(-60)}{ .mfi      add             rAD_T2 = 0x1d0, rAD_TB1      nop.f           0      nop.i           0};;{ .mmb      ldfe            fA6 = [rAD_T1],16      ldfe            fA5 = [rAD_T2],16(p6)  br.cond.spnt    SINH_VERY_SMALL           // Branch if |x| < 2^(-60)};;{ .mmi      ldfe            fA4 = [rAD_T1],16      ldfe            fA3 = [rAD_T2],16      nop.i           0};;{ .mmi      ldfe            fA2 = [rAD_T1]      ldfe            fA1 = [rAD_T2]      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fX3 = fNormX, fXsq, f0      nop.i           0}{ .mfi      nop.m           0      fma.s1          fX4 = fXsq, fXsq, f0      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fA65 = fXsq, fA6, fA5      nop.i           0}{ .mfi      nop.m           0      fma.s1          fA43 = fXsq, fA4, fA3      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fA21 = fXsq, fA2, fA1      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fA6543 = fX4, fA65, fA43      nop.i           0};;{ .mfi      nop.m           0      fma.s1          fA654321 = fX4, fA6543, fA21      nop.i           0};;// Dummy multiply to generate inexact{ .mfi      nop.m           0      fmpy.s0         fTmp = fA6, fA6      nop.i           0}{ .mfb      nop.m           0      fma.d.s0        f8 = fA654321, fX3, fNormX      br.ret.sptk     b0                // Exit if 2^-60 < |x| < 0.25};;SINH_VERY_SMALL:// Here if 0 < |x| < 2^-60// Compute result by x + sgn(x)*x^2 to get properly rounded result.pred.rel "mutex",p7,p8{ .mfi      nop.m           0(p7)  fnma.d.s0       f8 = fNormX, fNormX, fNormX // If x<0 result ~ x-x^2      nop.i           0}{ .mfb      nop.m           0(p8)  fma.d.s0        f8 = fNormX, fNormX, fNormX // If x>0 result ~ x+x^2      br.ret.sptk     b0                          // Exit if |x| < 2^-60};;SINH_POSSIBLE_OVERFLOW:// Here if fMAX_DBL_NORM_ARG < |x| < fMIN_DBL_OFLOW_ARG// This cannot happen if input is a double, only if input higher precision.// Overflow is a possibility, not a certainty.// Recompute result using status field 2 with user's rounding mode,// and wre set.  If result is larger than largest double, then we have// overflow{ .mfi      mov             rGt_ln  = 0x103ff // Exponent for largest dbl + 1 ulp      fsetc.s2        0x7F,0x42         // Get user's round mode, set wre      nop.i           0};;{ .mfi      setf.exp        fGt_pln = rGt_ln  // Create largest double + 1 ulp      fma.d.s2        fWre_urm_f8 = fS, fP, fS    // Result with wre set      nop.i           0};;{ .mfi      nop.m           0      fsetc.s2        0x7F,0x40                   // Turn off wre in sf2      nop.i           0};;{ .mfi      nop.m           0      fcmp.ge.s1      p6, p0 =  fWre_urm_f8, fGt_pln // Test for overflow      nop.i           0};;{ .mfb      nop.m           0      nop.f           0(p6)  br.cond.spnt    SINH_CERTAIN_OVERFLOW // Branch if overflow};;{ .mfb      nop.m           0      fma.d.s0        f8 = fS, fP, fS      br.ret.sptk     b0                     // Exit if really no overflow};;SINH_CERTAIN_OVERFLOW:{ .mfi      sub             rTmp = rExp_mask, r0, 1      fcmp.lt.s1      p6, p7 = fNormX, f0    // Test for x < 0      nop.i           0};;{ .mmf      alloc           r32=ar.pfs,1,4,4,0      setf.exp        fTmp = rTmp      fmerge.s        FR_X = f8,f8};;{ .mfi      mov             GR_Parameter_TAG = 127(p6)  fnma.d.s0       FR_RESULT = fTmp, fTmp, f0    // Set I,O and -INF result      nop.i           0}{ .mfb      nop.m           0(p7)  fma.d.s0        FR_RESULT = fTmp, fTmp, f0    // Set I,O and +INF result      br.cond.sptk    __libm_error_region};;// Here if x unormSINH_UNORM:{ .mfb      getf.exp        rSignexp_x = fNormX    // Must recompute if x unorm      fcmp.eq.s0      p6, p0 = f8, f0        // Set D flag      br.cond.sptk    SINH_COMMON};;GLOBAL_IEEE754_END(sinh)LOCAL_LIBM_ENTRY(__libm_error_region).prologue{ .mfi        add   GR_Parameter_Y=-32,sp             // Parameter 2 value        nop.f 0.save   ar.pfs,GR_SAVE_PFS        mov  GR_SAVE_PFS=ar.pfs                 // Save ar.pfs}{ .mfi.fframe 64        add sp=-64,sp                           // Create new stack        nop.f 0        mov GR_SAVE_GP=gp                       // Save gp};;{ .mmi        stfd [GR_Parameter_Y] = FR_Y,16         // STORE Parameter 2 on stack        add GR_Parameter_X = 16,sp              // Parameter 1 address.save   b0, GR_SAVE_B0        mov GR_SAVE_B0=b0                       // Save b0};;.body{ .mib        stfd [GR_Parameter_X] = FR_X            // STORE Parameter 1 on stack        add   GR_Parameter_RESULT = 0,GR_Parameter_Y  // Parameter 3 address        nop.b 0}{ .mib        stfd [GR_Parameter_Y] = FR_RESULT       // STORE Parameter 3 on stack        add   GR_Parameter_Y = -16,GR_Parameter_Y        br.call.sptk b0=__libm_error_support#   // Call error handling function};;{ .mmi        add   GR_Parameter_RESULT = 48,sp        nop.m 0        nop.i 0};;{ .mmi        ldfd  f8 = [GR_Parameter_RESULT]       // Get return result off stack.restore sp        add   sp = 64,sp                       // Restore stack pointer        mov   b0 = GR_SAVE_B0                  // Restore return address};;{ .mib        mov   gp = GR_SAVE_GP                  // Restore gp        mov   ar.pfs = GR_SAVE_PFS             // Restore ar.pfs        br.ret.sptk     b0                     // Return};;LOCAL_LIBM_END(__libm_error_region).type   __libm_error_support#,@function.global __libm_error_support#