开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 查看源码
_sparc_mul.s
资源名称：leda.tar.gz　[点击查看]
上传用户：gzelex
上传日期：2007-01-07
资源大小：707k
文件大小：5k
源码类别：

数值算法/人工智能

开发平台：

MultiPlatform

_sparc_mul.s：源码内容
  2. !------------------------------------------------------------------------------
  3. ! assembler code for unsigned 32 bit multiplication
  4. ! taken from 
  5. !
  6. !            "The SPARC Architecture Manual"  (Appendix E.2)
  7. !
  8. !------------------------------------------------------------------------------
  11. !------------------------------------------------------------------------------
  12. ! extern "C" word Multiply_Words(word a, word b, word* high);
  13. ! (returns lower 32 bit of product, stores higher 32 bit in *high)
  14. !------------------------------------------------------------------------------
  16. .seg "text"
  18. .global _Multiply_Words
  19. _Multiply_Words:
  20. or      %o0, %o1, %o4   !logical or of multiplier and miltiplicand
  21. mov     %o0, %y         !multiplier to Y register 
  22. andncc  %o4, 0xfff, %o5 !mask out lower 12 bits
  23. be      mul_shortway    !can do it the short way
  24. andcc   %g0, %g0, %o4   !zero the partial product; clear N and V cond.
  26. ! long  multiply
  27. mulscc  %o4, %o1, %o4   !first of 33 iterations
  28. mulscc  %o4, %o1, %o4
  29. mulscc  %o4, %o1, %o4
  30. mulscc  %o4, %o1, %o4
  31. mulscc  %o4, %o1, %o4
  32. mulscc  %o4, %o1, %o4
  33. mulscc  %o4, %o1, %o4
  34. mulscc  %o4, %o1, %o4
  35. mulscc  %o4, %o1, %o4
  36. mulscc  %o4, %o1, %o4
  37. mulscc  %o4, %o1, %o4
  38. mulscc  %o4, %o1, %o4
  39. mulscc  %o4, %o1, %o4
  40. mulscc  %o4, %o1, %o4
  41. mulscc  %o4, %o1, %o4
  42. mulscc  %o4, %o1, %o4
  43. mulscc  %o4, %o1, %o4
  44. mulscc  %o4, %o1, %o4
  45. mulscc  %o4, %o1, %o4
  46. mulscc  %o4, %o1, %o4
  47. mulscc  %o4, %o1, %o4
  48. mulscc  %o4, %o1, %o4
  49. mulscc  %o4, %o1, %o4
  50. mulscc  %o4, %o1, %o4
  51. mulscc  %o4, %o1, %o4
  52. mulscc  %o4, %o1, %o4
  53. mulscc  %o4, %o1, %o4
  54. mulscc  %o4, %o1, %o4
  55. mulscc  %o4, %o1, %o4
  56. mulscc  %o4, %o1, %o4
  57. mulscc  %o4, %o1, %o4
  58. mulscc  %o4, %o1, %o4  ! 32nd iteration;
  59. mulscc  %o4, %g0, %o4  ! last iteration only shits
  60. tst     %o1
  61. bge     1f
  62. nop
  63. add     %o4, %o0, %o4
  64. 1:
  65. rd      %y, %o0 
  66. addcc   %o4, %g0, %o1  ! set zero bit appropriately
  67. st      %o1, [%o2]     ! return high bits
  68. retl                   ! leaf-routine return
  69. nop
  71. ! short multiply
  72. mul_shortway:
  73. mulscc  %o4, %o1, %o4 ! first iteration of 13
  74. mulscc  %o4, %o1, %o4
  75. mulscc  %o4, %o1, %o4
  76. mulscc  %o4, %o1, %o4
  77. mulscc  %o4, %o1, %o4
  78. mulscc  %o4, %o1, %o4
  79. mulscc  %o4, %o1, %o4
  80. mulscc  %o4, %o1, %o4
  81. mulscc  %o4, %o1, %o4
  82. mulscc  %o4, %o1, %o4
  83. mulscc  %o4, %o1, %o4
  84. mulscc  %o4, %o1, %o4  ! 12nd iteration
  85. mulscc  %o4, %g0, %o4  ! last iteration only shits
  86. rd      %y, %o5
  87. sll     %o4, 12, %o4   ! left shift partial product by 12 bits
  88. srl     %o5, 20, %o5   ! right shift product by 12 bits
  89. or      %o5, %o4, %o0  ! merge for true product
  90. addcc   %g0, %g0, %o1
  91. st      %o1, [%o2]     ! return high bits
  92. retl                   ! leaf routine return
  93. nop
  96. !------------------------------------------------------------------------------
  97. ! extern "C" word Mult_Inner_Loop(word *p, word *a, word* a_stop, word B)
  98. !------------------------------------------------------------------------------
  99. ! %o0 p
  100. ! %o1 a
  101. ! %o2 a_stop
  102. ! %o3 B
  103. !
  104. .global    _Mult_Inner_Loop
  105. _Mult_Inner_Loop:
  106. cmp        %o1,%o2       ! if a >= a_stop return 
  107. bge        IMULT_L1
  108. mov        0, %g5        ! carry = 0
  109. IMULT_L0:                 ! Loop: while (a < a_stop)
  110. mov        %o3, %y       ! load B to Y register
  111. andcc      %g0, %g0, %g3
  112. ld         [%o1], %g2    ! A = *a
  113. ld         [%o0], %g4    ! P = *p
  114. !
  115. ! compute 64 bit product A*B = (high,low)
  116. !
  117. mulscc     %g3, %g2, %g3 
  118. mulscc     %g3, %g2, %g3
  119. mulscc     %g3, %g2, %g3
  120. mulscc     %g3, %g2, %g3
  121. mulscc     %g3, %g2, %g3
  122. mulscc     %g3, %g2, %g3
  123. mulscc     %g3, %g2, %g3
  124. mulscc     %g3, %g2, %g3
  125. mulscc     %g3, %g2, %g3
  126. mulscc     %g3, %g2, %g3
  127. mulscc     %g3, %g2, %g3
  128. mulscc     %g3, %g2, %g3
  129. mulscc     %g3, %g2, %g3
  130. mulscc     %g3, %g2, %g3
  131. mulscc     %g3, %g2, %g3
  132. mulscc     %g3, %g2, %g3
  133. mulscc     %g3, %g2, %g3
  134. mulscc     %g3, %g2, %g3
  135. mulscc     %g3, %g2, %g3
  136. mulscc     %g3, %g2, %g3
  137. mulscc     %g3, %g2, %g3
  138. mulscc     %g3, %g2, %g3
  139. mulscc     %g3, %g2, %g3
  140. mulscc     %g3, %g2, %g3
  141. mulscc     %g3, %g2, %g3
  142. mulscc     %g3, %g2, %g3
  143. mulscc     %g3, %g2, %g3
  144. mulscc     %g3, %g2, %g3
  145. mulscc     %g3, %g2, %g3
  146. mulscc     %g3, %g2, %g3
  147. mulscc     %g3, %g2, %g3
  148. mulscc     %g3, %g2, %g3
  149. mulscc     %g3, %g0, %g3
  150. tst        %g2            ! sign correction
  151. bge        IMULT_L2
  152. rd         %y, %g2        ! low
  153. add        %g3, %o3, %g3  ! sign correction
  154. IMULT_L2:
  155. addcc      %g5, %g4, %g4  ! P += carry
  156. addx       %g0, %g0, %g5  ! new carry from addition
  157. addcc      %g2, %g4, %g2  ! low += P
  158. addx       %g3, %g5, %g5  ! carry += high
  159. st         %g2, [%o0]     ! *p = low
  160. inc        4,   %o1       ! a++
  161. cmp        %o1, %o2       ! if (a < a_stop) goto IMULT_L0
  162. bl         IMULT_L0       ! loop
  163. inc        4,   %o0       ! p++
  164. IMULT_L1:
  165. st         %g5, [%o0]     ! *p = carry
  166. retl
  167. mov        %g5, %o0       ! return carry