Return to tanh87.s CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm / i387

Annotation of coherent/b/lib/libm/i387/tanh87.s, revision 1.1.1.1

1.1       root        1: //////////
                      2: / /usr/src/libm/i387/tanh87.s
                      3: / i387 mathematics library
                      4: / tanh(d)
                      5: //////////
                      6: 
                      7: RASIZE =       4
                      8: 
                      9:        .globl  tanh
                     10:        .globl  _exp
                     11:        .globl  _tstcc
                     12: 
                     13: //////////
                     14: / double
                     15: / tanh(d)
                     16: / double d;
                     17: /
                     18: / tanh(d) = (exp(d) - exp(-d)) / (exp(d) + exp(-d))
                     19: /        = (1 - exp(-2 * d)) / (1 + exp(-2 * d))       (used for d > 0)
                     20: /        = - (1 - exp(2 * d)) / (1 + exp(2 * d))       (used for d < 0)
                     21: //////////
                     22: 
                     23: d      =       RASIZE          / arg offset
                     24: 
                     25: tanh:
                     26:        fldl    d(%esp)         / Load argument d.
                     27: 
                     28:        call    _tstcc          / Check sign of d.
                     29:        je      ?1              / d is zero, just return it.
                     30:        jb      ?0              / Jump if negative.
                     31:        fchs                    / Change positive to negative.
                     32: 
                     33: ?0:
                     34:        pushf                   / Save the flags.
                     35:        fadd    %st, %st        / 2 * d
                     36:        call    _exp            / exp(2 * d)
                     37:        fld1                    / 1, exp(2 * d)
                     38:        fsub    %st, %st(1)     / 1 - exp(2 * d), exp(2 * d)
                     39:        fld1                    / 1, 1 - exp(2 * d), exp(2 * d)
                     40:        faddp   %st(2), %st     / 1 - exp(2 * d), 1 + exp(2 * d)
                     41:        fdiv                    / (1 -  exp(2 * d)) / (1 + exp(2 * d))
                     42:        popf                    / Restore the flags.
                     43:        ja      ?1              / d was positive, done.
                     44:        fchs                    / d was negative, flip result sign.
                     45: 
                     46: ?1:
                     47:        ret
                     48: 
                     49: / end of tanh87.s

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.