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

数学计算

开发平台:

Unix_Linux

  1. # VAX __gmpn_add_n -- Add two limb vectors of the same length > 0 and store
  2. # sum in a third limb vector.
  3. # Copyright 1999, 2000 Free Software Foundation, Inc.
  4. # This file is part of the GNU MP Library.
  5. # The GNU MP Library is free software; you can redistribute it and/or modify
  6. # it under the terms of the GNU Lesser General Public License as published by
  7. # the Free Software Foundation; either version 3 of the License, or (at your
  8. # option) any later version.
  9. # The GNU MP Library is distributed in the hope that it will be useful, but
  10. # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  11. # or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  12. # License for more details.
  13. # You should have received a copy of the GNU Lesser General Public License
  14. # along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  15. # INPUT PARAMETERS
  16. # res_ptr (sp + 4)
  17. # s1_ptr (sp + 8)
  18. # s2_ptr (sp + 12)
  19. # size (sp + 16)
  20. .text
  21. .align 1
  22. .globl ___gmpn_add_n
  23. ___gmpn_add_n:
  24. .word 0x0
  25. movl 16(ap),r0
  26. movl 12(ap),r1
  27. movl 8(ap),r2
  28. movl 4(ap),r3
  29. mnegl r0,r5
  30. addl2 $3,r0
  31. ashl $-2,r0,r0 # unroll loop count
  32. bicl2 $-4,r5 # mask out low 2 bits
  33. movaq (r5)[r5],r5 # 9x
  34. jmp Loop(r5)
  35. Loop: movl (r2)+,r4
  36. adwc (r1)+,r4
  37. movl r4,(r3)+
  38. movl (r2)+,r4
  39. adwc (r1)+,r4
  40. movl r4,(r3)+
  41. movl (r2)+,r4
  42. adwc (r1)+,r4
  43. movl r4,(r3)+
  44. movl (r2)+,r4
  45. adwc (r1)+,r4
  46. movl r4,(r3)+
  47. sobgtr r0,Loop
  48. adwc r0,r0
  49. ret