X-Git-Url: http://git.mutantstargoat.com/user/nuclear/?p=csgray;a=blobdiff_plain;f=src%2Fcsgray.c;h=c9d2ed2e820336bb96d7f86b396ecef8975d5cab;hp=53b01468d021feb2a52c32087b742fc0dcb08d7f;hb=aab39eb873290ce4b268fe4ae303def6bf932bc0;hpb=3dd5ba3e29ccff62d188bb4e9e8d23d1879a6024

diff --git a/src/csgray.c b/src/csgray.c
index 53b0146..c9d2ed2 100644
--- a/src/csgray.c
+++ b/src/csgray.c
@@ -3,6 +3,7 @@
 #include <string.h>
 #include <math.h>
 #include <float.h>
+#include <treestore.h>
 #include "csgimpl.h"
 #include "matrix.h"
 #include "geom.h"
@@ -12,6 +13,7 @@ static int ray_trace(struct ray *ray, float *col);
 static void shade(float *col, struct ray *ray, struct hit *hit);
 static void background(float *col, struct ray *ray);
 static int find_intersection(struct ray *ray, struct hit *best);
+static csg_object *load_object(struct ts_node *node);
 
 static float ambient[3];
 static struct camera cam;
@@ -23,7 +25,7 @@ int csg_init(void)
 	oblist = 0;
 	plights = 0;
 
-	csg_ambient(0.05, 0.05, 0.05);
+	csg_ambient(0, 0, 0);
 	csg_view(0, 0, 5, 0, 0, 0);
 	csg_fov(50);
 
@@ -42,12 +44,30 @@ void csg_destroy(void)
 
 void csg_view(float x, float y, float z, float tx, float ty, float tz)
 {
+	float dir[3];
+	float len;
+
 	cam.x = x;
 	cam.y = y;
 	cam.z = z;
 	cam.tx = tx;
 	cam.ty = ty;
 	cam.tz = tz;
+
+	dir[0] = tx - x;
+	dir[1] = ty - y;
+	dir[2] = tz - z;
+	len = sqrt(dir[0] * dir[0] + dir[1] * dir[1] + dir[2] * dir[2]);
+
+	if(1.0f - fabs(ty - y) / len < 1e-6f) {
+		cam.ux = cam.uy = 0.0f;
+		cam.uz = -1.0f;
+	} else {
+		cam.ux = cam.uz = 0.0f;
+		cam.uy = 1.0f;
+	}
+
+	mat4_lookat(cam.xform, x, y, z, tx, ty, tz, cam.ux, cam.uy, cam.uz);
 }
 
 void csg_fov(float fov)
@@ -58,12 +78,49 @@ void csg_fov(float fov)
 
 int csg_load(const char *fname)
 {
-	return 0;	/* TODO */
+	struct ts_node *root = 0, *c;
+	csg_object *o;
+
+	if(!(root = ts_load(fname))) {
+		fprintf(stderr, "failed to open %s\n", fname);
+		return -1;
+	}
+	if(strcmp(root->name, "csgray_scene") != 0) {
+		fprintf(stderr, "invalid scene file: %s\n", fname);
+		goto err;
+	}
+
+	c = root->child_list;
+	while(c) {
+		if(strcmp(c->name, "viewer") == 0) {
+			static float def_pos[] = {0, 0, 5};
+			static float def_targ[] = {0, 0, 0};
+
+			float *p = ts_get_attr_vec(c, "position", def_pos);
+			float *t = ts_get_attr_vec(c, "target", def_targ);
+
+			csg_view(p[0], p[1], p[2], t[0], t[1], t[2]);
+			csg_fov(ts_get_attr_num(c, "fov", 50.0f));
+
+		} else if((o = load_object(c))) {
+			csg_add_object(o);
+		}
+		c = c->next;
+	}
+
+	ts_free_tree(root);
+	return 0;
+
+err:
+	if(root) {
+		ts_free_tree(root);
+	}
+	return -1;
 }
 
 int csg_save(const char *fname)
 {
-	return 0;	/* TODO */
+	return -1;	/* TODO */
 }
 
 void csg_add_object(csg_object *o)
