smallint.c 2.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475
  1. /*
  2. * This file is part of the MicroPython project, http://micropython.org/
  3. *
  4. * The MIT License (MIT)
  5. *
  6. * Copyright (c) 2013, 2014 Damien P. George
  7. *
  8. * Permission is hereby granted, free of charge, to any person obtaining a copy
  9. * of this software and associated documentation files (the "Software"), to deal
  10. * in the Software without restriction, including without limitation the rights
  11. * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
  12. * copies of the Software, and to permit persons to whom the Software is
  13. * furnished to do so, subject to the following conditions:
  14. *
  15. * The above copyright notice and this permission notice shall be included in
  16. * all copies or substantial portions of the Software.
  17. *
  18. * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  19. * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  20. * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  21. * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  22. * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  23. * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
  24. * THE SOFTWARE.
  25. */
  26. #include "py/smallint.h"
  27. bool mp_small_int_mul_overflow(mp_int_t x, mp_int_t y) {
  28. // Check for multiply overflow; see CERT INT32-C
  29. if (x > 0) { // x is positive
  30. if (y > 0) { // x and y are positive
  31. if (x > (MP_SMALL_INT_MAX / y)) {
  32. return true;
  33. }
  34. } else { // x positive, y nonpositive
  35. if (y < (MP_SMALL_INT_MIN / x)) {
  36. return true;
  37. }
  38. } // x positive, y nonpositive
  39. } else { // x is nonpositive
  40. if (y > 0) { // x is nonpositive, y is positive
  41. if (x < (MP_SMALL_INT_MIN / y)) {
  42. return true;
  43. }
  44. } else { // x and y are nonpositive
  45. if (x != 0 && y < (MP_SMALL_INT_MAX / x)) {
  46. return true;
  47. }
  48. } // End if x and y are nonpositive
  49. } // End if x is nonpositive
  50. return false;
  51. }
  52. mp_int_t mp_small_int_modulo(mp_int_t dividend, mp_int_t divisor) {
  53. // Python specs require that mod has same sign as second operand
  54. dividend %= divisor;
  55. if ((dividend < 0 && divisor > 0) || (dividend > 0 && divisor < 0)) {
  56. dividend += divisor;
  57. }
  58. return dividend;
  59. }
  60. mp_int_t mp_small_int_floor_divide(mp_int_t num, mp_int_t denom) {
  61. if (num >= 0) {
  62. if (denom < 0) {
  63. num += -denom - 1;
  64. }
  65. } else {
  66. if (denom >= 0) {
  67. num += -denom + 1;
  68. }
  69. }
  70. return num / denom;
  71. }