Quelle  submul_1.asm   Sprache: Masm

dnl  HP-PA 7100/7200 mpn_submul_1 -- Multiply a limb vector with a limb and
dnl  subtract the result from a second limb vector.

dnl  Copyright 1995, 2000-2003 Free Software Foundation, Inc.

dnl  This file is part of the GNU MP Library.
dnl
dnl  The GNU MP Library is free software; you can redistribute it and/or modify
dnl  it under the terms of either:
dnl
dnl    * the GNU Lesser General Public License as published by the Free
dnl      Software Foundation; either version 3 of the License, or (at your
dnl      option) any later version.
dnl
dnl  or
dnl
dnl    * the GNU General Public License as published by the Free Software
dnl      Foundation; either version 2 of the License, or (at your option) any
dnl      later version.
dnl
dnl  or both in parallel, as here.
dnl
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 General Public License
dnl  for more details.
dnl
dnl  You should have received copies of the GNU General Public License and the
dnl  GNU Lesser General Public License along with the GNU MP Library.  If not,
dnl  see https://www.gnu.org/licenses/.

include(`../config.m4')

C INPUT PARAMETERS
define(`res_ptr',`%r26')
define(`s1_ptr',`%r25')
define(`size_param',`%r24')
define(`s2_limb',`%r23')

define(`cylimb',`%r28')
define(`s0',`%r19')
define(`s1',`%r20')
define(`s2',`%r3')
define(`s3',`%r4')
define(`lo0',`%r21')
define(`lo1',`%r5')
define(`lo2',`%r6')
define(`lo3',`%r7')
define(`hi0',`%r22')
define(`hi1',`%r23')    C safe to reuse
define(`hi2',`%r29')
define(`hi3',`%r1')

ASM_START()
PROLOGUE(mpn_submul_1)
C .callinfo frame=128,no_calls

 ldo 128(%r30),%r30
 stws s2_limb,-16(%r30)
 add  %r0,%r0,cylimb   C clear cy and cylimb
 addib,< -4,size_param,L(few_limbs)
 fldws -16(%r30),%fr31R

 ldo -112(%r30),%r31
 stw %r3,-96(%r30)
 stw %r4,-92(%r30)
 stw %r5,-88(%r30)
 stw %r6,-84(%r30)
 stw %r7,-80(%r30)

 bb,>=,n  s1_ptr,29,L(0)

 fldws,ma 4(s1_ptr),%fr4
 ldws  0(res_ptr),s0
 xmpyu  %fr4,%fr31R,%fr5
 fstds  %fr5,-16(%r31)
 ldws -16(%r31),cylimb
 ldws -12(%r31),lo0
 sub  s0,lo0,s0
 add  s0,lo0,%r0   C invert cy
 addib,< -1,size_param,L(few_limbs)
 stws,ma  s0,4(res_ptr)

C start software pipeline ----------------------------------------------------
LDEF(0)
 fldds,ma 8(s1_ptr),%fr4
 fldds,ma 8(s1_ptr),%fr8

 xmpyu  %fr4L,%fr31R,%fr5
 xmpyu  %fr4R,%fr31R,%fr6
 xmpyu  %fr8L,%fr31R,%fr9
 xmpyu  %fr8R,%fr31R,%fr10

 fstds  %fr5,-16(%r31)
 fstds  %fr6,-8(%r31)
 fstds  %fr9,0(%r31)
 fstds  %fr10,8(%r31)

 ldws   -16(%r31),hi0
 ldws   -12(%r31),lo0
 ldws -8(%r31),hi1
 ldws -4(%r31),lo1
 ldws  0(%r31),hi2
 ldws  4(%r31),lo2
 ldws  8(%r31),hi3
 ldws 12(%r31),lo3

 addc  lo0,cylimb,lo0
 addc  lo1,hi0,lo1
 addc  lo2,hi1,lo2
 addc  lo3,hi2,lo3

 addib,<  -4,size_param,L(end)
 addc  %r0,hi3,cylimb   C propagate carry into cylimb
C main loop ------------------------------------------------------------------
LDEF(loop)
 fldds,ma 8(s1_ptr),%fr4
 fldds,ma 8(s1_ptr),%fr8

 ldws  0(res_ptr),s0
 xmpyu  %fr4L,%fr31R,%fr5
 ldws  4(res_ptr),s1
 xmpyu  %fr4R,%fr31R,%fr6
 ldws  8(res_ptr),s2
 xmpyu  %fr8L,%fr31R,%fr9
 ldws 12(res_ptr),s3
 xmpyu  %fr8R,%fr31R,%fr10

 fstds  %fr5,-16(%r31)
 sub  s0,lo0,s0
 fstds  %fr6,-8(%r31)
 subb  s1,lo1,s1
 fstds  %fr9,0(%r31)
 subb  s2,lo2,s2
 fstds  %fr10,8(%r31)
 subb  s3,lo3,s3
 subb  %r0,%r0,lo0   C these two insns ...
 add  lo0,lo0,%r0   C ... just invert cy

 ldws   -16(%r31),hi0
 ldws   -12(%r31),lo0
 ldws -8(%r31),hi1
 ldws -4(%r31),lo1
 ldws  0(%r31),hi2
 ldws  4(%r31),lo2
 ldws  8(%r31),hi3
 ldws 12(%r31),lo3

 addc  lo0,cylimb,lo0
 stws,ma  s0,4(res_ptr)
 addc  lo1,hi0,lo1
 stws,ma  s1,4(res_ptr)
 addc  lo2,hi1,lo2
 stws,ma  s2,4(res_ptr)
 addc  lo3,hi2,lo3
 stws,ma  s3,4(res_ptr)

 addib,>= -4,size_param,L(loop)
 addc  %r0,hi3,cylimb   C propagate carry into cylimb
C finish software pipeline ---------------------------------------------------
LDEF(end)
 ldws  0(res_ptr),s0
 ldws  4(res_ptr),s1
 ldws  8(res_ptr),s2
 ldws 12(res_ptr),s3

 sub  s0,lo0,s0
 stws,ma  s0,4(res_ptr)
 subb  s1,lo1,s1
 stws,ma  s1,4(res_ptr)
 subb  s2,lo2,s2
 stws,ma  s2,4(res_ptr)
 subb  s3,lo3,s3
 stws,ma  s3,4(res_ptr)
 subb  %r0,%r0,lo0   C these two insns ...
 add  lo0,lo0,%r0   C ... invert cy

C restore callee-saves registers ---------------------------------------------
 ldw -96(%r30),%r3
 ldw -92(%r30),%r4
 ldw -88(%r30),%r5
 ldw -84(%r30),%r6
 ldw -80(%r30),%r7

LDEF(few_limbs)
 addib,=,n 4,size_param,L(ret)

LDEF(loop2)
 fldws,ma 4(s1_ptr),%fr4
 ldws  0(res_ptr),s0
 xmpyu  %fr4,%fr31R,%fr5
 fstds  %fr5,-16(%r30)
 ldws -16(%r30),hi0
 ldws -12(%r30),lo0
 addc  lo0,cylimb,lo0
 addc  %r0,hi0,cylimb
 sub  s0,lo0,s0
 add  s0,lo0,%r0   C invert cy
 stws,ma  s0,4(res_ptr)
 addib,<> -1,size_param,L(loop2)
 nop

LDEF(ret)
 addc  %r0,cylimb,cylimb
 bv  0(%r2)
 ldo  -128(%r30),%r30
EPILOGUE(mpn_submul_1)