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

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

数值算法/人工智能

开发平台：

MultiPlatform

_seg_set.c：源码内容
  1. /*******************************************************************************
  2. +
  3. +  LEDA-R  3.2.3
  4. +
  5. +  _seg_set.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. #include <LEDA/segment_set.h>
  15. segment  SegmentSet::key(seg_item it)
  16. { segment s(x0(it),y(it),x1(it),y(it)); 
  17.   return s.rotate(point(0,0),alpha);
  18.  }
  20. seg_item SegmentSet::insert(segment s, GenPtr p)
  21. { segment t = s.rotate(point(0,0),-alpha);
  22.   double x0 = t.xcoord1();
  23.   double x1 = t.xcoord2();
  24.   double  y = t.ycoord1();
  25.   return segment_tree<double,double,GenPtr>::insert(x0,x1,y,p);
  26. }
  28. seg_item SegmentSet::lookup(segment s)
  29. { segment t = s.rotate(point(0,0),-alpha);
  30.   double x0 = t.xcoord1();
  31.   double x1 = t.xcoord2();
  32.   double  y = t.ycoord1();
  33.   return segment_tree<double,double,GenPtr>::lookup(x0,x1,y);
  34. }
  36. void     SegmentSet::del(segment s)
  37. { segment t = s.rotate(point(0,0),-alpha);
  38.   double x0 = t.xcoord1();
  39.   double x1 = t.xcoord2();
  40.   double  y = t.ycoord1();
  41.   (void)segment_tree<double,double,GenPtr>::del(x0,x1,y); 
  42. }
  44. list<seg_item>  SegmentSet::intersection(segment q)
  45. { list<seg_tree_item> L;
  46.   segment t = q.rotate(point(0,0),-alpha);
  47.   double  x = t.xcoord1();
  48.   double y0 = t.ycoord1();
  49.   double y1 = t.ycoord2();
  50.   L = segment_tree<double,double,GenPtr>::query(x,y0,y1); 
  51.   return *((list<seg_item>*)&L);
  52. }
  54. list<seg_item>  SegmentSet::intersection(line q)
  55. { list<seg_tree_item> L;
  56.   line t = q.rotate(-alpha);
  57.   double  x = t.x_proj(0);
  58.   L = segment_tree<double,double,GenPtr>::y_infinity_query(x); 
  59.   return *((list<seg_item>*)&L);
  60. }