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

数学计算

开发平台:

Unix_Linux

  1. dnl  x86-64 mpn_rshift optimized for "Core 2".
  2. dnl  Copyright 2007, 2009 Free Software Foundation, Inc.
  3. dnl
  4. dnl  This file is part of the GNU MP Library.
  5. dnl
  6. dnl  The GNU MP Library is free software; you can redistribute it and/or
  7. dnl  modify it under the terms of the GNU Lesser General Public License as
  8. dnl  published by the Free Software Foundation; either version 3 of the
  9. dnl  License, or (at your option) any later version.
  10. dnl
  11. dnl  The GNU MP Library is distributed in the hope that it will be useful,
  12. dnl  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13. dnl  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14. dnl  Lesser General Public License for more details.
  15. dnl
  16. dnl  You should have received a copy of the GNU Lesser General Public License
  17. dnl  along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  18. include(`../config.m4')
  19. C      cycles/limb
  20. C K8,K9:  4.25
  21. C K10:  4.25
  22. C P4: 14.7
  23. C P6 core2:  1.27
  24. C P6 corei7:  1.5
  25. C INPUT PARAMETERS
  26. define(`rp', `%rdi')
  27. define(`up', `%rsi')
  28. define(`n', `%rdx')
  29. define(`cnt', `%cl')
  30. ASM_START()
  31. TEXT
  32. ALIGN(16)
  33. PROLOGUE(mpn_rshift)
  34. mov %edx, %eax
  35. and $3, %eax
  36. jne L(nb00)
  37. L(b00): C n = 4, 8, 12, ...
  38. mov (up), %r10
  39. mov 8(up), %r11
  40. xor %eax, %eax
  41. shrd %cl, %r10, %rax
  42. mov 16(up), %r8
  43. lea 8(up), up
  44. lea -24(rp), rp
  45. sub $4, n
  46. jmp L(00)
  47. L(nb00):C n = 1, 5, 9, ...
  48. cmp $2, %eax
  49. jae L(nb01)
  50. L(b01): mov (up), %r9
  51. xor %eax, %eax
  52. shrd %cl, %r9, %rax
  53. sub $2, n
  54. jb L(le1)
  55. mov 8(up), %r10
  56. mov 16(up), %r11
  57. lea 16(up), up
  58. lea -16(rp), rp
  59. jmp L(01)
  60. L(le1): shr %cl, %r9
  61. mov %r9, (rp)
  62. ret
  63. L(nb01):C n = 2, 6, 10, ...
  64. jne L(b11)
  65. L(b10): mov (up), %r8
  66. mov 8(up), %r9
  67. xor %eax, %eax
  68. shrd %cl, %r8, %rax
  69. sub $3, n
  70. jb L(le2)
  71. mov 16(up), %r10
  72. lea 24(up), up
  73. lea -8(rp), rp
  74. jmp L(10)
  75. L(le2): shrd %cl, %r9, %r8
  76. mov %r8, (rp)
  77. shr %cl, %r9
  78. mov %r9, 8(rp)
  79. ret
  80. ALIGN(16)
  81. L(b11): C n = 3, 7, 11, ...
  82. mov (up), %r11
  83. mov 8(up), %r8
  84. xor %eax, %eax
  85. shrd %cl, %r11, %rax
  86. mov 16(up), %r9
  87. lea 32(up), up
  88. sub $4, n
  89. jb L(end)
  90. ALIGN(16)
  91. L(top): shrd %cl, %r8, %r11
  92. mov -8(up), %r10
  93. mov %r11, (rp)
  94. L(10): shrd %cl, %r9, %r8
  95. mov (up), %r11
  96. mov %r8, 8(rp)
  97. L(01): shrd %cl, %r10, %r9
  98. mov 8(up), %r8
  99. mov %r9, 16(rp)
  100. L(00): shrd %cl, %r11, %r10
  101. mov 16(up), %r9
  102. mov %r10, 24(rp)
  103. add $32, up
  104. lea 32(rp), rp
  105. sub $4, n
  106. jnc L(top)
  107. L(end): shrd %cl, %r8, %r11
  108. mov %r11, (rp)
  109. shrd %cl, %r9, %r8
  110. mov %r8, 8(rp)
  111. shr %cl, %r9
  112. mov %r9, 16(rp)
  113. ret
  114. EPILOGUE()