home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / ledar34.zip / leda-r-3_4_tar / LEDA-3.4 / prog / window / hilbert.c < prev    next >
Encoding:
C/C++ Source or Header  |  1996-09-03  |  1.3 KB  |  107 lines
  1. #include <LEDA/window.h>
  2.  
  3. window W;
  4.  
  5. double x, y, dx, dy;
  6.  
  7. void A(int);
  8. void B(int);
  9. void C(int);
  10. void D(int);
  11.  
  12.  
  13. void plot(double new_x, double new_y)
  14. { W.draw_segment(x,y,new_x,new_y);
  15.   x = new_x;
  16.   y = new_y;
  17.  }
  18.  
  19.  
  20. void A(int i)
  21. { 
  22.   if (i > 0)
  23.   { D(i-1); plot(x-dx,y);
  24.     A(i-1); plot(x,y-dy);
  25.     A(i-1); plot(x+dx,y); 
  26.     B(i-1);
  27.    }
  28.  }
  29.  
  30. void B(int i)
  31. { 
  32.   if (i > 0)
  33.   { C(i-1); plot(x,y+dy);
  34.     B(i-1); plot(x+dx,y);
  35.     B(i-1); plot(x,y-dy); 
  36.     A(i-1);
  37.    }
  38.  }
  39.  
  40.  
  41.  
  42. void C(int i)
  43. { 
  44.   if (i > 0)
  45.   { B(i-1); plot(x+dx,y);
  46.     C(i-1); plot(x,y+dy);
  47.     C(i-1); plot(x-dx,y); 
  48.     D(i-1);
  49.    }
  50.  }
  51.  
  52. void D(int i)
  53. { 
  54.   if (i > 0)
  55.   { A(i-1); plot(x,y-dy);
  56.     D(i-1); plot(x-dx,y);
  57.     D(i-1); plot(x,y+dy); 
  58.     C(i-1);
  59.    }
  60.  }
  61.  
  62.  
  63. int n = 5;
  64.  
  65. void hilbert()
  66. {
  67.    double lx = W.xmax() - W.xmin();
  68.    double ly = W.ymax() - W.ymin();
  69.  
  70.    double x0 = W.xmin() + 0.98*lx;
  71.    double y0 = W.ymin() + 0.98*ly;
  72.  
  73.    dx = 0.96 * lx/(1 << n);
  74.    dy = 0.96 * ly/(1 << n);
  75.  
  76.    x = x0;
  77.    y = y0;
  78.  
  79.    A(n);
  80.  
  81.    W.draw_segment(x0,y0,x0+dx,y0);
  82.    W.draw_segment(x0+dx,y0,x0+dx,y);
  83.    W.draw_segment(x0+dx,y,x,y);
  84.  
  85.   }
  86.  
  87.  
  88. int main()
  89. {   
  90.  
  91.  W.open();
  92.  
  93.  panel P("hilbert curve");
  94.  P.int_item("n = ",n,1,10);
  95.  
  96.  W.set_redraw(hilbert);
  97.  
  98.  for(;;)
  99.  { P.open();
  100.    W.clear();
  101.    hilbert();
  102.    W.read_mouse();
  103.   }
  104.  
  105.   return 0;
  106. }
  107.