nuclear@1: /*
nuclear@1: libvmath - a vector math library
nuclear@1: Copyright (C) 2004-2011 John Tsiombikas <nuclear@member.fsf.org>
nuclear@1: 
nuclear@1: This program is free software: you can redistribute it and/or modify
nuclear@1: it under the terms of the GNU Lesser General Public License as published
nuclear@1: by the Free Software Foundation, either version 3 of the License, or
nuclear@1: (at your option) any later version.
nuclear@1: 
nuclear@1: This program is distributed in the hope that it will be useful,
nuclear@1: but WITHOUT ANY WARRANTY; without even the implied warranty of
nuclear@1: MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
nuclear@1: GNU Lesser General Public License for more details.
nuclear@1: 
nuclear@1: You should have received a copy of the GNU Lesser General Public License
nuclear@1: along with this program.  If not, see <http://www.gnu.org/licenses/>.
nuclear@1: */
nuclear@1: 
nuclear@1: #ifdef __cplusplus
nuclear@1: extern "C" {
nuclear@1: #endif	/* __cplusplus */
nuclear@1: 
nuclear@1: static inline ray_t ray_cons(vec3_t origin, vec3_t dir)
nuclear@1: {
nuclear@1: 	ray_t r;
nuclear@1: 	r.origin = origin;
nuclear@1: 	r.dir = dir;
nuclear@1: 	return r;
nuclear@1: }
nuclear@1: 
nuclear@1: #ifdef __cplusplus
nuclear@1: }
nuclear@1: 
nuclear@1: inline Ray reflect_ray(const Ray &inray, const Vector3 &norm)
nuclear@1: {
nuclear@1: 	Ray ray = inray;
nuclear@1: 	ray.iter--;
nuclear@1: 	ray.dir = ray.dir.reflection(norm);
nuclear@1: 	return ray;
nuclear@1: }
nuclear@1: 
nuclear@1: inline Ray refract_ray(const Ray &inray, const Vector3 &norm, scalar_t mat_ior, bool entering, scalar_t ray_mag)
nuclear@1: {
nuclear@1: 	Ray ray = inray;
nuclear@1: 	ray.iter--;
nuclear@1: 
nuclear@1: 	scalar_t ior = ray.calc_ior(entering, mat_ior);
nuclear@1: 	
nuclear@1: 	if(entering) {
nuclear@1: 		ray.enter(mat_ior);
nuclear@1: 	} else {
nuclear@1: 		ray.leave();
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(ray_mag < 0.0) {
nuclear@1: 		ray_mag = ray.dir.length();
nuclear@1: 	}
nuclear@1: 	ray.dir = (ray.dir / ray_mag).refraction(norm, ior) * ray_mag;
nuclear@1: 
nuclear@1: 	/* check TIR */
nuclear@1: 	if(dot_product(ray.dir, norm) > 0.0) {
nuclear@1: 		if(entering) {
nuclear@1: 			ray.leave();
nuclear@1: 		} else {
nuclear@1: 			ray.enter(mat_ior);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	return ray;
nuclear@1: }
nuclear@1: #endif	/* __cplusplus */