Return to cosh87.m CVS log

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

Annotation of coherent/b/lib/libm/i8087/cosh87.m, revision 1.1

1.1     ! root        1: //////////
        !             2: / libm 8087
        !             3: / cosh(d), sinh(d)
        !             4: //////////
        !             5: 
        !             6: #include "larges.h"
        !             7: #include "ifno8087.h"
        !             8: 
        !             9:        .globl  cosh_
        !            10:        .globl  sinh_
        !            11:        .globl  exp
        !            12: 
        !            13: //////////
        !            14: / double
        !            15: / cosh(d)
        !            16: / double d;
        !            17: /
        !            18: / double
        !            19: / sinh(d)
        !            20: / double d;
        !            21: /
        !            22: / cosh(d) = (exp(d) + 1/exp(d))/2
        !            23: / sinh(d) = (exp(d) - 1/exp(d))/2
        !            24: //////////
        !            25: 
        !            26: d      =       RASIZE          / arg offset
        !            27: 
        !            28: cosh_:
        !            29:        ifno8087(_cosh_)
        !            30:        clc                     / Clear Carry to indicate cosh.
        !            31:        jmp     0f              / Go to common code.
        !            32: 
        !            33: sinh_:
        !            34:        ifno8087(_sinh_)
        !            35:        stc                     / Set Carry to indicate sinh.
        !            36: 
        !            37: 0:
        !            38:        mov     bx, sp
        !            39:        fdld    Pss d(bx)       / Load argument d.
        !            40: 
        !            41:        pushf                   / Save the Carry flag.
        !            42:        Gcall   exp             / exp(d)
        !            43:        fld     st              / exp(d), exp(d)
        !            44:        fld1                    / 1, exp(d), exp(d)
        !            45:        fdiv    st,st1          / 1/exp(d), exp(d)
        !            46:        popf                    / Restore the flags.
        !            47:        jnc     1f              / Jump if cosh.
        !            48:        fchs                    / -1/exp(d), exp(d) for cosh
        !            49: 
        !            50: 1:     fadd                    / exp(d) +- 1/exp(d)
        !            51:        fidiv   cs:two          / (exp(d) +- 1/exp(d)) / 2
        !            52:        Gret
        !            53: 
        !            54: //////////
        !            55: / data
        !            56: //////////
        !            57: two:   .word   2

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