aors_n.asm revision 1.1.1.3
1dnl Intel mpn_add_n/mpn_sub_n optimised for Conroe, Nehalem. 2 3dnl Copyright 2006, 2007, 2011-2013 Free Software Foundation, Inc. 4 5dnl This file is part of the GNU MP Library. 6dnl 7dnl The GNU MP Library is free software; you can redistribute it and/or modify 8dnl it under the terms of either: 9dnl 10dnl * the GNU Lesser General Public License as published by the Free 11dnl Software Foundation; either version 3 of the License, or (at your 12dnl option) any later version. 13dnl 14dnl or 15dnl 16dnl * the GNU General Public License as published by the Free Software 17dnl Foundation; either version 2 of the License, or (at your option) any 18dnl later version. 19dnl 20dnl or both in parallel, as here. 21dnl 22dnl The GNU MP Library is distributed in the hope that it will be useful, but 23dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25dnl for more details. 26dnl 27dnl You should have received copies of the GNU General Public License and the 28dnl GNU Lesser General Public License along with the GNU MP Library. If not, 29dnl see https://www.gnu.org/licenses/. 30 31include(`../config.m4') 32 33 34C cycles/limb 35C AMD K8,K9 2 36C AMD K10 2 37C Intel P4 10 38C Intel core2 2 39C Intel NHM 2 40C Intel SBR 2 41C Intel atom 9 42C VIA nano 3 43 44C INPUT PARAMETERS 45define(`rp', `%rdi') 46define(`up', `%rsi') 47define(`vp', `%rdx') 48define(`n', `%rcx') 49define(`cy', `%r8') 50 51ifdef(`OPERATION_add_n', ` 52 define(ADCSBB, adc) 53 define(func, mpn_add_n) 54 define(func_nc, mpn_add_nc)') 55ifdef(`OPERATION_sub_n', ` 56 define(ADCSBB, sbb) 57 define(func, mpn_sub_n) 58 define(func_nc, mpn_sub_nc)') 59 60MULFUNC_PROLOGUE(mpn_add_n mpn_add_nc mpn_sub_n mpn_sub_nc) 61 62ABI_SUPPORT(DOS64) 63ABI_SUPPORT(STD64) 64 65ASM_START() 66 TEXT 67 ALIGN(16) 68PROLOGUE(func) 69 FUNC_ENTRY(4) 70 xor %r8, %r8 71L(start): 72 mov (up), %r10 73 mov (vp), %r11 74 75 lea (up,n,8), up 76 lea (vp,n,8), vp 77 lea (rp,n,8), rp 78 mov R32(n), R32(%rax) 79 neg n 80 and $3, R32(%rax) 81 je L(b00) 82 add %rax, n C clear low rcx bits for jrcxz 83 cmp $2, R32(%rax) 84 jl L(b01) 85 je L(b10) 86 87L(b11): neg %r8 C set cy 88 jmp L(e11) 89 90L(b00): neg %r8 C set cy 91 mov %r10, %r8 92 mov %r11, %r9 93 lea 4(n), n 94 jmp L(e00) 95 96 nop 97 nop 98 nop 99L(b01): neg %r8 C set cy 100 jmp L(top) 101 102L(b10): neg %r8 C set cy 103 mov %r10, %r8 104 mov %r11, %r9 105 jmp L(e10) 106 107L(end): ADCSBB %r11, %r10 108 mov %r10, -8(rp) 109 mov R32(%rcx), R32(%rax) C clear eax, ecx contains 0 110 adc R32(%rax), R32(%rax) 111 FUNC_EXIT() 112 ret 113 114 ALIGN(16) 115L(top): jrcxz L(end) 116 mov (up,n,8), %r8 117 mov (vp,n,8), %r9 118 lea 4(n), n 119 ADCSBB %r11, %r10 120 mov %r10, -40(rp,n,8) 121L(e00): mov -24(up,n,8), %r10 122 mov -24(vp,n,8), %r11 123 ADCSBB %r9, %r8 124 mov %r8, -32(rp,n,8) 125L(e11): mov -16(up,n,8), %r8 126 mov -16(vp,n,8), %r9 127 ADCSBB %r11, %r10 128 mov %r10, -24(rp,n,8) 129L(e10): mov -8(up,n,8), %r10 130 mov -8(vp,n,8), %r11 131 ADCSBB %r9, %r8 132 mov %r8, -16(rp,n,8) 133 jmp L(top) 134EPILOGUE() 135 136PROLOGUE(func_nc) 137 FUNC_ENTRY(4) 138IFDOS(` mov 56(%rsp), %r8 ') 139 jmp L(start) 140EPILOGUE() 141 142