3d-raytrace.html 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491
  1. <!DOCTYPE html>
  2. <head>
  3. <!--
  4. Copyright (C) 2007 Apple Inc. All rights reserved.
  5. Redistribution and use in source and binary forms, with or without
  6. modification, are permitted provided that the following conditions
  7. are met:
  8. 1. Redistributions of source code must retain the above copyright
  9. notice, this list of conditions and the following disclaimer.
  10. 2. Redistributions in binary form must reproduce the above copyright
  11. notice, this list of conditions and the following disclaimer in the
  12. documentation and/or other materials provided with the distribution.
  13. THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY
  14. EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  15. IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
  16. PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL APPLE COMPUTER, INC. OR
  17. CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
  18. EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
  19. PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
  20. PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
  21. OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  22. (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
  23. OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  24. -->
  25. <title>SunSpider 3d-raytrace</title>
  26. </head>
  27. <body>
  28. <h3>3d-raytrace</h3>
  29. <div id="console">
  30. </div>
  31. <script>
  32. var _sunSpiderStartDate = new Date();
  33. /*
  34. * Copyright (C) 2007 Apple Inc. All rights reserved.
  35. *
  36. * Redistribution and use in source and binary forms, with or without
  37. * modification, are permitted provided that the following conditions
  38. * are met:
  39. * 1. Redistributions of source code must retain the above copyright
  40. * notice, this list of conditions and the following disclaimer.
  41. * 2. Redistributions in binary form must reproduce the above copyright
  42. * notice, this list of conditions and the following disclaimer in the
  43. * documentation and/or other materials provided with the distribution.
  44. *
  45. * THIS SOFTWARE IS PROVIDED BY APPLE COMPUTER, INC. ``AS IS'' AND ANY
  46. * EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  47. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
  48. * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL APPLE COMPUTER, INC. OR
  49. * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
  50. * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
  51. * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
  52. * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
  53. * OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  54. * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
  55. * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  56. */
  57. function createVector(x,y,z) {
  58. return new Array(x,y,z);
  59. }
  60. function sqrLengthVector(self) {
  61. return self[0] * self[0] + self[1] * self[1] + self[2] * self[2];
  62. }
  63. function lengthVector(self) {
  64. return Math.sqrt(self[0] * self[0] + self[1] * self[1] + self[2] * self[2]);
  65. }
  66. function addVector(self, v) {
  67. self[0] += v[0];
  68. self[1] += v[1];
  69. self[2] += v[2];
  70. return self;
  71. }
  72. function subVector(self, v) {
  73. self[0] -= v[0];
  74. self[1] -= v[1];
  75. self[2] -= v[2];
  76. return self;
  77. }
  78. function scaleVector(self, scale) {
  79. self[0] *= scale;
  80. self[1] *= scale;
  81. self[2] *= scale;
  82. return self;
  83. }
  84. function normaliseVector(self) {
  85. var len = Math.sqrt(self[0] * self[0] + self[1] * self[1] + self[2] * self[2]);
  86. self[0] /= len;
  87. self[1] /= len;
  88. self[2] /= len;
  89. return self;
  90. }
  91. function add(v1, v2) {
  92. return new Array(v1[0] + v2[0], v1[1] + v2[1], v1[2] + v2[2]);
  93. }
  94. function sub(v1, v2) {
  95. return new Array(v1[0] - v2[0], v1[1] - v2[1], v1[2] - v2[2]);
  96. }
  97. function scalev(v1, v2) {
  98. return new Array(v1[0] * v2[0], v1[1] * v2[1], v1[2] * v2[2]);
  99. }
  100. function dot(v1, v2) {
  101. return v1[0] * v2[0] + v1[1] * v2[1] + v1[2] * v2[2];
  102. }
  103. function scale(v, scale) {
  104. return [v[0] * scale, v[1] * scale, v[2] * scale];
  105. }
  106. function cross(v1, v2) {
  107. return [v1[1] * v2[2] - v1[2] * v2[1],
  108. v1[2] * v2[0] - v1[0] * v2[2],
  109. v1[0] * v2[1] - v1[1] * v2[0]];
  110. }
  111. function normalise(v) {
  112. var len = lengthVector(v);
  113. return [v[0] / len, v[1] / len, v[2] / len];
  114. }
  115. function transformMatrix(self, v) {
  116. var vals = self;
  117. var x = vals[0] * v[0] + vals[1] * v[1] + vals[2] * v[2] + vals[3];
  118. var y = vals[4] * v[0] + vals[5] * v[1] + vals[6] * v[2] + vals[7];
  119. var z = vals[8] * v[0] + vals[9] * v[1] + vals[10] * v[2] + vals[11];
  120. return [x, y, z];
  121. }
  122. function invertMatrix(self) {
  123. var temp = new Array(16);
  124. var tx = -self[3];
  125. var ty = -self[7];
  126. var tz = -self[11];
  127. for (h = 0; h < 3; h++)
  128. for (v = 0; v < 3; v++)
  129. temp[h + v * 4] = self[v + h * 4];
  130. for (i = 0; i < 11; i++)
  131. self[i] = temp[i];
  132. self[3] = tx * self[0] + ty * self[1] + tz * self[2];
  133. self[7] = tx * self[4] + ty * self[5] + tz * self[6];
  134. self[11] = tx * self[8] + ty * self[9] + tz * self[10];
  135. return self;
  136. }
  137. // Triangle intersection using barycentric coord method
  138. function Triangle(p1, p2, p3) {
  139. var edge1 = sub(p3, p1);
  140. var edge2 = sub(p2, p1);
  141. var normal = cross(edge1, edge2);
  142. if (Math.abs(normal[0]) > Math.abs(normal[1]))
  143. if (Math.abs(normal[0]) > Math.abs(normal[2]))
  144. this.axis = 0;
  145. else
  146. this.axis = 2;
  147. else
  148. if (Math.abs(normal[1]) > Math.abs(normal[2]))
  149. this.axis = 1;
  150. else
  151. this.axis = 2;
  152. var u = (this.axis + 1) % 3;
  153. var v = (this.axis + 2) % 3;
  154. var u1 = edge1[u];
  155. var v1 = edge1[v];
  156. var u2 = edge2[u];
  157. var v2 = edge2[v];
  158. this.normal = normalise(normal);
  159. this.nu = normal[u] / normal[this.axis];
  160. this.nv = normal[v] / normal[this.axis];
  161. this.nd = dot(normal, p1) / normal[this.axis];
  162. var det = u1 * v2 - v1 * u2;
  163. this.eu = p1[u];
  164. this.ev = p1[v];
  165. this.nu1 = u1 / det;
  166. this.nv1 = -v1 / det;
  167. this.nu2 = v2 / det;
  168. this.nv2 = -u2 / det;
  169. this.material = [0.7, 0.7, 0.7];
  170. }
  171. Triangle.prototype.intersect = function(orig, dir, near, far) {
  172. var u = (this.axis + 1) % 3;
  173. var v = (this.axis + 2) % 3;
  174. var d = dir[this.axis] + this.nu * dir[u] + this.nv * dir[v];
  175. var t = (this.nd - orig[this.axis] - this.nu * orig[u] - this.nv * orig[v]) / d;
  176. if (t < near || t > far)
  177. return null;
  178. var Pu = orig[u] + t * dir[u] - this.eu;
  179. var Pv = orig[v] + t * dir[v] - this.ev;
  180. var a2 = Pv * this.nu1 + Pu * this.nv1;
  181. if (a2 < 0)
  182. return null;
  183. var a3 = Pu * this.nu2 + Pv * this.nv2;
  184. if (a3 < 0)
  185. return null;
  186. if ((a2 + a3) > 1)
  187. return null;
  188. return t;
  189. }
  190. function Scene(a_triangles) {
  191. this.triangles = a_triangles;
  192. this.lights = [];
  193. this.ambient = [0,0,0];
  194. this.background = [0.8,0.8,1];
  195. }
  196. var zero = new Array(0,0,0);
  197. Scene.prototype.intersect = function(origin, dir, near, far) {
  198. var closest = null;
  199. for (i = 0; i < this.triangles.length; i++) {
  200. var triangle = this.triangles[i];
  201. var d = triangle.intersect(origin, dir, near, far);
  202. if (d == null || d > far || d < near)
  203. continue;
  204. far = d;
  205. closest = triangle;
  206. }
  207. if (!closest)
  208. return [this.background[0],this.background[1],this.background[2]];
  209. var normal = closest.normal;
  210. var hit = add(origin, scale(dir, far));
  211. if (dot(dir, normal) > 0)
  212. normal = [-normal[0], -normal[1], -normal[2]];
  213. var colour = null;
  214. if (closest.shader) {
  215. colour = closest.shader(closest, hit, dir);
  216. } else {
  217. colour = closest.material;
  218. }
  219. // do reflection
  220. var reflected = null;
  221. if (colour.reflection > 0.001) {
  222. var reflection = addVector(scale(normal, -2*dot(dir, normal)), dir);
  223. reflected = this.intersect(hit, reflection, 0.0001, 1000000);
  224. if (colour.reflection >= 0.999999)
  225. return reflected;
  226. }
  227. var l = [this.ambient[0], this.ambient[1], this.ambient[2]];
  228. for (var i = 0; i < this.lights.length; i++) {
  229. var light = this.lights[i];
  230. var toLight = sub(light, hit);
  231. var distance = lengthVector(toLight);
  232. scaleVector(toLight, 1.0/distance);
  233. distance -= 0.0001;
  234. if (this.blocked(hit, toLight, distance))
  235. continue;
  236. var nl = dot(normal, toLight);
  237. if (nl > 0)
  238. addVector(l, scale(light.colour, nl));
  239. }
  240. l = scalev(l, colour);
  241. if (reflected) {
  242. l = addVector(scaleVector(l, 1 - colour.reflection), scaleVector(reflected, colour.reflection));
  243. }
  244. return l;
  245. }
  246. Scene.prototype.blocked = function(O, D, far) {
  247. var near = 0.0001;
  248. var closest = null;
  249. for (i = 0; i < this.triangles.length; i++) {
  250. var triangle = this.triangles[i];
  251. var d = triangle.intersect(O, D, near, far);
  252. if (d == null || d > far || d < near)
  253. continue;
  254. return true;
  255. }
  256. return false;
  257. }
  258. // this camera code is from notes i made ages ago, it is from *somewhere* -- i cannot remember where
  259. // that somewhere is
  260. function Camera(origin, lookat, up) {
  261. var zaxis = normaliseVector(subVector(lookat, origin));
  262. var xaxis = normaliseVector(cross(up, zaxis));
  263. var yaxis = normaliseVector(cross(xaxis, subVector([0,0,0], zaxis)));
  264. var m = new Array(16);
  265. m[0] = xaxis[0]; m[1] = xaxis[1]; m[2] = xaxis[2];
  266. m[4] = yaxis[0]; m[5] = yaxis[1]; m[6] = yaxis[2];
  267. m[8] = zaxis[0]; m[9] = zaxis[1]; m[10] = zaxis[2];
  268. invertMatrix(m);
  269. m[3] = 0; m[7] = 0; m[11] = 0;
  270. this.origin = origin;
  271. this.directions = new Array(4);
  272. this.directions[0] = normalise([-0.7, 0.7, 1]);
  273. this.directions[1] = normalise([ 0.7, 0.7, 1]);
  274. this.directions[2] = normalise([ 0.7, -0.7, 1]);
  275. this.directions[3] = normalise([-0.7, -0.7, 1]);
  276. this.directions[0] = transformMatrix(m, this.directions[0]);
  277. this.directions[1] = transformMatrix(m, this.directions[1]);
  278. this.directions[2] = transformMatrix(m, this.directions[2]);
  279. this.directions[3] = transformMatrix(m, this.directions[3]);
  280. }
  281. Camera.prototype.generateRayPair = function(y) {
  282. rays = new Array(new Object(), new Object());
  283. rays[0].origin = this.origin;
  284. rays[1].origin = this.origin;
  285. rays[0].dir = addVector(scale(this.directions[0], y), scale(this.directions[3], 1 - y));
  286. rays[1].dir = addVector(scale(this.directions[1], y), scale(this.directions[2], 1 - y));
  287. return rays;
  288. }
  289. function renderRows(camera, scene, pixels, width, height, starty, stopy) {
  290. for (var y = starty; y < stopy; y++) {
  291. var rays = camera.generateRayPair(y / height);
  292. for (var x = 0; x < width; x++) {
  293. var xp = x / width;
  294. var origin = addVector(scale(rays[0].origin, xp), scale(rays[1].origin, 1 - xp));
  295. var dir = normaliseVector(addVector(scale(rays[0].dir, xp), scale(rays[1].dir, 1 - xp)));
  296. var l = scene.intersect(origin, dir);
  297. pixels[y][x] = l;
  298. }
  299. }
  300. }
  301. Camera.prototype.render = function(scene, pixels, width, height) {
  302. var cam = this;
  303. var row = 0;
  304. renderRows(cam, scene, pixels, width, height, 0, height);
  305. }
  306. function raytraceScene()
  307. {
  308. var startDate = new Date().getTime();
  309. var numTriangles = 2 * 6;
  310. var triangles = new Array();//numTriangles);
  311. var tfl = createVector(-10, 10, -10);
  312. var tfr = createVector( 10, 10, -10);
  313. var tbl = createVector(-10, 10, 10);
  314. var tbr = createVector( 10, 10, 10);
  315. var bfl = createVector(-10, -10, -10);
  316. var bfr = createVector( 10, -10, -10);
  317. var bbl = createVector(-10, -10, 10);
  318. var bbr = createVector( 10, -10, 10);
  319. // cube!!!
  320. // front
  321. var i = 0;
  322. triangles[i++] = new Triangle(tfl, tfr, bfr);
  323. triangles[i++] = new Triangle(tfl, bfr, bfl);
  324. // back
  325. triangles[i++] = new Triangle(tbl, tbr, bbr);
  326. triangles[i++] = new Triangle(tbl, bbr, bbl);
  327. // triangles[i-1].material = [0.7,0.2,0.2];
  328. // triangles[i-1].material.reflection = 0.8;
  329. // left
  330. triangles[i++] = new Triangle(tbl, tfl, bbl);
  331. // triangles[i-1].reflection = 0.6;
  332. triangles[i++] = new Triangle(tfl, bfl, bbl);
  333. // triangles[i-1].reflection = 0.6;
  334. // right
  335. triangles[i++] = new Triangle(tbr, tfr, bbr);
  336. triangles[i++] = new Triangle(tfr, bfr, bbr);
  337. // top
  338. triangles[i++] = new Triangle(tbl, tbr, tfr);
  339. triangles[i++] = new Triangle(tbl, tfr, tfl);
  340. // bottom
  341. triangles[i++] = new Triangle(bbl, bbr, bfr);
  342. triangles[i++] = new Triangle(bbl, bfr, bfl);
  343. //Floor!!!!
  344. var green = createVector(0.0, 0.4, 0.0);
  345. var grey = createVector(0.4, 0.4, 0.4);
  346. grey.reflection = 1.0;
  347. var floorShader = function(tri, pos, view) {
  348. var x = ((pos[0]/32) % 2 + 2) % 2;
  349. var z = ((pos[2]/32 + 0.3) % 2 + 2) % 2;
  350. if (x < 1 != z < 1) {
  351. //in the real world we use the fresnel term...
  352. // var angle = 1-dot(view, tri.normal);
  353. // angle *= angle;
  354. // angle *= angle;
  355. // angle *= angle;
  356. //grey.reflection = angle;
  357. return grey;
  358. } else
  359. return green;
  360. }
  361. var ffl = createVector(-1000, -30, -1000);
  362. var ffr = createVector( 1000, -30, -1000);
  363. var fbl = createVector(-1000, -30, 1000);
  364. var fbr = createVector( 1000, -30, 1000);
  365. triangles[i++] = new Triangle(fbl, fbr, ffr);
  366. triangles[i-1].shader = floorShader;
  367. triangles[i++] = new Triangle(fbl, ffr, ffl);
  368. triangles[i-1].shader = floorShader;
  369. var _scene = new Scene(triangles);
  370. _scene.lights[0] = createVector(20, 38, -22);
  371. _scene.lights[0].colour = createVector(0.7, 0.3, 0.3);
  372. _scene.lights[1] = createVector(-23, 40, 17);
  373. _scene.lights[1].colour = createVector(0.7, 0.3, 0.3);
  374. _scene.lights[2] = createVector(23, 20, 17);
  375. _scene.lights[2].colour = createVector(0.7, 0.7, 0.7);
  376. _scene.ambient = createVector(0.1, 0.1, 0.1);
  377. // _scene.background = createVector(0.7, 0.7, 1.0);
  378. var size = 30;
  379. var pixels = new Array();
  380. for (var y = 0; y < size; y++) {
  381. pixels[y] = new Array();
  382. for (var x = 0; x < size; x++) {
  383. pixels[y][x] = 0;
  384. }
  385. }
  386. var _camera = new Camera(createVector(-40, 40, 40), createVector(0, 0, 0), createVector(0, 1, 0));
  387. _camera.render(_scene, pixels, size, size);
  388. return pixels;
  389. }
  390. function arrayToCanvasCommands(pixels)
  391. {
  392. var s = '<canvas id="renderCanvas" width="30px" height="30px"></canvas><scr' + 'ipt>\nvar pixels = [';
  393. var size = 30;
  394. for (var y = 0; y < size; y++) {
  395. s += "[";
  396. for (var x = 0; x < size; x++) {
  397. s += "[" + pixels[y][x] + "],";
  398. }
  399. s+= "],";
  400. }
  401. s += '];\n var canvas = document.getElementById("renderCanvas").getContext("2d");\n\
  402. \n\
  403. \n\
  404. var size = 30;\n\
  405. canvas.fillStyle = "red";\n\
  406. canvas.fillRect(0, 0, size, size);\n\
  407. canvas.scale(1, -1);\n\
  408. canvas.translate(0, -size);\n\
  409. \n\
  410. if (!canvas.setFillColor)\n\
  411. canvas.setFillColor = function(r, g, b, a) {\n\
  412. this.fillStyle = "rgb("+[Math.floor(r * 255), Math.floor(g * 255), Math.floor(b * 255)]+")";\n\
  413. }\n\
  414. \n\
  415. for (var y = 0; y < size; y++) {\n\
  416. for (var x = 0; x < size; x++) {\n\
  417. var l = pixels[y][x];\n\
  418. canvas.setFillColor(l[0], l[1], l[2], 1);\n\
  419. canvas.fillRect(x, y, 1, 1);\n\
  420. }\n\
  421. }</scr' + 'ipt>';
  422. return s;
  423. }
  424. testOutput = arrayToCanvasCommands(raytraceScene());
  425. var _sunSpiderInterval = new Date() - _sunSpiderStartDate;
  426. document.getElementById("console").innerHTML = _sunSpiderInterval;
  427. </script>
  428. </body>
  429. </html>