sparc64.h
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:8k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. /* UltraSPARC 64 support macros.
  2.    THE FUNCTIONS IN THIS FILE ARE FOR INTERNAL USE ONLY.  THEY'RE ALMOST
  3.    CERTAIN TO BE SUBJECT TO INCOMPATIBLE CHANGES OR DISAPPEAR COMPLETELY IN
  4.    FUTURE GNU MP RELEASES.
  5. Copyright 2003 Free Software Foundation, Inc.
  6. This file is part of the GNU MP Library.
  7. The GNU MP Library is free software; you can redistribute it and/or modify
  8. it under the terms of the GNU Lesser General Public License as published by
  9. the Free Software Foundation; either version 3 of the License, or (at your
  10. option) any later version.
  11. The GNU MP Library is distributed in the hope that it will be useful, but
  12. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  13. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  14. License for more details.
  15. You should have received a copy of the GNU Lesser General Public License
  16. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  17. #define LOW32(x)   ((x) & 0xFFFFFFFF)
  18. #define HIGH32(x)  ((x) >> 32)
  19. /* Halfword number i in src is accessed as src[i+HALF_ENDIAN_ADJ(i)].
  20.    Plain src[i] would be incorrect in big endian, HALF_ENDIAN_ADJ has the
  21.    effect of swapping the two halves in this case.  */
  22. #if HAVE_LIMB_BIG_ENDIAN
  23. #define HALF_ENDIAN_ADJ(i)  (1 - (((i) & 1) << 1))   /* +1 even, -1 odd */
  24. #endif
  25. #if HAVE_LIMB_LITTLE_ENDIAN
  26. #define HALF_ENDIAN_ADJ(i)  0                        /* no adjust */
  27. #endif
  28. #ifndef HALF_ENDIAN_ADJ
  29. Error, error, unknown limb endianness;
  30. #endif
  31. /* umul_ppmm_lowequal sets h to the high limb of q*d, assuming the low limb
  32.    of that product is equal to l.  dh and dl are the 32-bit halves of d.
  33.    |-----high----||----low-----|
  34.    +------+------+
  35.    |             |                 ph = qh * dh
  36.    +------+------+
  37.           +------+------+
  38.           |             |          pm1 = ql * dh
  39.           +------+------+
  40.           +------+------+
  41.           |             |          pm2 = qh * dl
  42.           +------+------+
  43.                  +------+------+
  44.                  |             |   pl = ql * dl (not calculated)
  45.                  +------+------+
  46.    Knowing that the low 64 bits is equal to l means that LOW(pm1) + LOW(pm2)
  47.    + HIGH(pl) == HIGH(l).  The only thing we need from those product parts
  48.    is whether they produce a carry into the high.
  49.    pm_l = LOW(pm1)+LOW(pm2) is done to contribute its carry, then the only
  50.    time there's a further carry from LOW(pm_l)+HIGH(pl) is if LOW(pm_l) >
  51.    HIGH(l).  pl is never actually calculated.  */
  52. #define umul_ppmm_lowequal(h, q, d, dh, dl, l)  
  53.   do {                                          
  54.     mp_limb_t  ql, qh, ph, pm1, pm2, pm_l;      
  55.     ASSERT (dh == HIGH32(d));                   
  56.     ASSERT (dl == LOW32(d));                    
  57.     ASSERT (q*d == l);                          
  58.                                                 
  59.     ql = LOW32 (q);                             
  60.     qh = HIGH32 (q);                            
  61.                                                 
  62.     pm1 = ql * dh;                              
  63.     pm2 = qh * dl;                              
  64.     ph  = qh * dh;                              
  65.                                                 
  66.     pm_l = LOW32 (pm1) + LOW32 (pm2);           
  67.                                                 
  68.     (h) = ph + HIGH32 (pm1) + HIGH32 (pm2)      
  69.       + HIGH32 (pm_l) + ((pm_l << 32) > l);     
  70.                                                 
  71.     ASSERT_HIGH_PRODUCT (h, q, d);              
  72.   } while (0)
  73. /* Set h to the high of q*d, assuming the low limb of that product is equal
  74.    to l, and that d fits in 32-bits.
  75.    |-----high----||----low-----|
  76.           +------+------+
  77.           |             |          pm = qh * dl
  78.           +------+------+
  79.                  +------+------+
  80.                  |             |   pl = ql * dl (not calculated)
  81.                  +------+------+
  82.    Knowing that LOW(pm) + HIGH(pl) == HIGH(l) (mod 2^32) means that the only
  83.    time there's a carry from that sum is when LOW(pm) > HIGH(l).  There's no
  84.    need to calculate pl to determine this.  */
  85. #define umul_ppmm_half_lowequal(h, q, d, l)     
  86.   do {                                          
  87.     mp_limb_t pm;                               
  88.     ASSERT (q*d == l);                          
  89.     ASSERT (HIGH32(d) == 0);                    
  90.                                                 
  91.     pm = HIGH32(q) * d;                         
  92.     (h) = HIGH32(pm) + ((pm << 32) > l);        
  93.     ASSERT_HIGH_PRODUCT (h, q, d);              
  94.   } while (0)
  95. /* check that h is the high limb of x*y */
  96. #if WANT_ASSERT
  97. #define ASSERT_HIGH_PRODUCT(h, x, y)    
  98.   do {                                  
  99.     mp_limb_t  want_h, dummy;           
  100.     umul_ppmm (want_h, dummy, x, y);    
  101.     ASSERT (h == want_h);               
  102.   } while (0)
  103. #else
  104. #define ASSERT_HIGH_PRODUCT(h, q, d)    
  105.   do { } while (0)
  106. #endif
  107. /* Count the leading zeros on a limb, but assuming it fits in 32 bits.
  108.    The count returned will be in the range 32 to 63.
  109.    This is the 32-bit generic C count_leading_zeros from longlong.h. */
  110. #define count_leading_zeros_32(count, x)                                      
  111.   do {                                                                        
  112.     mp_limb_t  __xr = (x);                                                    
  113.     unsigned   __a;                                                           
  114.     ASSERT ((x) != 0);                                                        
  115.     ASSERT ((x) <= CNST_LIMB(0xFFFFFFFF));                                    
  116.     __a = __xr < ((UWtype) 1 << 16) ? (__xr < ((UWtype) 1 << 8) ? 1 : 8 + 1)  
  117.       : (__xr < ((UWtype) 1 << 24)  ? 16 + 1 : 24 + 1);                       
  118.                                                                               
  119.     (count) = W_TYPE_SIZE + 1 - __a - __clz_tab[__xr >> __a];                 
  120.   } while (0)
  121. /* Set inv to a 32-bit inverse floor((b*(b-d)-1) / d), knowing that d fits
  122.    32 bits and is normalized (high bit set).  */
  123. #define invert_half_limb(inv, d)                
  124.   do {                                          
  125.     mp_limb_t  _n;                              
  126.     ASSERT ((d) <= 0xFFFFFFFF);                 
  127.     ASSERT ((d) & 0x80000000);                  
  128.     _n = (((mp_limb_t) -(d)) << 32) - 1;        
  129.     (inv) = (mp_limb_t) (unsigned) (_n / (d));  
  130.   } while (0)
  131. /* Divide nh:nl by d, setting q to the quotient and r to the remainder.
  132.    q, r, nh and nl are 32-bits each, d_limb is 32-bits but in an mp_limb_t,
  133.    dinv_limb is similarly a 32-bit inverse but in an mp_limb_t.  */
  134. #define udiv_qrnnd_half_preinv(q, r, nh, nl, d_limb, dinv_limb)         
  135.   do {                                                                  
  136.     unsigned   _n2, _n10, _n1, _nadj, _q11n, _xh, _r, _q;               
  137.     mp_limb_t  _n, _x;                                                  
  138.     ASSERT (d_limb <= 0xFFFFFFFF);                                      
  139.     ASSERT (dinv_limb <= 0xFFFFFFFF);                                   
  140.     ASSERT (d_limb & 0x80000000);                                       
  141.     ASSERT (nh < d_limb);                                               
  142.     _n10 = (nl);                                                        
  143.     _n2 = (nh);                                                         
  144.     _n1 = (int) _n10 >> 31;                                             
  145.     _nadj = _n10 + (_n1 & d_limb);                                      
  146.     _x = dinv_limb * (_n2 - _n1) + _nadj;                               
  147.     _q11n = ~(_n2 + HIGH32 (_x));             /* -q1-1 */               
  148.     _n = ((mp_limb_t) _n2 << 32) + _n10;                                
  149.     _x = _n + d_limb * _q11n;                 /* n-q1*d-d */            
  150.     _xh = HIGH32 (_x) - d_limb;               /* high(n-q1*d-d) */      
  151.     ASSERT (_xh == 0 || _xh == ~0);                                     
  152.     _r = _x + (d_limb & _xh);                 /* addback */             
  153.     _q = _xh - _q11n;                         /* q1+1-addback */        
  154.     ASSERT (_r < d_limb);                                               
  155.     ASSERT (d_limb * _q + _r == _n);                                    
  156.     (r) = _r;                                                           
  157.     (q) = _q;                                                           
  158.   } while (0)