nuclear@28: #include "ray.h"
nuclear@28: #include "vector.h"
nuclear@28: 
nuclear@28: ray_t ray_transform(ray_t r, mat4_t xform)
nuclear@28: {
nuclear@28: 	mat4_t upper;
nuclear@28: 	vec3_t dir;
nuclear@28: 
nuclear@28: 	m4_copy(upper, xform);
nuclear@28: 	upper[0][3] = upper[1][3] = upper[2][3] = upper[3][0] = upper[3][1] = upper[3][2] = 0.0;
nuclear@28: 	upper[3][3] = 1.0;
nuclear@28: 
nuclear@28: 	dir = v3_sub(r.dir, r.origin);
nuclear@28: 	dir = v3_transform(dir, upper);
nuclear@28: 	r.origin = v3_transform(r.origin, xform);
nuclear@28: 	r.dir = v3_add(dir, r.origin);
nuclear@28: 	return r;
nuclear@28: }