Home | History | Annotate | Line # | Download | only in builtins
      1      1.1  joerg /*===-- mulodi4.c - Implement __mulodi4 -----------------------------------===
      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 __mulodi4 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 /* Returns: a * b */
     18      1.1  joerg 
     19      1.1  joerg /* Effects: sets *overflow to 1  if a * b overflows */
     20      1.1  joerg 
     21  1.1.1.2  joerg COMPILER_RT_ABI di_int
     22      1.1  joerg __mulodi4(di_int a, di_int b, int* overflow)
     23      1.1  joerg {
     24      1.1  joerg     const int N = (int)(sizeof(di_int) * CHAR_BIT);
     25      1.1  joerg     const di_int MIN = (di_int)1 << (N-1);
     26      1.1  joerg     const di_int MAX = ~MIN;
     27      1.1  joerg     *overflow = 0;
     28      1.1  joerg     di_int result = a * b;
     29      1.1  joerg     if (a == MIN)
     30      1.1  joerg     {
     31      1.1  joerg         if (b != 0 && b != 1)
     32      1.1  joerg 	    *overflow = 1;
     33      1.1  joerg 	return result;
     34      1.1  joerg     }
     35      1.1  joerg     if (b == MIN)
     36      1.1  joerg     {
     37      1.1  joerg         if (a != 0 && a != 1)
     38      1.1  joerg 	    *overflow = 1;
     39      1.1  joerg         return result;
     40      1.1  joerg     }
     41      1.1  joerg     di_int sa = a >> (N - 1);
     42      1.1  joerg     di_int abs_a = (a ^ sa) - sa;
     43      1.1  joerg     di_int sb = b >> (N - 1);
     44      1.1  joerg     di_int abs_b = (b ^ sb) - sb;
     45      1.1  joerg     if (abs_a < 2 || abs_b < 2)
     46      1.1  joerg         return result;
     47      1.1  joerg     if (sa == sb)
     48      1.1  joerg     {
     49      1.1  joerg         if (abs_a > MAX / abs_b)
     50      1.1  joerg             *overflow = 1;
     51      1.1  joerg     }
     52      1.1  joerg     else
     53      1.1  joerg     {
     54      1.1  joerg         if (abs_a > MIN / -abs_b)
     55      1.1  joerg             *overflow = 1;
     56      1.1  joerg     }
     57      1.1  joerg     return result;
     58      1.1  joerg }
     59