1/* AMD64 __mpn_addmul_1 -- Multiply a limb vector with a limb and add 2 the result to a second limb vector. 3 Copyright (C) 2004 Free Software Foundation, Inc. 4 This file is part of the GNU MP Library. 5 6 The GNU MP Library is free software; you can redistribute it and/or modify 7 it under the terms of the GNU Lesser General Public License as published by 8 the Free Software Foundation; either version 2.1 of the License, or (at your 9 option) any later version. 10 11 The GNU MP Library is distributed in the hope that it will be useful, but 12 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 13 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 14 License for more details. 15 16 You should have received a copy of the GNU Lesser General Public License 17 along with the GNU MP Library; see the file COPYING.LIB. If not, write to 18 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, 19 MA 02111-1307, USA. */ 20 21#include "sysdep.h" 22#include "asm-syntax.h" 23 24 .text 25ENTRY (__mpn_addmul_1) 26 movq %rdx, %r11 27 leaq (%rsi,%rdx,8), %rsi 28 leaq (%rdi,%rdx,8), %rdi 29 negq %r11 30 xorl %r8d, %r8d 31 xorl %r10d, %r10d 32 .p2align 2 33L(loop): 34 movq (%rsi,%r11,8), %rax 35 mulq %rcx 36 addq (%rdi,%r11,8), %rax 37 adcq %r10, %rdx 38 addq %r8, %rax 39 movq %r10, %r8 40 movq %rax, (%rdi,%r11,8) 41 adcq %rdx, %r8 42 incq %r11 43 jne L(loop) 44 movq %r8, %rax 45 ret 46END (__mpn_addmul_1) 47