Category

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

Mathimatics-Numerical algorithms

Development Platform:

MultiPlatform

interval_set.c：Code Content
  1. #include <LEDA/plane.h>
  2. #include <LEDA/interval_set.h>
  3. #include <LEDA/window.h>
  8. main()
  9. { 
  10.   interval_set<segment>  I;
  12.   window W;
  14.   W.init(-1000,1000,-1000);
  15.   W.set_line_width(1);
  17.   panel P("interval_set");
  19.   P.text_item("Use the left mouse button to insert a set of segments, ");
  20.   P.text_item("terminate the input by clicking the right button. Now, ");
  21.   P.text_item("you can define query slabs by giving two vertical lines");
  22.   P.text_item("with the left mouse button. For each slab, the program ");
  23.   P.text_item("computes and shows all segments intersecting it. Leave ");
  24.   P.text_item("the program by clicking the right mouse button.        ");
  26.   P.button("continue");
  28.   P.open();
  31.   segment s;
  33.   while (W >> s)
  34.   { W << s;
  35.     if (s.xcoord1() <= s.xcoord2())
  36.       I.insert(s.xcoord1(),s.xcoord2(),s);
  37.     else
  38.       I.insert(s.xcoord2(),s.xcoord1(),s);
  39.    }
  41.   W.set_mode(xor_mode);
  43.   list<is_item> L;
  44.   is_item it;
  45.   double    x0 = W.xmin();
  46.   double    x1 = W.xmin();
  48.   point   p;
  50.   while (W >> p)
  51.   { 
  52.     W.draw_vline(x0);
  53.     W.draw_vline(x1);
  55.     W.set_line_width(4);
  56.     forall(it,L) W << I.inf(it);
  57.     W.set_line_width(1);
  58.     forall(it,L) W << I.inf(it);
  60.     x0 = p.xcoord();
  61.     W.draw_vline(x0);
  63.     W >> p;
  64.     x1 = p.xcoord();
  65.     W.draw_vline(x1);
  67.     L = I.intersection(x0,x1);
  69.     forall(it,L) W << I.inf(it);
  70.     W.set_line_width(4);
  71.     forall(it,L) W << I.inf(it);
  72.     W.set_line_width(1);
  74.     forall(it,L) cout << I.inf(it) << "n";
  75.     newline;
  77.   }
  79.   return 0;
  80. }