Return to cos87.m CVS log

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

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

1.1     ! root        1: //////////
        !             2: / libm 8087
        !             3: / cos(d), sin(d)
        !             4: //////////
        !             5: 
        !             6: #include "larges.h"
        !             7: #include "ifno8087.h"
        !             8: 
        !             9:        .globl  cos_
        !            10:        .globl  sin_
        !            11:        .globl  tan
        !            12: 
        !            13: //////////
        !            14: / double
        !            15: / cos(angle)
        !            16: / double angle;
        !            17: /
        !            18: / double
        !            19: / sin(angle)
        !            20: / double angle;
        !            21: //////////
        !            22: 
        !            23: d      =       RASIZE          / arg offset
        !            24: 
        !            25: cos_:
        !            26:        ifno8087(_cos_)
        !            27:        clc                     / Clear Carry for cos.
        !            28:        jmp     0f
        !            29: 
        !            30: sin_:
        !            31:        ifno8087(_sin_)
        !            32:        stc                     / Set Carry for sin.
        !            33: 
        !            34: 0:     mov     bx, sp
        !            35:        fdld    Pss d(bx)       / Load argument d.
        !            36: 
        !            37:        pushf                   / Save the flags.
        !            38:        Gcall   tan             / Get x, y.
        !            39:        fld     st              / x, x, y
        !            40:        fmul    st,st           / x*x, x, y
        !            41:        fld     st2             / y, x*x, x, y
        !            42:        fmul    st,st           / y*y, x*x, x, y
        !            43:        fadd                    / x*x + y*y, x, y
        !            44:        fsqrt                   / r = sqrt(x*x + y*y), x, y
        !            45:        popf                    / Restore the flags.
        !            46:        jc      1f              / Carry set for sin.
        !            47:        fdiv                    / cos(d) = x / r, y.
        !            48:        fstp    st1             / Pop the stack, return result.
        !            49:        Gret
        !            50: 
        !            51: 1:     fstp    st1             / r, y.
        !            52:        fdiv                    / sin(d) = y / r
        !            53:        Gret

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