add_n.asm
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:4k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. dnl  SPARC v9 mpn_add_n -- Add two limb vectors of the same length > 0 and
  2. dnl  store sum in a third limb vector.
  3. dnl  Copyright 2001, 2002, 2003 Free Software Foundation, Inc.
  4. dnl  This file is part of the GNU MP Library.
  5. dnl  The GNU MP Library is free software; you can redistribute it and/or modify
  6. dnl  it under the terms of the GNU Lesser General Public License as published
  7. dnl  by the Free Software Foundation; either version 3 of the License, or (at
  8. dnl  your option) any later version.
  9. dnl  The GNU MP Library is distributed in the hope that it will be useful, but
  10. dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  11. dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  12. dnl  License for more details.
  13. dnl  You should have received a copy of the GNU Lesser General Public License
  14. dnl  along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  15. include(`../config.m4')
  16. C    cycles/limb
  17. C UltraSPARC 1&2:     4
  18. C UltraSPARC 3:       4.5
  19. C Compute carry-out from the most significant bits of u,v, and r, where
  20. C r=u+v+carry_in, using logic operations.
  21. C This code runs at 4 cycles/limb on UltraSPARC 1 and 2.  It has a 4 insn
  22. C recurrency, and the UltraSPARC 1 and 2 the IE units are 100% saturated.
  23. C Therefore, it seems futile to try to optimize this any further...
  24. C INPUT PARAMETERS
  25. define(`rp',`%i0')
  26. define(`up',`%i1')
  27. define(`vp',`%i2')
  28. define(`n',`%i3')
  29. define(`u0',`%l0')
  30. define(`u1',`%l2')
  31. define(`u2',`%l4')
  32. define(`u3',`%l6')
  33. define(`v0',`%l1')
  34. define(`v1',`%l3')
  35. define(`v2',`%l5')
  36. define(`v3',`%l7')
  37. define(`cy',`%i4')
  38. define(`fanop',`fitod %f0,%f2') dnl  A quasi nop running in the FA pipe
  39. define(`fmnop',`fmuld %f0,%f0,%f4') dnl  A quasi nop running in the FM pipe
  40. ASM_START()
  41. REGISTER(%g2,#scratch)
  42. REGISTER(%g3,#scratch)
  43. PROLOGUE(mpn_add_n)
  44. save %sp,-160,%sp
  45. fitod %f0,%f0 C make sure f0 contains small, quiet number
  46. subcc n,4,%g0
  47. bl,pn %icc,.Loop0
  48. mov 0,cy
  49. ldx [up+0],u0
  50. ldx [vp+0],v0
  51. add up,32,up
  52. ldx [up-24],u1
  53. ldx [vp+8],v1
  54. add vp,32,vp
  55. ldx [up-16],u2
  56. ldx [vp-16],v2
  57. ldx [up-8],u3
  58. ldx [vp-8],v3
  59. subcc n,8,n
  60. add u0,v0,%g1 C main add
  61. add %g1,cy,%g4 C carry add
  62. or u0,v0,%g2
  63. bl,pn %icc,.Lend4567
  64. fanop
  65. b,a .Loop
  66. .align 16
  67. C START MAIN LOOP
  68. .Loop: andn %g2,%g4,%g2
  69. and u0,v0,%g3
  70. ldx [up+0],u0
  71. fanop
  72. C --
  73. or %g3,%g2,%g2
  74. ldx [vp+0],v0
  75. add up,32,up
  76. fanop
  77. C --
  78. srlx %g2,63,cy
  79. add u1,v1,%g1
  80. stx %g4,[rp+0]
  81. fanop
  82. C --
  83. add %g1,cy,%g4
  84. or u1,v1,%g2
  85. fmnop
  86. fanop
  87. C --
  88. andn %g2,%g4,%g2
  89. and u1,v1,%g3
  90. ldx [up-24],u1
  91. fanop
  92. C --
  93. or %g3,%g2,%g2
  94. ldx [vp+8],v1
  95. add vp,32,vp
  96. fanop
  97. C --
  98. srlx %g2,63,cy
  99. add u2,v2,%g1
  100. stx %g4,[rp+8]
  101. fanop
  102. C --
  103. add %g1,cy,%g4
  104. or u2,v2,%g2
  105. fmnop
  106. fanop
  107. C --
  108. andn %g2,%g4,%g2
  109. and u2,v2,%g3
  110. ldx [up-16],u2
  111. fanop
  112. C --
  113. or %g3,%g2,%g2
  114. ldx [vp-16],v2
  115. add rp,32,rp
  116. fanop
  117. C --
  118. srlx %g2,63,cy
  119. add u3,v3,%g1
  120. stx %g4,[rp-16]
  121. fanop
  122. C --
  123. add %g1,cy,%g4
  124. or u3,v3,%g2
  125. fmnop
  126. fanop
  127. C --
  128. andn %g2,%g4,%g2
  129. and u3,v3,%g3
  130. ldx [up-8],u3
  131. fanop
  132. C --
  133. or %g3,%g2,%g2
  134. subcc n,4,n
  135. ldx [vp-8],v3
  136. fanop
  137. C --
  138. srlx %g2,63,cy
  139. add u0,v0,%g1
  140. stx %g4,[rp-8]
  141. fanop
  142. C --
  143. add %g1,cy,%g4
  144. or u0,v0,%g2
  145. bge,pt %icc,.Loop
  146. fanop
  147. C END MAIN LOOP
  148. .Lend4567:
  149. andn %g2,%g4,%g2
  150. and u0,v0,%g3
  151. or %g3,%g2,%g2
  152. srlx %g2,63,cy
  153. add u1,v1,%g1
  154. stx %g4,[rp+0]
  155. add %g1,cy,%g4
  156. or u1,v1,%g2
  157. andn %g2,%g4,%g2
  158. and u1,v1,%g3
  159. or %g3,%g2,%g2
  160. srlx %g2,63,cy
  161. add u2,v2,%g1
  162. stx %g4,[rp+8]
  163. add %g1,cy,%g4
  164. or u2,v2,%g2
  165. andn %g2,%g4,%g2
  166. and u2,v2,%g3
  167. or %g3,%g2,%g2
  168. add rp,32,rp
  169. srlx %g2,63,cy
  170. add u3,v3,%g1
  171. stx %g4,[rp-16]
  172. add %g1,cy,%g4
  173. or u3,v3,%g2
  174. andn %g2,%g4,%g2
  175. and u3,v3,%g3
  176. or %g3,%g2,%g2
  177. srlx %g2,63,cy
  178. stx %g4,[rp-8]
  179. addcc n,4,n
  180. bz,pn %icc,.Lret
  181. fanop
  182. .Loop0: ldx [up],u0
  183. add up,8,up
  184. ldx [vp],v0
  185. add vp,8,vp
  186. add rp,8,rp
  187. subcc n,1,n
  188. add u0,v0,%g1
  189. or u0,v0,%g2
  190. add %g1,cy,%g4
  191. and u0,v0,%g3
  192. andn %g2,%g4,%g2
  193. stx %g4,[rp-8]
  194. or %g3,%g2,%g2
  195. bnz,pt %icc,.Loop0
  196. srlx %g2,63,cy
  197. .Lret: mov cy,%i0
  198. ret
  199. restore
  200. EPILOGUE(mpn_add_n)