开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 查看源码
_sgnu_sub.c
资源名称：leda.tar.gz　[点击查看]
上传用户：gzelex
上传日期：2007-01-07
资源大小：707k
文件大小：2k
源码类别：

数值算法/人工智能

开发平台：

MultiPlatform

_sgnu_sub.c：源码内容
  1. /*******************************************************************************
  2. +
  3. +  LEDA-R  3.2.3
  4. +
  5. +  _sgnu_sub.c
  6. +
  7. +  Copyright (c) 1995  by  Max-Planck-Institut fuer Informatik
  8. +  Im Stadtwald, 66123 Saarbruecken, Germany     
  9. +  All rights reserved.
  10. + 
  11. *******************************************************************************/
  14. #define MAX_WORD     0xFFFFFFFF
  16. typedef unsigned long  word;
  17. typedef unsigned int sz_t;
  21. extern sz_t School_Sub(word *a, sz_t a_used, word *b, sz_t b_used, word* diff)
  22. {
  23.   /* compute diff = a - b    (a > b) */
  25.   word* d_stop = diff + a_used;
  26.   int n = b_used / 16;
  27.   int r = b_used % 16; 
  29.   word  aa;
  30.   word  bb;
  32. #define SUB_LOOP_BODY {
  33. aa = *a++;
  34. bb = *b++;
  35. asm volatile ("subxcc  %2,%1,%0" : "=r"(aa) : "r"(bb), "0"(aa) );
  36. *diff++ = aa; }
  40.   asm volatile ("addcc  %%g0,%%g0,%0" : "=r"(bb) : );
  42.   switch (r) 
  43.    { case 15: SUB_LOOP_BODY;
  44.      case 14: SUB_LOOP_BODY;
  45.      case 13: SUB_LOOP_BODY;
  46.      case 12: SUB_LOOP_BODY;
  47.      case 11: SUB_LOOP_BODY;
  48.      case 10: SUB_LOOP_BODY;
  49.      case  9: SUB_LOOP_BODY;
  50.      case  8: SUB_LOOP_BODY;
  51.      case  7: SUB_LOOP_BODY;
  52.      case  6: SUB_LOOP_BODY;
  53.      case  5: SUB_LOOP_BODY;
  54.      case  4: SUB_LOOP_BODY;
  55.      case  3: SUB_LOOP_BODY;
  56.      case  2: SUB_LOOP_BODY;
  57.      case  1: SUB_LOOP_BODY;
  58.      case  0: asm volatile ("addx  %%g0,%%g0,%0" : "=r"(bb) : );
  59.     }
  62.   while (n--)
  63.   { asm volatile ("addcc 0xffffffff,%1,%0" : "=r"(bb) : "0"(bb) );
  64.     SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; 
  65.     SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; 
  66.     SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; 
  67.     SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY; SUB_LOOP_BODY;
  68.     asm volatile ("addx  %%g0,%%g0,%0" : "=r"(bb) : );
  69.    }
  71.   if (diff != a) /* Copy_Vector(diff,a,a_used-b_used); */
  72.       for(b=diff; b != d_stop; b++) *b = *a++;
  74.   if (bb)
  75.     while (--*diff == MAX_WORD) diff++;
  77.   while (*--d_stop == 0) a_used--;
  79.   return a_used;
  80. }