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

数学计算

开发平台:

Unix_Linux

  1. dnl  HP-PA 2.0 32-bit mpn_add_n -- Add two limb vectors of the same length > 0
  2. dnl  and store sum in a third limb vector.
  3. dnl  Copyright 1997, 1998, 2000, 2001, 2002 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 INPUT PARAMETERS
  17. C res_ptr gr26
  18. C s1_ptr gr25
  19. C s2_ptr gr24
  20. C size gr23
  21. C This runs at 2 cycles/limb on PA8000.
  22. ASM_START()
  23. PROLOGUE(mpn_add_n)
  24. sub %r0,%r23,%r22
  25. zdep %r22,30,3,%r28 C r28 = 2 * (-n & 7)
  26. zdep %r22,29,3,%r22 C r22 = 4 * (-n & 7)
  27. sub %r25,%r22,%r25 C offset s1_ptr
  28. sub %r24,%r22,%r24 C offset s2_ptr
  29. sub %r26,%r22,%r26 C offset res_ptr
  30. blr %r28,%r0 C branch into loop
  31. add %r0,%r0,%r0 C reset carry
  32. LDEF(loop)
  33. ldw 0(%r25),%r20
  34. ldw 0(%r24),%r31
  35. addc %r20,%r31,%r20
  36. stw %r20,0(%r26)
  37. LDEF(7)
  38. ldw 4(%r25),%r21
  39. ldw 4(%r24),%r19
  40. addc %r21,%r19,%r21
  41. stw %r21,4(%r26)
  42. LDEF(6)
  43. ldw 8(%r25),%r20
  44. ldw 8(%r24),%r31
  45. addc %r20,%r31,%r20
  46. stw %r20,8(%r26)
  47. LDEF(5)
  48. ldw 12(%r25),%r21
  49. ldw 12(%r24),%r19
  50. addc %r21,%r19,%r21
  51. stw %r21,12(%r26)
  52. LDEF(4)
  53. ldw 16(%r25),%r20
  54. ldw 16(%r24),%r31
  55. addc %r20,%r31,%r20
  56. stw %r20,16(%r26)
  57. LDEF(3)
  58. ldw 20(%r25),%r21
  59. ldw 20(%r24),%r19
  60. addc %r21,%r19,%r21
  61. stw %r21,20(%r26)
  62. LDEF(2)
  63. ldw 24(%r25),%r20
  64. ldw 24(%r24),%r31
  65. addc %r20,%r31,%r20
  66. stw %r20,24(%r26)
  67. LDEF(1)
  68. ldw 28(%r25),%r21
  69. ldo 32(%r25),%r25
  70. ldw 28(%r24),%r19
  71. addc %r21,%r19,%r21
  72. stw %r21,28(%r26)
  73. ldo 32(%r24),%r24
  74. addib,> -8,%r23,L(loop)
  75. ldo 32(%r26),%r26
  76. bv (%r2)
  77. addc %r0,%r0,%r28
  78. EPILOGUE()