Return to sincos.s CVS log

Up to [CSRG BSD Unix] / 43BSDReno / lib / libm / vax

Annotation of 43BSDReno/lib/libm/vax/sincos.s, revision 1.1.1.1

1.1       root        1: # Copyright (c) 1985 Regents of the University of California.
                      2: # All rights reserved.
                      3: #
                      4: # Redistribution and use in source and binary forms are permitted
                      5: # provided that the above copyright notice and this paragraph are
                      6: # duplicated in all such forms and that any documentation,
                      7: # advertising materials, and other materials related to such
                      8: # distribution and use acknowledge that the software was developed
                      9: # by the University of California, Berkeley.  The name of the
                     10: # University may not be used to endorse or promote products derived
                     11: # from this software without specific prior written permission.
                     12: # THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
                     13: # IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
                     14: # WARRANTIES OF MERCHANTIBILITY AND FITNESS FOR A PARTICULAR PURPOSE.
                     15: #
                     16: # All recipients should regard themselves as participants in an ongoing
                     17: # research project and hence should feel obligated to report their
                     18: # experiences (good or bad) with these elementary function codes, using
                     19: # the sendbug(8) program, to the authors.
                     20: #
                     21: #      @(#)sincos.s    5.3 (Berkeley) 6/30/88
                     22: #
                     23:        .data
                     24:        .align  2
                     25: _sccsid:
                     26: .asciz "@(#)sincos.s   1.1 (Berkeley) 8/21/85; 5.3 (ucb.elefunt) 6/30/88"
                     27: 
                     28: #  This is the implementation of Peter Tang's double precision  
                     29: #  sine and cosine for the VAX using Bob Corbett's argument reduction.
                     30: #  
                     31: #  Notes:
                     32: #       under 1,024,000 random arguments testing on [0,2*pi] 
                     33: #       sin() observed maximum error = 0.814 ulps
                     34: #       cos() observed maximum error = 0.792 ulps
                     35: #
                     36: # double sin(arg)
                     37: # double arg;
                     38: # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
                     39: # S. McDonald, April 4,  1985
                     40: #
                     41:        .globl  _sin
                     42:        .text
                     43:        .align  1
                     44: 
                     45: _sin:  .word   0xffc           # save r2-r11
                     46:        movq    4(ap),r0
                     47:        bicw3   $0x807f,r0,r2
                     48:        beql    1f              # if x is zero or reserved operand then return x
                     49: #
                     50: # Save the PSL's IV & FU bits on the stack.
                     51: #
                     52:        movpsl  r2
                     53:        bicw3   $0xff9f,r2,-(sp)
                     54: #
                     55: # Clear the IV & FU bits.
                     56: #
                     57:        bicpsw  $0x0060
                     58: #
                     59: #  Entered by  sine    ; save  0  in  r4 .
                     60: #
                     61:        jsb     libm$argred
                     62:        movl    $0,r4
                     63:        jsb     libm$sincos
                     64:        bispsw  (sp)+
                     65: 1:     ret
                     66: 
                     67: #
                     68: # double cos(arg)
                     69: # double arg;
                     70: # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
                     71: # S. McDonald, April 4,  1985
                     72: #
                     73:        .globl  _cos
                     74:        .text
                     75:        .align  1
                     76: 
                     77: _cos:  .word   0xffc           # save r2-r11
                     78:        movq    4(ap),r0
                     79:        bicw3   $0x7f,r0,r2
                     80:        cmpw    $0x8000,r2
                     81:        beql    1f              # if x is reserved operand then return x
                     82: #
                     83: # Save the PSL's IV & FU bits on the stack.
                     84: #
                     85:        movpsl  r2
                     86:        bicw3   $0xff9f,r2,-(sp)
                     87: #
                     88: # Clear the IV & FU bits.
                     89: #
                     90:        bicpsw  $0x0060
                     91: #
                     92: #  Entered by  cosine  ; save  1  in  r4 .
                     93: #
                     94:        jsb     libm$argred
                     95:        movl    $1,r4
                     96:        jsb     libm$sincos
                     97:        bispsw  (sp)+
                     98: 1:     ret

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