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

数学计算

开发平台:

Unix_Linux

  1. dnl  PowerPC-32 mpn_addlsh1_n -- rp[] = up[] + (vp[] << 1)
  2. dnl  Copyright 2003, 2005, 2007 Free Software Foundation, Inc.
  3. dnl  This file is part of the GNU MP Library.
  4. dnl  The GNU MP Library is free software; you can redistribute it and/or modify
  5. dnl  it under the terms of the GNU Lesser General Public License as published
  6. dnl  by the Free Software Foundation; either version 3 of the License, or (at
  7. dnl  your option) any later version.
  8. dnl  The GNU MP Library is distributed in the hope that it will be useful, but
  9. dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  10. dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  11. dnl  License for more details.
  12. dnl  You should have received a copy of the GNU Lesser General Public License
  13. dnl  along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  14. include(`../config.m4')
  15. C                cycles/limb
  16. C 603e:            ?
  17. C 604e:            4.0
  18. C 75x (G3):        5.0
  19. C 7400,7410 (G4):  5.0
  20. C 744x,745x (G4+): 5.0
  21. C power4/ppc970:   4.25
  22. C power5:          5.0
  23. C INPUT PARAMETERS
  24. C rp r3
  25. C up r4
  26. C vp r5
  27. C n r6
  28. define(`rp',`r3')
  29. define(`up',`r4')
  30. define(`vp',`r5')
  31. define(`s0',`r6')
  32. define(`s1',`r7')
  33. define(`u0',`r8')
  34. define(`v0',`r10')
  35. define(`v1',`r11')
  36. ASM_START()
  37. PROLOGUE(mpn_addlsh1_n)
  38. mtctr r6 C copy n in ctr
  39. addic r31, r31, 0 C clear cy
  40. lwz v0, 0(vp) C load v limb
  41. lwz u0, 0(up) C load u limb
  42. addi up, up, -4 C update up
  43. addi rp, rp, -4 C update rp
  44. slwi s1, v0, 1
  45. bdz L(end) C If done, skip loop
  46. L(loop):
  47. lwz v1, 4(vp) C load v limb
  48. adde s1, s1, u0 C add limbs with cy, set cy
  49. srwi s0, v0, 31 C shift down previous v limb
  50. stw s1, 4(rp) C store result limb
  51. lwzu u0, 8(up) C load u limb and update up
  52. rlwimi s0, v1, 1, 0,30 C left shift v limb and merge with prev v limb
  53. bdz L(exit) C decrement ctr and exit if done
  54. lwzu v0, 8(vp) C load v limb and update vp
  55. adde s0, s0, u0 C add limbs with cy, set cy
  56. srwi s1, v1, 31 C shift down previous v limb
  57. stwu s0, 8(rp) C store result limb and update rp
  58. lwz u0, 4(up) C load u limb
  59. rlwimi s1, v0, 1, 0,30 C left shift v limb and merge with prev v limb
  60. bdnz L(loop) C decrement ctr and loop back
  61. L(end): adde r7, s1, u0
  62. srwi r4, v0, 31
  63. stw r7, 4(rp) C store last result limb
  64. addze r3, r4
  65. blr
  66. L(exit):
  67. adde r7, s0, u0
  68. srwi r4, v1, 31
  69. stw r7, 8(rp) C store last result limb
  70. addze r3, r4
  71. blr
  72. EPILOGUE()