nuclear@28: #ifndef VMATH_RAY_H_
nuclear@28: #define VMATH_RAY_H_
nuclear@28: 
nuclear@28: #include "matrix.h"
nuclear@28: #include "vector.h"
nuclear@28: 
nuclear@28: typedef struct {
nuclear@28: 	vec3_t origin, dir;
nuclear@28: } ray_t;
nuclear@28: 
nuclear@28: #ifdef __cplusplus
nuclear@28: extern "C" {
nuclear@28: #endif	/* __cplusplus */
nuclear@28: 
nuclear@28: static inline ray_t ray_cons(vec3_t origin, vec3_t dir);
nuclear@28: ray_t ray_transform(ray_t r, mat4_t m);
nuclear@28: 
nuclear@28: #ifdef __cplusplus
nuclear@28: }	/* __cplusplus */
nuclear@28: 
nuclear@28: #include <stack>
nuclear@28: 
nuclear@28: class Ray {
nuclear@28: private:
nuclear@28: 	std::stack<scalar_t> ior_stack;
nuclear@28: 
nuclear@28: public:
nuclear@28: 	/* enviornmental index of refraction, normally 1.0 */
nuclear@28: 	static scalar_t env_ior;
nuclear@28: 	
nuclear@28: 	Vector3 origin, dir;
nuclear@28: 	scalar_t energy;
nuclear@28: 	int iter;
nuclear@28: 	scalar_t ior;
nuclear@28: 	long time;
nuclear@28: 
nuclear@28: 	Ray();
nuclear@28: 	Ray(const Vector3 &origin, const Vector3 &dir);
nuclear@28: 
nuclear@28: 	void transform(const Matrix4x4 &xform);
nuclear@28: 	Ray transformed(const Matrix4x4 &xform) const;
nuclear@28: 
nuclear@28: 	void enter(scalar_t new_ior);
nuclear@28: 	void leave();
nuclear@28: 
nuclear@28: 	scalar_t calc_ior(bool entering, scalar_t mat_ior = 1.0) const;
nuclear@28: };
nuclear@28: 
nuclear@28: inline Ray reflect_ray(const Ray &inray, const Vector3 &norm);
nuclear@28: inline Ray refract_ray(const Ray &inray, const Vector3 &norm, scalar_t ior, bool entering, scalar_t ray_mag = -1.0);
nuclear@28: #endif	/* __cplusplus */
nuclear@28: 
nuclear@28: #include "ray.inl"
nuclear@28: 
nuclear@28: #endif	/* VMATH_RAY_H_ */