Return to sin.s CVS log

Up to [Research Unix] / researchv10dc / libnm

Annotation of researchv10dc/libnm/sin.s, revision 1.1

1.1     ! root        1: # double float sin(),cos()
        !             2: #coefficients are from Hart & Cheney  SIN3370 (18.80D)
        !             3: .text
        !             4: .align 1
        !             5: .globl _sin
        !             6: .globl _cos
        !             7: _sin:
        !             8:        .word 0x03c0
        !             9:        clrl    r9      
        !            10:        movd    4(ap),r0
        !            11:        jgeq    range
        !            12:        movl    $2,r9
        !            13:        jbr     negarg
        !            14: 
        !            15: .align 1
        !            16: _cos:
        !            17:        .word   0x0fc0
        !            18:        movl    $1,r9
        !            19:        movd    4(ap),r0
        !            20:        jgeq    range
        !            21: negarg:
        !            22:        mnegd   r0,r0
        !            23: 
        !            24: range:
        !            25:        emodd   twoopi,twoopix,r0,r8,r6
        !            26:        addl2   r9,r8
        !            27:        bicb2   $~03,r8
        !            28:        caseb   r8,$0,$3
        !            29: bse:   .word   qda-bse,qdb-bse,qdc-bse,qdd-bse
        !            30: qdb:   subd3   r6,$0d1.0e+0,r6
        !            31:        jbr     qda
        !            32: qdc:   mnegd   r6,r6
        !            33:        jbr     qda
        !            34: qdd:   subd2   $0d1.0e+0,r6
        !            35: qda:
        !            36:        muld3   r6,r6,r8
        !            37:        polyd   r8,$4,pcoef
        !            38:        muld2   r0,r6
        !            39:        polyd   r8,$4,qcoef
        !            40:        divd3   r0,r6,r0
        !            41:        ret
        !            42: .data
        !            43: .align 2
        !            44: twoopi:.double 0d6.36619772367581340000e-01
        !            45: qcoef:.double  0d1.00000000000000000000e+00
        !            46:  .double       0d1.32653490878613630000e+02
        !            47:  .double       0d9.46309610153820810000e+03
        !            48:  .double       0d4.08179225234329970000e+05
        !            49:  .double       0d8.64455865292253430000e+06
        !            50: pcoef:.double  0d1.45968840666576870000e+02
        !            51:  .double       0d-1.38472724998245280000e+04
        !            52:  .double       0d4.40103053537526640000e+05
        !            53:  .double       0d-4.94290810090284410000e+06
        !            54:  .double       0d1.35788409787737560000e+07
        !            55: twoopix: .byte 0x29

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