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

数学计算

开发平台:

Unix_Linux

  1. dnl  mc68020 mpn_rshift -- mpn right shift.
  2. dnl  Copyright 1996, 1999, 2000, 2001, 2002, 2003 Free Software Foundation,
  3. dnl  Inc.
  4. dnl
  5. dnl  This file is part of the GNU MP Library.
  6. dnl
  7. dnl  The GNU MP Library is free software; you can redistribute it and/or
  8. dnl  modify it under the terms of the GNU Lesser General Public License as
  9. dnl  published by the Free Software Foundation; either version 3 of the
  10. dnl  License, or (at your option) any later version.
  11. dnl
  12. dnl  The GNU MP Library is distributed in the hope that it will be useful,
  13. dnl  but WITHOUT ANY WARRANTY; without even the implied warranty of
  14. dnl  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15. dnl  Lesser General Public License for more details.
  16. dnl
  17. dnl  You should have received a copy of the GNU Lesser General Public License
  18. dnl  along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  19. include(`../config.m4')
  20. C           cycles/limb
  21. C        shift==1  shift>1
  22. C 68040:    9         12
  23. C mp_limb_t mpn_rshift (mp_ptr res_ptr, mp_srcptr s_ptr, mp_size_t s_size,
  24. C                       unsigned cnt);
  25. C
  26. C The "cnt" parameter is either 16 bits or 32 bits depending on
  27. C SIZEOF_UNSIGNED (see ABI notes in mpn/m68k/README).  The value is of
  28. C course only 1 to 31.  When loaded as 16 bits there's garbage in the upper
  29. C half, hence the use of cmpw.  The shift instructions take the their count
  30. C modulo 64, so the upper part doesn't matter to them either.
  31. C
  32. C INPUT PARAMETERS
  33. C res_ptr (sp + 4)
  34. C s_ptr (sp + 8)
  35. C s_size (sp + 12)
  36. C cnt (sp + 16)
  37. define(res_ptr, `a1')
  38. define(s_ptr,   `a0')
  39. define(s_size,  `d6')
  40. define(cnt,     `d4')
  41. ifdef(`SIZEOF_UNSIGNED',,
  42. `m4_error(`SIZEOF_UNSIGNED not defined, should be in config.m4
  43. ')')
  44. PROLOGUE(mpn_rshift)
  45. C Save used registers on the stack.
  46. moveml d2-d6/a2, M(-,sp)
  47. C Copy the arguments to registers.
  48. movel M(sp,28), res_ptr
  49. movel M(sp,32), s_ptr
  50. movel M(sp,36), s_size
  51. ifelse(SIZEOF_UNSIGNED,2,
  52. ` movew M(sp,40), cnt',
  53. ` movel M(sp,40), cnt')
  54. moveql #1, d5
  55. cmpw d5, cnt
  56. bne L(Lnormal)
  57. cmpl res_ptr, s_ptr
  58. bls L(Lspecial) C jump if res_ptr >= s_ptr
  59. ifelse(scale_available_p,1,`
  60. lea M(res_ptr,s_size,l,4), a2
  61. ',`
  62. movel s_size, d0
  63. asll #2, d0
  64. lea M(res_ptr,d0,l), a2
  65. ')
  66. cmpl s_ptr, a2
  67. bls L(Lspecial) C jump if s_ptr >= res_ptr + s_size
  68. L(Lnormal:)
  69. moveql #32, d5
  70. subl cnt, d5
  71. movel M(s_ptr,+), d2
  72. movel d2, d0
  73. lsll d5, d0 C compute carry limb
  74. lsrl cnt, d2
  75. movel d2, d1
  76. subql #1, s_size
  77. beq L(Lend)
  78. lsrl #1, s_size
  79. bcs L(L1)
  80. subql #1, s_size
  81. L(Loop:)
  82. movel M(s_ptr,+), d2
  83. movel d2, d3
  84. lsll d5, d3
  85. orl d3, d1
  86. movel d1, M(res_ptr,+)
  87. lsrl cnt, d2
  88. L(L1:)
  89. movel M(s_ptr,+), d1
  90. movel d1, d3
  91. lsll d5, d3
  92. orl d3, d2
  93. movel d2, M(res_ptr,+)
  94. lsrl cnt, d1
  95. dbf s_size, L(Loop)
  96. subl #0x10000, s_size
  97. bcc L(Loop)
  98. L(Lend:)
  99. movel d1, M(res_ptr) C store most significant limb
  100. C Restore used registers from stack frame.
  101. moveml M(sp,+), d2-d6/a2
  102. rts
  103. C We loop from most significant end of the arrays, which is only permissable
  104. C if the source and destination don't overlap, since the function is
  105. C documented to work for overlapping source and destination.
  106. L(Lspecial:)
  107. ifelse(scale_available_p,1,`
  108. lea M(s_ptr,s_size,l,4), s_ptr
  109. lea M(res_ptr,s_size,l,4), res_ptr
  110. ',`
  111. movel s_size, d0
  112. asll #2, d0
  113. addl d0, s_ptr
  114. addl d0, res_ptr
  115. ')
  116. clrl d0 C initialize carry
  117. eorw #1, s_size
  118. lsrl #1, s_size
  119. bcc L(LL1)
  120. subql #1, s_size
  121. L(LLoop:)
  122. movel M(-,s_ptr), d2
  123. roxrl #1, d2
  124. movel d2, M(-,res_ptr)
  125. L(LL1:)
  126. movel M(-,s_ptr), d2
  127. roxrl #1, d2
  128. movel d2, M(-,res_ptr)
  129. dbf s_size, L(LLoop)
  130. roxrl #1, d0 C save cy in msb
  131. subl #0x10000, s_size
  132. bcs L(LLend)
  133. addl d0, d0 C restore cy
  134. bra L(LLoop)
  135. L(LLend:)
  136. C Restore used registers from stack frame.
  137. moveml M(sp,+), d2-d6/a2
  138. rts
  139. EPILOGUE(mpn_rshift)