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

数学计算

开发平台:

Unix_Linux

  1. dnl  PowerPC-64 mpn_lshift -- rp[] = up[] << cnt
  2. dnl  Copyright 2003, 2005 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 POWER3/PPC630:     1.5
  17. C POWER4/PPC970:     3.0
  18. C INPUT PARAMETERS
  19. define(`rp',`r3')
  20. define(`up',`r4')
  21. define(`n',`r5')
  22. define(`cnt',`r6')
  23. define(`tnc',`r5')
  24. define(`v0',`r0')
  25. define(`v1',`r7')
  26. define(`u0',`r8')
  27. define(`u1',`r9')
  28. define(`h0',`r10')
  29. define(`h1',`r11')
  30. ASM_START()
  31. PROLOGUE(mpn_lshift)
  32. ifdef(`HAVE_ABI_mode32',
  33. ` rldicl r7, r5, 0, 32 C zero extend n
  34. mtctr r7', C copy n to count register
  35. ` mtctr n') C copy n to count register
  36. ifdef(`HAVE_ABI_mode32',
  37. ` rldic r0, n, 3, 32', C byte count corresponding to n
  38. ` rldicr r0, n, 3, 60') C byte count corresponding to n
  39. add rp, rp, r0 C rp = rp + n
  40. add up, up, r0 C up = up + n
  41. addi rp, rp, 8 C rp now points 16 beyond end
  42. addi up, up, -8 C up now points to last limb
  43. subfic tnc, cnt, 64 C reverse shift count
  44. ld u0, 0(up)
  45. sld h0, u0, cnt
  46. srd r12, u0, tnc C return value
  47. bdz L(1) C jump for n = 1
  48. ld u1, -8(up)
  49. bdz L(2) C jump for n = 2
  50. ldu u0, -16(up)
  51. bdz L(end) C jump for n = 3
  52. L(oop): srd v1, u1, tnc
  53. sld h1, u1, cnt
  54. ld u1, -8(up)
  55. or h0, v1, h0
  56. stdu h0, -16(rp)
  57. bdz L(exit)
  58. srd v0, u0, tnc
  59. sld h0, u0, cnt
  60. ldu u0, -16(up)
  61. or h1, v0, h1
  62. std h1, -8(rp)
  63. bdnz L(oop)
  64. L(end): srd v1, u1, tnc
  65. sld h1, u1, cnt
  66. or h0, v1, h0
  67. stdu h0, -16(rp)
  68. srd v0, u0, tnc
  69. sld h0, u0, cnt
  70. or h1, v0, h1
  71. std h1, -8(rp)
  72. L(1): std h0, -16(rp)
  73. ifdef(`HAVE_ABI_mode32',
  74. ` srdi r3, r12, 32
  75. mr r4, r12
  76. ',` mr r3, r12
  77. ')
  78. blr
  79. L(exit): srd v0, u0, tnc
  80. sld h0, u0, cnt
  81. or h1, v0, h1
  82. std h1, -8(rp)
  83. L(2): srd v1, u1, tnc
  84. sld h1, u1, cnt
  85. or h0, v1, h0
  86. stdu h0, -16(rp)
  87. std h1, -8(rp)
  88. ifdef(`HAVE_ABI_mode32',
  89. ` srdi r3, r12, 32
  90. mr r4, r12
  91. ',` mr r3, r12
  92. ')
  93. blr
  94. EPILOGUE()