Return to sin.s CVS log

Up to [CSRG BSD Unix] / 42BSD / usr.lib / libnm

Annotation of 42BSD/usr.lib/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:        bispsw  $0xe0
        !            10:        clrl    r9      
        !            11:        movd    4(ap),r0
        !            12:        jgeq    range
        !            13:        movl    $2,r9
        !            14:        jbr     negarg
        !            15: 
        !            16: .align 1
        !            17: _cos:
        !            18:        .word   0x0fc0
        !            19:        bispsw  $0xe0
        !            20:        movl    $1,r9
        !            21:        movd    4(ap),r0
        !            22:        jgeq    range
        !            23: negarg:
        !            24:        mnegd   r0,r0
        !            25: 
        !            26: range:
        !            27:        emodd   twoopi,twoopix,r0,r8,r6
        !            28:        addl2   r9,r8
        !            29:        bicb2   $~03,r8
        !            30:        caseb   r8,$0,$3
        !            31: bse:   .word   qda-bse,qdb-bse,qdc-bse,qdd-bse
        !            32: qdb:   subd3   r6,$0d1.0e+0,r6
        !            33:        jbr     qda
        !            34: qdc:   mnegd   r6,r6
        !            35:        jbr     qda
        !            36: qdd:   subd2   $0d1.0e+0,r6
        !            37: qda:
        !            38:        muld3   r6,r6,r8
        !            39:        polyd   r8,$4,pcoef
        !            40:        muld2   r0,r6
        !            41:        polyd   r8,$4,qcoef
        !            42:        divd3   r0,r6,r0
        !            43:        ret
        !            44: .data
        !            45: .align 2
        !            46: twoopi:.double 0d6.36619772367581340000e-01
        !            47: qcoef:.double  0d1.00000000000000000000e+00
        !            48:  .double       0d1.32653490878613630000e+02
        !            49:  .double       0d9.46309610153820810000e+03
        !            50:  .double       0d4.08179225234329970000e+05
        !            51:  .double       0d8.64455865292253430000e+06
        !            52: pcoef:.double  0d1.45968840666576870000e+02
        !            53:  .double       0d-1.38472724998245280000e+04
        !            54:  .double       0d4.40103053537526640000e+05
        !            55:  .double       0d-4.94290810090284410000e+06
        !            56:  .double       0d1.35788409787737560000e+07
        !            57: twoopix: .byte 0x29

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