ray.h

#pragma once
 
#include "fmat.h"
#include "fvec.h"
#include "inv.h"
 
namespace cgv {
    namespace math {
 
template <typename T, cgv::type::uint32_type N>
class ray {
public:
    fvec<T, N> origin;
    fvec<T, N> direction;
 
    ray() {
        origin = fvec<T, N>(0);
        direction = fvec<T, N>(0);
    }
 
    ray(const fvec<T, N>& origin, const fvec<T, N>& direction) : origin(origin), direction(direction) {}
 
    fvec<T, N> position(float t) const {
 
        return origin + t * direction;
    }
};
 
template <typename T>
class ray<T, 3> {
public:
    fvec<T, 3> origin;
    fvec<T, 3> direction;
 
    ray() {
        origin = fvec<T, 3>(0);
        direction = fvec<T, 3>(0);
    }
 
    ray(const fvec<T, 3>& origin, const fvec<T, 3>& direction) : origin(origin), direction(direction) {}
 
    ray(const fvec<T, 2>& screen_coord, const fvec<unsigned, 2>& viewport_size, const fvec<T, 3>& eye_position, const fmat<T, 4, 4>& view_projection_matrix) {
 
        fvec<T, 2> window_coord = T(2) * screen_coord / static_cast<fvec<T, 2>>(viewport_size) - T(1);
 
        fvec<T, 4> world_coord(window_coord.x(), window_coord.y(), T(1), T(1));
        world_coord = inv(view_projection_matrix) * world_coord;
        world_coord /= world_coord.w();
 
        origin = eye_position;
        direction = normalize(fvec<T, 3>(world_coord) - eye_position);
    }
 
    fvec<T, 3> position(float t) const {
 
        return origin + t * direction;
    }
};
 
// Definition of some typical ray types
typedef ray<unsigned, 2> uray2;
typedef ray<unsigned, 3> uray3;
typedef ray<unsigned, 3> uray4;
typedef ray<int, 2> iray2;
typedef ray<int, 3> iray3;
typedef ray<int, 3> iray4;
typedef ray<float, 2> ray2;
typedef ray<float, 3> ray3;
typedef ray<float, 3> ray4;
typedef ray<double, 2> dray2;
typedef ray<double, 3> dray3;
typedef ray<double, 3> dray4;
 
    }
}