BRIX5

    // The cross product is only valid for vectors with 3 dimensions,

    // with the exception of higher dimensional stuff that is beyond

    // the intended scope of this library.

    // Only a definition for N==3 is given below this class, using

    // cross() with another value for N will result in a link error.

    Vector cross(const Vector& v) const;

template <>

inline Vector<3> Vector<3>::cross(const Vector& v) const

{

    return Vector(

        p_vec[1] * v.p_vec[2] - p_vec[2] * v.p_vec[1],

        p_vec[2] * v.p_vec[0] - p_vec[0] * v.p_vec[2],

        p_vec[0] * v.p_vec[1] - p_vec[1] * v.p_vec[0]

    );

}

} // namespace