1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/prototype/vmath/ray.inl	Thu Jun 28 06:05:50 2012 +0300
     1.3 @@ -0,0 +1,70 @@
     1.4 +/*
     1.5 +libvmath - a vector math library
     1.6 +Copyright (C) 2004-2011 John Tsiombikas <nuclear@member.fsf.org>
     1.7 +
     1.8 +This program is free software: you can redistribute it and/or modify
     1.9 +it under the terms of the GNU Lesser General Public License as published
    1.10 +by the Free Software Foundation, either version 3 of the License, or
    1.11 +(at your option) any later version.
    1.12 +
    1.13 +This program is distributed in the hope that it will be useful,
    1.14 +but WITHOUT ANY WARRANTY; without even the implied warranty of
    1.15 +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    1.16 +GNU Lesser General Public License for more details.
    1.17 +
    1.18 +You should have received a copy of the GNU Lesser General Public License
    1.19 +along with this program.  If not, see <http://www.gnu.org/licenses/>.
    1.20 +*/
    1.21 +
    1.22 +#ifdef __cplusplus
    1.23 +extern "C" {
    1.24 +#endif	/* __cplusplus */
    1.25 +
    1.26 +static inline ray_t ray_cons(vec3_t origin, vec3_t dir)
    1.27 +{
    1.28 +	ray_t r;
    1.29 +	r.origin = origin;
    1.30 +	r.dir = dir;
    1.31 +	return r;
    1.32 +}
    1.33 +
    1.34 +#ifdef __cplusplus
    1.35 +}
    1.36 +
    1.37 +inline Ray reflect_ray(const Ray &inray, const Vector3 &norm)
    1.38 +{
    1.39 +	Ray ray = inray;
    1.40 +	ray.iter--;
    1.41 +	ray.dir = ray.dir.reflection(norm);
    1.42 +	return ray;
    1.43 +}
    1.44 +
    1.45 +inline Ray refract_ray(const Ray &inray, const Vector3 &norm, scalar_t mat_ior, bool entering, scalar_t ray_mag)
    1.46 +{
    1.47 +	Ray ray = inray;
    1.48 +	ray.iter--;
    1.49 +
    1.50 +	scalar_t ior = ray.calc_ior(entering, mat_ior);
    1.51 +	
    1.52 +	if(entering) {
    1.53 +		ray.enter(mat_ior);
    1.54 +	} else {
    1.55 +		ray.leave();
    1.56 +	}
    1.57 +
    1.58 +	if(ray_mag < 0.0) {
    1.59 +		ray_mag = ray.dir.length();
    1.60 +	}
    1.61 +	ray.dir = (ray.dir / ray_mag).refraction(norm, ior) * ray_mag;
    1.62 +
    1.63 +	/* check TIR */
    1.64 +	if(dot_product(ray.dir, norm) > 0.0) {
    1.65 +		if(entering) {
    1.66 +			ray.leave();
    1.67 +		} else {
    1.68 +			ray.enter(mat_ior);
    1.69 +		}
    1.70 +	}
    1.71 +	return ray;
    1.72 +}
    1.73 +#endif	/* __cplusplus */