bn_mp_reduce_2k_setup.c 656 B

1234567891011121314151617181920212223242526272829303132
  1. #include "tommath_private.h"
  2. #ifdef BN_MP_REDUCE_2K_SETUP_C
  3. /* LibTomMath, multiple-precision integer library -- Tom St Denis */
  4. /* SPDX-License-Identifier: Unlicense */
  5. /* determines the setup value */
  6. mp_err mp_reduce_2k_setup(const mp_int *a, mp_digit *d)
  7. {
  8. mp_err err;
  9. mp_int tmp;
  10. int p;
  11. if ((err = mp_init(&tmp)) != MP_OKAY) {
  12. return err;
  13. }
  14. p = mp_count_bits(a);
  15. if ((err = mp_2expt(&tmp, p)) != MP_OKAY) {
  16. mp_clear(&tmp);
  17. return err;
  18. }
  19. if ((err = s_mp_sub(&tmp, a, &tmp)) != MP_OKAY) {
  20. mp_clear(&tmp);
  21. return err;
  22. }
  23. *d = tmp.dp[0];
  24. mp_clear(&tmp);
  25. return MP_OKAY;
  26. }
  27. #endif