Return to sin.s CVS log

Up to [CSRG BSD Unix] / 43BSD / usr.lib / liboldnm

Annotation of 43BSD/usr.lib/liboldnm/sin.s, revision 1.1.1.1

1.1       root        1: #
                      2: # Copyright (c) 1980 Regents of the University of California.
                      3: # All rights reserved.  The Berkeley software License Agreement
                      4: # specifies the terms and conditions for redistribution.
                      5: #
                      6: #      @(#)sin.s       5.1 (Berkeley) 5/8/85
                      7: # 
                      8: #
                      9: # double float sin(),cos()
                     10: #coefficients are from Hart & Cheney  SIN3370 (18.80D)
                     11: .text
                     12: .align 1
                     13: .globl _sin
                     14: .globl _cos
                     15: _sin:
                     16:        .word 0x03c0
                     17:        bispsw  $0xe0
                     18:        clrl    r9      
                     19:        movd    4(ap),r0
                     20:        jgeq    range
                     21:        movl    $2,r9
                     22:        jbr     negarg
                     23: 
                     24: .align 1
                     25: _cos:
                     26:        .word   0x0fc0
                     27:        bispsw  $0xe0
                     28:        movl    $1,r9
                     29:        movd    4(ap),r0
                     30:        jgeq    range
                     31: negarg:
                     32:        mnegd   r0,r0
                     33: 
                     34: range:
                     35:        emodd   twoopi,twoopix,r0,r8,r6
                     36:        addl2   r9,r8
                     37:        bicb2   $~03,r8
                     38:        caseb   r8,$0,$3
                     39: bse:   .word   qda-bse,qdb-bse,qdc-bse,qdd-bse
                     40: qdb:   subd3   r6,$0d1.0e+0,r6
                     41:        jbr     qda
                     42: qdc:   mnegd   r6,r6
                     43:        jbr     qda
                     44: qdd:   subd2   $0d1.0e+0,r6
                     45: qda:
                     46:        muld3   r6,r6,r8
                     47:        polyd   r8,$4,pcoef
                     48:        muld2   r0,r6
                     49:        polyd   r8,$4,qcoef
                     50:        divd3   r0,r6,r0
                     51:        ret
                     52: .data
                     53: .align 2
                     54: twoopi:.double 0d6.36619772367581340000e-01
                     55: qcoef:.double  0d1.00000000000000000000e+00
                     56:  .double       0d1.32653490878613630000e+02
                     57:  .double       0d9.46309610153820810000e+03
                     58:  .double       0d4.08179225234329970000e+05
                     59:  .double       0d8.64455865292253430000e+06
                     60: pcoef:.double  0d1.45968840666576870000e+02
                     61:  .double       0d-1.38472724998245280000e+04
                     62:  .double       0d4.40103053537526640000e+05
                     63:  .double       0d-4.94290810090284410000e+06
                     64:  .double       0d1.35788409787737560000e+07
                     65: twoopix: .byte 0x29

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