---

vector.cpp

//FILE:         vector.cpp
//AUTHOR:       Nathan Cournia <nathan@cournia.com>

#include "vector.h"

const math::vector3 math::vector3::ZERO( 0.0, 0.0, 0.0 );
const math::vector3 math::vector3::X_AXIS( 1.0, 0.0, 0.0 );
const math::vector3 math::vector3::Y_AXIS( 0.0, 1.0, 0.0 );
const math::vector3 math::vector3::Z_AXIS( 0.0, 0.0, 1.0 );

math::vector3 
operator* ( real_t lhs, const math::vector3& rhs )
{
        return math::vector3( (rhs.m_data[ 0 ] * lhs), 
                (rhs.m_data[ 1 ] * lhs),
                (rhs.m_data[ 2 ] * lhs)
        );
}

std::ostream& 
operator<< ( std::ostream& o, const math::vector3& vec )
{
        return o << vec.m_data[ 0 ] << " " << vec.m_data[ 1 ] 
                << " " << vec.m_data[ 2 ];
}


math::vector3 
math::cross( const math::vector3& u, const math::vector3& v )
{
        return math::vector3( 
                u.m_data[ 1 ] * v.m_data[ 2 ] - u.m_data[ 2 ] * v.m_data[ 1 ],
                u.m_data[ 2 ] * v.m_data[ 0 ] - u.m_data[ 0 ] * v.m_data[ 2 ],
                u.m_data[ 0 ] * v.m_data[ 1 ] - u.m_data[ 1 ] * v.m_data[ 0 ] 
        );
}

void 
math::add( math::vector3& dest, const math::vector3& lhs, 
        const math::vector3& rhs )
{
        dest.m_data[ 0 ] = lhs.m_data[ 0 ] + rhs.m_data[ 0 ];
        dest.m_data[ 1 ] = lhs.m_data[ 1 ] + rhs.m_data[ 1 ];
        dest.m_data[ 2 ] = lhs.m_data[ 2 ] + rhs.m_data[ 2 ];
}

void 
math::subtract( math::vector3& dest, const math::vector3& lhs, 
        const math::vector3& rhs )
{
        dest.m_data[ 0 ] = lhs.m_data[ 0 ] - rhs.m_data[ 0 ];
        dest.m_data[ 1 ] = lhs.m_data[ 1 ] - rhs.m_data[ 1 ];
        dest.m_data[ 2 ] = lhs.m_data[ 2 ] - rhs.m_data[ 2 ];
}

void 
math::multiply( math::vector3& dest, const math::vector3& lhs, 
        const math::vector3& rhs )
{
        dest.m_data[ 0 ] = lhs.m_data[ 0 ] * rhs.m_data[ 0 ];
        dest.m_data[ 1 ] = lhs.m_data[ 1 ] * rhs.m_data[ 1 ];
        dest.m_data[ 2 ] = lhs.m_data[ 2 ] * rhs.m_data[ 2 ];
}

void 
math::divide( math::vector3& dest, const math::vector3& lhs, 
        const math::vector3& rhs )
{
        dest.m_data[ 0 ] = lhs.m_data[ 0 ] / rhs.m_data[ 0 ];
        dest.m_data[ 1 ] = lhs.m_data[ 1 ] / rhs.m_data[ 1 ];
        dest.m_data[ 2 ] = lhs.m_data[ 2 ] / rhs.m_data[ 2 ];
}

void 
math::multiply( math::vector3& dest, const math::vector3& lhs, real_t rhs )
{
        dest.m_data[ 0 ] = lhs.m_data[ 0 ] * rhs;
        dest.m_data[ 1 ] = lhs.m_data[ 1 ] * rhs;
        dest.m_data[ 2 ] = lhs.m_data[ 2 ] * rhs;
}

void 
math::multiply( math::vector3& dest, real_t lhs, const math::vector3& rhs )
{
        dest.m_data[ 0 ] = rhs.m_data[ 0 ] * lhs;
        dest.m_data[ 1 ] = rhs.m_data[ 1 ] * lhs;
        dest.m_data[ 2 ] = rhs.m_data[ 2 ] * lhs;
}

real_t 
math::dot( const math::vector3& lhs, const math::vector3& rhs )
{
        return ( 
                (lhs.m_data[ 0 ] * rhs.m_data[ 0 ]) +
                (lhs.m_data[ 1 ] * rhs.m_data[ 1 ]) +
                (lhs.m_data[ 2 ] * rhs.m_data[ 2 ])
        );
}

---