bn_mp_sqr.c 881 B

12345678910111213141516171819202122232425262728
  1. #include "tommath_private.h"
  2. #ifdef BN_MP_SQR_C
  3. /* LibTomMath, multiple-precision integer library -- Tom St Denis */
  4. /* SPDX-License-Identifier: Unlicense */
  5. /* computes b = a*a */
  6. mp_err mp_sqr(const mp_int *a, mp_int *b)
  7. {
  8. mp_err err;
  9. if (MP_HAS(S_MP_TOOM_SQR) && /* use Toom-Cook? */
  10. (a->used >= MP_TOOM_SQR_CUTOFF)) {
  11. err = s_mp_toom_sqr(a, b);
  12. } else if (MP_HAS(S_MP_KARATSUBA_SQR) && /* Karatsuba? */
  13. (a->used >= MP_KARATSUBA_SQR_CUTOFF)) {
  14. err = s_mp_karatsuba_sqr(a, b);
  15. } else if (MP_HAS(S_MP_SQR_FAST) && /* can we use the fast comba multiplier? */
  16. (((a->used * 2) + 1) < MP_WARRAY) &&
  17. (a->used < (MP_MAXFAST / 2))) {
  18. err = s_mp_sqr_fast(a, b);
  19. } else if (MP_HAS(S_MP_SQR)) {
  20. err = s_mp_sqr(a, b);
  21. } else {
  22. err = MP_VAL;
  23. }
  24. b->sign = MP_ZPOS;
  25. return err;
  26. }
  27. #endif