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