geometry.impl

#include <math.h>

template <unsigned int n>
const double Vector<n>::NORMAL_EPSILON = 1e-6;

template <unsigned int n>
const double Object<n>::INTERSECT_EPSILON = 1e-6;

/*
 * Implementation of Point class.
 */

template <unsigned int n>
Point<n>::Point() {
    assert (n > 0);
    data(n,0) = 1.0;
}

template <unsigned int n>
Point<n>::Point (const double data[n]) {
    assert (n > 0);
    for (unsigned int i = 0; i < n; i ++)
        this->data(i,0) = data[i];
    this->data(n,0) = 1.0;
}

template <unsigned int n>
std::ostream& operator<< (std::ostream &os, const Point<n> &p) {
    os << "p<";
    for (unsigned int i = 0; i < n; i ++) {
        if (i > 0)
            os << ',';
        os << p[i];
    }
    os << '>';
    return os;
}

/*
 * Implementation of Vector class.
 */

template <unsigned int n>
Vector<n>::Vector (const double data[n]) {
    assert (n > 0);
    for (unsigned int i = 0; i < n; i ++)
        this->data(i,0) = data[i];
}

template <unsigned int n>
double Vector<n>::getLength() const {
    double lensq = 0.0;
    for (unsigned int i = 0; i < n; i ++)
        lensq += data(i,0)*data(i,0);
    return sqrt(lensq);
}

template <unsigned int n>
void Vector<n>::normalize() {
    double len = getLength();
    assert (len > 0.0);
    for (unsigned int i = 0; i < n; i ++)
        data(i,0) /= len;
}

template <unsigned int n>
void Vector<n>::translateGL() const {
    assert (n <= 3);
    if (n == 1) {
        glTranslated (data(0,0), 0.0,       0.0);
    } else if (n == 2) {
        glTranslated (data(0,0), data(1,0), 0.0);
    } else {
        glTranslated (data(0,0), data(1,0), data(2,0));
    }
}

template <unsigned int n>
Vector<n> Vector<n>::operator-() const {
    Vector<n> v;
    for (unsigned int i = 0; i < n; i ++)
        v.data(i,0) = -data(i,0);
    return v;
}

template <unsigned int n>
double operator* (const Vector<n> &a, const Vector<n> &b) {
    double dot = 0.0;
    for (unsigned int i = 0; i < n; i ++)
        dot += a[i] * b[i];
    return dot;
}

template <unsigned int n>
std::ostream& operator<< (std::ostream &os, const Vector<n> &p) {
    os << "v<";
    for (unsigned int i = 0; i < n; i ++) {
        if (i > 0)
            os << ',';
        os << p[i];
    }
    os << '>';
    return os;
}

template <unsigned int n>
Vector<n> operator- (const Point<n> &a, const Point<n> &b) {
    Vector<n> v;
    for (unsigned int i = 0; i < n; i ++)
        v.data(i,0) = a[i] - b[i];
    return v;
}

/*
 * Implementation of Ray class.
 */

template <unsigned int n>
Ray<n>::Ray() {
    assert (n > 0);
    Matrix<n+1,1> vdata;
    for (unsigned int i = 0; i < n; i ++)
        vdata(i,0) = 1.0 / static_cast<double>(n);
    dir.setMatrix (vdata);
}

template <unsigned int n>
std::ostream& operator<< (std::ostream &os, const Ray<n> &r) {
    os << "ray(" << r.getStart() << ',' << r.getDirection() << ")";
    return os;
}

---