From 11da511c784eca003deb90c23570f0873954e0de Mon Sep 17 00:00:00 2001 From: Duncan Wilkie Date: Sat, 18 Nov 2023 06:11:09 -0600 Subject: Initial commit. --- gmp-6.3.0/mpn/sh/add_n.asm | 59 +++++++++++++++++++++++++++++++++++ gmp-6.3.0/mpn/sh/sh2/addmul_1.asm | 65 +++++++++++++++++++++++++++++++++++++++ gmp-6.3.0/mpn/sh/sh2/mul_1.asm | 62 +++++++++++++++++++++++++++++++++++++ gmp-6.3.0/mpn/sh/sh2/submul_1.asm | 65 +++++++++++++++++++++++++++++++++++++++ gmp-6.3.0/mpn/sh/sub_n.asm | 59 +++++++++++++++++++++++++++++++++++ 5 files changed, 310 insertions(+) create mode 100644 gmp-6.3.0/mpn/sh/add_n.asm create mode 100644 gmp-6.3.0/mpn/sh/sh2/addmul_1.asm create mode 100644 gmp-6.3.0/mpn/sh/sh2/mul_1.asm create mode 100644 gmp-6.3.0/mpn/sh/sh2/submul_1.asm create mode 100644 gmp-6.3.0/mpn/sh/sub_n.asm (limited to 'gmp-6.3.0/mpn/sh') diff --git a/gmp-6.3.0/mpn/sh/add_n.asm b/gmp-6.3.0/mpn/sh/add_n.asm new file mode 100644 index 0000000..79d17d0 --- /dev/null +++ b/gmp-6.3.0/mpn/sh/add_n.asm @@ -0,0 +1,59 @@ +dnl SH mpn_add_n -- Add two limb vectors of the same length > 0 and store sum +dnl in a third limb vector. + +dnl Copyright 1995, 1997, 2000, 2011 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 +C rp r4 +C up r5 +C vp r6 +C n r7 + +changecom(blah) C disable # to make all C comments below work + +ASM_START() +PROLOGUE(mpn_add_n) + mov #0,r3 C clear cy save reg + +L(top): mov.l @r5+,r1 + mov.l @r6+,r2 + shlr r3 C restore cy + addc r2,r1 + movt r3 C save cy + mov.l r1,@r4 + dt r7 + bf.s L(top) + add #4,r4 + + rts + mov r3,r0 C return carry-out from most significant limb +EPILOGUE() diff --git a/gmp-6.3.0/mpn/sh/sh2/addmul_1.asm b/gmp-6.3.0/mpn/sh/sh2/addmul_1.asm new file mode 100644 index 0000000..c914b29 --- /dev/null +++ b/gmp-6.3.0/mpn/sh/sh2/addmul_1.asm @@ -0,0 +1,65 @@ +dnl SH2 mpn_addmul_1 -- Multiply a limb vector with a limb and add the result +dnl to a second limb vector. + +dnl Copyright 1995, 2000, 2011 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 +C res_ptr r4 +C s1_ptr r5 +C size r6 +C s2_limb r7 + +changecom(blah) C disable # to make all C comments below work + +ASM_START() +PROLOGUE(mpn_addmul_1) + mov #0,r2 C cy_limb = 0 + mov #0,r0 C Keep r0 = 0 for entire loop + clrt + +L(top): mov.l @r5+,r3 + dmulu.l r3,r7 + sts macl,r1 + addc r2,r1 C lo_prod += old cy_limb + sts mach,r2 C new cy_limb = hi_prod + mov.l @r4,r3 + addc r0,r2 C cy_limb += T, T = 0 + addc r3,r1 + addc r0,r2 C cy_limb += T, T = 0 + dt r6 + mov.l r1,@r4 + bf.s L(top) + add #4,r4 + + rts + mov r2,r0 +EPILOGUE() diff --git a/gmp-6.3.0/mpn/sh/sh2/mul_1.asm b/gmp-6.3.0/mpn/sh/sh2/mul_1.asm new file mode 100644 index 0000000..83548a6 --- /dev/null +++ b/gmp-6.3.0/mpn/sh/sh2/mul_1.asm @@ -0,0 +1,62 @@ +dnl SH2 mpn_mul_1 -- Multiply a limb vector with a limb and store the result +dnl in a second limb vector. + +dnl Copyright 1995, 2000, 2011 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 +C res_ptr r4 +C s1_ptr r5 +C size r6 +C s2_limb r7 + +changecom(blah) C disable # to make all C comments below work + +ASM_START() +PROLOGUE(mpn_mul_1) + mov #0,r2 C cy_limb = 0 + mov #0,r0 C Keep r0 = 0 for entire loop + clrt + +L(top): mov.l @r5+,r3 + dmulu.l r3,r7 + sts macl,r1 + addc r2,r1 + sts mach,r2 + addc r0,r2 C propagate carry to cy_limb (dt clobbers T) + dt r6 + mov.l r1,@r4 + bf.s L(top) + add #4,r4 + + rts + mov r2,r0 +EPILOGUE() diff --git a/gmp-6.3.0/mpn/sh/sh2/submul_1.asm b/gmp-6.3.0/mpn/sh/sh2/submul_1.asm new file mode 100644 index 0000000..bef2abd --- /dev/null +++ b/gmp-6.3.0/mpn/sh/sh2/submul_1.asm @@ -0,0 +1,65 @@ +dnl SH2 mpn_submul_1 -- Multiply a limb vector with a limb and subtract the +dnl result from a second limb vector. + +dnl Copyright 1995, 2000, 2011 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 +C res_ptr r4 +C s1_ptr r5 +C size r6 +C s2_limb r7 + +changecom(blah) C disable # to make all C comments below work + +ASM_START() +PROLOGUE(mpn_submul_1) + mov #0,r2 C cy_limb = 0 + mov #0,r0 C Keep r0 = 0 for entire loop + clrt + +L(top): mov.l @r5+,r3 + dmulu.l r3,r7 + sts macl,r1 + addc r2,r1 C lo_prod += old cy_limb + sts mach,r2 C new cy_limb = hi_prod + mov.l @r4,r3 + addc r0,r2 C cy_limb += T, T = 0 + subc r1,r3 + addc r0,r2 C cy_limb += T, T = 0 + dt r6 + mov.l r3,@r4 + bf.s L(top) + add #4,r4 + + rts + mov r2,r0 +EPILOGUE() diff --git a/gmp-6.3.0/mpn/sh/sub_n.asm b/gmp-6.3.0/mpn/sh/sub_n.asm new file mode 100644 index 0000000..465bc80 --- /dev/null +++ b/gmp-6.3.0/mpn/sh/sub_n.asm @@ -0,0 +1,59 @@ +dnl SH mpn_sub_n -- Subtract two limb vectors of the same length > 0 and store +dnl difference in a third limb vector. + +dnl Copyright 1995, 1997, 2000, 2011 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 +C rp r4 +C up r5 +C vp r6 +C n r7 + +changecom(blah) C disable # to make all C comments below work + +ASM_START() +PROLOGUE(mpn_sub_n) + mov #0,r3 C clear cy save reg + +L(top): mov.l @r5+,r1 + mov.l @r6+,r2 + shlr r3 C restore cy + subc r2,r1 + movt r3 C save cy + mov.l r1,@r4 + dt r7 + bf.s L(top) + add #4,r4 + + rts + mov r3,r0 C return carry-out from most significant limb +EPILOGUE() -- cgit v1.2.3