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

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

数值算法/人工智能

开发平台：

MultiPlatform

leda.c：源码内容
  1. #include <LEDA/graph.h>
  2. #include <LEDA/stream.h>
  3. #include <LEDA/array.h>
  4. #include <LEDA/window.h>
  5. #include <math.h>
  8. typedef GRAPH<vector,int> polyhedron;
  10. void rotate(float alpha1,float alpha2, vector& p)
  11. { 
  12.   // rotate 3d-point p about the origin
  13.   // by alpha2 in yz-plane and alpha1 in xy-plane
  15.     double R  = hypot(p[1],p[2]);
  16.     double phi = asin(p[1]/R);
  17.   
  18.     if (p[2] < 0) phi = LEDA_PI - phi;
  19.   
  20.     p[1]  = ((R != 0) ? R*sin(phi+alpha2) : 0);
  21.     p[2]  = ((R != 0) ? R*cos(phi+alpha2) : 0);
  23.     R   = hypot(p[0],p[2]);
  24.     phi = asin(p[0]/R);
  26.     if (p[2] < 0) phi = LEDA_PI - phi;
  27.   
  28.     p[0]  = ((R != 0) ? R*sin(phi+alpha1) : 0);
  29.     p[2]  = ((R != 0) ? R*cos(phi+alpha1) : 0);
  30. }
  33. inline point project(vector p)   // project p into xy-plane
  34. { return point(p[0],p[1]); }
  37. void draw_poly(window& W, polyhedron& poly, vector trans)
  38. { edge e;
  39.   forall_edges(e,poly) 
  40.   { point a = project(poly[source(e)] + trans );
  41.     point b = project(poly[target(e)] + trans );
  42.     W.draw_segment(a,b,blue);
  43.    }
  44.  }
  47. void make_side(polyhedron& poly, node* L, float z)
  48. {
  49.   L[0]= poly.new_node(vector(-70,-20,z));
  50.   L[1]= poly.new_node(vector(-40,-20,z));
  51.   L[2]= poly.new_node(vector(-40,-10,z));
  52.   L[3]= poly.new_node(vector(-60,-10,z));
  53.   L[4]= poly.new_node(vector(-60,+20,z));
  54.   L[5]= poly.new_node(vector(-70,+20,z));
  56.   poly.new_edge(L[0],L[5]);
  57.   int i;
  58.   for(i = 1; i<=5; i++) poly.new_edge(L[i], L[i-1]);
  60.   L[6] = poly.new_node(vector(-30,-20,z));
  61.   L[7] = poly.new_node(vector( -5,-20,z));
  62.   L[8] = poly.new_node(vector( -5,-10,z));
  63.   L[9] = poly.new_node(vector(-20,-10,z));
  64.   L[10]= poly.new_node(vector(-20, -5,z));
  65.   L[11]= poly.new_node(vector( -5, -5,z));
  66.   L[12]= poly.new_node(vector( -5, +5,z));
  67.   L[13]= poly.new_node(vector(-20, +5,z));
  68.   L[14]= poly.new_node(vector(-20,+10,z));
  69.   L[15]= poly.new_node(vector( -5,+10,z));
  70.   L[16]= poly.new_node(vector( -5,+20,z));
  71.   L[17]= poly.new_node(vector(-30,+20,z));
  73.   poly.new_edge(L[6],L[17]);
  74.   for(i = 7; i<=17; i++) poly.new_edge(L[i], L[i-1]);
  76.   L[18]= poly.new_node(vector( 5,-20,z));
  77.   L[19]= poly.new_node(vector(20,-20,z));
  78.   L[20]= poly.new_node(vector(35,-10,z));
  79.   L[21]= poly.new_node(vector(35,+10,z));
  80.   L[22]= poly.new_node(vector(20,+20,z));
  81.   L[23]= poly.new_node(vector( 5,+20,z));
  83.   poly.new_edge(L[18],L[23]);
  84.   for(i = 19; i<=23; i++) poly.new_edge(L[i], L[i-1]);
  86.   L[24]= poly.new_node(vector(15,-10,z));
  87.   L[25]= poly.new_node(vector(20,-10,z));
  88.   L[26]= poly.new_node(vector(25, -5,z));
  89.   L[27]= poly.new_node(vector(25, +5,z));
  90.   L[28]= poly.new_node(vector(20,+10,z));
  91.   L[29]= poly.new_node(vector(15,+10,z));
  93.   poly.new_edge(L[24],L[29]);
  94.   for(i = 25; i<=29; i++) poly.new_edge(L[i], L[i-1]);
  96.   L[30]= poly.new_node(vector(40,-20,z));
  97.   L[31]= poly.new_node(vector(50,-20,z));
  98.   L[32]= poly.new_node(vector(55, -5,z));
  99.   L[33]= poly.new_node(vector(65, -5,z));
  100.   L[34]= poly.new_node(vector(70,-20,z));
  101.   L[35]= poly.new_node(vector(80,-20,z));
  102.   L[36]= poly.new_node(vector(65, 20,z));
  103.   L[37]= poly.new_node(vector(55, 20,z));
  105.   poly.new_edge(L[30],L[37]);
  106.   for(i = 31; i<=37; i++) poly.new_edge(L[i], L[i-1]);
  108.   L[38]= poly.new_node(vector(55, 0,z));
  109.   L[39]= poly.new_node(vector(65, 0,z));
  110.   L[40]= poly.new_node(vector(60,15,z));
  111.   poly.new_edge(L[38],L[40]);
  112.   poly.new_edge(L[39],L[38]);
  113.   poly.new_edge(L[40],L[39]);
  115. }
  119. void make_logo(polyhedron& poly)
  120. { 
  121.   node L[41];
  122.   node R[41];
  124.   make_side(poly,L,-5);
  125.   make_side(poly,R,+5);
  127.   for(int i = 0; i<41; i++) poly.new_edge(L[i],R[i]);
  129. }
  133.   
  135. main()
  136. { 
  137.   int w,h;
  139.   char* map = Read_Leda_Bitmap("../../incl/LEDA/bitmaps/leda.lbm",w,h);
  141.   window W(w,h,window::center,window::center);
  142.   W.set_frame_label(
  143.               "The LEDA Platform for Combinatorial and Geometric Computing");
  144.   W.init(-100,100,-100);
  145.   W.set_show_coordinates(false);
  146.   W.insert_bitmap(w,h,map);
  147.   W.set_node_width(4);
  148.   W.set_mode(xor_mode);
  150.   node v;
  152.   int speed = 30;
  154.   // define a polyhedron in 3d space
  156.   polyhedron poly;
  158.   make_logo(poly);
  160.   //W.draw_point(0,0);
  162.   forall_nodes(v,poly) rotate(0.3,0.5,poly[v]);
  164.   polyhedron poly0 = poly;
  166.   vector dir(2);   // direction and speed of rotation
  168.   vector trans0(0,0,0);
  169.   vector trans(0,0,0);
  171.   draw_poly(W,poly,trans);
  172.   
  173.   for(;;)
  174.   { dir[0] += rand_int(-10000,10000)/1000000.0;
  175.     dir[1] += rand_int(-10000,10000)/1000000.0;
  176.     dir = dir.norm()*(speed/500.0);
  178.     //trans[0] += rand_int(-1000,1000)/5000.0;
  179.     //trans[1] += rand_int(-1000,1000)/5000.0;
  180.     //trans[1] += rand_int(-1000,1000)/5000.0;
  182.     if (W.get_button() && W.confirm("quit leda demo ?")) break;
  184.     forall_nodes(v,poly) rotate(dir[0],dir[1],poly[v]);
  185.  
  186.     draw_poly(W,poly,trans);   // draw new position
  187.     draw_poly(W,poly0,trans0);  // erase old position  (xor_mode !!)
  189.     //trans0 = trans;
  191.     // poly0 = poly;
  192.     node w = poly.first_node();
  193.     forall_nodes(v,poly0) 
  194.     { poly0[v] = poly[w];
  195.       w = poly.succ_node(w);
  196.      }
  197.   }
  199.  return 0;
  200. }