mulodi4.c 1.5 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364
  1. /*===-- mulodi4.c - Implement __mulodi4 -----------------------------------===
  2. *
  3. * The LLVM Compiler Infrastructure
  4. *
  5. * This file is dual licensed under the MIT and the University of Illinois Open
  6. * Source Licenses. See LICENSE.TXT for details.
  7. *
  8. * ===----------------------------------------------------------------------===
  9. *
  10. * This file implements __mulodi4 for the compiler_rt library.
  11. *
  12. * ===----------------------------------------------------------------------===
  13. */
  14. #if 0
  15. #include "int_lib.h"
  16. #else
  17. #define COMPILER_RT_ABI
  18. #define di_int int64_t
  19. #include "torint.h"
  20. #endif
  21. /* Returns: a * b */
  22. /* Effects: sets *overflow to 1 if a * b overflows */
  23. COMPILER_RT_ABI di_int
  24. __mulodi4(di_int a, di_int b, int* overflow)
  25. {
  26. const int N = (int)(sizeof(di_int) * CHAR_BIT);
  27. const di_int MIN = (di_int)1 << (N-1);
  28. const di_int MAX = ~MIN;
  29. *overflow = 0;
  30. di_int result = a * b;
  31. if (a == MIN)
  32. {
  33. if (b != 0 && b != 1)
  34. *overflow = 1;
  35. return result;
  36. }
  37. if (b == MIN)
  38. {
  39. if (a != 0 && a != 1)
  40. *overflow = 1;
  41. return result;
  42. }
  43. di_int sa = a >> (N - 1);
  44. di_int abs_a = (a ^ sa) - sa;
  45. di_int sb = b >> (N - 1);
  46. di_int abs_b = (b ^ sb) - sb;
  47. if (abs_a < 2 || abs_b < 2)
  48. return result;
  49. if (sa == sb)
  50. {
  51. if (abs_a > MAX / abs_b)
  52. *overflow = 1;
  53. }
  54. else
  55. {
  56. if (abs_a > MIN / -abs_b)
  57. *overflow = 1;
  58. }
  59. return result;
  60. }