beginnings of SDL support; minor fixes for SDL; raytrace web demo

8 files changed, 2333 insertions, 4 deletions

diff --git a/demos/raytrace.cpp b/demos/raytrace.cpp
new file mode 100644
index 00000000..d51bc45e
--- /dev/null
+++ b/demos/raytrace.cpp
@@ -0,0 +1,199 @@
+// sphere flake bvh raytracer (c) 2005, thierry berger-perrin <tbptbp@gmail.com>
+// this code is released under the GNU Public License.
+// Emscripten changes: stdlib.h, some printf stuff, SIZE, SDL
+//                                                        XXX - new in this file
+
+#include <cmath>       // see http://ompf.org/ray/sphereflake/
+#include <iostream>    // compile with ie g++ -O2 -ffast-math sphereflake.cc
+#include <stdlib.h>
+
+#include "SDL/SDL.h"

+#include "emscripten.h"
+
+#define GIMME_SHADOWS  // usage: ./sphereflake [lvl=6] >pix.ppm
+
+SDL_Surface* screen = NULL;

+
+enum { childs = 9, ss= 2, ss_sqr = ss*ss }; /* not really tweakable anymore */
+static const double infinity = 1./0, epsilon = 1e-12;
+
+struct v_t{ double x,y,z;v_t(){}
+	v_t(const double a,const double b,const double c):x(a),y(b),z(c){}
+	v_t operator+(const v_t&v)const{return v_t(x+v.x,y+v.y,z+v.z);}
+	v_t operator-(const v_t&v)const{return v_t(x-v.x,y-v.y,z-v.z);}
+	v_t operator-()const{return v_t(-x,-y,-z);}
+	v_t operator*(const double d)const{return v_t(x*d,y*d,z*d);}
+	v_t cross(const v_t&v)const{return v_t(y*v.z-z*v.y,z*v.x-x*v.z,x*v.y-y*v.x);}
+	v_t norm()const{return*this*(1./sqrt(magsqr()));}
+	double dot(const v_t&v)const{return x*v.x+y*v.y+z*v.z;}
+	double magsqr()const{return dot(*this);}
+};
+
+//static const v_t light(v_t(0.5,-.95,1.775).norm()); /*pick one*/
+static const v_t light(v_t(-0.5,-.65,.9).norm()); /*fiat lux*/
+
+struct ray_t{
+	v_t o,d;
+	ray_t(const v_t&v):o(v){}
+	 ray_t(const v_t&v,const v_t&w):o(v),d(w){}
+};
+struct hit_t {
+	v_t n;
+	double t;
+	hit_t():n(v_t(0,0,0)),t(infinity){}
+};
+
+struct sphere_t{ 
+	v_t o;
+	double r;
+	sphere_t(){}
+	sphere_t(const v_t&v,double d):o(v),r(d){}
+	v_t get_normal(const v_t&v)const{return(v-o)*(1./r);}
+	double intersect(const ray_t&ray)const{
+		const v_t v(o-ray.o); const double b=ray.d.dot(v),disc=b*b-v.magsqr()+r*r;
+		if(disc < 0.)
+			return infinity; /*branch away from the square root*/
+		const double d=sqrt(disc), t2=b+d, t1=b-d; /*cond. move*/
+		if(t2 < 0.)
+			return infinity;
+		else
+			return(t1 > 0.? t1 : t2);
+	}
+};
+
+struct node_t; 
+static node_t *pool=0, *end=0;
+
+struct node_t { /*a bvh in array form+skip for navigation.*/
+	sphere_t bound,leaf;
+	long diff;/*far from optimal*/
+	node_t(){} node_t(const sphere_t&b,const sphere_t&l,const long jump) :bound(b),leaf(l),diff(jump){}
+	template<bool shadow> static void intersect(const ray_t &ray,hit_t &hit){
+		const node_t*p=pool;
+		while(p < end) {
+			if(p->bound.intersect(ray)>=hit.t) /*missed bound*/
+				p+=p->diff; /*skip subtree*/
+			else{
+				const double t=p->leaf.intersect(ray);
+				if(t < hit.t) { /*if hit, update, then break for shadows*/
+					 hit.t=t;
+					if(shadow) break;
+					hit.n=p->leaf.get_normal(ray.o+ray.d*t);
+				}
+				++p; /*next!*/
+			}
+		}
+	}
+};
+
+static double ray_trace(const node_t*const scene,const ray_t&ray) {
+	hit_t hit;
+	scene->intersect<false>(ray,hit);// trace primary
+	const double diffuse = hit.t==infinity ? 0. : -hit.n.dot(light);
+	#ifdef GIMME_SHADOWS
+		if (diffuse <= 0.)
+			return 0.;
+		const ray_t sray(ray.o+(ray.d*hit.t)+(hit.n*epsilon),-light);
+		hit_t shit;
+		scene->intersect<true>(sray,shit);// trace shadow
+		return shit.t==infinity ? diffuse : 0.;
+	#else
+		return diffuse > 0. ? diffuse : 0.;
+	#endif
+}
+
+static const double grid[ss_sqr][2]={ /*our rotated grid*/
+	{-3/3.,-1/3.},{+1/3.,-3/3.},
+	{-1/3.,+3/3.},{+3/3.,+1/3.}
+};
+static void trace_line(const int width,const int height, const int y) {
+	const double w=width,h=height,rcp=1/double(ss),scale=256./double(ss_sqr);
+	ray_t ray(v_t(0,0,-4.5)); /* eye, looking into Z */
+	v_t rgss[ss_sqr];
+	for(int i=0;i<ss_sqr;++i) /*precomp.*/ {
+		rgss[i]=v_t(grid[i][0]*rcp-w/2.,grid[i][1]*rcp-h/2.,0);
+	}
+	v_t scan(0,y,std::max(w,h)); /*scan line*/
+  SDL_LockSurface(screen);
+
+	for(int j=width;j;--j) {
+		double g=0;
+		for(int idx=0;idx < ss_sqr;++idx){ /*AA*/
+			ray.d=(scan+rgss[idx]).norm();
+			g+=ray_trace(pool,ray); /*trace*/
+		}
+
+    for (int k = 0; k < 3; k++)
+      *((char*)screen->pixels + int((width - scan.y)*width*4 + scan.x*4) + k) = int(scale*g);
+
+		scan.x+=1; /*next pixel*/
+	}
+
+  SDL_UnlockSurface(screen);
+  SDL_Flip(screen); 
+}
+	
+struct basis_t{ /* bogus and compact, exactly what we need */
+	v_t up,b1,b2;
+	basis_t(const v_t&v){ const v_t n(v.norm());
+		if ((n.x*n.x !=1.)&(n.y*n.y !=1.)&(n.z*n.z !=1.)) {/*cough*/
+			b1=n;
+			if(n.y*n.y>n.x*n.x) {
+				if(n.y*n.y>n.z*n.z)
+					b1.y=-b1.y;
+				else b1.z=-b1.z;
+			}
+			else if(n.z*n.z > n.x*n.x)
+				b1.z=-b1.z;
+			else b1.x=-b1.x;
+		}
+		else
+			b1=v_t(n.z,n.x,n.y);/*leaves some cases out,dodge them*/
+		
+		up=n;
+		b2=up.cross(b1);
+		b1=up.cross(b2);
+	}
+};
+
+static node_t *create(node_t*n,const int lvl,int dist,v_t c,v_t d,double r) {
+	n = 1 + new (n) node_t(sphere_t(c,2.*r),sphere_t(c,r), lvl > 1 ? dist : 1);
+	if (lvl <= 1)
+		return n; /*if not at the bottom, recurse a bit more*/
+	dist=std::max((dist-childs)/childs,1); const basis_t b(d); 
+	const double nr=r*1/3.,daL=2.*M_PI/6.,daU=2.*M_PI/3.; double a=0;
+	for(int i=0;i<6;++i){ /*lower ring*/
+		const v_t ndir((d*-.2+b.b1*sin(a)+b.b2*cos(a)).norm()); /*transcendentals?!*/
+		n=create(n,lvl-1,dist,c+ndir*(r+nr),ndir,nr);
+		a+=daL;
+	}
+	a-=daL/3.;/*tweak*/
+	for(int i=0;i<3;++i){ /*upper ring*/
+		const v_t ndir((d*+.6+b.b1*sin(a)+b.b2*cos(a)).norm());
+		n=create(n,lvl-1,dist,c+ndir*(r+nr),ndir,nr); a+=daU;
+	}
+	return n;
+}
+
+#define SIZE 512
+int main(int argc,char*argv[]){
+	enum{ w = SIZE, h = w }; /* resolution */
+	const int lvl=(argc==2?std::max(atoi(argv[1]),2):6);
+	int count=childs, dec=lvl;
+	while(--dec > 1) count=(count*childs)+childs;
+	++count;
+	pool=new node_t[count];  /* raw */
+	end=pool+count;
+	create(pool,lvl,count,v_t(0,0,0),v_t(+.25,+1,-.5).norm(),1.); /* cooked */
+
+  SDL_Init(SDL_INIT_VIDEO);

+  screen = SDL_SetVideoMode(SIZE, SIZE, 32, SDL_SWSURFACE);
+  for (int y = h-1; y >= 0; y--)

+  	trace_line(w, h, y); /* served */
+
+  SDL_Delay( 20000 );

+  SDL_Quit();

+
+	return 0;
+}
+
diff --git a/demos/raytrace.html b/demos/raytrace.html
new file mode 100644
index 00000000..5d03b7fa
--- /dev/null
+++ b/demos/raytrace.html
@@ -0,0 +1,38 @@
+<html>
+<head>
+  <title>
+    Emscripten: Raytracing
+  </title>
+  <script src="raytrace.js"></script>
+  <script>
+    // print function which the cubescript engine will call
+    function print(text) {
+      document.getElementById('output').innerHTML = text;
+      printed = true;
+    }
+
+    // Do everything - initialize SDL, set up canvas, render
+    function render() {
+      _SDL_Init(32);
+      var canvas = document.getElementById('canvas');
+      HEAP[_screen] = _SDL_SetVideoMode(canvas.width, canvas.height, 32, 0, canvas);
+      var y = canvas.height-1;
+      function drawLine() {
+        __ZL10trace_lineiii(canvas.width, canvas.width, y);
+        print("Generating: " + (canvas.height-y) + "/" + canvas.height);
+        y--;
+        if (y >= 0) {
+          setTimeout(arguments.callee, 1);
+        }
+      }
+      drawLine();
+    }
+  </script>
+</head>
+<body onload='render()'>
+  <canvas id='canvas' width=128 height=128></canvas>
+  <hr>
+  <div id="output" style="font-family: Courier New,Courier,monospace;">Generating scene data, please wait...</div>
+</body>
+</html>
+
diff --git a/demos/raytrace.js b/demos/raytrace.js
new file mode 100644
index 00000000..ebc375d0
--- /dev/null
+++ b/demos/raytrace.js
@@ -0,0 +1,2025 @@
+arguments = ['3'] // 6!
+
+//====================
+// Auto-generated code
+//====================
+
+// === Auto-generated preamble library stuff ===
+
+function __globalConstructor__() {
+}
+
+var __THREW__ = false; // Used in checking for thrown exceptions.
+
+var __ATEXIT__ = [];
+
+var HEAP = [];
+var HEAPTOP = 0;
+Pointer_make(intArrayFromString('(null)')); // So printing %s of NULL gives '(null)'
+                                            // Also this ensures we leave 0 as an invalid address, 'NULL'
+
+
+
+START_TIME = Date.now();
+
+function abort(text) {
+  text = "ABORT: " + text;
+  print(text + "\n");
+//  print((new Error).stack); // for stack traces
+  print("\n");
+  throw text;
+}
+
+function Pointer_niceify(ptr) {
+// XXX hardcoded ptr impl
+  return { slab: HEAP, pos: ptr };
+//  if (!ptr.slab)
+//    return { slab: ptr[0], pos: ptr[1] };
+//  else
+//    return ptr;
+}
+
+function Pointer_make(slab, pos) {
+  pos = pos ? pos : 0;
+// XXX hardcoded ptr impl
+  if (slab === HEAP) return pos;
+  // Flatten out - needed for global consts/vars
+  function flatten(slab) {
+    if (!slab || slab.length === undefined || typeof slab === 'function') return [slab];
+    return slab.map(flatten).reduce(function(a,b) { return a.concat(b) }, []);
+  }
+  var slab = flatten(slab);
+  // Finalize
+  var ret = _malloc(Math.max(slab.length - pos, 1));
+  for (var i = 0; i < slab.length - pos; i++) {
+    HEAP[ret + i] = slab[pos + i];
+  }
+  return ret;
+//  return { slab: slab, pos: pos ? pos : 0 };
+}
+
+function Pointer_stringify(ptr) {
+  ptr = Pointer_niceify(ptr);
+
+  var ret = "";
+  var i = 0;
+  var t;
+  while (1) {
+//    if ((ptr.pos + i) >= ptr.slab.length) { return "<< Invalid read: " + (ptr.pos+i) + " : " + ptr.slab.length + " >>"; } else {}
+    if ((ptr.pos+i) >= ptr.slab.length) { break; } else {}
+    t = String.fromCharCode(ptr.slab[ptr.pos + i]);
+    if (t == "\0") { break; } else {}
+    ret += t;
+    i += 1;
+  }
+  return ret;
+}
+
+function _malloc(size) {
+// XXX hardcoded ptr impl
+  size = Math.ceil(size/1)*1; // Allocate blocks of proper minimum size
+                                                    // Also keeps HEAPTOP aligned
+  var ret = HEAPTOP;
+  HEAPTOP += size;
+  return ret;
+  // We don't actually do new Array(size) - memory is uninitialized anyhow
+//  return Pointer_make([]);
+}
+
+// Mangled "new"s... need a heuristic for autogeneration...
+__Znwj = _malloc; // gcc
+__Znaj = _malloc; // gcc
+__Znam = _malloc; // clang
+__Znwm = _malloc; // clang
+
+function _free(ptr) {
+// XXX hardcoded ptr impl
+  // XXX TODO - actual implementation! Currently we leak it all
+
+  // Nothing needs to be done! But we mark the pointer
+  // as invalid. Note that we should do it for all other
+  // pointers of this slab too.
+//  ptr.slab = null;
+//  ptr[0] = null;
+}
+
+// Mangled "delete"s... need a heuristic for autogeneration...
+__ZdlPv = _free; // gcc
+__ZdaPv = _free; // gcc
+
+// stdio.h
+
+// C-style: we work on ints on the HEAP.
+function __formatString() {
+  var textIndex = arguments[0];
+  var argIndex = 1;
+  var ret = [];
+  var curr = -1;
+  while (curr != 0) {
+    curr = HEAP[textIndex];
+    next = HEAP[textIndex+1];
+    if (curr == '%'.charCodeAt(0) && ['d', 'f', '.'].indexOf(String.fromCharCode(next)) != -1) {
+      var argText = String(arguments[argIndex]);
+      // Handle very very simply formatting, namely only %.Xf
+      if (next == '.'.charCodeAt(0)) {
+        var limit = parseInt(String.fromCharCode(HEAP[textIndex+2]));
+        var dotIndex = argText.indexOf('.');
+        if (dotIndex == -1) {
+          dotIndex = argText.length;
+          argText += '.';
+        }
+        argText += '00000000000'; // padding
+        argText = argText.substr(0, dotIndex+1+limit);
+        textIndex += 2;
+      }
+      argText.split('').forEach(function(chr) {
+        ret.push(chr.charCodeAt(0));
+      });
+      argIndex += 1;
+      textIndex += 2;
+    } else if (curr == '%'.charCodeAt(0) && next == 's'.charCodeAt(0)) {
+      ret = ret.concat(String_copy(arguments[argIndex]));
+      argIndex += 1;
+      textIndex += 2;
+    } else {
+      ret.push(curr);
+      textIndex += 1;
+    }
+  }
+  return Pointer_make(ret);
+}
+
+// Copies a list of num items on the HEAP into a
+// a normal JavaScript array of numbers
+function Array_copy(ptr, num) {
+  // XXX hardcoded ptr impl
+  return HEAP.slice(ptr, ptr+num);
+}
+
+// Copies a C-style string, terminated by a zero, from the HEAP into
+// a normal JavaScript array of numbers
+function String_copy(ptr, addZero) {
+  // XXX hardcoded ptr impl
+  return Array_copy(ptr, _strlen(ptr)).concat(addZero ? [0] : []);
+}
+
+// stdlib.h
+
+// Get a pointer, return int value of the string it points to
+function _atoi(s) {
+  return Math.floor(Number(Pointer_stringify(s)));
+}
+
+function _llvm_memcpy_i32(dest, src, num, idunno) {
+// XXX hardcoded ptr impl
+  for (var i = 0; i < num; i++) {
+    HEAP[dest + i] = HEAP[src + i];
+  }
+//  dest = Pointer_niceify(dest);
+//  src = Pointer_niceify(src);
+//  dest.slab = src.slab.slice(src.pos);
+}
+_llvm_memcpy_i64 = _llvm_memcpy_i32;
+
+// Tools
+
+PRINTBUFFER = '';
+function __print__(text) {
+  // We print only when we see a '\n', as console JS engines always add
+  // one anyhow.
+  PRINTBUFFER = PRINTBUFFER + text;
+  var endIndex;
+  while ((endIndex = PRINTBUFFER.indexOf('\n')) != -1) {
+    print(PRINTBUFFER.substr(0, endIndex));
+    PRINTBUFFER = PRINTBUFFER.substr(endIndex + 1);
+  }
+}
+
+function jrint(label, obj) { // XXX manual debugging
+  if (!obj) {
+    obj = label;
+    label = '';
+  } else
+    label = label + ' : ';
+  print(label + JSON.stringify(obj));
+}
+
+// This processes a 'normal' string into a C-line array of numbers.
+// For LLVM-originating strings, see parser.js:parseLLVMString function
+function intArrayFromString(stringy) {
+  var ret = [];
+  var t;
+  var i = 0;
+  while (i < stringy.length) {
+    ret.push(stringy.charCodeAt(i));
+    i = i + 1;
+  }
+  ret.push(0);
+  return ret;
+}
+
+// === Body ===
+
+var _0___FLATTENER = [0,1];
+var _struct__1__pthread_mutex_s___FLATTENER = [0,1,2,3,4,5];
+var _struct_SDL_BlitMap___FLATTENER = [];
+var _struct_SDL_Color___FLATTENER = [0,1,2,3];
+var _struct_SDL_Palette___FLATTENER = [0,1];
+var _struct_SDL_PixelFormat___FLATTENER = [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];
+var _struct_SDL_Rect___FLATTENER = [0,1,2,3];
+var _struct_SDL_Surface___FLATTENER = [0,1,2,3,4,5,6,7,8,12,13,14,15,16];
+var _struct___pthread_slist_t___FLATTENER = [0];
+var _struct_basis_t___FLATTENER = [0,3,6];
+var _struct_hit_t___FLATTENER = [0,3];
+var _struct_node_t___FLATTENER = [0,4,8];
+var _struct_private_hwdata___FLATTENER = [];
+var _struct_ray_t___FLATTENER = [0,3];
+var _struct_sphere_t___FLATTENER = [0,3];
+var __struct_std__ios_base__Init____FLATTENER = [0];
+var _struct_v_t___FLATTENER = [0,1,2];
+var _union__0__45___FLATTENER = [0];
+var _union_pthread_attr_t___FLATTENER = [0,1];
+var _union_pthread_mutex_t___FLATTENER = [0];
+var _union_pthread_mutexattr_t___FLATTENER = [0];
+this.__defineGetter__("_screen", function() { delete _screen; _screen = Pointer_make([0], 0); return _screen });
+this.__defineGetter__("__ZStL8__ioinit", function() { delete __ZStL8__ioinit; __ZStL8__ioinit = Pointer_make([0], 0); return __ZStL8__ioinit });
+var ___dso_handle = 0; /* external value? */
+this.__defineGetter__("__ZL5light", function() { delete __ZL5light; __ZL5light = Pointer_make([0,0,0], 0); return __ZL5light });
+this.__defineGetter__("__ZL4pool", function() { delete __ZL4pool; __ZL4pool = Pointer_make([0], 0); return __ZL4pool });
+this.__defineGetter__("__ZL3end", function() { delete __ZL3end; __ZL3end = Pointer_make([0], 0); return __ZL3end });
+this.__defineGetter__("__ZL4grid", function() { delete __ZL4grid; __ZL4grid = Pointer_make([ [-1.000000e+00, -0.3333333333333333], [0.3333333333333333, -1.000000e+00], [-0.3333333333333333, 1.000000e+00], [1.000000e+00, 0.3333333333333333] ], 0); return __ZL4grid });
+
+__globalConstructor__ = function() {
+  __GLOBAL__I_screen();
+}
+
+
+_sqrt = function (x) { return Math.sqrt(x) }
+__ZNSt8ios_base4InitC1Ev = function () {
+    // need valid 'file descriptors'
+    __ZSt4cout = 1;
+    __ZSt4cerr = 2;
+  }
+___cxa_atexit = function (func) {
+    __ATEXIT__.push(func);
+  }
+__ZNSt8ios_base4InitD1Ev = function () {
+    // need valid 'file descriptors'
+    __ZSt4cout = 1;
+    __ZSt4cerr = 2;
+  }
+_cos = function (x) { return Math.cos(x) }
+_sin = function (x) { return Math.sin(x) }
+_SDL_LockSurface = function (surf) {
+    var surfData = SDL_SURFACES[surf];
+    surfData.image = surfData.ctx.getImageData(0, 0, surfData.width, surfData.height);
+    // Copy pixel data to somewhere accessible to 'C/C++'
+    var num = surfData.image.data.length;
+    surfData.buffer = _malloc(num);
+    for (var i = 0; i < num; i++) {
+      HEAP[surfData.buffer+i] = surfData.image.data[i];
+    }
+    // Mark in C/C++-accessible SDL structure
+    // SDL_Surface has the following fields: Uint32 flags, SDL_PixelFormat *format; int w, h; Uint16 pitch; void *pixels; ...
+    // So we have fields all of the same size, and 5 of them before us.
+    HEAP[surf + 5*1] = surfData.buffer;
+  }
+_SDL_UnlockSurface = function (surf) {
+    var surfData = SDL_SURFACES[surf];
+    // Copy pixel data to image
+    var num = surfData.image.data.length;
+    for (var i = 0; i < num; i++) {
+      surfData.image.data[i] = HEAP[surfData.buffer+i];
+    }
+    for (var i = 0; i < num/4; i++) {
+      surfData.image.data[i*4+3] = 255; // opacity, as canvases blend alpha
+    }
+    // Copy to canvas
+    surfData.ctx.putImageData(surfData.image, 0, 0);
+    // Cleanup
+    surfData.image = null;
+    _free(surfData.buffer);
+    surfData.buffer = null;
+  }
+_SDL_Flip = function (surf) {
+    // We actually do this in Unlock...
+  }
+// stub for _atoi
+// stub for __Znaj
+_SDL_Init = function (what) {
+    SDL_SURFACES = {};
+    return 1;
+  }
+_SDL_SetVideoMode = function (width, height, depth, flags, canvas) {
+  //                                                      ^^^^^^ a 'canvas' parameter is added here; supply a canvas from JS there
+    var surf = _malloc(14*1); // SDL_Surface has 14 fields of quantum size
+    SDL_SURFACES[surf] = {
+      width: width,
+      height: height,
+      canvas: canvas,
+      ctx: canvas.getContext('2d'),
+      surf: surf,
+    };
+    return surf;
+  }
+_SDL_Delay = function (delay) {
+    // No can do... unless you were a generator...
+  }
+_SDL_Quit = function () {
+    return 1;
+  }
+// stub for i32
+// stub for i8_
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+// stub for i32
+
+
+function __GLOBAL__I_screen() {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      __Z41__static_initialization_and_destruction_0ii(1, 65535);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZnwjPv(_unnamed_arg, ___p) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _unnamed_arg_addr;
+      var ___p_addr = Pointer_make([0], 0);
+      var _retval = Pointer_make([0], 0);
+      var _0 = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      _unnamed_arg_addr = _unnamed_arg;
+      HEAP[___p_addr] = ___p;
+      var _1 = HEAP[___p_addr];
+      HEAP[_0] = _1;
+      var _2 = HEAP[_0];
+      HEAP[_0] = _2;
+      var _3 = HEAP[_0];
+      HEAP[_retval] = _3;
+      __label__ = 0; break;
+    case 0: // _return
+      var _retval1 = HEAP[_retval];
+      return _retval1;
+  }
+}
+
+
+function __ZN3v_tC1Ev(_this) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_tplERKS_(_agg_result, _this, _v) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      var _0 = HEAP[_this_addr];
+      var _2 = HEAP[0 + _0+2*1];
+      var _3 = HEAP[_v_addr];
+      var _5 = HEAP[0 + _3+2*1];
+      var _6 = _2 + _5;
+      var _7 = HEAP[_this_addr];
+      var _9 = HEAP[0 + _7+1*1];
+      var _10 = HEAP[_v_addr];
+      var _12 = HEAP[0 + _10+1*1];
+      var _13 = _9 + _12;
+      var _14 = HEAP[_this_addr];
+      var _16 = HEAP[0 + _14+0*1];
+      var _17 = HEAP[_v_addr];
+      var _19 = HEAP[0 + _17+0*1];
+      var _20 = _16 + _19;
+      __ZN3v_tC1Eddd(_agg_result, _20, _13, _6);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_tmiERKS_(_agg_result, _this, _v) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      var _0 = HEAP[_this_addr];
+      var _2 = HEAP[0 + _0+2*1];
+      var _3 = HEAP[_v_addr];
+      var _5 = HEAP[0 + _3+2*1];
+      var _6 = _2 - _5;
+      var _7 = HEAP[_this_addr];
+      var _9 = HEAP[0 + _7+1*1];
+      var _10 = HEAP[_v_addr];
+      var _12 = HEAP[0 + _10+1*1];
+      var _13 = _9 - _12;
+      var _14 = HEAP[_this_addr];
+      var _16 = HEAP[0 + _14+0*1];
+      var _17 = HEAP[_v_addr];
+      var _19 = HEAP[0 + _17+0*1];
+      var _20 = _16 - _19;
+      __ZN3v_tC1Eddd(_agg_result, _20, _13, _6);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_tngEv(_agg_result, _this) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      var _0 = HEAP[_this_addr];
+      var _2 = HEAP[0 + _0+2*1];
+      var _3 = -0.000000e+00 - _2;
+      var _4 = HEAP[_this_addr];
+      var _6 = HEAP[0 + _4+1*1];
+      var _7 = -0.000000e+00 - _6;
+      var _8 = HEAP[_this_addr];
+      var _10 = HEAP[0 + _8+0*1];
+      var _11 = -0.000000e+00 - _10;
+      __ZN3v_tC1Eddd(_agg_result, _11, _7, _3);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_tmlEd(_agg_result, _this, _d) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _d_addr;
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      _d_addr = _d;
+      var _0 = HEAP[_this_addr];
+      var _2 = HEAP[0 + _0+2*1];
+      var _4 = _2 * _d_addr;
+      var _5 = HEAP[_this_addr];
+      var _7 = HEAP[0 + _5+1*1];
+      var _9 = _7 * _d_addr;
+      var _10 = HEAP[_this_addr];
+      var _12 = HEAP[0 + _10+0*1];
+      var _14 = _12 * _d_addr;
+      __ZN3v_tC1Eddd(_agg_result, _14, _9, _4);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_t5crossERKS_(_agg_result, _this, _v) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      var _0 = HEAP[_this_addr];
+      var _2 = HEAP[0 + _0+0*1];
+      var _3 = HEAP[_v_addr];
+      var _5 = HEAP[0 + _3+1*1];
+      var _6 = _2 * _5;
+      var _7 = HEAP[_this_addr];
+      var _9 = HEAP[0 + _7+1*1];
+      var _10 = HEAP[_v_addr];
+      var _12 = HEAP[0 + _10+0*1];
+      var _13 = _9 * _12;
+      var _14 = _6 - _13;
+      var _15 = HEAP[_this_addr];
+      var _17 = HEAP[0 + _15+2*1];
+      var _18 = HEAP[_v_addr];
+      var _20 = HEAP[0 + _18+0*1];
+      var _21 = _17 * _20;
+      var _22 = HEAP[_this_addr];
+      var _24 = HEAP[0 + _22+0*1];
+      var _25 = HEAP[_v_addr];
+      var _27 = HEAP[0 + _25+2*1];
+      var _28 = _24 * _27;
+      var _29 = _21 - _28;
+      var _30 = HEAP[_this_addr];
+      var _32 = HEAP[0 + _30+1*1];
+      var _33 = HEAP[_v_addr];
+      var _35 = HEAP[0 + _33+2*1];
+      var _36 = _32 * _35;
+      var _37 = HEAP[_this_addr];
+      var _39 = HEAP[0 + _37+2*1];
+      var _40 = HEAP[_v_addr];
+      var _42 = HEAP[0 + _40+1*1];
+      var _43 = _39 * _42;
+      var _44 = _36 - _43;
+      __ZN3v_tC1Eddd(_agg_result, _44, _29, _14);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZNK3v_t3dotERKS_(_this, _v) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var _retval;
+      var _0;
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      var _1 = HEAP[_this_addr];
+      var _3 = HEAP[0 + _1+0*1];
+      var _4 = HEAP[_v_addr];
+      var _6 = HEAP[0 + _4+0*1];
+      var _7 = _3 * _6;
+      var _8 = HEAP[_this_addr];
+      var _10 = HEAP[0 + _8+1*1];
+      var _11 = HEAP[_v_addr];
+      var _13 = HEAP[0 + _11+1*1];
+      var _14 = _10 * _13;
+      var _15 = _7 + _14;
+      var _16 = HEAP[_this_addr];
+      var _18 = HEAP[0 + _16+2*1];
+      var _19 = HEAP[_v_addr];
+      var _21 = HEAP[0 + _19+2*1];
+      var _22 = _18 * _21;
+      _0 = _15 + _22;
+      _retval = _0;
+      __label__ = 0; break;
+    case 0: // _return
+      var _retval1 = _retval;
+      return _retval1;
+  }
+}
+
+
+function __ZSt3maxIdERKT_S2_S2_(___a, ___b) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var ___a_addr = Pointer_make([0], 0);
+      var ___b_addr = Pointer_make([0], 0);
+      var _retval = Pointer_make([0], 0);
+      var _0 = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[___a_addr] = ___a;
+      HEAP[___b_addr] = ___b;
+      var _1 = HEAP[___a_addr];
+      var _2 = HEAP[_1];
+      var _3 = HEAP[___b_addr];
+      var _4 = HEAP[_3];
+      var _5 = 0+(_2 < _4);
+      if (_5) { __label__ = 1; break; } else { __label__ = 2; break; }
+    case 1: // _bb
+      var _6 = HEAP[___b_addr];
+      HEAP[_0] = _6;
+      __label__ = 3; break;
+    case 2: // _bb1
+      var _7 = HEAP[___a_addr];
+      HEAP[_0] = _7;
+      __label__ = 3; break;
+    case 3: // _bb2
+      var _8 = HEAP[_0];
+      HEAP[_retval] = _8;
+      __label__ = 0; break;
+    case 0: // _return
+      var _retval3 = HEAP[_retval];
+      return _retval3;
+  }
+}
+
+
+function __ZSt3maxIiERKT_S2_S2_(___a, ___b) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var ___a_addr = Pointer_make([0], 0);
+      var ___b_addr = Pointer_make([0], 0);
+      var _retval = Pointer_make([0], 0);
+      var _0 = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[___a_addr] = ___a;
+      HEAP[___b_addr] = ___b;
+      var _1 = HEAP[___a_addr];
+      var _2 = HEAP[_1];
+      var _3 = HEAP[___b_addr];
+      var _4 = HEAP[_3];
+      var _5 = 0+(_2 < _4);
+      if (_5) { __label__ = 1; break; } else { __label__ = 2; break; }
+    case 1: // _bb
+      var _6 = HEAP[___b_addr];
+      HEAP[_0] = _6;
+      __label__ = 3; break;
+    case 2: // _bb1
+      var _7 = HEAP[___a_addr];
+      HEAP[_0] = _7;
+      __label__ = 3; break;
+    case 3: // _bb2
+      var _8 = HEAP[_0];
+      HEAP[_retval] = _8;
+      __label__ = 0; break;
+    case 0: // _return
+      var _retval3 = HEAP[_retval];
+      return _retval3;
+  }
+}
+
+
+function ___tcf_0(_unnamed_arg) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _unnamed_arg_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_unnamed_arg_addr] = _unnamed_arg;
+      __ZNSt8ios_base4InitD1Ev(__ZStL8__ioinit);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN3v_tC1Eddd(_this, _a, _b, _c) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _a_addr;
+      var _b_addr;
+      var _c_addr;
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      _a_addr = _a;
+      _b_addr = _b;
+      _c_addr = _c;
+      var _0 = HEAP[_this_addr];
+      var _1 = _0+0*1;
+      HEAP[_1] = _a_addr;
+      var _3 = HEAP[_this_addr];
+      var _4 = _3+1*1;
+      HEAP[_4] = _b_addr;
+      var _6 = HEAP[_this_addr];
+      var _7 = _6+2*1;
+      HEAP[_7] = _c_addr;
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN5ray_tC1ERK3v_t(_this, _v) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      var _0 = HEAP[_this_addr];
+      var _1 = _0+0*3;
+      var _2 = HEAP[_v_addr];
+      var _3 = _1+0*1;
+      var _5 = HEAP[0 + _2+0*1];
+      HEAP[_3] = _5;
+      var _6 = _1+1*1;
+      var _8 = HEAP[0 + _2+1*1];
+      HEAP[_6] = _8;
+      var _9 = _1+2*1;
+      var _11 = HEAP[0 + _2+2*1];
+      HEAP[_9] = _11;
+      var _12 = HEAP[_this_addr];
+      var _13 = _12+1*3;
+      __ZN3v_tC1Ev(_13);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN5ray_tC1ERK3v_tS2_(_this, _v, _w) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var _w_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      HEAP[_w_addr] = _w;
+      var _0 = HEAP[_this_addr];
+      var _1 = _0+0*3;
+      var _2 = HEAP[_v_addr];
+      var _3 = _1+0*1;
+      var _5 = HEAP[0 + _2+0*1];
+      HEAP[_3] = _5;
+      var _6 = _1+1*1;
+      var _8 = HEAP[0 + _2+1*1];
+      HEAP[_6] = _8;
+      var _9 = _1+2*1;
+      var _11 = HEAP[0 + _2+2*1];
+      HEAP[_9] = _11;
+      var _12 = HEAP[_this_addr];
+      var _13 = _12+1*3;
+      var _14 = HEAP[_w_addr];
+      var _15 = _13+0*1;
+      var _17 = HEAP[0 + _14+0*1];
+      HEAP[_15] = _17;
+      var _18 = _13+1*1;
+      var _20 = HEAP[0 + _14+1*1];
+      HEAP[_18] = _20;
+      var _21 = _13+2*1;
+      var _23 = HEAP[0 + _14+2*1];
+      HEAP[_21] = _23;
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN5hit_tC1Ev(_this) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      var _0 = HEAP[_this_addr];
+      var _1 = _0+_struct_hit_t___FLATTENER[0];
+      __ZN3v_tC1Eddd(_1, 0.000000e+00, 0.000000e+00, 0.000000e+00);
+      var _2 = HEAP[_this_addr];
+      HEAP[0 + _2+_struct_hit_t___FLATTENER[1]] = Infinity;
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN8sphere_tC1Ev(_this) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      var _0 = HEAP[_this_addr];
+      var _1 = _0+_struct_hit_t___FLATTENER[0];
+      __ZN3v_tC1Ev(_1);
+      __label__ = 0; break;
+    case 0: // _return
+      return;
+  }
+}
+
+
+function __ZN8sphere_tC1ERK3v_td(_this, _v, _d) {
+  var __label__ = 18; /* _entry */
+  while(1) switch(__label__) {
+    case 18: // _entry
+      var _this_addr = Pointer_make([0], 0);
+      var _v_addr = Pointer_make([0], 0);
+      var _d_addr;
+      var __alloca_point_ = 0;
+      HEAP[_this_addr] = _this;
+      HEAP[_v_addr] = _v;
+      _d_addr = _d;
+