home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / diverses / leda / src / plane / _line.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-11-15  |  2.8 KB  |  122 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  2.1.1                                                 11-15-1991
  4. +
  5. +
  6. +  _line.c
  7. +
  8. +
  9. +  Copyright (c) 1991  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #include <LEDA/plane.h>
  17. #include <math.h>
  18.  
  19. static const double eps = 1e-10;
  20.  
  21. //------------------------------------------------------------------------------
  22. // lines 
  23. //------------------------------------------------------------------------------
  24.  
  25.  
  26. line_rep::line_rep(segment s)  
  27. { seg = s; 
  28.   count = 1; 
  29.  }
  30.  
  31. line_rep::line_rep() { count = 1; }
  32.  
  33. line::line()                    { ptr = new line_rep; }
  34. line::line(int)                 { ptr = new line_rep; }
  35. line::line(ent x)               { ptr= (line_rep*)x; ptr->count++; }
  36. line::line(line& l)             { ptr = l.ptr; ptr->count++; }
  37. line::line(segment s)           { ptr = new line_rep(s); }
  38. line::line(point x, point y)    { ptr = new line_rep(segment(x,y)); }
  39. line::line(point p, double alpha) { ptr = new line_rep(segment(p,alpha,1)); }
  40.   
  41. void line::clear()              { if (--(ptr->count)==0) delete ptr; }
  42.  
  43. double line::y_proj(double x)         
  44. { return slope() * x + y_abs(); }
  45.  
  46. double line::x_proj(double y)         
  47. { return (vertical()) ? ptr->seg.xcoord1() : (y-y_abs())/slope(); }
  48.  
  49. bool line::contains(point p)   
  50. { if (vertical())
  51.       return p.xcoord() == ptr->seg.xcoord1(); 
  52.   else
  53.       return p.ycoord() == y_proj(p.xcoord());
  54.  }
  55.  
  56. bool line::contains(segment s) 
  57. { point p = s.start(); 
  58.   point q = s.end();
  59.   if (contains(p))
  60.     if (contains(q)) return true; 
  61.   return false;
  62.  }
  63.  
  64.  
  65. ostream& operator<<(ostream& out, const line& l) 
  66. { return out << l.ptr->seg; }
  67.  
  68. istream& operator>>(istream& in, line& l)  
  69. { segment s; 
  70.   in >> s; 
  71.   l = line(s); 
  72.   return in; 
  73.  }
  74.  
  75. line& line::operator=(const line& l)
  76. { l.ptr->count++;
  77.   if (--ptr->count == 0)  delete ptr;
  78.   ptr = l.ptr;
  79.   return *this;
  80. }
  81.  
  82. bool line::intersection(line L, point& inter)
  83. {
  84.   point p = L.ptr->seg.ptr->start;
  85.   segment s = perpendicular(p);
  86.  
  87.   double l = s.length();
  88.  
  89.   if (l==0)
  90.   { inter = p;
  91.     return true;
  92.    }
  93.  
  94.   double fi = L.angle(s);
  95.   double cosfi = cos(fi);
  96.  
  97.   if (cosfi==0) return false;
  98.  
  99.   inter = p.translate(L.angle(),l/cosfi);
  100.   return true;
  101.  
  102.  }
  103.  
  104. bool line::intersection(segment s, point& inter)
  105. { if (intersection(line(s),inter))
  106.   { double d1  = inter.distance(s.ptr->start);
  107.     double d2  = inter.distance(s.ptr->end);
  108.     double l   = s.length();
  109.     if ( d1<=l && d2 <=l) return true;
  110.    }
  111.    return false;
  112. }
  113.  
  114.  
  115. segment line::perpendicular(point q)
  116. { segment s = ptr->seg;
  117.   point r = q.translate(s.angle()+M_PI_2,-distance(q));
  118.   return segment(q,r);
  119.  }
  120.  
  121.  
  122.