nuclear@2: /*
nuclear@2: libvmath - a vector math library
nuclear@2: Copyright (C) 2004-2011 John Tsiombikas <nuclear@member.fsf.org>
nuclear@2: 
nuclear@2: This program is free software: you can redistribute it and/or modify
nuclear@2: it under the terms of the GNU Lesser General Public License as published
nuclear@2: by the Free Software Foundation, either version 3 of the License, or
nuclear@2: (at your option) any later version.
nuclear@2: 
nuclear@2: This program is distributed in the hope that it will be useful,
nuclear@2: but WITHOUT ANY WARRANTY; without even the implied warranty of
nuclear@2: MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
nuclear@2: GNU Lesser General Public License for more details.
nuclear@2: 
nuclear@2: You should have received a copy of the GNU Lesser General Public License
nuclear@2: along with this program.  If not, see <http://www.gnu.org/licenses/>.
nuclear@2: */
nuclear@2: 
nuclear@2: #ifndef VMATH_RAY_H_
nuclear@2: #define VMATH_RAY_H_
nuclear@2: 
nuclear@2: #include "matrix.h"
nuclear@2: #include "vector.h"
nuclear@2: 
nuclear@2: typedef struct {
nuclear@2: 	vec3_t origin, dir;
nuclear@2: } ray_t;
nuclear@2: 
nuclear@2: #ifdef __cplusplus
nuclear@2: extern "C" {
nuclear@2: #endif	/* __cplusplus */
nuclear@2: 
nuclear@2: static inline ray_t ray_cons(vec3_t origin, vec3_t dir);
nuclear@2: ray_t ray_transform(ray_t r, mat4_t m);
nuclear@2: 
nuclear@2: #ifdef __cplusplus
nuclear@2: }	/* __cplusplus */
nuclear@2: 
nuclear@2: #include <stack>
nuclear@2: 
nuclear@2: class Ray {
nuclear@2: private:
nuclear@2: 	std::stack<scalar_t> ior_stack;
nuclear@2: 
nuclear@2: public:
nuclear@2: 	/* enviornmental index of refraction, normally 1.0 */
nuclear@2: 	static scalar_t env_ior;
nuclear@2: 	
nuclear@2: 	Vector3 origin, dir;
nuclear@2: 	scalar_t energy;
nuclear@2: 	int iter;
nuclear@2: 	scalar_t ior;
nuclear@2: 	long time;
nuclear@2: 
nuclear@2: 	Ray();
nuclear@2: 	Ray(const Vector3 &origin, const Vector3 &dir);
nuclear@2: 
nuclear@2: 	void transform(const Matrix4x4 &xform);
nuclear@2: 	Ray transformed(const Matrix4x4 &xform) const;
nuclear@2: 
nuclear@2: 	void enter(scalar_t new_ior);
nuclear@2: 	void leave();
nuclear@2: 
nuclear@2: 	scalar_t calc_ior(bool entering, scalar_t mat_ior = 1.0) const;
nuclear@2: };
nuclear@2: 
nuclear@2: inline Ray reflect_ray(const Ray &inray, const Vector3 &norm);
nuclear@2: inline Ray refract_ray(const Ray &inray, const Vector3 &norm, scalar_t ior, bool entering, scalar_t ray_mag = -1.0);
nuclear@2: #endif	/* __cplusplus */
nuclear@2: 
nuclear@2: #include "ray.inl"
nuclear@2: 
nuclear@2: #endif	/* VMATH_RAY_H_ */