# method: tanh(arg)=sinh(arg)/cosh(arg); if |arg|<0.5 sinh(9 is computed from
# a polynomial approx otherwise from exp().
# Only 1 call to exp() is made.
# J F Jarvis, August 17, 1978
divd3 r8,$0d1.0e+0,-(sp) # exp(-|arg|)
muld3 r6,r6,r0 # |arg|<0.5, sinh() from poly approx
t4: # r0,1 = 2*sinh(|arg|)
addd3 (sp),r8,r2 # r2,3 = cosh(|arg|);
divd2 r2,r0 # r0,1 = tanh(|arg|)
xorw2 r10,r0 # r0,1 = tanh(arg)
.double 0d0.251726188251e-7
.double 0d0.275569807356154e-5
.double 0d0.1984127027907999e-3
.double 0d0.833333333307759961e-2
.double 0d0.16666666666667212324e+0
.double 0d0.99999999999999998116e+0