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