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

数学计算

开发平台:

Unix_Linux

  1. # VAX __gmpn_mul_1 -- Multiply a limb vector with a limb and store
  2. # the result in a second limb vector.
  3. # Copyright 1992, 1994, 1996, 2000 Free Software Foundation, Inc.
  4. # This file is part of the GNU MP Library.
  5. # The GNU MP Library is free software; you can redistribute it and/or modify
  6. # it under the terms of the GNU Lesser General Public License as published by
  7. # the Free Software Foundation; either version 3 of the License, or (at your
  8. # option) any later version.
  9. # The GNU MP Library is distributed in the hope that it will be useful, but
  10. # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  11. # or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  12. # License for more details.
  13. # You should have received a copy of the GNU Lesser General Public License
  14. # along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.
  15. # INPUT PARAMETERS
  16. # res_ptr (sp + 4)
  17. # s1_ptr (sp + 8)
  18. # size (sp + 12)
  19. # s2_limb (sp + 16)
  20. .text
  21. .align 1
  22. .globl ___gmpn_mul_1
  23. ___gmpn_mul_1:
  24. .word 0xfc0
  25. movl 12(ap),r4
  26. movl 8(ap),r8
  27. movl 4(ap),r9
  28. movl 16(ap),r6
  29. jlss s2_big
  30. # One might want to combine the addl2 and the store below, but that
  31. # is actually just slower according to my timing tests.  (VAX 3600)
  32. clrl r3
  33. incl r4
  34. ashl $-1,r4,r7
  35. jlbc r4,L1
  36. clrl r11
  37. # Loop for S2_LIMB < 0x80000000
  38. Loop1: movl (r8)+,r1
  39. jlss L1n0
  40. emul r1,r6,$0,r2
  41. addl2 r11,r2
  42. adwc $0,r3
  43. movl r2,(r9)+
  44. L1: movl (r8)+,r1
  45. jlss L1n1
  46. L1p1: emul r1,r6,$0,r10
  47. addl2 r3,r10
  48. adwc $0,r11
  49. movl r10,(r9)+
  50. sobgtr r7,Loop1
  51. movl r11,r0
  52. ret
  53. L1n0: emul r1,r6,$0,r2
  54. addl2 r11,r2
  55. adwc r6,r3
  56. movl r2,(r9)+
  57. movl (r8)+,r1
  58. jgeq L1p1
  59. L1n1: emul r1,r6,$0,r10
  60. addl2 r3,r10
  61. adwc r6,r11
  62. movl r10,(r9)+
  63. sobgtr r7,Loop1
  64. movl r11,r0
  65. ret
  66. s2_big: clrl r3
  67. incl r4
  68. ashl $-1,r4,r7
  69. jlbc r4,L2
  70. clrl r11
  71. # Loop for S2_LIMB >= 0x80000000
  72. Loop2: movl (r8)+,r1
  73. jlss L2n0
  74. emul r1,r6,$0,r2
  75. addl2 r11,r2
  76. adwc r1,r3
  77. movl r2,(r9)+
  78. L2: movl (r8)+,r1
  79. jlss L2n1
  80. L2p1: emul r1,r6,$0,r10
  81. addl2 r3,r10
  82. adwc r1,r11
  83. movl r10,(r9)+
  84. sobgtr r7,Loop2
  85. movl r11,r0
  86. ret
  87. L2n0: emul r1,r6,$0,r2
  88. addl2 r1,r3
  89. addl2 r11,r2
  90. adwc r6,r3
  91. movl r2,(r9)+
  92. movl (r8)+,r1
  93. jgeq L2p1
  94. L2n1: emul r1,r6,$0,r10
  95. addl2 r1,r11
  96. addl2 r3,r10
  97. adwc r6,r11
  98. movl r10,(r9)+
  99. sobgtr r7,Loop2
  100. movl r11,r0
  101. ret