mitkVis2Vector4.h

/*=========================================================================

  Program:   3DMed
  Date:      $Date: 2014-02-25 18:30:00 +0800 $
  Version:   $Version: 4.6.0 $
  Copyright: MIPG, Institute of Automation, Chinese Academy of Sciences

=========================================================================*/

#ifndef _mitkVis2Vector4_h
#define _mitkVis2Vector4_h

#include "mitkVisualization2Includes.h"

#include <math.h>

class MITK_VISUALIZATION2_API mitkVis2Vector4
{
public:
    float ele[4];
    mitkVis2Vector4();  
    mitkVis2Vector4(const mitkVis2Vector4 &v);
    mitkVis2Vector4(float x, float y, float z, float w = 1.0f);
    mitkVis2Vector4(float srcVector[4]);
    operator float* () { return ele; }  
    operator const float* () const { return ele; }

    mitkVis2Vector4& operator = (const mitkVis2Vector4 &a)
    {
        ele[0] = a.ele[0];
        ele[1] = a.ele[1];
        ele[2] = a.ele[2];
        ele[3] = a.ele[3];
        return *this; 
    }
    mitkVis2Vector4& operator *= (float);
    mitkVis2Vector4& operator += (const mitkVis2Vector4 &);
    mitkVis2Vector4& operator -= (const mitkVis2Vector4 &);

    void Normalize();
    float Length();
};

inline float operator * (const mitkVis2Vector4& A, const mitkVis2Vector4& B) 
{
    return A.ele[0]*B.ele[0] + A.ele[1]*B.ele[1] + A.ele[2]*B.ele[2];
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator % (const mitkVis2Vector4& A, const mitkVis2Vector4& B) 
{
    mitkVis2Vector4 Res;
    Res.ele[0] = A.ele[1] * B.ele[2] - A.ele[2] * B.ele[1]; 
    Res.ele[1] = A.ele[2] * B.ele[0] - A.ele[0] * B.ele[2]; 
    Res.ele[2] = A.ele[0] * B.ele[1] - A.ele[1] * B.ele[0]; 
    Res.ele[3] = 0.0f;
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator * (const mitkVis2Vector4 &V, float s) 
{
    mitkVis2Vector4 Res;
    Res.ele[0] = V.ele[0] * s;
    Res.ele[1] = V.ele[1] * s;
    Res.ele[2] = V.ele[2] * s;
    Res.ele[3] = V.ele[3] * s;  
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator * (float s, const mitkVis2Vector4 &V) 
{
    mitkVis2Vector4 Res;
    Res.ele[0] = V.ele[0] * s;
    Res.ele[1] = V.ele[1] * s;
    Res.ele[2] = V.ele[2] * s;
    Res.ele[3] = V.ele[3] * s;  
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4& mitkVis2Vector4::operator *= (float s) 
{
    ele[0] *= s;
    ele[1] *= s;
    ele[2] *= s;
    ele[3] *= s;    
    return *this;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator + (const mitkVis2Vector4& A, const mitkVis2Vector4& B) 
{
    mitkVis2Vector4 Res;
    Res.ele[0] = A.ele[0] + B.ele[0];
    Res.ele[1] = A.ele[1] + B.ele[1];
    Res.ele[2] = A.ele[2] + B.ele[2];
    Res.ele[3] = A.ele[3] + B.ele[3];
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator - (const mitkVis2Vector4& A, const mitkVis2Vector4& B) 
{
    mitkVis2Vector4 Res;
    Res.ele[0] = A.ele[0] - B.ele[0];
    Res.ele[1] = A.ele[1] - B.ele[1];
    Res.ele[2] = A.ele[2] - B.ele[2];
    Res.ele[3] = A.ele[3] - B.ele[3];
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 operator - (const mitkVis2Vector4& A)
{
    mitkVis2Vector4 Res;
    Res.ele[0] = - A.ele[0];
    Res.ele[1] = - A.ele[1];
    Res.ele[2] = - A.ele[2];
    Res.ele[3] = - A.ele[3];
    return Res;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 & mitkVis2Vector4::operator += (const mitkVis2Vector4 &B) 
{
    ele[0] += B.ele[0];
    ele[1] += B.ele[1];
    ele[2] += B.ele[2];
    ele[3] += B.ele[3];
    return *this;
}
//----------------------------------------------------------------------------
inline mitkVis2Vector4 & mitkVis2Vector4::operator -= (const mitkVis2Vector4 &B) 
{
    ele[0] -= B.ele[0];
    ele[1] -= B.ele[1];
    ele[2] -= B.ele[2];
    ele[3] -= B.ele[3];
    return *this;
}
//----------------------------------------------------------------------------
inline void mitkVis2Vector4::Normalize() 
{
    float vecLenInv = 1.0f / sqrtf(ele[0]*ele[0] + ele[1]*ele[1] + ele[2]*ele[2]);
    ele[0] *= vecLenInv;
    ele[1] *= vecLenInv;
    ele[2] *= vecLenInv;
}
//----------------------------------------------------------------------------
inline float mitkVis2Vector4::Length()
{
    return sqrtf(ele[0]*ele[0] + ele[1]*ele[1] + ele[2]*ele[2]);
}


#endif

---