bn_mp_lshd.c 1.6 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667
  1. #include <tommath.h>
  2. #ifdef BN_MP_LSHD_C
  3. /* LibTomMath, multiple-precision integer library -- Tom St Denis
  4. *
  5. * LibTomMath is a library that provides multiple-precision
  6. * integer arithmetic as well as number theoretic functionality.
  7. *
  8. * The library was designed directly after the MPI library by
  9. * Michael Fromberger but has been written from scratch with
  10. * additional optimizations in place.
  11. *
  12. * The library is free for all purposes without any express
  13. * guarantee it works.
  14. *
  15. * Tom St Denis, tomstdenis@gmail.com, http://math.libtomcrypt.com
  16. */
  17. /* shift left a certain amount of digits */
  18. int mp_lshd (mp_int * a, int b)
  19. {
  20. int x, res;
  21. /* if its less than zero return */
  22. if (b <= 0) {
  23. return MP_OKAY;
  24. }
  25. /* grow to fit the new digits */
  26. if (a->alloc < a->used + b) {
  27. if ((res = mp_grow (a, a->used + b)) != MP_OKAY) {
  28. return res;
  29. }
  30. }
  31. {
  32. register mp_digit *top, *bottom;
  33. /* increment the used by the shift amount then copy upwards */
  34. a->used += b;
  35. /* top */
  36. top = a->dp + a->used - 1;
  37. /* base */
  38. bottom = a->dp + a->used - 1 - b;
  39. /* much like mp_rshd this is implemented using a sliding window
  40. * except the window goes the otherway around. Copying from
  41. * the bottom to the top. see bn_mp_rshd.c for more info.
  42. */
  43. for (x = a->used - 1; x >= b; x--) {
  44. *top-- = *bottom--;
  45. }
  46. /* zero the lower digits */
  47. top = a->dp;
  48. for (x = 0; x < b; x++) {
  49. *top++ = 0;
  50. }
  51. }
  52. return MP_OKAY;
  53. }
  54. #endif
  55. /* $Source: /cvs/libtom/libtommath/bn_mp_lshd.c,v $ */
  56. /* $Revision: 1.3 $ */
  57. /* $Date: 2006/03/31 14:18:44 $ */