1dnl Alpha mpn_mul_1 -- Multiply a limb vector with a limb and store 2dnl the result in a second limb vector. 3 4dnl Copyright 1992, 1994, 1995, 2000, 2002 Free Software Foundation, Inc. 5 6dnl This file is part of the GNU MP Library. 7 8dnl The GNU MP Library is free software; you can redistribute it and/or modify 9dnl it under the terms of the GNU Lesser General Public License as published 10dnl by the Free Software Foundation; either version 3 of the License, or (at 11dnl your option) any later version. 12 13dnl The GNU MP Library is distributed in the hope that it will be useful, but 14dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 15dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 16dnl License for more details. 17 18dnl You should have received a copy of the GNU Lesser General Public License 19dnl along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. 20 21include(`../config.m4') 22 23C cycles/limb 24C EV4: 42 25C EV5: 18 26C EV6: 7 27 28C INPUT PARAMETERS 29C rp r16 30C up r17 31C n r18 32C vl r19 33C cl r20 34 35 36ASM_START() 37PROLOGUE(mpn_mul_1c) 38 ldq r2,0(r17) C r2 = s1_limb 39 lda r18,-1(r18) C size-- 40 mulq r2,r19,r3 C r3 = prod_low 41 umulh r2,r19,r4 C r4 = prod_high 42 beq r18,$Le1c C jump if size was == 1 43 ldq r2,8(r17) C r2 = s1_limb 44 lda r18,-1(r18) C size-- 45 addq r3,r20,r3 C r3 = cy_limb + cl 46 stq r3,0(r16) 47 cmpult r3,r20,r0 C r0 = carry from (cy_limb + cl) 48 bne r18,$Loop C jump if size was == 2 49 br r31,$Le2 50$Le1c: addq r3,r20,r3 C r3 = cy_limb + cl 51 cmpult r3,r20,r0 C r0 = carry from (cy_limb + cl) 52$Le1: stq r3,0(r16) 53 addq r4,r0,r0 54 ret r31,(r26),1 55EPILOGUE(mpn_mul_1c) 56 57PROLOGUE(mpn_mul_1) 58 ldq r2,0(r17) C r2 = s1_limb 59 lda r18,-1(r18) C size-- 60 mulq r2,r19,r3 C r3 = prod_low 61 bic r31,r31,r0 C clear cy_limb 62 umulh r2,r19,r4 C r4 = prod_high 63 beq r18,$Le1 C jump if size was == 1 64 ldq r2,8(r17) C r2 = s1_limb 65 lda r18,-1(r18) C size-- 66 stq r3,0(r16) 67 beq r18,$Le2 C jump if size was == 2 68 69 ALIGN(8) 70$Loop: mulq r2,r19,r3 C r3 = prod_low 71 addq r4,r0,r0 C cy_limb = cy_limb + 'cy' 72 lda r18,-1(r18) C size-- 73 umulh r2,r19,r4 C r4 = prod_high 74 ldq r2,16(r17) C r2 = s1_limb 75 lda r17,8(r17) C s1_ptr++ 76 addq r3,r0,r3 C r3 = cy_limb + prod_low 77 stq r3,8(r16) 78 cmpult r3,r0,r0 C r0 = carry from (cy_limb + prod_low) 79 lda r16,8(r16) C res_ptr++ 80 bne r18,$Loop 81 82$Le2: mulq r2,r19,r3 C r3 = prod_low 83 addq r4,r0,r0 C cy_limb = cy_limb + 'cy' 84 umulh r2,r19,r4 C r4 = prod_high 85 addq r3,r0,r3 C r3 = cy_limb + prod_low 86 cmpult r3,r0,r0 C r0 = carry from (cy_limb + prod_low) 87 stq r3,8(r16) 88 addq r4,r0,r0 C cy_limb = prod_high + cy 89 ret r31,(r26),1 90EPILOGUE(mpn_mul_1) 91ASM_END() 92