12345678910111213141516171819202122232425262728 |
- #include "tommath_private.h"
- #ifdef BN_MP_SQR_C
- /* LibTomMath, multiple-precision integer library -- Tom St Denis */
- /* SPDX-License-Identifier: Unlicense */
- /* computes b = a*a */
- mp_err mp_sqr(const mp_int *a, mp_int *b)
- {
- mp_err err;
- if (MP_HAS(S_MP_TOOM_SQR) && /* use Toom-Cook? */
- (a->used >= MP_TOOM_SQR_CUTOFF)) {
- err = s_mp_toom_sqr(a, b);
- } else if (MP_HAS(S_MP_KARATSUBA_SQR) && /* Karatsuba? */
- (a->used >= MP_KARATSUBA_SQR_CUTOFF)) {
- err = s_mp_karatsuba_sqr(a, b);
- } else if (MP_HAS(S_MP_SQR_FAST) && /* can we use the fast comba multiplier? */
- (((a->used * 2) + 1) < MP_WARRAY) &&
- (a->used < (MP_MAXFAST / 2))) {
- err = s_mp_sqr_fast(a, b);
- } else if (MP_HAS(S_MP_SQR)) {
- err = s_mp_sqr(a, b);
- } else {
- err = MP_VAL;
- }
- b->sign = MP_ZPOS;
- return err;
- }
- #endif
|