- Code/Resource
- Windows Develop
- Linux-Unix program
- Internet-Socket-Network
- Web Server
- Browser Client
- Ftp Server
- Ftp Client
- Browser Plugins
- Proxy Server
- Email Server
- Email Client
- WEB Mail
- Firewall-Security
- Telnet Server
- Telnet Client
- ICQ-IM-Chat
- Search Engine
- Sniffer Package capture
- Remote Control
- xml-soap-webservice
- P2P
- WEB(ASP,PHP,...)
- TCP/IP Stack
- SNMP
- Grid Computing
- SilverLight
- DNS
- Cluster Service
- Network Security
- Communication-Mobile
- Game Program
- Editor
- Multimedia program
- Graph program
- Compiler program
- Compress-Decompress algrithms
- Crypt_Decrypt algrithms
- Mathimatics-Numerical algorithms
- MultiLanguage
- Disk/Storage
- Java Develop
- assembly language
- Applications
- Other systems
- Database system
- Embeded-SCM Develop
- FlashMX/Flex
- source in ebook
- Delphi VCL
- OS Develop
- MiddleWare
- MPI
- MacOS develop
- LabView
- ELanguage
- Software/Tools
- E-Books
- Artical/Document
addmul_1.asm
Package: gmp-5.0.1.tar.gz [view]
Upload User: qaz666999
Upload Date: 2022-08-06
Package Size: 2570k
Code Size: 2k
Category:
Algorithm
Development Platform:
Unix_Linux
- dnl Alpha mpn_addmul_1 -- Multiply a limb vector with a limb and add the
- dnl result to a second limb vector.
- dnl Copyright 1992, 1994, 1995, 2000, 2002 Free Software Foundation, Inc.
- dnl This file is part of the GNU MP Library.
- dnl The GNU MP Library is free software; you can redistribute it and/or modify
- dnl it under the terms of the GNU Lesser General Public License as published
- dnl by the Free Software Foundation; either version 3 of the License, or (at
- dnl your option) any later version.
- dnl The GNU MP Library is distributed in the hope that it will be useful, but
- dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
- dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
- dnl License for more details.
- dnl You should have received a copy of the GNU Lesser General Public License
- dnl along with the GNU MP Library. If not, see http://www.gnu.org/licenses/.
- include(`../config.m4')
- C cycles/limb
- C EV4: 42
- C EV5: 18
- C EV6: 7
- C INPUT PARAMETERS
- C rp r16
- C up r17
- C n r18
- C vl r19
- ASM_START()
- PROLOGUE(mpn_addmul_1)
- ldq r2,0(r17) C r2 = s1_limb
- addq r17,8,r17 C s1_ptr++
- subq r18,1,r18 C size--
- mulq r2,r19,r3 C r3 = prod_low
- ldq r5,0(r16) C r5 = *res_ptr
- umulh r2,r19,r0 C r0 = prod_high
- beq r18,$Lend1 C jump if size was == 1
- ldq r2,0(r17) C r2 = s1_limb
- addq r17,8,r17 C s1_ptr++
- subq r18,1,r18 C size--
- addq r5,r3,r3
- cmpult r3,r5,r4
- stq r3,0(r16)
- addq r16,8,r16 C res_ptr++
- beq r18,$Lend2 C jump if size was == 2
- ALIGN(8)
- $Loop: mulq r2,r19,r3 C r3 = prod_low
- ldq r5,0(r16) C r5 = *res_ptr
- addq r4,r0,r0 C cy_limb = cy_limb + 'cy'
- subq r18,1,r18 C size--
- umulh r2,r19,r4 C r4 = cy_limb
- ldq r2,0(r17) C r2 = s1_limb
- addq r17,8,r17 C s1_ptr++
- addq r3,r0,r3 C r3 = cy_limb + prod_low
- cmpult r3,r0,r0 C r0 = carry from (cy_limb + prod_low)
- addq r5,r3,r3
- cmpult r3,r5,r5
- stq r3,0(r16)
- addq r16,8,r16 C res_ptr++
- addq r5,r0,r0 C combine carries
- bne r18,$Loop
- $Lend2: mulq r2,r19,r3 C r3 = prod_low
- ldq r5,0(r16) C r5 = *res_ptr
- addq r4,r0,r0 C cy_limb = cy_limb + 'cy'
- umulh r2,r19,r4 C r4 = cy_limb
- addq r3,r0,r3 C r3 = cy_limb + prod_low
- cmpult r3,r0,r0 C r0 = carry from (cy_limb + prod_low)
- addq r5,r3,r3
- cmpult r3,r5,r5
- stq r3,0(r16)
- addq r5,r0,r0 C combine carries
- addq r4,r0,r0 C cy_limb = prod_high + cy
- ret r31,(r26),1
- $Lend1: addq r5,r3,r3
- cmpult r3,r5,r5
- stq r3,0(r16)
- addq r0,r5,r0
- ret r31,(r26),1
- EPILOGUE(mpn_addmul_1)
- ASM_END()