summaryrefslogtreecommitdiff
path: root/gmp/mpn/generic/sqr.c
diff options
context:
space:
mode:
Diffstat (limited to 'gmp/mpn/generic/sqr.c')
-rw-r--r--gmp/mpn/generic/sqr.c99
1 files changed, 0 insertions, 99 deletions
diff --git a/gmp/mpn/generic/sqr.c b/gmp/mpn/generic/sqr.c
deleted file mode 100644
index 3743761f78..0000000000
--- a/gmp/mpn/generic/sqr.c
+++ /dev/null
@@ -1,99 +0,0 @@
-/* mpn_sqr -- square natural numbers.
-
-Copyright 1991, 1993, 1994, 1996-2003, 2005, 2008, 2009 Free Software
-Foundation, Inc.
-
-This file is part of the GNU MP Library.
-
-The GNU MP Library is free software; you can redistribute it and/or modify
-it under the terms of either:
-
- * the GNU Lesser General Public License as published by the Free
- Software Foundation; either version 3 of the License, or (at your
- option) any later version.
-
-or
-
- * the GNU General Public License as published by the Free Software
- Foundation; either version 2 of the License, or (at your option) any
- later version.
-
-or both in parallel, as here.
-
-The GNU MP Library is distributed in the hope that it will be useful, but
-WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
-or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
-for more details.
-
-You should have received copies of the GNU General Public License and the
-GNU Lesser General Public License along with the GNU MP Library. If not,
-see https://www.gnu.org/licenses/. */
-
-#include "gmp.h"
-#include "gmp-impl.h"
-#include "longlong.h"
-
-void
-mpn_sqr (mp_ptr p, mp_srcptr a, mp_size_t n)
-{
- ASSERT (n >= 1);
- ASSERT (! MPN_OVERLAP_P (p, 2 * n, a, n));
-
- if (BELOW_THRESHOLD (n, SQR_BASECASE_THRESHOLD))
- { /* mul_basecase is faster than sqr_basecase on small sizes sometimes */
- mpn_mul_basecase (p, a, n, a, n);
- }
- else if (BELOW_THRESHOLD (n, SQR_TOOM2_THRESHOLD))
- {
- mpn_sqr_basecase (p, a, n);
- }
- else if (BELOW_THRESHOLD (n, SQR_TOOM3_THRESHOLD))
- {
- /* Allocate workspace of fixed size on stack: fast! */
- mp_limb_t ws[mpn_toom2_sqr_itch (SQR_TOOM3_THRESHOLD_LIMIT-1)];
- ASSERT (SQR_TOOM3_THRESHOLD <= SQR_TOOM3_THRESHOLD_LIMIT);
- mpn_toom2_sqr (p, a, n, ws);
- }
- else if (BELOW_THRESHOLD (n, SQR_TOOM4_THRESHOLD))
- {
- mp_ptr ws;
- TMP_SDECL;
- TMP_SMARK;
- ws = TMP_SALLOC_LIMBS (mpn_toom3_sqr_itch (n));
- mpn_toom3_sqr (p, a, n, ws);
- TMP_SFREE;
- }
- else if (BELOW_THRESHOLD (n, SQR_TOOM6_THRESHOLD))
- {
- mp_ptr ws;
- TMP_SDECL;
- TMP_SMARK;
- ws = TMP_SALLOC_LIMBS (mpn_toom4_sqr_itch (n));
- mpn_toom4_sqr (p, a, n, ws);
- TMP_SFREE;
- }
- else if (BELOW_THRESHOLD (n, SQR_TOOM8_THRESHOLD))
- {
- mp_ptr ws;
- TMP_SDECL;
- TMP_SMARK;
- ws = TMP_SALLOC_LIMBS (mpn_toom6_sqr_itch (n));
- mpn_toom6_sqr (p, a, n, ws);
- TMP_SFREE;
- }
- else if (BELOW_THRESHOLD (n, SQR_FFT_THRESHOLD))
- {
- mp_ptr ws;
- TMP_DECL;
- TMP_MARK;
- ws = TMP_ALLOC_LIMBS (mpn_toom8_sqr_itch (n));
- mpn_toom8_sqr (p, a, n, ws);
- TMP_FREE;
- }
- else
- {
- /* The current FFT code allocates its own space. That should probably
- change. */
- mpn_fft_mul (p, a, n, a, n);
- }
-}