Category

Development Platform
Home > SourceCode/Document > Mathimatics-Numerical algorithms > View Code
_sgnu_add.c
Package: leda.tar.gz　[view]
Upload User: gzelex
Upload Date: 2007-01-07
Package Size: 707k
Code Size: 2k
Category:

Mathimatics-Numerical algorithms

Development Platform:

MultiPlatform

_sgnu_add.c：Code Content
  1. /*******************************************************************************
  2. +
  3. +  LEDA-R  3.2.3
  4. +
  5. +  _sgnu_add.c
  6. +
  7. +  Copyright (c) 1995  by  Max-Planck-Institut fuer Informatik
  8. +  Im Stadtwald, 66123 Saarbruecken, Germany     
  9. +  All rights reserved.
  10. + 
  11. *******************************************************************************/
  13. #define MAX_WORD     0xFFFFFFFF
  15. typedef unsigned long  word;
  16. typedef unsigned int sz_t;
  18. extern sz_t School_Add(word *a, sz_t a_used, word *b, sz_t b_used, word* sum)
  19. {
  20.   /* compute sum = a + b   (a_used >= b_used) */
  22.   int n = b_used / 16;
  23.   word* a_stop = a + a_used;
  24.   word* s_last = sum + a_used;
  25.   word aa;
  26.   word bb;  /* is used for b and to save carry (asm) */
  28. #define ADD_LOOP_BODY {
  29. aa = *a++; bb = *b++;
  30. asm volatile ("addxcc  %2, %1,  %0" : "=r"(aa)    : "0"(aa), "r"(bb) );
  31. *sum++ = aa; }
  33. asm volatile ("addcc  %%g0,%%g0,%0" : "=r"(bb) :  ); 
  35.   switch (b_used%16) {
  36.       case 15: ADD_LOOP_BODY;
  37.       case 14: ADD_LOOP_BODY;
  38.       case 13: ADD_LOOP_BODY;
  39.       case 12: ADD_LOOP_BODY;
  40.       case 11: ADD_LOOP_BODY;
  41.       case 10: ADD_LOOP_BODY;
  42.       case  9: ADD_LOOP_BODY;
  43.       case  8: ADD_LOOP_BODY;
  44.       case  7: ADD_LOOP_BODY;
  45.       case  6: ADD_LOOP_BODY;
  46.       case  5: ADD_LOOP_BODY;
  47.       case  4: ADD_LOOP_BODY;
  48.       case  3: ADD_LOOP_BODY;
  49.       case  2: ADD_LOOP_BODY;
  50.       case  1: ADD_LOOP_BODY;
  51.       case  0: asm volatile ("addx  %%g0,%%g0,%0" : "=r"(bb) :  ); 
  52.      }
  55.   while (n--)
  56.   { asm volatile ("addcc 0xffffffff,%1,%0" : "=r"(bb) : "0"(bb) );
  57.     ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; 
  58.     ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; 
  59.     ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; 
  60.     ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY; ADD_LOOP_BODY;
  61.     asm volatile ("addx  %%g0,%%g0,%0" : "=r"(bb) : );
  62.    }
  64.   if (sum != a) 
  65.   { /* Copy_Vector(sum,a,a_used-b_used); */
  66.     for(b =sum; b < s_last; b++) *b = *a++;
  67.     *s_last = 0;
  68.    }
  70.   if (bb) /* propagate carry */
  71.   { while (++*sum == 0) sum++;
  72.     if (sum == s_last) a_used++;
  73.    }
  75.   return a_used;
  76. }