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

数学计算

开发平台:

Unix_Linux

  1. dnl  IBM POWER mpn_mul_1 -- Multiply a limb vector with a limb and store the
  2. dnl  result in a second limb vector.
  3. dnl  Copyright 1992, 1994, 1999, 2000, 2001 Free Software 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. dnl  INPUT PARAMETERS
  16. dnl  res_ptr r3
  17. dnl  s1_ptr r4
  18. dnl  size r5
  19. dnl  s2_limb r6
  20. dnl  The POWER architecture has no unsigned 32x32->64 bit multiplication
  21. dnl  instruction.  To obtain that operation, we have to use the 32x32->64
  22. dnl  signed multiplication instruction, and add the appropriate compensation to
  23. dnl  the high limb of the result.  We add the multiplicand if the multiplier
  24. dnl  has its most significant bit set, and we add the multiplier if the
  25. dnl  multiplicand has its most significant bit set.  We need to preserve the
  26. dnl  carry flag between each iteration, so we have to compute the compensation
  27. dnl  carefully (the natural, srai+and doesn't work).  Since all POWER can
  28. dnl  branch in zero cycles, we use conditional branches for the compensation.
  29. include(`../config.m4')
  30. ASM_START()
  31. PROLOGUE(mpn_mul_1)
  32. cal 3,-4(3)
  33. l 0,0(4)
  34. cmpi 0,6,0
  35. mtctr 5
  36. mul 9,0,6
  37. srai 7,0,31
  38. and 7,7,6
  39. mfmq 8
  40. ai 0,0,0 C reset carry
  41. cax 9,9,7
  42. blt Lneg
  43. Lpos: bdz Lend
  44. Lploop: lu 0,4(4)
  45. stu 8,4(3)
  46. cmpi 0,0,0
  47. mul 10,0,6
  48. mfmq 0
  49. ae 8,0,9
  50. bge Lp0
  51. cax 10,10,6 C adjust high limb for negative limb from s1
  52. Lp0: bdz Lend0
  53. lu 0,4(4)
  54. stu 8,4(3)
  55. cmpi 0,0,0
  56. mul 9,0,6
  57. mfmq 0
  58. ae 8,0,10
  59. bge Lp1
  60. cax 9,9,6 C adjust high limb for negative limb from s1
  61. Lp1: bdn Lploop
  62. b Lend
  63. Lneg: cax 9,9,0
  64. bdz Lend
  65. Lnloop: lu 0,4(4)
  66. stu 8,4(3)
  67. cmpi 0,0,0
  68. mul 10,0,6
  69. cax 10,10,0 C adjust high limb for negative s2_limb
  70. mfmq 0
  71. ae 8,0,9
  72. bge Ln0
  73. cax 10,10,6 C adjust high limb for negative limb from s1
  74. Ln0: bdz Lend0
  75. lu 0,4(4)
  76. stu 8,4(3)
  77. cmpi 0,0,0
  78. mul 9,0,6
  79. cax 9,9,0 C adjust high limb for negative s2_limb
  80. mfmq 0
  81. ae 8,0,10
  82. bge Ln1
  83. cax 9,9,6 C adjust high limb for negative limb from s1
  84. Ln1: bdn Lnloop
  85. b Lend
  86. Lend0: cal 9,0(10)
  87. Lend: st 8,4(3)
  88. aze 3,9
  89. br
  90. EPILOGUE(mpn_mul_1)