home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / IRIT / IRITS.ZIP / CONTOURS.IRT < prev    next >
Encoding:
Text File  |  1990-05-05  |  1.1 KB  |  52 lines
  1. #
  2. # Intersection of cone and a cylinder:
  3. # Try this one with resolution equal 20 - slower, but much nicer!
  4. #
  5. #            Created by Gershon Elber,    Jan. 89
  6. #
  7.  
  8. save_mat = view_mat;
  9. view_mat = view_mat * scale(vector(0.2, 0.2, 0.2));
  10. save_res = resolution;
  11. resolution = 8;
  12.  
  13. cone1 = cone(vector(0,0,-1),vector(0,0,4),2);
  14. cylin1= cylin(vector(0,3,1),vector(0,-6,0),0.7);
  15.  
  16. a = cone1+cylin1;
  17. free(cylin1);
  18. free(cone1);
  19. beep(500,300);
  20. interact(list(a), false);
  21.  
  22. intercrv = true;
  23. #
  24. # Currently variables can not be introduced in a loop (except the iteration
  25. # vairable), so we prepare all variables in advance.
  26. #
  27. p = plane(vector(0.0, 0.0, 1.0), vector( 0.0, 0.0, -0.9), 6.0);
  28. c = cntrs = a * p;
  29. view(list(cntrs), false);
  30. for ( (i = -0.7), 0.2, 2.5,
  31.         (
  32.           (p = plane(vector(0.0, 0.0, 1.0), vector( 0.0, 0.0, i), 6.0)):
  33.       (c = a * p):
  34.       view(list(c), false):
  35.       (cntrs = cntrs + c)
  36.         )
  37.     );
  38. intercrv = false;
  39.  
  40. gdump("contours", cntrs);
  41. mdump("contours", view_mat);
  42. free(a);
  43. free(p);
  44. free(c);
  45. free(cntrs);
  46.  
  47. resolution = save_res;
  48. view_mat = save_mat;
  49.  
  50.  
  51.  
  52.