home *** CD-ROM | disk | FTP | other *** search
/ Mega A/V / mega_av.zip / mega_av / GRAPHUTL / POVSCN.ZIP / INCLUDE.ZIP / SHAPESQ.INC < prev    next >
Text File  |  1992-07-03  |  13KB  |  346 lines
  1. // Persistence Of Vision Raytracer 
  2. // Standard include file.
  3.  
  4. // Quartic shapes include file
  5. //
  6. // Several cubic and quartic shape definitions 
  7. // by Alexander Enzmann 
  8.  
  9. /* In the following descriptions, multiplication of two terms is
  10.   shown as the two terms next to each other (i.e. x y, rather than
  11.   x*y.  The expression c(n, m) is the binomial coefficient, n!/m!(n-m)!. */
  12.  
  13. /* Bicorn
  14.   This curve looks like the top part of a paraboloid, bounded
  15.   from below by another paraboloid.  The basic equation is:
  16.      y^2 - (x^2 + z^2) y^2 - (x^2 + z^2 + 2 y - 1)^2 = 0.  */
  17. #declare Bicorn =
  18. quartic {
  19.    < 1.0   0.0   0.0   0.0  1.0   0.0   4.0   2.0   0.0 -2.0
  20.    0.0   0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0  0.0
  21.    0.0   0.0   0.0   1.0  0.0   3.0   0.0   4.0   0.0 -4.0
  22.    1.0   0.0  -2.0   0.0  1.0 >
  23. }
  24.  
  25. /* Crossed Trough
  26.   This is a surface with four pieces that sweep up from the x-z plane.
  27.   The equation is: y = x^2 z^2.  */
  28. #declare Crossed_Trough =
  29. quartic {
  30.    < 0.0   0.0   0.0   0.0  0.0   0.0   0.0   4.0   0.0  0.0
  31.    0.0   0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0  0.0
  32.    0.0   0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0 -1.0
  33.    0.0   0.0   0.0   0.0  0.0 >
  34. }
  35.  
  36. /* a drop coming out of water? This is a curve formed by using the equation
  37.   y = 1/2 x^2 (x + 1) as the radius of a cylinder having the x-axis as
  38.   its central axis. The final form of the equation is:
  39.      y^2 + z^2 = 0.5 (x^3 + x^2) */
  40. #declare Cubic_Cylinder =
  41. quartic {
  42.    < 0.0   0.0   0.0   -0.5  0.0   0.0   0.0   0.0   0.0 -0.5
  43.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  44.    0.0   0.0   0.0    0.0  0.0   1.0   0.0   0.0   0.0  0.0
  45.    0.0   0.0   1.0    0.0  0.0 >
  46. }
  47.  
  48. /* a cubic saddle. The equation is: z = x^3 - y^3. */
  49. #declare Cubic_Saddle_1 =
  50. quartic {
  51.    <  0.0   0.0   0.0    1.0   0.0   0.0   0.0   0.0   0.0  0.0
  52.    0.0   0.0   0.0    0.0   0.0   0.0   0.0   0.0   0.0  0.0
  53.    0.0   0.0  -1.0    0.0   0.0   0.0   0.0   0.0   0.0  0.0
  54.    0.0   0.0   0.0   -1.0   0.0 >
  55. }
  56.  
  57. /* Variant of a devil's curve in 3-space.  This figure has a top and
  58.   bottom part that are very similar to a hyperboloid of one sheet,
  59.   however the central region is pinched in the middle leaving two
  60.   teardrop shaped holes. The equation is:
  61.      x^4 + 2 x^2 z^2 - 0.36 x^2 - y^4 + 0.25 y^2 + z^4 = 0.  */
  62. #declare Devils_Curve =
  63. quartic {
  64.    <-1.0   0.0   0.0    0.0  0.0   0.0   0.0  -2.0   0.0  0.36
  65.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  66.    1.0   0.0   0.0    0.0  0.0  -0.25  0.0   0.0   0.0  0.0
  67.    -1.0   0.0   0.0    0.0  0.0 >
  68. }
  69.  
  70. /* Folium
  71.   This is a folium rotated about the x-axis.  The formula is:
  72.      2 x^2 - 3 x y^2 - 3 x z^2 + y^2 + z^2 = 0. */
  73. #declare Folium =
  74. quartic {
  75.    < 0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  2.0
  76.    0.0   0.0  -3.0    0.0  0.0   0.0   0.0  -3.0   0.0  0.0
  77.    0.0   0.0   0.0    0.0  0.0   1.0   0.0   0.0   0.0  0.0
  78.    0.0   0.0   1.0    0.0  0.0 >
  79. }
  80.  
  81. /* Glob - sort of like basic teardrop shape. The equation is:
  82.    y^2 + z^2 = 0.5 x^5 + 0.5 x^4. */
  83. #declare Glob_5 =
  84. poly { 5
  85.    <-0.5   0.0   0.0  -0.5   0.0   0.0   0.0   0.0   0.0  0.0
  86.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  87.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  88.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  89.    0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0  0.0
  90.    0.0   0.0   0.0   1.0   0.0   0.0 >
  91. }
  92.  
  93. /* Variant of a lemniscate - the two lobes are much more teardrop-like. */
  94. #declare Twin_Glob =
  95. poly { 6
  96.    < 4.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0 -4.0
  97.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  98.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  99.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  100.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  101.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  102.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  103.    1.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  104.    0.0   1.0   0.0   0.0 >
  105. }
  106.  
  107. /*  Approximation to the helix z = arctan(y/x).
  108.  
  109.    The helix can be approximated with an algebraic equation (kept to the
  110.    range of a quartic) with the following steps:
  111.  
  112.       tan(z) = y/x   =>  sin(z)/cos(z) = y/x   =>
  113.  
  114.    (1) x sin(z) - y cos(z) = 0
  115.  
  116.    Using the taylor expansions for sin, cos about z = 0,
  117.  
  118.       sin(z) = z - z^3/3! + z^5/5! - ...
  119.       cos(z) = 1 - z^2/2! + z^6/6! - ...
  120.  
  121.    Throwing out the high order terms, the expression (1) can be written as:
  122.  
  123.       x (z - z^3/6) - y (1 + z^2/2) = 0, or
  124.  
  125.   (2) -1/6 x z^3 + x z + 1/2 y z^2 - y = 0
  126.  
  127.   This helix (2) turns 90 degrees in the range 0 <= z <= sqrt(2)/2.  By using
  128.   scale <2 2 2>, the helix defined below turns 90 degrees in the range
  129.   0 <= z <= sqrt(2) = 1.4042.
  130. */
  131. #declare Helix =
  132. object {
  133.    intersection {
  134.       quadric { Cylinder_Z
  135.          scale <2.0 2.0 2.0>
  136.          texture { color Clear }
  137.       }
  138.       plane { <0.0 0.0 1.0> 1.4142 texture { color Clear } }
  139.       plane { <0.0 0.0 -1.0> 0.0 texture { color Clear } }
  140.       quartic {
  141.          <  0.0   0.0   0.0    0.0  0.0   0.0   0.0    0.0   0.0  0.0
  142.          0.0   0.0   0.0    0.0  0.0   0.0  -0.1666 0.0   1.0  0.0
  143.          0.0   0.0   0.0    0.0  0.0   0.0   0.0    0.5   0.0 -1.0
  144.          0.0   0.0   0.0    0.0  0.0 >
  145.       }
  146.    }
  147. }
  148.  
  149. /* This is an alternate Helix, using clipped_by instead of csg intersection. */
  150. #declare Helix_1 =
  151. object {
  152.    quartic {
  153.       < 0.0   0.0   0.0    0.0  0.0   0.0   0.0    0.0   0.0  0.0
  154.       0.0   0.0   0.0    0.0  0.0   0.0  -0.1666 0.0   1.0  0.0
  155.       0.0   0.0   0.0    0.0  0.0   0.0   0.0    0.5   0.0 -1.0
  156.       0.0   0.0   0.0    0.0  0.0 >
  157.    }
  158.    clipped_by {
  159.       quadric { Cylinder_Z scale <2.0 2.0 2.0> }
  160.       plane { <0.0 0.0 1.0> 1.4142 }
  161.       plane { <0.0 0.0 -1.0> 0.0 }
  162.    }
  163.    bounded_by {
  164.       intersection {
  165.          quadric { Cylinder_Z scale <2.0 2.0 2.0> }
  166.          plane { <0.0 0.0 1.0> 1.4142 }
  167.          plane { <0.0 0.0 -1.0> 0.0 }
  168.       }
  169.    }
  170. }
  171.  
  172. /* Hyperbolic Torus having major radius sqrt(40), minor radius sqrt(12).
  173.   This figure is generated by sweeping a circle along the arms of a
  174.   hyperbola.  The equation is:
  175.  
  176.      x^4 + 2 x^2 y^2 - 2 x^2 z^2 - 104 x^2 + y^4 - 2 y^2 z^2 +
  177.      56 y^2 + z^4 + 104 z^2 + 784 = 0.
  178.  
  179.   See the description for the torus below. */
  180. #declare Hyperbolic_Torus_40_12 =
  181. quartic {
  182.    < 1.0   0.0   0.0    0.0     2.0   0.0   0.0  -2.0   0.0 -104.0
  183.    0.0   0.0   0.0    0.0     0.0   0.0   0.0   0.0   0.0    0.0
  184.    1.0   0.0   0.0   -2.0     0.0  56.0   0.0   0.0   0.0    0.0
  185.    1.0   0.0  104.0   0.0   784.0 >
  186. }
  187.  
  188. /* Lemniscate of Gerono
  189.   This figure looks like two teardrops with their pointed ends connected.
  190.   It is formed by rotating the Lemniscate of Gerono about the x-axis.
  191.   The formula is:
  192.      x^4 - x^2 + y^2 + z^2 = 0. */
  193. #declare Lemniscate =
  194. quartic {
  195.    < 1.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0 -1.0
  196.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0  0.0
  197.    0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0  0.0
  198.    0.0   0.0   1.0   0.0   0.0 >
  199. }
  200.  
  201. /* This is a figure with a bumpy sheet on one side and something that
  202.   looks like a paraboloid (but with an internal bubble).  The formula
  203.   is:
  204.      (x^2 + y^2 + a c x)^2 - (x^2 + y^2)(c - a x)^2.
  205.  
  206.    -99*x^4+40*x^3-98*x^2*y^2-98*x^2*z^2+99*x^2+40*x*y^2+40*x*z^2+y^4+2*y^2*z^2
  207.    -y^2+z^4-z^2
  208.  
  209. */
  210. #declare Quartic_Loop_1 =
  211. quartic {
  212.    <99.0   0.0   0.0 -40.0  98.0   0.0   0.0  98.0   0.0 -99.0
  213.    0.0   0.0 -40.0   0.0   0.0   0.0   0.0 -40.0   0.0   0.0
  214.    -1.0   0.0   0.0  -2.0   0.0   1.0   0.0   0.0   0.0   0.0
  215.    -1.0   0.0   1.0   0.0   0.0 >
  216. }
  217.  
  218. /* Monkey Saddle
  219.   This surface has three parts that sweep up and three down.  This gives
  220.   a saddle that has a place for two legs and a tail... The equation is:
  221.  
  222.      z = c (x^3 - 3 x y^2).
  223.  
  224.   The value c gives a vertical scale to the surface - the smaller the
  225.   value of c, the flatter the surface will be (near the origin). */
  226. #declare Monkey_Saddle =
  227. quartic {
  228.    < 0.0   0.0   0.0   1.0  0.0   0.0   0.0   0.0   0.0  0.0
  229.    0.0   0.0  -3.0   0.0  0.0   0.0   0.0   0.0   0.0  0.0
  230.    0.0   0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0  0.0
  231.    0.0   0.0   0.0  -1.0  0.0 >
  232. }
  233.  
  234. /* Parabolic Torus having major radius sqrt(40), minor radius sqrt(12).
  235.   This figure is generated by sweeping a circle along the arms of a
  236.   parabola.  The equation is:
  237.  
  238.      x^4 + 2 x^2 y^2 - 2 x^2 z - 104 x^2 + y^4 - 2 y^2 z +
  239.      56 y^2 + z^2 + 104 z + 784 = 0.
  240.  
  241.   See the description for the torus below. */
  242. #declare Parabolic_Torus_40_12 =
  243. quartic {
  244.    < 1.0   0.0   0.0    0.0     2.0   0.0   0.0   0.0  -2.0 -104.0
  245.    0.0   0.0   0.0    0.0     0.0   0.0   0.0   0.0   0.0    0.0
  246.    1.0   0.0   0.0    0.0    -2.0  56.0   0.0   0.0   0.0    0.0
  247.    0.0   0.0   1.0  104.0   784.0 >
  248. }
  249.  
  250. /* Piriform
  251.   This figure looks like a hersheys kiss. It is formed by sweeping
  252.   a Piriform about the x-axis.  a basic form of the equation is:
  253.      (x^4 - x^3) + y^2 + z^2 = 0.
  254. */
  255. #declare Piriform =
  256. quartic {
  257.    < 4.0   0.0   0.0   -4.0  0.0   0.0   0.0   0.0   0.0  0.0
  258.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  259.    0.0   0.0   0.0    0.0  0.0   1.0   0.0   0.0   0.0  0.0
  260.    0.0   0.0   1.0    0.0  0.0 >
  261. }
  262.  
  263. /* n-Roll Mill
  264.   This curve in the plane looks like several hyperbolas with their
  265.   bumps arranged about the origin.  The general formula is:
  266.  
  267.      x^n - c(n,2) x^(n-2) y^2 + c(n,4) x^(n-4) y^4 - ... = a
  268.  
  269.   When rendering in 3-Space, the resulting figure looks like a
  270.   cylinder with indented sides.
  271. */
  272.  
  273. /* Quartic parabola - a 4th degree polynomial (has two bumps at the bottom)
  274.   that has been swept around the z axis. The equation is:
  275.      0.1 x^4 - x^2 - y^2 - z^2 + 0.9 = 0. */
  276. #declare Quartic_Paraboloid =
  277. quartic {
  278.    < 0.1   0.0   0.0  0.0   0.0   0.0   0.0   0.0   0.0  -1.0
  279.    0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0   0.0   0.0
  280.    0.0   0.0   0.0  0.0   0.0   0.0   0.0   0.0   0.0  -1.0
  281.    0.0   0.0  -1.0  0.0   0.9 >
  282. }
  283.  
  284. /* Quartic Cylinder - a Space Needle?  */
  285. #declare Quartic_Cylinder =
  286. quartic {
  287.    < 0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.01
  288.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0
  289.    0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.0   0.0
  290.    0.0   0.0   0.01  0.0  -0.01 >
  291. }
  292.  
  293. /* Steiners quartic surface */
  294. #declare Steiner_Surface =
  295. quartic {
  296.    < 0.0   0.0   0.0  0.0  1.0   0.0   0.0   1.0   0.0   0.0
  297.    0.0   0.0   0.0  0.0  1.0   0.0   0.0   0.0   0.0   0.0
  298.    0.0   0.0   0.0  1.0  0.0   0.0   0.0   0.0   0.0   0.0
  299.    0.0   0.0   0.0  0.0  0.0 >
  300. }
  301.  
  302. /* Torus having major radius sqrt(40), minor radius sqrt(12) */
  303. #declare Torus_40_12 =
  304. quartic {
  305.    < 1.0   0.0   0.0    0.0     2.0   0.0   0.0   2.0   0.0 -104.0
  306.    0.0   0.0   0.0    0.0     0.0   0.0   0.0   0.0   0.0    0.0
  307.    1.0   0.0   0.0    2.0     0.0  56.0   0.0   0.0   0.0    0.0
  308.    1.0   0.0 -104.0   0.0   784.0 >
  309. }
  310.  
  311. /* Witch of Agnesi */
  312. #declare Witch_Hat =
  313. quartic {
  314.    <  0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0
  315.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0
  316.    0.0   0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.04
  317.    0.0   0.0   0.0   0.0   0.04 >
  318. }
  319.  
  320. /* very rough approximation to the sin-wave surface z = sin(2 pi x y).
  321.   In order to get an approximation good to 7 decimals at a distance of
  322.   1 from the origin would require a polynomial of degree around 60.  This
  323.   would require around 200k coefficients. For best results, scale by
  324.   something like <1 1 0.2>. */
  325. #declare Sinsurf =
  326. poly { 6
  327.    < 0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  328.    -1116.226 0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  329.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  330.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  331.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0 18.8496
  332.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  333.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  334.    0.0   0.0   0.0    0.0  0.0   0.0   0.0   0.0   0.0  0.0
  335.    0.0   0.0  -1.0    0.0 >
  336. }
  337.  
  338. /* Empty quartic equation.  Ready to be filled with numbers...
  339.    quartic
  340.     {< 0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0
  341.        0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0
  342.        0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0
  343.        0.0   0.0   0.0   0.0   0.0 >
  344.    }
  345. */
  346.