source: GTP/branches/IllumWPdeliver2008dec/IlluminationWP/demos/Standalone/Hierarchical Systems Demo [OpenGL]/RESOURCES/include/My3DGraphRes/Vector.h @ 3255

//************************************************************************* //
// 3D Vector osztály
//
// Szirmay-Kalos Laszlo, 2002. November.
//************************************************************************* //
#ifndef VECTOR_H
#define VECTOR_H

#include <includes.h>

//===============================================================
class Vector {
//===============================================================
public:
    float x, y, z,w;   // a Descartes koordináták
        Vector( ) { x = y = z = 0.0; }
        Vector( float x0, float y0, float z0, float w0 = 1.0 ) { x = x0; y = y0; z = z0;w=w0; }

        Vector operator+( const Vector& v ) { // két vektor összege
                float X = x + v.x, Y = y + v.y, Z = z + v.z, W=w+v.w;
                return Vector(X, Y, Z,W);
        }
        Vector operator-( const Vector& v ) { 
                float X = x - v.x, Y = y - v.y, Z = z - v.z,W=w-v.w;
                return Vector(X, Y, Z,W);
        }       
        Vector operator*( float f ) {         // vektor és szám szorzata
                return Vector( x * f, y * f, z * f ,w*f);
        }
        float operator*( const Vector& v ) {  // két vektor skaláris szorzata
                return (x * v.x + y * v.y + z * v.z);
        }
        Vector operator%( const Vector& v ) { // két vektor vektoriális szorzata
                float X = y * v.z - z * v.y, Y = z * v.x - x * v.z, Z = x * v.y - y * v.x;
                return Vector(X, Y, Z);
        }
        float Length( ) {                     // vektor abszolút értéke
                return (float)sqrt( x * x + y * y + z * z );
        }
        void operator+=( const Vector& v ) {  // vektor összeadás
                x += v.x; y += v.y; z += v.z;w+=v.w; 
        }
        void operator-=( const Vector& v ) {  // vektor különbség
                x -= v.x; y -= v.y; z -= v.z;w-=v.w; 
        }
        void operator*=( float f ) {              // vektor és szám szorzata
                x *= f; y *= f; z *= f;w*=f; 
        }
        Vector operator/( float f ) {             // vektor osztva egy számmal
                return Vector( x/f, y/f, z/f ,w/f); 
        }
        Vector Normalize( ) {                                     // vektor normalizálása
                float l = Length( );
                if ( l < 0.000001f) { x = 1; y = 0; z = 0; }
                else { x /= l; y /= l; z /= l;
                return *this;}
        }
        Vector UnitVector( ) {                            // egy vektorral párhuzamos egységvektor
                Vector r = * this;
                r.Normalize();
                return r;
        }
        Vector Rotate( Vector& axis, float angle ) {    // vektor forgatása egy tengely körül
                Vector iv = this -> UnitVector();
                Vector jv = axis.UnitVector() % this -> UnitVector();
                float radian = angle * M_PI/180;
                return (iv * cos(radian) +  jv * sin(radian));
        }
        Vector RotateX(float angle)
        {
                float radian = angle * M_PI/180;
                Vector ret;
                ret.x=x;
                ret.y=y*cos(radian)-z*sin(radian);
                ret.z=y*sin(radian)+z*cos(radian);
                return ret;
        }
        Vector RotateY(float angle)
        {
                float radian = angle * M_PI/180;
                Vector ret;
                ret.y=y;
                ret.x=x*cos(radian)+z*sin(radian);
                ret.z=-x*sin(radian)+z*cos(radian);
                return ret;
        }
        Vector RotateZ(float angle)
        {
                float radian = angle * M_PI/180;
                Vector ret;
                ret.z=z;
                ret.x=x*cos(radian)-y*sin(radian);
                ret.y=x*sin(radian)+y*cos(radian);
                return ret;
        }
        float * GetArray() { return &x; }
        float * GetArrayf() { return &x;}

        float& X() { return x; }
        float& Y() { return y; }
        float& Z() { return z; }
        float& W() { return w; }
};

//--------------------------------------------
class Matrix {
//--------------------------------------------
 public:
        float m[4][4];
        Matrix( ) { }
        Matrix( float d1, float d2, float d3 ) {
                Clear();
                m[0][0] = d1; m[1][1] = d2; m[2][2] = d3;
        }
        void Clear( ) { memset( &m[0][0], 0, sizeof( m ) ); }
        
        Vector operator*( const Vector& v ) { 
                return Vector(m[0][0] * v.x + m[0][1] * v.y + m[0][2] * v.z,
                                  m[1][0] * v.x + m[1][1] * v.y + m[1][2] * v.z,
                                          m[2][0] * v.x + m[2][1] * v.y + m[2][2] * v.z); 
        }

        Matrix operator*( const Matrix& mat ) {
                Matrix result;
                for( int i = 0; i < 3; i++ )
                        for( int j = 0; j < 3; j++ ) {
                                result.m[i][j] = 0;
                                for( int k = 0; k < 3; k++ ) result.m[i][j] += m[i][k] * mat.m[k][j];
                        }
                return result;
        }

        Matrix Transpose( ) {
                Matrix result;
                for( int i = 0; i < 3; i++ )
                        for( int j = 0; j < 3; j++ ) result.m[j][i] = m[i][j];
                return result;
        }
        Vector Transform( Vector& v ) {
                return Vector( v.x * m[0][0] + v.y * m[1][0] + v.z * m[2][0] + m[3][0],
                                   v.x * m[0][1] + v.y * m[1][1] + v.z * m[2][1] + m[3][1],
                                       v.x * m[0][2] + v.y * m[1][2] + v.z * m[2][2] + m[3][2]);
        }
        float * GetArray() { return &m[0][0]; }
};


#endif