@@ -198,7 +255,7 @@ csg_object *csg_plane(float x, float y, float z, float nx, float ny, float nz)
 	o->plane.ny = ny;
 	o->plane.nz = nz;
 	o->plane.d = x * nx + y * ny + z * nz;
-	return 0;
+	return o;
 }
 
 csg_object *csg_box(float x, float y, float z, float xsz, float ysz, float zsz)
@@ -273,6 +330,40 @@ void csg_opacity(csg_object *o, float p)
 	o->ob.opacity = p;
 }
 
+void csg_metallic(csg_object *o, int m)
+{
+	o->ob.metallic = m;
+}
+
+void csg_reset_xform(csg_object *o)
+{
+	mat4_identity(o->ob.xform);
+	mat4_identity(o->ob.inv_xform);
+}
+
+void csg_translate(csg_object *o, float x, float y, float z)
+{
+	mat4_translate(o->ob.xform, x, y, z);
+	mat4_pre_translate(o->ob.inv_xform, -x, -y, -z);
+}
+
+void csg_rotate(csg_object *o, float angle, float x, float y, float z)
+{
+	mat4_rotate(o->ob.xform, angle, x, y, z);
+	mat4_pre_rotate(o->ob.inv_xform, -angle, x, y, z);
+}
+
+void csg_scale(csg_object *o, float x, float y, float z)
+{
+	mat4_scale(o->ob.xform, x, y, z);
+	mat4_pre_scale(o->ob.inv_xform, 1.0f / x, 1.0f / y, 1.0f / z);
+}
+
+void csg_lookat(csg_object *o, float x, float y, float z, float tx, float ty, float tz, float ux, float uy, float uz)
+{
+	mat4_lookat(o->ob.xform, x, y, z, tx, ty, tz, ux, uy, uz);
+	mat4_inv_lookat(o->ob.inv_xform, x, y, z, tx, ty, tz, ux, uy, uz);
+}
 
 void csg_render_pixel(int x, int y, int width, int height, float aspect, float *color)
 {
@@ -297,14 +388,15 @@ void csg_render_image(float *pixels, int width, int height)
 
 static void calc_primary_ray(struct ray *ray, int x, int y, int w, int h, float aspect)
 {
-	/* TODO */
 	ray->dx = aspect * ((float)x / (float)w * 2.0f - 1.0f);
 	ray->dy = 1.0f - (float)y / (float)h * 2.0f;
 	ray->dz = -1.0f / tan(cam.fov * 0.5f);
 
-	ray->x = cam.x;
-	ray->y = cam.y;
-	ray->z = cam.z;
+	ray->x = 0;
+	ray->y = 0;
+	ray->z = 0;
+
+	xform_ray(ray, cam.xform);
 }
 
 static int ray_trace(struct ray *ray, float *col)
@@ -326,10 +418,10 @@ static int ray_trace(struct ray *ray, float *col)
 
 static void shade(float *col, struct ray *ray, struct hit *hit)
 {
-	float ndotl, len;
+	float ndotl, ndoth, len, falloff, spec;
 	csg_object *o, *lt = plights;
 	float dcol[3], scol[3] = {0};
-	float ldir[3];
+	float ldir[3], lcol[3], hdir[3];
 	struct ray sray;
 	struct hit tmphit;
 
@@ -350,24 +442,57 @@ static void shade(float *col, struct ray *ray, struct hit *hit)
 		sray.dy = ldir[1];
 		sray.dz = ldir[2];
 
-		if(!find_intersection(&sray, &tmphit) || tmphit.t < 1.0f) {
+		if(!find_intersection(&sray, &tmphit) || tmphit.t < 0.000001 || tmphit.t > 1.0f) {
 			if((len = sqrt(ldir[0] * ldir[0] + ldir[1] * ldir[1] + ldir[2] * ldir[2])) != 0.0f) {
 				float s = 1.0f / len;
 				ldir[0] *= s;
 				ldir[1] *= s;
 				ldir[2] *= s;
 			}
+			falloff = 1.0f / (len * len);
+
+			lcol[0] = lt->ob.emr * falloff;
+			lcol[1] = lt->ob.emg * falloff;
+			lcol[2] = lt->ob.emb * falloff;
 
 			if((ndotl = hit->nx * ldir[0] + hit->ny * ldir[1] + hit->nz * ldir[2]) < 0.0f) {
 				ndotl = 0.0f;
 			}
 
-			dcol[0] += o->ob.r * lt->ob.emr * ndotl;
-			dcol[1] += o->ob.g * lt->ob.emg * ndotl;
-			dcol[2] += o->ob.b * lt->ob.emb * ndotl;
+			dcol[0] += o->ob.r * lcol[0] * ndotl;
+			dcol[1] += o->ob.g * lcol[1] * ndotl;
+			dcol[2] += o->ob.b * lcol[2] * ndotl;
+
+			if(o->ob.roughness < 1.0f) {
+				float gloss = 1.0f - o->ob.roughness;
+
+				hdir[0] = ldir[0] - ray->dx;
+				hdir[1] = ldir[1] - ray->dy;
+				hdir[2] = ldir[2] - ray->dz;
+				if((len = sqrt(hdir[0] * hdir[0] + hdir[1] * hdir[1] + hdir[2] * hdir[2])) != 0.0f) {
+					float s = 1.0f / len;
+					hdir[0] *= s;
+					hdir[1] *= s;
+					hdir[2] *= s;
+				}
+
+				if((ndoth = hit->nx * hdir[0] + hit->ny * hdir[1] + hit->nz * hdir[2]) < 0.0f) {
+					ndoth = 0.0f;
+				}
+				spec = gloss * pow(ndoth, 100.0f * gloss);
+
+				if(o->ob.metallic) {
+					lcol[0] *= o->ob.r;
+					lcol[1] *= o->ob.g;
+					lcol[2] *= o->ob.b;
+				}
+				scol[0] += lcol[0] * spec;
+				scol[1] += lcol[1] * spec;
+				scol[2] += lcol[2] * spec;
+			}
 		}
 
-		lt = lt->ob.next;
+		lt = lt->ob.plt_next;
 	}
 
 	col[0] = dcol[0] + scol[0];
@@ -382,16 +507,32 @@ static void background(float *col, struct ray *ray)
 
 static int find_intersection(struct ray *ray, struct hit *best)
 {
+	int idx = 0;
 	csg_object *o;
-	struct hit *hit;
+	struct hinterv *hit, *it;
 
 	best->t = FLT_MAX;
 	best->o = 0;
 
 	o = oblist;
 	while(o) {
-		if((hit = ray_intersect(ray, o)) && hit->t < best->t) {
-			*best = *hit;
+		if((hit = ray_intersect(ray, o))) {
+			it = hit;
+			while(it) {
+				if(it->end[0].t > 1e-6) {
+					idx = 0;
+					break;
+				}
+				if(it->end[1].t > 1e-6) {
+					idx = 1;
+					break;
+				}
+				it = it->next;
+			}
+
+			if(it && it->end[idx].t < best->t) {
+				*best = it->end[idx];
+			}
 		}
 		free_hit_list(hit);
 		o = o->ob.next;
@@ -399,3 +540,122 @@ static int find_intersection(struct ray *ray, struct hit *best)
 
 	return best->o != 0;
 }
+
+static csg_object *load_object(struct ts_node *node)
+{
+	float *avec;
+	struct ts_node *c;
+	csg_object *sub, *o = 0, *olist = 0, *otail = 0;
+	int num_subobj = 0, is_csgop = 0;
+
+	if(strcmp(node->name, "null") == 0) {
+		if(!(o = csg_null(0, 0, 0))) {
+			goto err;
+		}
+
+	} else if(strcmp(node->name, "sphere") == 0) {
+		float rad = ts_get_attr_num(node, "radius", 1.0f);
+		if(!(o = csg_sphere(0, 0, 0, rad))) {
+			goto err;
+		}
+
+	} else if(strcmp(node->name, "plane") == 0) {
+		static float def_norm[] = {0, 1, 0};
+		float *norm = ts_get_attr_vec(node, "normal", def_norm);
+		if(!(o = csg_plane(0, 0, 0, norm[0], norm[1], norm[2]))) {
+			goto err;
+		}
+
+	} else if(strcmp(node->name, "box") == 0) {
+		static float def_sz[] = {1, 1, 1};
+		float *sz = ts_get_attr_vec(node, "size", def_sz);
+		if(!(o = csg_box(0, 0, 0, sz[0], sz[1], sz[2]))) {
+			goto err;
+		}
+
+	} else if(strcmp(node->name, "union") == 0) {
+		if(!(o = csg_union(0, 0))) {
+			goto err;
+		}
+		is_csgop = 1;
+
+	} else if(strcmp(node->name, "intersect") == 0) {
+		if(!(o = csg_intersection(0, 0))) {
+			goto err;
+		}
+		is_csgop = 1;
+
+	} else if(strcmp(node->name, "subtract") == 0) {
+		if(!(o = csg_subtraction(0, 0))) {
+			goto err;
+		}
+		is_csgop = 1;
+
+	} else {
+		goto err;
+	}
+
+	if(is_csgop) {
+		c = node->child_list;
+		while(c) {
+			if((sub = load_object(c))) {
+				if(olist) {
+					otail->ob.next = sub;
+					otail = sub;
+				} else {
+					olist = otail = sub;
+				}
+				++num_subobj;
+			}
+			c = c->next;
+		}
+
+		if(num_subobj != 2) {
+			goto err;
+		}
+		o->un.a = olist;
+		o->un.b = olist->ob.next;
+		olist->ob.next = 0;
+	}
+
+	if((avec = ts_get_attr_vec(node, "position", 0))) {
+		csg_translate(o, avec[0], avec[1], avec[2]);
+	}
+	if((avec = ts_get_attr_vec(node, "rotaxis", 0))) {
+		csg_rotate(o, ts_get_attr_num(node, "rotangle", 0.0f), avec[0], avec[1], avec[2]);
+	}
+	if((avec = ts_get_attr_vec(node, "scaling", 0))) {
+		csg_scale(o, avec[0], avec[1], avec[2]);
+	}
+	if((avec = ts_get_attr_vec(node, "target", 0))) {
+		/* don't move this before position */
+		float def_up[] = {0, 1, 0};
+		float *up = ts_get_attr_vec(node, "up", def_up);
+		float x = o->ob.xform[12];
+		float y = o->ob.xform[13];
+		float z = o->ob.xform[14];
+		csg_lookat(o, x, y, z, avec[0], avec[1], avec[2], up[0], up[1], up[2]);
+	}
+
+	if((avec = ts_get_attr_vec(node, "color", 0))) {
+		csg_color(o, avec[0], avec[1], avec[2]);
+	}
+	if((avec = ts_get_attr_vec(node, "emission", 0))) {
+		csg_emission(o, avec[0], avec[1], avec[2]);
+	}
+
+	csg_roughness(o, ts_get_attr_num(node, "roughness", o->ob.roughness));
+	csg_opacity(o, ts_get_attr_num(node, "opacity", o->ob.opacity));
+	csg_metallic(o, ts_get_attr_int(node, "metallic", o->ob.metallic));
+
+	return o;
+
+err:
+	csg_free_object(o);
+	while(olist) {
+		o = olist;
+		olist = olist->ob.next;
+		csg_free_object(o);
+	}
+	return 0;
+}