home *** CD-ROM | disk | FTP | other *** search
/ InfoMagic Source Code 1993 July / THE_SOURCE_CODE_CD_ROM.iso / bsd_srcs / lib / libm / vax / sincos.s < prev    next >
Encoding:
Text File  |  1991-04-12  |  3.2 KB  |  110 lines
  1. # Copyright (c) 1985 Regents of the University of California.
  2. # All rights reserved.
  3. #
  4. # Redistribution and use in source and binary forms, with or without
  5. # modification, are permitted provided that the following conditions
  6. # are met:
  7. # 1. Redistributions of source code must retain the above copyright
  8. #    notice, this list of conditions and the following disclaimer.
  9. # 2. Redistributions in binary form must reproduce the above copyright
  10. #    notice, this list of conditions and the following disclaimer in the
  11. #    documentation and/or other materials provided with the distribution.
  12. # 3. All advertising materials mentioning features or use of this software
  13. #    must display the following acknowledgement:
  14. #    This product includes software developed by the University of
  15. #    California, Berkeley and its contributors.
  16. # 4. Neither the name of the University nor the names of its contributors
  17. #    may be used to endorse or promote products derived from this software
  18. #    without specific prior written permission.
  19. #
  20. # THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  21. # ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  22. # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  23. # ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  24. # FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  25. # DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  26. # OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  27. # HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  28. # LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  29. # OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  30. # SUCH DAMAGE.
  31. #
  32. #    @(#)sincos.s    5.4 (Berkeley) 10/9/90
  33. #
  34.     .data
  35.     .align    2
  36. _sccsid:
  37. .asciz    "@(#)sincos.s    1.1 (Berkeley) 8/21/85; 5.4 (ucb.elefunt) 10/9/90"
  38.  
  39. #  This is the implementation of Peter Tang's double precision  
  40. #  sine and cosine for the VAX using Bob Corbett's argument reduction.
  41. #  
  42. #  Notes:
  43. #       under 1,024,000 random arguments testing on [0,2*pi] 
  44. #       sin() observed maximum error = 0.814 ulps
  45. #       cos() observed maximum error = 0.792 ulps
  46. #
  47. # double sin(arg)
  48. # double arg;
  49. # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
  50. # S. McDonald, April 4,  1985
  51. #
  52.     .globl    _sin
  53.     .text
  54.     .align    1
  55.  
  56. _sin:    .word    0xffc        # save r2-r11
  57.     movq    4(ap),r0
  58.     bicw3    $0x807f,r0,r2
  59.     beql    1f        # if x is zero or reserved operand then return x
  60. #
  61. # Save the PSL's IV & FU bits on the stack.
  62. #
  63.     movpsl    r2
  64.     bicw3    $0xff9f,r2,-(sp)
  65. #
  66. # Clear the IV & FU bits.
  67. #
  68.     bicpsw    $0x0060
  69. #
  70. #  Entered by  sine    ; save  0  in  r4 .
  71. #
  72.     jsb    libm$argred
  73.     movl    $0,r4
  74.     jsb    libm$sincos
  75.     bispsw    (sp)+
  76. 1:    ret
  77.  
  78. #
  79. # double cos(arg)
  80. # double arg;
  81. # method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
  82. # S. McDonald, April 4,  1985
  83. #
  84.     .globl    _cos
  85.     .text
  86.     .align    1
  87.  
  88. _cos:    .word    0xffc        # save r2-r11
  89.     movq    4(ap),r0
  90.     bicw3    $0x7f,r0,r2
  91.     cmpw    $0x8000,r2
  92.     beql    1f        # if x is reserved operand then return x
  93. #
  94. # Save the PSL's IV & FU bits on the stack.
  95. #
  96.     movpsl    r2
  97.     bicw3    $0xff9f,r2,-(sp)
  98. #
  99. # Clear the IV & FU bits.
  100. #
  101.     bicpsw    $0x0060
  102. #
  103. #  Entered by  cosine  ; save  1  in  r4 .
  104. #
  105.     jsb    libm$argred
  106.     movl    $1,r4
  107.     jsb    libm$sincos
  108.     bispsw    (sp)+
  109. 1:    ret
  110.