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

数学计算

开发平台:

Unix_Linux

  1. dnl  PowerPC-32 mpn_lshift -- Shift a number left.
  2. dnl  Copyright 1995, 1998, 2000, 2002, 2003, 2004, 2005 Free Software
  3. dnl  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                cycles/limb
  17. C 603e:            ?
  18. C 604e:            3.0
  19. C 75x (G3):        3.0
  20. C 7400,7410 (G4):  3.0
  21. C 7445,7455 (G4+): 2.5
  22. C 7447,7457 (G4+): 2.25
  23. C power4/ppc970:   2.5
  24. C power5:          2.5
  25. C INPUT PARAMETERS
  26. C rp r3
  27. C up r4
  28. C n r5
  29. C cnt r6
  30. ASM_START()
  31. PROLOGUE(mpn_lshift)
  32. cmpwi cr0, r5, 12 C more than 12 limbs?
  33. slwi r0, r5, 2
  34. add r4, r4, r0 C make r4 point at end of s1
  35. add r7, r3, r0 C make r7 point at end of res
  36. bgt L(BIG) C branch if more than 12 limbs
  37. mtctr r5 C copy size into CTR
  38. subfic r8, r6, 32
  39. lwzu r11, -4(r4) C load first s1 limb
  40. srw r3, r11, r8 C compute function return value
  41. bdz L(end1)
  42. L(oop): lwzu r10, -4(r4)
  43. slw r9, r11, r6
  44. srw r12, r10, r8
  45. or r9, r9, r12
  46. stwu r9, -4(r7)
  47. bdz L(end2)
  48. lwzu r11, -4(r4)
  49. slw r9, r10, r6
  50. srw r12, r11, r8
  51. or r9, r9, r12
  52. stwu r9, -4(r7)
  53. bdnz L(oop)
  54. L(end1):
  55. slw r0, r11, r6
  56. stw r0, -4(r7)
  57. blr
  58. L(end2):
  59. slw r0, r10, r6
  60. stw r0, -4(r7)
  61. blr
  62. L(BIG):
  63. stmw r24, -32(r1) C save registers we are supposed to preserve
  64. lwzu r9, -4(r4)
  65. subfic r8, r6, 32
  66. srw r3, r9, r8 C compute function return value
  67. slw r0, r9, r6
  68. addi r5, r5, -1
  69. andi. r10, r5, 3 C count for spill loop
  70. beq L(e)
  71. mtctr r10
  72. lwzu r28, -4(r4)
  73. bdz L(xe0)
  74. L(loop0):
  75. slw r12, r28, r6
  76. srw r24, r28, r8
  77. lwzu r28, -4(r4)
  78. or r24, r0, r24
  79. stwu r24, -4(r7)
  80. mr r0, r12
  81. bdnz L(loop0) C taken at most once!
  82. L(xe0): slw r12, r28, r6
  83. srw r24, r28, r8
  84. or r24, r0, r24
  85. stwu r24, -4(r7)
  86. mr r0, r12
  87. L(e): srwi r5, r5, 2 C count for unrolled loop
  88. addi r5, r5, -1
  89. mtctr r5
  90. lwz r28, -4(r4)
  91. lwz r29, -8(r4)
  92. lwz r30, -12(r4)
  93. lwzu r31, -16(r4)
  94. L(loopU):
  95. slw r9, r28, r6
  96. srw r24, r28, r8
  97. lwz r28, -4(r4)
  98. slw r10, r29, r6
  99. srw r25, r29, r8
  100. lwz r29, -8(r4)
  101. slw r11, r30, r6
  102. srw r26, r30, r8
  103. lwz r30, -12(r4)
  104. slw r12, r31, r6
  105. srw r27, r31, r8
  106. lwzu r31, -16(r4)
  107. or r24, r0, r24
  108. stw r24, -4(r7)
  109. or r25, r9, r25
  110. stw r25, -8(r7)
  111. or r26, r10, r26
  112. stw r26, -12(r7)
  113. or r27, r11, r27
  114. stwu r27, -16(r7)
  115. mr r0, r12
  116. bdnz L(loopU)
  117. slw r9, r28, r6
  118. srw r24, r28, r8
  119. slw r10, r29, r6
  120. srw r25, r29, r8
  121. slw r11, r30, r6
  122. srw r26, r30, r8
  123. slw r12, r31, r6
  124. srw r27, r31, r8
  125. or r24, r0, r24
  126. stw r24, -4(r7)
  127. or r25, r9, r25
  128. stw r25, -8(r7)
  129. or r26, r10, r26
  130. stw r26, -12(r7)
  131. or r27, r11, r27
  132. stw r27, -16(r7)
  133. stw r12, -20(r7)
  134. lmw r24, -32(r1) C restore registers
  135. blr
  136. EPILOGUE(mpn_lshift)