ray.h

#pragma once
 
#include "vec.h"
#include "inv.h"
 
namespace cgv {
namespace math {
 
template<class T>
struct ray {
    vec<T> origin;
    vec<T> direction;
 
    ray(unsigned int n = 3) {
        origin.resize(n);
        origin.fill(T(0));
        direction.resize(n);
        direction.fill(T(0));
    }
 
    ray(const vec<T>& origin, const vec<T>& direction) : origin(origin), direction(direction) {
        assert(origin.dim() == directin.dim());
    }
 
    vec<T> position(T t) const {
        return origin + t * direction;
    }
};
 
typedef ray<float> rayn;
typedef ray<double> drayn;
 
} // namespace math
} // namespace cgv