Return to sincos.s CVS log

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

Annotation of 43BSD/usr.lib/libm/VAX/sincos.s, revision 1.1.1.1

1.1       root        1: # 
                      2: # Copyright (c) 1985 Regents of the University of California.
                      3: # 
                      4: # Use and reproduction of this software are granted  in  accordance  with
                      5: # the terms and conditions specified in  the  Berkeley  Software  License
                      6: # Agreement (in particular, this entails acknowledgement of the programs'
                      7: # source, and inclusion of this notice) with the additional understanding
                      8: # that  all  recipients  should regard themselves as participants  in  an
                      9: # ongoing  research  project and hence should  feel  obligated  to report
                     10: # their  experiences (good or bad) with these elementary function  codes,
                     11: # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
                     12: #
                     13: 
                     14: # @(#)sincos.s 1.1 (Berkeley) 8/21/85
                     15: 
                     16: #  This is the implementation of Peter Tang's double precision  
                     17: #  sine and cosine for the VAX using Bob Corbett's argument reduction.
                     18: #  
                     19: #  Notes:
                     20: #       under 1,024,000 random arguments testing on [0,2*pi] 
                     21: #       sin() observed maximum error = 0.814 ulps
                     22: #       cos() observed maximum error = 0.792 ulps
                     23: #
                     24: # double sin(arg)
                     25: # double arg;
                     26: # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
                     27: # S. McDonald, April 4,  1985
                     28: #
                     29:        .globl  _sin
                     30:        .text
                     31:        .align  1
                     32: 
                     33: _sin:  .word   0xffc           # save r2-r11
                     34:        movq    4(ap),r0
                     35:        bicw3   $0x807f,r0,r2
                     36:        beql    1f              # if x is zero or reserved operand then return x
                     37: #
                     38: # Save the PSL's IV & FU bits on the stack.
                     39: #
                     40:        movpsl  r2
                     41:        bicw3   $0xff9f,r2,-(sp)
                     42: #
                     43: # Clear the IV & FU bits.
                     44: #
                     45:        bicpsw  $0x0060
                     46: #
                     47: #  Entered by  sine    ; save  0  in  r4 .
                     48: #
                     49:        jsb     libm$argred
                     50:        movl    $0,r4
                     51:        jsb     libm$sincos
                     52:        bispsw  (sp)+
                     53: 1:     ret
                     54: 
                     55: #
                     56: # double cos(arg)
                     57: # double arg;
                     58: # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
                     59: # S. McDonald, April 4,  1985
                     60: #
                     61:        .globl  _cos
                     62:        .text
                     63:        .align  1
                     64: 
                     65: _cos:  .word   0xffc           # save r2-r11
                     66:        movq    4(ap),r0
                     67:        bicw3   $0x7f,r0,r2
                     68:        cmpw    $0x8000,r2
                     69:        beql    1f              # if x is reserved operand then return x
                     70: #
                     71: # Save the PSL's IV & FU bits on the stack.
                     72: #
                     73:        movpsl  r2
                     74:        bicw3   $0xff9f,r2,-(sp)
                     75: #
                     76: # Clear the IV & FU bits.
                     77: #
                     78:        bicpsw  $0x0060
                     79: #
                     80: #  Entered by  cosine  ; save  1  in  r4 .
                     81: #
                     82:        jsb     libm$argred
                     83:        movl    $1,r4
                     84:        jsb     libm$sincos
                     85:        bispsw  (sp)+
                     86: 1:     ret

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