1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/libs/vmath/ray.cc	Sat Feb 01 19:58:19 2014 +0200
     1.3 @@ -0,0 +1,61 @@
     1.4 +#include "ray.h"
     1.5 +#include "vector.h"
     1.6 +
     1.7 +scalar_t Ray::env_ior = 1.0;
     1.8 +
     1.9 +Ray::Ray()
    1.10 +{
    1.11 +	ior = env_ior;
    1.12 +	energy = 1.0;
    1.13 +	time = 0;
    1.14 +	iter = 0;
    1.15 +}
    1.16 +
    1.17 +Ray::Ray(const Vector3 &origin, const Vector3 &dir)
    1.18 +{
    1.19 +	this->origin = origin;
    1.20 +	this->dir = dir;
    1.21 +	ior = env_ior;
    1.22 +	energy = 1.0;
    1.23 +	time = 0;
    1.24 +	iter = 0;
    1.25 +}
    1.26 +
    1.27 +void Ray::transform(const Matrix4x4 &xform)
    1.28 +{
    1.29 +	Matrix4x4 upper = xform;
    1.30 +	upper[0][3] = upper[1][3] = upper[2][3] = upper[3][0] = upper[3][1] = upper[3][2] = 0.0;
    1.31 +	upper[3][3] = 1.0;
    1.32 +
    1.33 +	dir.transform(upper);
    1.34 +	origin.transform(xform);
    1.35 +}
    1.36 +
    1.37 +Ray Ray::transformed(const Matrix4x4 &xform) const
    1.38 +{
    1.39 +	Ray foo = *this;
    1.40 +	foo.transform(xform);
    1.41 +	return foo;
    1.42 +}
    1.43 +
    1.44 +void Ray::enter(scalar_t new_ior)
    1.45 +{
    1.46 +	ior = new_ior;
    1.47 +}
    1.48 +
    1.49 +void Ray::leave()
    1.50 +{
    1.51 +}
    1.52 +
    1.53 +scalar_t Ray::calc_ior(bool entering, scalar_t mat_ior) const
    1.54 +{
    1.55 +	scalar_t from_ior = this->ior;
    1.56 +
    1.57 +	if(entering) {
    1.58 +		return from_ior / mat_ior;
    1.59 +	}
    1.60 +
    1.61 +	Ray tmp = *this;
    1.62 +	tmp.leave();
    1.63 +	return from_ior / tmp.ior;
    1.64 +}