1  /* { dg-require-effective-target int32plus } */
       2  
       3  typedef long unsigned int size_t;
       4  extern void *memset (void *__s, int __c, size_t __n) __attribute__ ((__nothrow__)) __attribute__ ((__nonnull__ (1)));
       5  
       6  static void
       7  encode (words, low, hi)
       8       long *words;
       9       unsigned long low;
      10       long hi;
      11  {
      12    words[0] = ((low) & (((unsigned long) 1 << (sizeof(unsigned long) / 2)) - 1));
      13    words[1] = ((unsigned long) (low) >> sizeof(unsigned long) / 2);
      14    words[2] = ((hi) & (((unsigned long) 1 << (sizeof(unsigned long) / 2)) - 1));
      15    words[3] = ((unsigned long) (hi) >> sizeof(unsigned long) / 2);
      16  }
      17  
      18  static void
      19  decode (words, low, hi)
      20       long *words;
      21       unsigned long *low;
      22       long *hi;
      23  {
      24    *low = words[0] + words[1] * ((unsigned long) 1 << sizeof(unsigned long) / 2);
      25    *hi = words[2] + words[3] * ((unsigned long) 1 << sizeof(unsigned long) / 2);
      26  }
      27  
      28  int
      29  neg_double (l1, h1, lv, hv)
      30       unsigned long l1;
      31       long h1;
      32       unsigned long *lv;
      33       long *hv;
      34  {
      35    if (l1 == 0)
      36      {
      37        *lv = 0;
      38        *hv = - h1;
      39        return (*hv & h1) < 0;
      40      }
      41    else
      42      {
      43        *lv = -l1;
      44        *hv = ~h1;
      45        return 0;
      46      }
      47  }
      48  
      49  int
      50  add_double (l1, h1, l2, h2, lv, hv)
      51       unsigned long l1, l2;
      52       long h1, h2;
      53       unsigned long *lv;
      54       long *hv;
      55  {
      56    unsigned long l;
      57    long h;
      58  
      59    l = l1 + l2;
      60    h = h1 + h2 + (l < l1);
      61  
      62    *lv = l;
      63    *hv = h;
      64    return ((~((h1) ^ (h2)) & ((h1) ^ (h))) < 0);
      65  }
      66  
      67  int
      68  mul_double (l1, h1, l2, h2, lv, hv)
      69       unsigned long l1, l2;
      70       long h1, h2;
      71       unsigned long *lv;
      72       long *hv;
      73  {
      74    long arg1[4];
      75    long arg2[4];
      76    long prod[4 * 2];
      77    unsigned long carry;
      78    int i, j, k;
      79    unsigned long toplow, neglow;
      80    long tophigh, neghigh;
      81  
      82    encode (arg1, l1, h1);
      83    encode (arg2, l2, h2);
      84  
      85    memset ((char *) prod, 0, sizeof prod);
      86  
      87    for (i = 0; i < 4; i++)
      88      {
      89        carry = 0;
      90        for (j = 0; j < 4; j++)
      91  	{
      92  	  k = i + j;
      93  
      94  	  carry += arg1[i] * arg2[j];
      95  
      96  	  carry += prod[k];
      97  	  prod[k] = ((carry) & (((unsigned long) 1 << (sizeof(unsigned long) / 2)) - 1));
      98  	  carry = ((unsigned long) (carry) >> sizeof(unsigned long) / 2);
      99  	}
     100        prod[i + 4] = carry;
     101      }
     102  
     103    decode (prod, lv, hv);
     104  
     105  
     106  
     107    decode (prod + 4, &toplow, &tophigh);
     108    if (h1 < 0)
     109      {
     110        neg_double (l2, h2, &neglow, &neghigh);
     111        add_double (neglow, neghigh, toplow, tophigh, &toplow, &tophigh);
     112      }
     113    if (h2 < 0)
     114      {
     115        neg_double (l1, h1, &neglow, &neghigh);
     116        add_double (neglow, neghigh, toplow, tophigh, &toplow, &tophigh);
     117      }
     118    return (*hv < 0 ? ~(toplow & tophigh) : toplow | tophigh) != 0;
     119  }