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

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

数值算法/人工智能

开发平台：

MultiPlatform

_line.c：源码内容
  1. /*******************************************************************************
  2. +
  3. +  LEDA-R  3.2.3
  4. +
  5. +  _line.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. #include <LEDA/line.h>
  15. #include <math.h>
  17. //------------------------------------------------------------------------------
  18. // lines 
  19. //------------------------------------------------------------------------------
  22. line::line()
  23. { PTR = new line_rep; }
  25. line::line(const segment& s) 
  26. { PTR = new line_rep(s); }
  28. line::line(const point& x, const point& y)    
  29. { PTR = new line_rep(segment(x,y)); }
  31. line::line(const point& p, const vector& v) 
  32. { PTR = new line_rep(segment(p,v)); }
  34. line::line(const point& p, double alpha) 
  35. { PTR = new line_rep(segment(p,alpha,1)); }
  36.   
  39. bool line::contains(const point& p) const
  40. { return orientation(ptr()->seg,p) == 0; }
  42. bool line::contains(const segment& s) const
  43. { return contains(s.start()) && contains(s.end()); }
  46. ostream& operator<<(ostream& out, const line& l) 
  47. { return out << l.seg(); }
  49. istream& operator>>(istream& in, line& l)  
  50. { segment s; 
  51.   in >> s; 
  52.   l = line(s); 
  53.   return in; 
  54.  }
  57. bool line::intersection(const line& s, point& inter) const
  58. { return ptr()->seg.intersection_of_lines(s.ptr()->seg,inter) ; }
  62. bool line::intersection(const segment& s, point& inter) const
  63. { 
  64.   double orient1 = orientation(ptr()->seg,s.start());
  65.   double orient2 = orientation(ptr()->seg,s.end());
  67.   if ( orient1*orient2 <= 0) 
  68.      { ptr()->seg.intersection_of_lines(s,inter);
  69.        return true;
  70.       }
  71.   else
  72.      return false;
  73. }
  76. segment line::perpendicular(const point& q) const
  77. { segment s0 = ptr()->seg;
  78.   segment s1 = s0.rotate90();
  79.   point r;
  80.   s0.intersection_of_lines(s1,r);
  81.   return segment(q,r);
  82.  }
  85. line p_bisector(const point& p, const point& q)
  86. { double dx = q.xcoord() - p.xcoord();
  87.   double dy = q.ycoord() - p.ycoord();
  88.   double x1 = (p.xcoord() + q.xcoord())/2;
  89.   double y1 = (p.ycoord() + q.ycoord())/2;
  90.   return line(point(x1,y1),point(x1+dy,y1-dx));
  91. }