1 /* Cray PVP/IEEE mpn_sqr_basecase.
3 Copyright 2000, 2001 Free Software Foundation, Inc.
5 This file is part of the GNU MP Library.
7 The GNU MP Library is free software; you can redistribute it and/or modify
8 it under the terms of the GNU Lesser General Public License as published by
9 the Free Software Foundation; either version 3 of the License, or (at your
10 option) any later version.
12 The GNU MP Library is distributed in the hope that it will be useful, but
13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15 License for more details.
17 You should have received a copy of the GNU Lesser General Public License
18 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
20 /* This is just mpn_mul_basecase with trivial modifications. */
22 #include <intrinsics.h>
27 mpn_sqr_basecase (mp_ptr rp,
28 mp_srcptr up, mp_size_t un)
30 mp_limb_t cy[un + un];
32 mp_limb_t a, b, r, s0, s1, c0, c1;
36 for (i = 0; i < un + un; i++)
43 for (j = 0; j < un; j++)
51 c0 = ((a & r) | ((a | r) & ~s0)) >> 63;
55 for (i = 1; i < un; i++)
58 b = _int_mult_upper (up[i - 1], ul);
60 c0 = ((a & b) | ((a | b) & ~s0)) >> 63;
64 c1 = ((s0 & r) | ((s0 | r) & ~s1)) >> 63;
67 rp[j + un] = _int_mult_upper (up[un - 1], ul);
72 for (i = 1; i < un + un; i++)
81 /* If that second loop generated carry, handle that in scalar loop. */
85 for (i = 1; i < un + un; i++)